結果

問題 No.569 3 x N グリッドのパスの数
ユーザー WA_TLEWA_TLE
提出日時 2017-09-08 23:38:22
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 48 ms / 2,000 ms
コード長 2,397 bytes
コンパイル時間 1,374 ms
コンパイル使用メモリ 127,624 KB
最終ジャッジ日時 2025-01-05 02:50:43
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 60
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<deque>
#include<queue>
#include<vector>
#include<algorithm>
#include<iostream>
#include<set>
#include<cmath>
#include<tuple>
#include<string>
#include<chrono>
#include<functional>
#include<iterator>
#include<random>
#include<unordered_set>
#include<unordered_map>
#include<array>
#include<map>
#include<iomanip>
//<< setprecision(5)

using namespace std;
typedef long long int llint;
typedef long double lldo;
#define mp make_pair
#define mt make_tuple
#define pub push_back
#define puf push_front
#define pob pop_back
#define pof pop_front
#define fir first
#define sec second
#define res resize
#define ins insert
#define era erase
#define dme cout<<-1<<endl;return 0
const int mod=1000000007;
const int big=1e9;
const long double pai=3.141592653589793238462643383279502884197;
const long double eps=1e-7;

template <class T,class U>void mineq(T& a,U b){if(a>b){a=b;}}
template <class T,class U>void maxeq(T& a,U b){if(a<b){a=b;}}
llint gcd(llint a,llint b){if(a%b==0){return b;}else return gcd(b,a%b);}
llint lcm(llint a,llint b){return a/gcd(a,b)*b;}
int main(void){
	array<array<array<llint,12>,12>,62> rui;//★　行　列　累　乗　★
	array<array<llint,12>,12> m;
	llint n;cin>>n;
	llint h,i,j,k;
	//	0	1	2	3	4	5	6	7	8	9	10	11
	//	->				->	->	->		>	>	>	
	//		->			<	<		->	<	<		>
	//			->		>		<	<	->		<	<
	//				->		>	>	>		->	->	->
	//遷移を手書きします(きつい)
	for(i=0;i<12;i++){for(j=0;j<12;j++){m[i][j]=0;}}
	for(i=0;i<12;i++){m[i][i]=1;}
	for(i=0;i<4;i++){for(j=0;j<4;j++){m[i][j]=1;}}
	for(i=4;i<7;i++){m[i][i+1]=1;m[i+1][i]=1;}
	for(i=8;i<11;i++){m[i][i+1]=1;m[i+1][i]=1;}
	for(i=4;i<8;i++){m[0][i]=1;m[i][3]=1;}
	m[1][6]=1;m[1][7]=1;m[4][2]=1;m[5][2]=1;
	for(i=8;i<12;i++){m[3][i]=1;m[i][0]=1;}
	m[2][8]=1;m[2][9]=1;m[10][1]=1;m[11][1]=1;
	rui[0]=m;
	
	for(h=1;h<62;h++){
		for(i=0;i<12;i++){for(j=0;j<12;j++){rui[h][i][j]=0;}}
		for(i=0;i<12;i++){
			for(j=0;j<12;j++){
				for(k=0;k<12;k++){
					rui[h][i][k]+=rui[h-1][i][j]*rui[h-1][j][k];
					rui[h][i][k]%=mod;
				}
			}
		}
	}
	array<llint,12> gen;//答え
	for(i=0;i<8;i++){gen[i]=1;}
	for(i=8;i<12;i++){gen[i]=0;}
	for(h=0;h<62;h++){
		if((n& ((llint)1L<<(llint)h) )!=0){
			array<llint,12> iti;
			for(i=0;i<12;i++){iti[i]=0;}
			for(i=0;i<12;i++){
				for(j=0;j<12;j++){
					iti[j]+=gen[i]*rui[h][i][j];
					iti[j]%=mod;
				}
			}
			swap(iti,gen);
		}
	}
	cout<<gen[3]<<endl;
	return 0;
}
0