結果
問題 | No.673 カブトムシ |
ユーザー | ojisan_IT |
提出日時 | 2018-04-13 23:25:54 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,889 bytes |
コンパイル時間 | 1,712 ms |
コンパイル使用メモリ | 174,360 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-06-30 05:54:19 |
合計ジャッジ時間 | 2,424 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
6,812 KB |
testcase_01 | AC | 2 ms
6,944 KB |
testcase_02 | AC | 2 ms
6,944 KB |
testcase_03 | AC | 2 ms
6,940 KB |
testcase_04 | AC | 2 ms
6,940 KB |
testcase_05 | AC | 2 ms
6,940 KB |
testcase_06 | AC | 2 ms
6,944 KB |
testcase_07 | AC | 2 ms
6,940 KB |
testcase_08 | AC | 2 ms
6,940 KB |
testcase_09 | AC | 2 ms
6,944 KB |
testcase_10 | AC | 2 ms
6,940 KB |
testcase_11 | AC | 2 ms
6,940 KB |
testcase_12 | AC | 2 ms
6,940 KB |
testcase_13 | AC | 2 ms
6,944 KB |
testcase_14 | AC | 2 ms
6,940 KB |
testcase_15 | AC | 2 ms
6,944 KB |
testcase_16 | AC | 2 ms
6,944 KB |
testcase_17 | AC | 2 ms
6,940 KB |
ソースコード
#include<bits/stdc++.h> using namespace std; using ll = long long; template<class T> using vt = vector<T>; template<class T> using vvt = vector<vt<T>>; template<class T> using ttt = tuple<T,T>; using tii = tuple<int,int>; using vi = vector<int>; #define rep(i,n) for(int i=0;i<(int)(n);i++) #define pb push_back #define mt make_tuple #define ALL(a) (a).begin(),(a).end() #define FST first #define SEC second #define DEB cerr<<"!"<<endl #define SHOW(a,b) cerr<<(a)<<" "<<(b)<<endl #define DIV ll(1e9+7) const int INF = (INT_MAX/2); const ll LLINF = (LLONG_MAX/2); const double eps = 1e-8; //const double PI = M_PI; inline ll pow(ll x,ll n,ll m){ll r=1;while(n>0){if((n&1)==1)r=r*x%m;x=x*x%m;n>>=1;}return r%m;} inline ll lcm(ll d1, ll d2){return d1 / __gcd(d1, d2) * d2;} // IT 5000兆 欲しい /* Coding space */ using Array = vector<ll>; using Matrix = vector<Array>; // O( n ) Matrix identity(ll n) { Matrix A(n, Array(n)); for (int i = 0; i < n; ++i) A[i][i] = 1; return A; } // O( n^2 ) Array mul(const Matrix &A, const Array &x) { Array y(A.size()); for (int i = 0; i < (int)A.size(); ++i) for (int j = 0; j < (int)A[0].size(); ++j) y[i] = (A[i][j] * x[j]) % DIV; return y; } // O( n^3 ) Matrix mul(const Matrix &A, const Matrix &B) { Matrix C(A.size(), Array(B[0].size())); for (int i = 0; i < (int)C.size(); ++i) for (int j = 0; j < (int)C[i].size(); ++j) for (int k = 0; k < (int)A[i].size(); ++k) C[i][j] += A[i][k] * B[k][j] % DIV, C[i][j] %= DIV; return C; } // O( n^3 log e ) Matrix pow(const Matrix &A, ll e) { return e == 0 ? identity(A.size()) : e % 2 == 0 ? pow(mul(A, A), e/2) : mul(A, pow(A, e-1)); } int main(){ ll b,c,d; cin >> b >> c >> d; b %= DIV; c%= DIV; Matrix A = {{c,b},{0,1}}; A = pow(A,d); Array B = {b,1}; Array Ans = mul(A,B); ll ans = Ans[0]; cout << (ans * c) % DIV << endl; }