結果

問題 No.534 フィボナッチフィボナッチ数
ユーザー bal4ubal4u
提出日時 2019-04-19 17:31:54
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,121 bytes
コンパイル時間 207 ms
コンパイル使用メモリ 30,592 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-09-22 11:37:13
合計ジャッジ時間 1,310 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 42
権限があれば一括ダウンロードができます

ソースコード

diff # 

// yukicoder: No.534 フィボナッチフィボナッチ数
// 2019.4.19 bal4u

#include <stdio.h>

#define MOD0 1000000007
int MOD;

//// fibonacci数の高速計算
void multiply(int F[2][2], int M[2][2])
{
	long long f00, f01, f10, f11, x, y, z, w;
	f00 = F[0][0], f01 = F[0][1], f10 = F[1][0], f11 = F[1][1];
	x = (f00 * M[0][0] % MOD + f01 * M[1][0] % MOD);
	if (x >= MOD) x -= MOD;
	y = (f00 * M[0][1] % MOD + f01 * M[1][1] % MOD);
	if (y >= MOD) y -= MOD;
	z = (f10 * M[0][0] % MOD + f11 * M[1][0] % MOD);
	if (z >= MOD) z -= MOD;
	w = (f10 * M[0][1] % MOD + f11 * M[1][1] % MOD);
	if (w >= MOD) w -= MOD;
	F[0][0] = (int)x, F[0][1] = (int)y, F[1][0] = (int)z, F[1][1] = (int)w;

}

int M[2][2] = { {1,1},{1,0} };
void power(int F[2][2], long long n)
{
	if (n == 0 || n == 1) return;
	power(F, n >> 1);
	multiply(F, F);
	if (n & 1) multiply(F, M);
}

int fib(long long n)
{
	int F[2][2] = { {1,1},{1,0} };
	if (n == 0)	return 0;
	power(F, n - 1);
	return F[0][0];
}

int main()
{
	int ans;
	long long n;
	scanf("%lld", &n);
	MOD = (MOD0+1) << 1;
	ans = fib(n);
	MOD = MOD0;
	printf("%d\n", fib(ans));
	return 0;
}
0