結果
問題 | No.1050 Zero (Maximum) |
ユーザー |
![]() |
提出日時 | 2020-05-08 23:02:57 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 11 ms / 2,000 ms |
コード長 | 1,198 bytes |
コンパイル時間 | 1,928 ms |
コンパイル使用メモリ | 174,808 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-07-04 01:08:15 |
合計ジャッジ時間 | 2,258 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 15 |
ソースコード
#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){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]+A[i][k]*B[k][j])%MOD;}}return C;}mat pow(mat A, ll 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]*v[j]; res = res%MOD; }ret[i] = res%MOD;}return ret;}int main(){ll M, K;cin >> M >> K;vec v(M,0);v[0] = 1;mat m(M, vec(M, 0));rep(i,M){rep(j,M){m[(i+j)%M][i]++;m[(i*j)%M][i]++;}}mat a = pow(m,K);vec ret = mulv(a, v);cout << ret[0] << endl;return 0;}