結果
| 問題 | No.1105 Many Triplets | 
| コンテスト | |
| ユーザー |  | 
| 提出日時 | 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()
      | ^~~~
            
            ソースコード
#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":" ");
	}
}
            
            
            
        