結果

問題 No.534 フィボナッチフィボナッチ数
ユーザー snrnsidy
提出日時 2021-10-31 18:57:14
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 3 ms / 2,000 ms
コード長 1,185 bytes
コンパイル時間 2,617 ms
コンパイル使用メモリ 201,320 KB
最終ジャッジ日時 2025-01-25 10:04:14
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 42
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h> 
using namespace std; 
#define ll long long int 
long long int MOD;

vector <vector<long long int>> mul(vector <vector<long long int>> a, vector<vector<long long int>> b)
{
	vector <vector<long long int>>c(2, vector<long long int>(2));
	c[0][0] = 0;
	c[0][1] = 0;
	c[1][0] = 0;
	c[1][1] = 0;

	for (int i = 0; i < 2; i++)
	{
			for (int j = 0; j < 2; j++)
			{
				for (int k = 0; k < 2; k++)
				{
					c[i][j] += (a[i][k] * b[k][j]);
					c[i][j] %= MOD;
				}
			}
	}

	return c;
}


int main(void)
{
	vector <vector <long long int>> a(2, vector<long long int>(2));
	vector <vector <long long int>> b(2, vector<long long int>(2));
	long long int n;
	a[0][0] = 1;
	a[0][1] = 1;
	a[1][0] = 1;
	a[1][1] = 0;

	b[0][0] = 1;
	b[0][1] = 0;
	b[1][0] = 0;
	b[1][1] = 1;

	MOD = 2e9 + 16;
	cin >> n;

	while (n > 0)
	{
		if (n & 1)
		{
			b = mul(a, b);
		}
		a = mul(a, a);
		n >>= 1;
	}

	n = b[1][0];

	MOD = 1e9 + 7;
	a[0][0] = 1;
	a[0][1] = 1;
	a[1][0] = 1;
	a[1][1] = 0;

	b[0][0] = 1;
	b[0][1] = 0;
	b[1][0] = 0;
	b[1][1] = 1;

	while (n > 0)
	{
		if (n & 1)
		{
			b = mul(a, b);
		}
		a = mul(a, a);
		n >>= 1;
	}
			
	cout << b[1][0] << '\n';

	return 0;
}
0