結果
問題 | No.840 ほむほむほむら |
ユーザー |
![]() |
提出日時 | 2019-06-14 21:37:34 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 147 ms / 4,000 ms |
コード長 | 1,414 bytes |
コンパイル時間 | 1,611 ms |
コンパイル使用メモリ | 166,284 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-09 00:35:48 |
合計ジャッジ時間 | 3,639 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 25 |
ソースコード
#include <bits/stdc++.h>#define syosu(x) fixed<<setprecision(x)using namespace std;typedef long long ll;typedef unsigned long long ull;typedef pair<int,int> P;typedef pair<double,double> pdd;typedef pair<ll,ll> pll;typedef vector<int> vi;typedef vector<vi> vvi;typedef vector<double> vd;typedef vector<vd> vvd;typedef vector<ll> vl;typedef vector<vl> vvl;typedef vector<string> vs;typedef vector<P> vp;typedef vector<vp> vvp;typedef vector<pll> vpll;typedef pair<int,P> pip;typedef vector<pip> vip;const int inf=1<<30;const ll INF=1ll<<60;const double pi=acos(-1);const double eps=1e-8;const ll mod=998244353 ;const int dx[4]={-1,0,1,0},dy[4]={0,-1,0,1};vvl mul(vvl A,vvl B){vvl C(A.size(),vl(B[0].size()));for(int i=0;i<A.size();i++) for(int k=0;k<B.size();k++) for(int j=0;j<B[0].size();j++) (C[i][j]+=A[i][k]*B[k][j])%=mod;return C;}vvl pow(vvl A,ll n){vvl B(A.size(),vl(A.size()));for(int i=0;i<A.size();i++) B[i][i]=1;while(n>0){if(n&1) B=mul(B,A);A=mul(A,A);n>>=1;}return B;}int n,k;int main(){cin>>n>>k;int N=k*k*k;vvl a(N,vl(N));for(int i=0;i<N;i++) for(int j=0;j<3;j++){int I=i%k,J=i/k%k,K=i/(k*k);if(j==0) (I+=1)%=k;if(j==1) (J+=I)%=k;if(j==2) (K+=J)%=k;a[i][I+J*k+K*k*k]++;}a=pow(a,n);vvl b(N,vl(1));b[0][0]++;a=mul(a,b);ll res=0;for(int i=0;i<N;i++) if(i/(k*k)==0) (res+=a[i][0])%=mod;cout<<res<<endl;}