結果
問題 | No.344 ある無理数の累乗 |
ユーザー |
![]() |
提出日時 | 2016-02-12 23:33:06 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,229 bytes |
コンパイル時間 | 1,340 ms |
コンパイル使用メモリ | 165,876 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-09-22 05:00:19 |
合計ジャッジ時間 | 2,284 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 30 |
ソースコード
#include <bits/stdc++.h> using namespace std; typedef long long ll; typedef vector<int> vi; typedef vector<ll> vl; typedef complex<double> P; typedef pair<int,int> pii; #define REP(i,n) for(ll i=0;i<n;++i) #define REPR(i,n) for(ll i=1;i<n;++i) #define FOR(i,a,b) for(ll i=a;i<b;++i) #define DEBUG(x) cout<<#x<<": "<<x<<endl #define DEBUG_VEC(v) cout<<#v<<":";REP(i,v.size())cout<<" "<<v[i];cout<<endl #define ALL(a) (a).begin(),(a).end() #define MOD (ll)(1e9+7) #define ADD(a,b) a=((a)+(b))%MOD #define FIX(a) ((a)%MOD+MOD)%MOD // https://oeis.org/A080041 typedef vector<vl> mat; mat mul(const mat &a, const mat &b){ mat ret(2,vl(2)); REP(k,2)REP(i,2)REP(j,2)ret[i][j]+=a[i][k]*b[k][j]; REP(i,2)REP(j,2)ret[i][j] %= 1000; return ret; } int main(){ ll n; cin>>n; if(n==0){ cout<<1<<endl; }else if(n==1){ cout<<2<<endl; }else if(n==2){ cout<<7<<endl; }else{ mat I(2,vl(2)); I[0][0] = I[1][1] = 1; mat X(2,vl(2)); X[0][0] = X[0][1] = 2; X[1][0] = 1; mat ret = I; ll t = n-1; while(t){ if(t&1) ret=mul(ret,X); X=mul(X,X); t>>=1; } ll a_n = (ret[0][0]*2 + ret[0][1]*2)%1000; if(n%2==0) --a_n; cout<<a_n<<endl; } return 0; }