結果
問題 | No.8036 Restricted Lucas (Easy) |
ユーザー | fura |
提出日時 | 2020-06-05 12:19:57 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 60 ms / 2,000 ms |
コード長 | 2,025 bytes |
コンパイル時間 | 929 ms |
コンパイル使用メモリ | 71,092 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-12-15 07:17:25 |
合計ジャッジ時間 | 1,501 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 3 ms
6,820 KB |
testcase_01 | AC | 5 ms
6,820 KB |
testcase_02 | AC | 59 ms
6,816 KB |
testcase_03 | AC | 60 ms
6,820 KB |
testcase_04 | AC | 58 ms
6,816 KB |
testcase_05 | AC | 58 ms
6,816 KB |
testcase_06 | AC | 58 ms
6,816 KB |
ソースコード
#include <iostream> #define rep(i,n) for(int i=ZERO;i<(n);i=add(i,ONE)) using namespace std; const int ZERO=false; const int ONE=true; const int TWO=ONE<<ONE; int MOD; int add(int a,int b){ if(b==ZERO){ return a; } else if(b>ZERO){ return add(a^b,(a&b)<<ONE); } else{ if(a>=ZERO) return add(b,a); else return add(~add(add(~a,ONE),add(~b,ONE)),ONE); } } int sub(int a,int b){ return add(a,add(~b,ONE)); } int modadd(int a,int b){ if(a<ZERO) a=add(a,MOD); if(a>=MOD) a=sub(a,MOD); if(b<ZERO) b=add(b,MOD); if(b>=MOD) b=sub(b,MOD); int res=add(a,b); if(res<ZERO) res=add(res,MOD); if(res>=MOD) res=sub(res,MOD); return res; } int modmul(int a,int b){ return b==ZERO?ZERO:modadd(modmul(a,b>>ONE)<<ONE,b&ONE?a:ZERO); } void modmul(int A[TWO][TWO],int B[TWO][TWO],int C[TWO][TWO]){ int tmp[TWO][TWO]={}; rep(i,TWO) rep(j,TWO) rep(k,TWO) { tmp[i][j]=modadd(tmp[i][j],modmul(A[i][k],B[k][j])); } rep(i,TWO) rep(j,TWO) C[i][j]=tmp[i][j]; } void solve(){ int A[TWO][TWO]={{ZERO,ONE},{ONE,ONE}}; int B[TWO][TWO]={{ONE,ZERO},{ZERO,ONE}}; int k; cin>>k; for(;k>ZERO;k>>=ONE){ if(k&ONE) modmul(B,A,B); modmul(A,A,A); } cout<<modadd(modmul(B[ZERO][ZERO],TWO),B[ZERO][ONE])<<endl; } int main(){ MOD=add(MOD<<ONE,ONE); MOD=add(MOD<<ONE,ONE); MOD=add(MOD<<ONE,ONE); MOD=add(MOD<<ONE,ZERO); MOD=add(MOD<<ONE,ONE); MOD=add(MOD<<ONE,ONE); MOD=add(MOD<<ONE,ONE); MOD=add(MOD<<ONE,ZERO); MOD=add(MOD<<ONE,ZERO); MOD=add(MOD<<ONE,ONE); MOD=add(MOD<<ONE,ONE); MOD=add(MOD<<ONE,ZERO); MOD=add(MOD<<ONE,ONE); MOD=add(MOD<<ONE,ZERO); MOD=add(MOD<<ONE,ONE); MOD=add(MOD<<ONE,ONE); MOD=add(MOD<<ONE,ZERO); MOD=add(MOD<<ONE,ZERO); MOD=add(MOD<<ONE,ONE); MOD=add(MOD<<ONE,ZERO); MOD=add(MOD<<ONE,ONE); MOD=add(MOD<<ONE,ZERO); MOD=add(MOD<<ONE,ZERO); MOD=add(MOD<<ONE,ZERO); MOD=add(MOD<<ONE,ZERO); MOD=add(MOD<<ONE,ZERO); MOD=add(MOD<<ONE,ZERO); MOD=add(MOD<<ONE,ONE); MOD=add(MOD<<ONE,ONE); MOD=add(MOD<<ONE,ONE); int q; cin>>q; rep(_,q) solve(); return ZERO; }