记录编号 |
437938 |
评测结果 |
AAAAAAAAAA |
题目名称 |
斐波那契平方和 |
最终得分 |
100 |
用户昵称 |
HeHe |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
0.003 s |
提交时间 |
2017-08-14 21:35:57 |
内存使用 |
0.32 MiB |
显示代码纯文本
#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
#define MOD (1000000007)
typedef long long LL;
inline void Pow(LL N);
LL res[2][2] = {{1, 0},
{0, 1}};
LL fib[2] = {0, 1};
LL tmp[2][2];
LL q[2][2] = {{0, 1},
{1, 1}};
LL N;
int main() {
#ifndef LOCAL
freopen("fibsqr.in", "r", stdin);
freopen("fibsqr.out", "w", stdout);
#endif
scanf("%lld", &N);
Pow(N);
memset(tmp, 0x00, sizeof(tmp));
for(int i = 0; i < 2; ++i)
for(int j = 0; j < 2; ++j)
(tmp[0][i] += fib[j] * res[j][i]) %= MOD;
printf("%lld", tmp[0][0] * tmp[0][1] % MOD);
return 0;
}
inline void Pow(LL N) {
while(N) {
if(N & 1) {
memset(tmp, 0x00, sizeof(tmp));
for(char i = 0; i < 2; ++i)
for(char j = 0; j < 2; ++j)
for(char k = 0; k < 2; ++k)
(tmp[i][j] += res[i][k] * q[k][j]) %= MOD;
for(char i = 0; i < 2; ++i)
for(char j = 0; j < 2; ++j)
res[i][j] = tmp[i][j];
}
memset(tmp, 0x00, sizeof(tmp));
for(char i = 0; i < 2; ++i)
for(char j = 0; j < 2; ++j)
for(char k = 0; k < 2; ++k)
(tmp[i][j] += q[i][k] * q[k][j]) %= MOD;
for(char i = 0; i < 2; ++i)
for(char j = 0; j < 2; ++j)
q[i][j] = tmp[i][j];
N >>= 1;
}
return ;
}