結果
問題 | No.891 隣接3項間の漸化式 |
ユーザー |
![]() |
提出日時 | 2020-02-21 22:33:43 |
言語 | C++17(gcc12) (gcc 12.3.0 + boost 1.87.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 2,207 bytes |
コンパイル時間 | 3,467 ms |
コンパイル使用メモリ | 138,944 KB |
最終ジャッジ日時 | 2025-01-09 01:38:21 |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 31 TLE * 1 MLE * 7 |
ソースコード
#include <iostream> #include <algorithm> #include <vector> #include <queue> #include <cmath> #include <functional> #include <iomanip> using namespace std; using ll = long long; using pii = pair<int, int>; using vi = vector<int>; using vvi = vector<vi>; // #define rep(i,n) for(int i=0;i<(n);++i) #define all(f,c,...) (([&](decltype((c)) cccc) { return (f)(begin(cccc), end(cccc), ## __VA_ARGS__); })(c)) #define _overload3(_1,_2,_3,name,...) name #define _rep(i,n) for(int i=0;i<(n);++i) #define repi(i,a,b) for(int i=(a);i<(b);++i) #define rep(...) _overload3(__VA_ARGS__,repi,_rep,)(__VA_ARGS__) #define _fromto2(i,n) for(int i=0;i<=(n);++i) #define _fromto3(i,a,b) for(int i=(a);i<=(b);++i) #define fromto(...) _overload3(__VA_ARGS__,_fromto3,_fromto2,)(__VA_ARGS__) #define rrep(i,n) for(int i=(n);i>=0;--i) const int MOD = 1e9+7; const int INF = 1e9; const int dx[4]={1,0,-1,0}, dy[4]={0,1,0,-1}; template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return 1;}return 0;} template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return 1;}return 0;} // debug template<class T>ostream& operator<<(ostream& os,const vector<T>& vec){os<<"{";for(int i=0;i<vec.size();++i)os<<(i?", ":"")<<vec[i];os<<"}";return os;} template<class T1,class T2>ostream& operator<<(ostream& os,const pair<T1,T2>& rhs){os<<"("<<rhs.first<<", "<<rhs.second<<")";return os;} #ifdef LOCAL void debug() {cerr << "\n";} template<class First> void debug(const First& first) {cerr<<first<<"\n";} template<class First, class... Rest> void debug(const First& first, const Rest&... rest) {cerr<<first<<",";debug(rest...);} #else #define debug(...) 42 #endif struct Matrix { ll a, b, c, d; public: Matrix(ll a, ll b, ll c, ll d) : a(a % MOD), b(b % MOD), c(c % MOD), d(d % MOD) {} }; Matrix mul(Matrix A, Matrix B) { return Matrix(A.a * B.a + A.b * B.c, A.a * B.b + A.b * B.d, A.c * B.a + A.d * B.c, A.c * B.b + A.d * B.d); } Matrix pow(Matrix A, int p) { if (p == 1) return A; if (p % 2 == 0) return pow(mul(A, A), p / 2); else return mul(pow(mul(A, A), p / 2), A); } int main() { int a, b, n; cin >> a >> b >> n; Matrix m(0, 1, b, a); m = pow(m, n-1); cout << m.d << endl; return 0; }