結果

問題 No.720 行列のできるフィボナッチ数列道場 (2)
ユーザー kotatsugamekotatsugame
提出日時 2020-03-27 05:01:43
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,058 bytes
コンパイル時間 639 ms
コンパイル使用メモリ 67,316 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-08-30 16:05:07
合計ジャッジ時間 2,075 ms
ジャッジサーバーID
(参考情報)
judge13 / judge11
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,376 KB
testcase_01 AC 2 ms
4,380 KB
testcase_02 AC 1 ms
4,376 KB
testcase_03 AC 2 ms
4,380 KB
testcase_04 AC 1 ms
4,376 KB
testcase_05 AC 1 ms
4,380 KB
testcase_06 AC 2 ms
4,376 KB
testcase_07 AC 2 ms
4,380 KB
testcase_08 AC 2 ms
4,376 KB
testcase_09 AC 2 ms
4,380 KB
testcase_10 AC 2 ms
4,376 KB
testcase_11 AC 2 ms
4,376 KB
testcase_12 AC 1 ms
4,380 KB
testcase_13 AC 1 ms
4,376 KB
testcase_14 AC 1 ms
4,376 KB
testcase_15 AC 2 ms
4,376 KB
testcase_16 AC 2 ms
4,380 KB
testcase_17 AC 2 ms
4,376 KB
testcase_18 AC 2 ms
4,376 KB
testcase_19 AC 2 ms
4,376 KB
testcase_20 AC 1 ms
4,376 KB
testcase_21 AC 2 ms
4,380 KB
testcase_22 AC 2 ms
4,376 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:31:1: 警告: ISO C++ では型の無い ‘main’ の宣言を禁止しています [-Wreturn-type]
   31 | main()
      | ^~~~

ソースコード

diff #

#include<iostream>
#include<array>
using namespace std;
const long mod=1e9+7;
template<int N>
array<array<long,N>,N>pw(array<array<long,N>,N>A,long K)
{
	array<array<long,N>,N>ret;
	for(int i=0;i<N;i++)for(int j=0;j<N;j++)ret[i][j]=i==j?1:0;
	while(K)
	{
		if(K&1)
		{
			array<array<long,N>,N>tmp;
			for(int i=0;i<N;i++)for(int j=0;j<N;j++)tmp[i][j]=ret[i][j],ret[i][j]=0;
			for(int i=0;i<N;i++)for(int j=0;j<N;j++)
			{
				for(int k=0;k<N;k++)(ret[i][j]+=tmp[i][k]*A[k][j])%=mod;
			}
		}
		K>>=1;
		array<array<long,N>,N>tmp;
		for(int i=0;i<N;i++)for(int j=0;j<N;j++)tmp[i][j]=A[i][j],A[i][j]=0;
		for(int i=0;i<N;i++)for(int j=0;j<N;j++)
		{
			for(int k=0;k<N;k++)(A[i][j]+=tmp[i][k]*tmp[k][j])%=mod;
		}
	}
	return ret;
}
main()
{
	long N,M;
	cin>>N>>M;
	array<array<long,2>,2>D;
	D[0][0]=1;
	D[0][1]=1;
	D[1][0]=1;
	D[1][1]=0;
	D=pw<2>(D,M);
	array<array<long,3>,3>F;
	F[0][0]=D[0][0];
	F[0][1]=D[0][1];
	F[1][0]=D[1][0];
	F[1][1]=D[1][1];
	F[0][2]=F[1][2]=0;
	F[2][0]=D[0][0];
	F[2][1]=D[0][1];
	F[2][2]=1;
	F=pw<3>(F,N);
	cout<<F[2][1]<<endl;
}
0