結果
問題 | No.578 3 x N グリッド上のサイクルのサイズ(easy) |
ユーザー | 👑 kmjp |
提出日時 | 2017-10-15 01:55:41 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 4 ms / 2,000 ms |
コード長 | 1,789 bytes |
コンパイル時間 | 1,430 ms |
コンパイル使用メモリ | 148,312 KB |
実行使用メモリ | 4,384 KB |
最終ジャッジ日時 | 2023-08-11 02:41:34 |
合計ジャッジ時間 | 3,483 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge11 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
4,380 KB |
testcase_01 | AC | 2 ms
4,376 KB |
testcase_02 | AC | 1 ms
4,384 KB |
testcase_03 | AC | 1 ms
4,380 KB |
testcase_04 | AC | 1 ms
4,380 KB |
testcase_05 | AC | 2 ms
4,380 KB |
testcase_06 | AC | 2 ms
4,380 KB |
testcase_07 | AC | 2 ms
4,380 KB |
testcase_08 | AC | 2 ms
4,376 KB |
testcase_09 | AC | 2 ms
4,376 KB |
testcase_10 | AC | 2 ms
4,376 KB |
testcase_11 | AC | 2 ms
4,376 KB |
testcase_12 | AC | 2 ms
4,380 KB |
testcase_13 | AC | 2 ms
4,376 KB |
testcase_14 | AC | 2 ms
4,380 KB |
testcase_15 | AC | 2 ms
4,380 KB |
testcase_16 | AC | 2 ms
4,380 KB |
testcase_17 | AC | 2 ms
4,380 KB |
testcase_18 | AC | 2 ms
4,380 KB |
testcase_19 | AC | 2 ms
4,380 KB |
testcase_20 | AC | 2 ms
4,380 KB |
testcase_21 | AC | 3 ms
4,376 KB |
testcase_22 | AC | 2 ms
4,380 KB |
testcase_23 | AC | 2 ms
4,376 KB |
testcase_24 | AC | 3 ms
4,380 KB |
testcase_25 | AC | 2 ms
4,376 KB |
testcase_26 | AC | 2 ms
4,376 KB |
testcase_27 | AC | 2 ms
4,376 KB |
testcase_28 | AC | 3 ms
4,380 KB |
testcase_29 | AC | 3 ms
4,380 KB |
testcase_30 | AC | 2 ms
4,376 KB |
testcase_31 | AC | 2 ms
4,380 KB |
testcase_32 | AC | 2 ms
4,376 KB |
testcase_33 | AC | 3 ms
4,384 KB |
testcase_34 | AC | 3 ms
4,376 KB |
testcase_35 | AC | 3 ms
4,376 KB |
testcase_36 | AC | 3 ms
4,376 KB |
testcase_37 | AC | 3 ms
4,376 KB |
testcase_38 | AC | 3 ms
4,376 KB |
testcase_39 | AC | 4 ms
4,376 KB |
testcase_40 | AC | 3 ms
4,380 KB |
testcase_41 | AC | 3 ms
4,380 KB |
testcase_42 | AC | 3 ms
4,376 KB |
testcase_43 | AC | 3 ms
4,384 KB |
testcase_44 | AC | 3 ms
4,380 KB |
testcase_45 | AC | 3 ms
4,380 KB |
testcase_46 | AC | 4 ms
4,380 KB |
testcase_47 | AC | 3 ms
4,380 KB |
testcase_48 | AC | 4 ms
4,380 KB |
testcase_49 | AC | 3 ms
4,384 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; typedef signed long long ll; #undef _P #define _P(...) (void)printf(__VA_ARGS__) #define FOR(x,to) for(x=0;x<(to);x++) #define FORR(x,arr) for(auto& x:arr) #define ITR(x,c) for(__typeof(c.begin()) x=c.begin();x!=c.end();x++) #define ALL(a) (a.begin()),(a.end()) #define ZERO(a) memset(a,0,sizeof(a)) #define MINUS(a) memset(a,0xff,sizeof(a)) //------------------------------------------------------- ll mo=1000000007; const int MAT=10; struct Mat { ll v[MAT][MAT]; }; Mat mulmat(Mat& a,Mat& b,int n=MAT) { ll mo2=4*mo*mo; int x,y,z; Mat r; FOR(x,n) FOR(y,n) r.v[x][y]=0; FOR(x,n) FOR(z,n) FOR(y,n) { r.v[x][y] += a.v[x][z]*b.v[z][y]; if(r.v[x][y]>mo2) r.v[x][y] -= mo2; } FOR(x,n) FOR(y,n) r.v[x][y]%=mo; return r; } Mat powmat(ll p,Mat a,int n=MAT) { int i,x,y; Mat r; FOR(x,n) FOR(y,n) r.v[x][y]=0; FOR(i,n) r.v[i][i]=1; while(p) { if(p%2) r=mulmat(r,a,n); a=mulmat(a,a,n); p>>=1; } return r; } ll N; string m[10]={ "0.........", "1012...1..", "210112....", "3210212.1.", "1.1201....", "2.21101...", "1..2.101..", "...1...0..", "21....1.0.", ".1231212.0", }; void solve() { int i,j,k,l,r,x,y; string s; Mat A,AD; FOR(y,10) FOR(x,10) { if(m[y][x]=='.') { A.v[y][x]=AD.v[y][x]=0; } else { A.v[y][x]=1; AD.v[y][x]=m[y][x]-'0'; if(y>=1 && y<=6) AD.v[y][x]+=2; if(y>=7 && y<=8) AD.v[y][x]+=4; } } cin>>N; ll ret=0; FOR(i,N+1) { Mat L,R,C; L=powmat(i,A); R=powmat(N-i,A); C=mulmat(L,AD); C=mulmat(C,R); ret += C.v[9][0]; } cout<<ret%mo<<endl; } int main(int argc,char** argv){ string s;int i; if(argc==1) ios::sync_with_stdio(false), cin.tie(0); FOR(i,argc-1) s+=argv[i+1],s+='\n'; FOR(i,s.size()) ungetc(s[s.size()-1-i],stdin); solve(); return 0; }