結果

問題 No.2530 Yellow Cards
ユーザー okkuukenkenokkuukenken
提出日時 2023-11-03 23:04:37
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 379 ms / 2,000 ms
コード長 1,463 bytes
コンパイル時間 7,714 ms
コンパイル使用メモリ 349,924 KB
実行使用メモリ 237,784 KB
最終ジャッジ日時 2024-09-25 21:21:38
合計ジャッジ時間 11,601 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 41 ms
43,520 KB
testcase_03 AC 18 ms
19,968 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 43 ms
43,392 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 370 ms
236,672 KB
testcase_08 AC 379 ms
236,800 KB
testcase_09 AC 362 ms
236,800 KB
testcase_10 AC 365 ms
237,148 KB
testcase_11 AC 370 ms
237,784 KB
testcase_12 AC 362 ms
236,928 KB
testcase_13 AC 367 ms
236,288 KB
testcase_14 AC 331 ms
214,400 KB
testcase_15 AC 239 ms
161,408 KB
testcase_16 AC 191 ms
128,128 KB
testcase_17 AC 331 ms
215,424 KB
testcase_18 AC 78 ms
57,216 KB
testcase_19 AC 20 ms
17,024 KB
testcase_20 AC 93 ms
67,840 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
constexpr int dx[]={1,0,0,-1},dy[]={0,1,-1,0};
constexpr int Dx[]={1,1,0,-1,-1,-1,0,1},Dy[]={0,1,1,1,0,-1,-1,-1};
constexpr int mod=998244353,Mod=1e9+7,inf=Mod;
constexpr ll linf=(ll)inf*inf;
typedef pair<ll,int>P;
#define m_p make_pair
struct fastio{
	fastio(){
		cin.tie(0);
		ios::sync_with_stdio(0);
		cout<<fixed<<setprecision(20);
	};
}fio;
template<class T,class U>
bool chmax(T&a,const U&b){
	if(a<b){
		a=b;
		return 1;
	}
	return 0;
}
template<class T,class U>
bool chmin(T&a,const U&b){
	if(a>b){
		a=b;
		return 1;
	}
	return 0;
}
//*
#include<atcoder/all>
using namespace atcoder;
//*/
ll mypow(int x,ll n){
	if(n==0)
		return 1;
	if(n%2==1)
		return mypow(x,n-1)*x%mod;
	ll res=mypow(x,n/2);
	return res*res%mod;
}
ll inv(int x){
	return mypow(x,mod-2);
}
ll dp1[5001][5001],dp2[5001][5001];
int main(){
	int n,k;
	cin>>n>>k;
	int in=inv(n);
	dp1[0][0]=1;
	dp2[0][0]=n;
	for(int i=0;i<k;i++){
		for(int j=0;j<=min(n,i);j++){
			if(j>0)
				(dp1[i+1][j-1]+=dp1[i][j]*in%mod*j)%=mod;
			if(j<n)
				(dp1[i+1][j+1]+=dp1[i][j]*in%mod*(n-j))%=mod;
			if(j>0)
				(dp2[i+1][j-1]+=(dp2[i][j]+dp1[i][j])*in%mod*j)%=mod;
			if(j<n)
				(dp2[i+1][j+1]+=dp2[i][j]*in%mod*(n-j))%=mod;
		}
	}
	int ans=0;
	for(int i=0;i<=n;i++)
		(ans+=dp2[k][i])%=mod;
	cout<<ans<<endl;
}
0