結果
問題 | No.1050 Zero (Maximum) |
ユーザー | kzyKT |
提出日時 | 2020-05-08 22:08:14 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 12 ms / 2,000 ms |
コード長 | 1,581 bytes |
コンパイル時間 | 1,551 ms |
コンパイル使用メモリ | 171,700 KB |
実行使用メモリ | 6,948 KB |
最終ジャッジ日時 | 2024-07-04 00:46:11 |
合計ジャッジ時間 | 2,180 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 15 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define F first #define S second #define R cin>> #define Z class #define ll long long #define ln cout<<'\n' #define in(a) insert(a) #define pb(a) push_back(a) #define pd(a) printf("%.10f\n",a) #define mem(a) memset(a,0,sizeof(a)) #define all(c) (c).begin(),(c).end() #define iter(c) __typeof((c).begin()) #define rrep(i,n) for(ll i=(ll)(n)-1;i>=0;i--) #define REP(i,m,n) for(ll i=(ll)(m);i<(ll)(n);i++) #define rep(i,n) REP(i,0,n) #define tr(it,c) for(iter(c) it=(c).begin();it!=(c).end();it++) void pr(){ln;} template<Z A,Z...B>void pr(const A &a,const B&...b){cout<<a<<' ';pr(b...);} template<Z A>void PR(A a,ll n){rep(i,n){if(i)cout<<' ';cout<<a[i];}ln;} ll check(ll n,ll m,ll x,ll y){return x>=0&&x<n&&y>=0&&y<m;} const ll MAX=1e9+7,MAXL=1LL<<61,dx[8]={-1,0,1,0,-1,-1,1,1},dy[8]={0,1,0,-1,-1,1,1,-1}; typedef pair<ll,ll> P; typedef vector<ll> Vec; typedef vector<Vec> Mat; ll const MOD = 1e9+7; Mat operator * (Mat const& A, Mat const& B) { Mat ret(A.size(), Vec(B[0].size())); rep(i,A.size())rep(k,B.size())rep(j,B[0].size())ret[i][j]=(ret[i][j]+A[i][k]*B[k][j])%MOD; return ret; } Mat pow(Mat A, ll n) { Mat ret(A.size(),Vec(A.size())); rep(i,A.size())ret[i][i]=1; for(;n>0;n>>=1) { if(n&1)ret=ret*A; A=A*A; } return ret; } void Main() { ll n,k; cin >> n >> k; Mat a;a.resize(n,Vec(n,0)); rep(i,n) { rep(j,n) { ll x=(i+j)%n; a[i][x]++; ll y=(i*j)%n; a[i][y]++; } } a=pow(a,k); pr(a[0][0]); } int main(){ios::sync_with_stdio(0);cin.tie(0);Main();return 0;}