結果
問題 |
No.1521 Playing Musical Chairs Alone
|
ユーザー |
|
提出日時 | 2021-05-28 21:13:38 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 184 ms / 2,000 ms |
コード長 | 830 bytes |
コンパイル時間 | 2,258 ms |
コンパイル使用メモリ | 201,264 KB |
最終ジャッジ日時 | 2025-01-21 19:34:32 |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 23 |
ソースコード
#include<bits/stdc++.h> #define int long long using namespace std; using Mat=vector<vector<int>>; constexpr int mod=1000000007; Mat operator*(Mat a,Mat b){ Mat c(a.size(),vector<int>(b[0].size())); for(int i=0;i<a.size();i++)for(int j=0;j<b[0].size();j++){ for(int k=0;k<b.size();k++)(c[i][j]+=a[i][k]*b[k][j])%=mod; } return c; } Mat e(int N){ Mat a(N,vector<int>(N)); for(int i=0;i<N;i++)a[i][i]=1; return a; } Mat mat_pow(Mat a,int b){ if(!b)return e(a.size()); if(b&1)return mat_pow(a,b-1)*a; return mat_pow(a*a,b/2); } signed main(){ int N,K,L; cin>>N>>K>>L; Mat A(N,vector<int>(1)),B(N,vector<int>(N)); A[0][0]=1; for(int i=0;i<N;i++){ for(int j=1;j<=L;j++)B[(i+j)%N][i]++; } Mat C=mat_pow(B,K)*A; for(int i=0;i<N;i++)cout<<C[i][0]<<endl; }