結果
問題 | No.2529 Treasure Hunter |
ユーザー |
|
提出日時 | 2023-11-03 22:14:23 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 22 ms / 2,000 ms |
コード長 | 1,471 bytes |
コンパイル時間 | 794 ms |
コンパイル使用メモリ | 81,700 KB |
最終ジャッジ日時 | 2025-02-17 18:26:14 |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 22 |
ソースコード
#include<iostream> #include<vector> #include<algorithm> #include<atcoder/modint> using namespace std; using namespace atcoder; using ll=long long; using mint=modint998244353; #define all(v) v.begin(),v.end() #define rall(v) v.rbegin(),v.rend() template<class T> bool chmax(T &a, T b){if (a < b){a = b;return true;} else return false;} template<class T> bool chmin(T &a, T b){if (a > b){a = b;return true;} else return false;} const int MAX=2000010; const ll MOD=998244353; mint fac[MAX+1],finv[MAX+1],inv[MAX+1]; void set_fac(){ fac[0]=fac[1]=1; finv[0]=finv[1]=1; inv[1]=1; for (int i=2;i<=MAX;i++){ fac[i]=fac[i-1]*i; inv[i]=-inv[MOD%i]*(MOD/i); finv[i]=finv[i-1]*inv[i]; } } mint cmb(int n,int r){ if (r<0||n<r)return 0; return fac[n]*finv[r]*finv[n-r]; } int main(){ ios::sync_with_stdio(false); std::cin.tie(nullptr); int T; cin>>T; while(T--){ ll N,M; cin>>N>>M; vector dp(M+1,vector<mint>(3)); dp[0][0]=1; for(int i=0;i<M;i++){ dp[i+1][0]+=dp[i][0]+dp[i][1]+dp[i][2]; dp[i+1][1]+=dp[i][0]*N+dp[i][1]*(N-1)+dp[i][2]*(N-2); if(N>=4){ dp[i+1][2]+=dp[i][0]*(N*(N-1)/2-N); dp[i+1][2]+=dp[i][1]*(N*(N-1)/2-N-(N-3)); dp[i+1][2]+=dp[i][2]*(N*(N-1)/2-N-2*(N-3)+1); } } mint ans=dp[M][0]+dp[M][1]+dp[M][2]; cout<<ans.val()<<"\n"; } }