結果

問題 No.1105 Many Triplets
ユーザー kotatsugame
提出日時 2020-07-03 22:10:39
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 755 bytes
コンパイル時間 499 ms
コンパイル使用メモリ 66,144 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-17 02:29:57
合計ジャッジ時間 1,426 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 25
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:6:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type]
    6 | main()
      | ^~~~

ソースコード

diff # 

#include<iostream>
using namespace std;
const long mod=1e9+7;
long N,A,B,C;
long X[3][3],E[3][3],T[3][3];
main()
{
	cin>>N>>A>>B>>C;
	N--;
	X[0][0]=X[1][1]=X[2][2]=1;
	X[0][1]=X[1][2]=X[2][0]=mod-1;
	E[0][0]=E[1][1]=E[2][2]=1;
	while(N)
	{
		if(N&1)
		{
			for(int i=0;i<3;i++)for(int j=0;j<3;j++)
			{
				T[i][j]=E[i][j];
				E[i][j]=0;
			}
			for(int i=0;i<3;i++)for(int j=0;j<3;j++)
			{
				for(int k=0;k<3;k++)(E[i][j]+=T[i][k]*X[k][j])%=mod;
			}
		}
		N>>=1;
		for(int i=0;i<3;i++)for(int j=0;j<3;j++)
		{
			T[i][j]=X[i][j];
			X[i][j]=0;
		}
		for(int i=0;i<3;i++)for(int j=0;j<3;j++)
		{
			for(int k=0;k<3;k++)(X[i][j]+=T[i][k]*T[k][j])%=mod;
		}
	}
	for(int i=0;i<3;i++)
	{
		cout<<(E[i][0]*A+E[i][1]*B+E[i][2]*C)%mod<<(i==2?"\n":" ");
	}
}
0