結果
問題 | No.673 カブトムシ |
ユーザー |
|
提出日時 | 2018-04-13 22:50:50 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 4 ms / 2,000 ms |
コード長 | 1,039 bytes |
コンパイル時間 | 717 ms |
コンパイル使用メモリ | 75,352 KB |
最終ジャッジ日時 | 2025-01-05 10:05:44 |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 14 |
ソースコード
#include <iostream> #include <vector> using namespace std; using ll = long long; constexpr ll MOD = 1000000007; using Mat = vector<vector<ll>>; Mat mul(const Mat& m1, const Mat& m2) { Mat ans(2, vector<ll>(2, 0LL)); for (int i = 0; i < 2; i++) { for (int j = 0; j < 2; j++) { for (int k = 0; k < 2; k++) { (ans[i][j] += m1[i][k] * m2[k][j] % MOD) %= MOD; } } } return ans; } Mat unit() { Mat ans(2, vector<ll>(2, 0)); ans[0][0] = 1; ans[1][1] = 1; return ans; } Mat power(const Mat& m, const ll n) { if (n == 0) { return unit(); } if (n % 2 == 1) { return mul(power(m, n - 1), m); } else { const Mat pp = power(m, n / 2); return mul(pp, pp); } } int main() { ll B, C, D; cin >> B >> C >> D; B %= MOD, C %= MOD; Mat mat(2, vector<ll>(2, 0)); mat[0][0] = C; mat[0][1] = B * C % MOD; mat[1][1] = 1; const Mat ans = power(mat, D); cout << ans[0][1] << endl; return 0; }