結果
問題 | No.673 カブトムシ |
ユーザー | どらら |
提出日時 | 2018-04-14 00:02:08 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,232 bytes |
コンパイル時間 | 1,711 ms |
コンパイル使用メモリ | 176,140 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-06-30 05:55:15 |
合計ジャッジ時間 | 2,388 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
6,816 KB |
testcase_01 | AC | 2 ms
6,940 KB |
testcase_02 | AC | 2 ms
6,940 KB |
testcase_03 | AC | 2 ms
6,940 KB |
testcase_04 | AC | 1 ms
6,940 KB |
testcase_05 | AC | 2 ms
6,944 KB |
testcase_06 | AC | 2 ms
6,940 KB |
testcase_07 | AC | 2 ms
6,940 KB |
testcase_08 | AC | 2 ms
6,940 KB |
testcase_09 | AC | 2 ms
6,940 KB |
testcase_10 | AC | 2 ms
6,940 KB |
testcase_11 | AC | 2 ms
6,944 KB |
testcase_12 | AC | 2 ms
6,944 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; #define REP(i,a,n) for(int i=(a); i<(int)(n); i++) #define rep(i,n) REP(i,0,n) #define FOR(it,c) for(__typeof((c).begin()) it=(c).begin(); it!=(c).end(); ++it) #define ALLOF(c) (c).begin(), (c).end() typedef long long ll; typedef unsigned long long ull; typedef vector<ll> vec; typedef vector<vec> mat; const ll MOD = 1000000007; mat mul(mat &A, mat &B){ // A*B mat C(A.size(), vec(B[0].size())); rep(i,A.size()) rep(k,B.size()){ rep(j,B[0].size()){ C[i][j] = ((C[i][j]%MOD)+(A[i][k]%MOD)*(B[k][j]%MOD))%MOD; } } return C; } mat pow(mat A, ll n){ // A^n mat B(A.size(), vec(A.size())); rep(i,A.size()) B[i][i] = 1; while(n>0){ if(n&1) B=mul(B,A); A=mul(A,A); n>>=1; } return B; } vector<ll> mulv(mat &A, vector<ll> v){ vector<ll> ret = v; rep(i,v.size()){ ll res = 0; rep(j,v.size()){ res += ((A[i][j]%MOD)*(v[j]%MOD))%MOD; res = res%MOD; } ret[i] = res%MOD; } return ret; } int main(){ ll B, C, D; cin >> B >> C >> D; B %= MOD; C %= MOD; vec a0 = {(B*C)%MOD,0,1}; mat A = {{C,0,(B*C)%MOD},{1,0,0},{0,0,1}}; mat nA = pow(A, D-1); vec ret = mulv(nA, a0); cout << ret[0] << endl; return 0; }