結果
問題 | No.891 隣接3項間の漸化式 |
ユーザー |
|
提出日時 | 2019-09-23 21:26:09 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,155 bytes |
コンパイル時間 | 1,722 ms |
コンパイル使用メモリ | 173,264 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-09-19 04:34:11 |
合計ジャッジ時間 | 2,701 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 39 |
ソースコード
#include <bits/stdc++.h> using namespace std; typedef long long ll; typedef pair<int,int> pii; typedef pair<ll,ll> pll; #define LLINF 9223372036854775807 #define MOD ll(1e9+7) #define all(x) (x).begin(),(x).end() #define dbg(x) cerr<<#x<<": "<<x<<endl typedef vector<vector<ll>> matrix; matrix mul(matrix &a, matrix &b){ matrix ret(a.size(), vector<ll>(b[0].size())); for(int i = 0; i < a.size(); i++){ for(int k = 0; k < b.size(); k++){ for(int j = 0; j < b[0].size(); j++){ ret[i][j] = (ret[i][j] + a[i][k]*b[k][j])%MOD; } } } return ret; } matrix matpow(matrix a, ll n){ matrix ret(a.size(), vector<ll>(a.size())); for(int i = 0; i < a.size(); i++){ ret[i][i] = 1; } while(n>0){ if(n&1){ ret = mul(ret, a); } a = mul(a, a); n = n >> 1; } return ret; } int main(){ ll a, b, n; cin >> a >> b >> n; matrix m(2,vector<ll>(2)); m[0][0] = a; m[0][1] = b; m[1][0] = 1; m[1][1] = 0; matrix mpow = matpow(m, n); ll ans = mpow[1][0]; cout << ans%MOD << endl; return 0; }