結果
問題 | No.891 隣接3項間の漸化式 |
ユーザー |
|
提出日時 | 2019-09-20 21:40:46 |
言語 | C++11 (gcc 13.3.0) |
結果 |
MLE
|
実行時間 | - |
コード長 | 1,209 bytes |
コンパイル時間 | 1,353 ms |
コンパイル使用メモリ | 165,160 KB |
実行使用メモリ | 813,952 KB |
最終ジャッジ日時 | 2024-09-14 16:57:03 |
合計ジャッジ時間 | 3,830 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge6 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 10 MLE * 1 -- * 28 |
ソースコード
/*** author: yuji9511 ***/ #include <bits/stdc++.h> using namespace std; typedef long long ll; typedef pair<ll, ll> lpair; const ll MOD = 1e9 + 7; const ll INF = 1e18; #define rep(i,m,n) for(ll i = (m); i < (n); i++) #define rrep(i,m,n) for(ll i = (m); i >= (n); i--) #define print(x) cout << (x) << endl; #define print2(x,y) cout << (x) << " " << (y) << endl; #define printa(x,n) for(ll i = 0; i < n; i++){ cout << (x[i]) << " \n"[i==n-1];}; typedef vector<ll> vec; typedef vector<vec> mat; mat mult(mat &A, mat &B){ mat C(A.size(), vec(B[0].size(),0)); rep(i,0,A.size()){ rep(k,0,B.size()){ rep(j,0,B[0].size()){ C[i][j] += A[i][k] * B[k][j]; C[i][j] %= MOD; } } } return C; } mat pow_mat(mat A, ll n){ if(n == 1) return A; if(n % 2 == 0){ mat B = pow_mat(A, n/2); return mult(B,B); }else{ mat B = pow_mat(A, n-1); return mult(A,B); } } int main(){ cin.tie(0); ios::sync_with_stdio(false); ll a,b,n; cin >> a >> b >> n; mat A(2, vec(2,0)); A[0][0] = a; A[0][1] = b; A[1][0] = 1; A[1][1] = 0; mat ans = pow_mat(A, n-1); ll res = ans[0][0]; print(res); }