結果
問題 | 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";}}