結果

問題 No.2951 Similar to Mex
ユーザー nouka28nouka28
提出日時 2024-09-14 11:19:40
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
RE  
実行時間 -
コード長 1,600 bytes
コンパイル時間 5,270 ms
コンパイル使用メモリ 310,092 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-10-25 20:50:53
合計ジャッジ時間 10,365 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
6,816 KB
testcase_01 AC 4 ms
6,820 KB
testcase_02 AC 3 ms
6,820 KB
testcase_03 AC 3 ms
6,816 KB
testcase_04 AC 4 ms
6,820 KB
testcase_05 AC 3 ms
6,824 KB
testcase_06 AC 4 ms
6,820 KB
testcase_07 AC 4 ms
6,816 KB
testcase_08 AC 3 ms
6,820 KB
testcase_09 AC 3 ms
6,816 KB
testcase_10 AC 3 ms
6,820 KB
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
testcase_29 RE -
testcase_30 RE -
testcase_31 RE -
testcase_32 RE -
testcase_33 RE -
testcase_34 RE -
testcase_35 RE -
testcase_36 AC 3 ms
6,820 KB
testcase_37 AC 4 ms
6,820 KB
testcase_38 AC 4 ms
6,816 KB
testcase_39 AC 7 ms
6,820 KB
testcase_40 RE -
testcase_41 RE -
testcase_42 RE -
testcase_43 RE -
testcase_44 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;

#include<atcoder/all>
using namespace atcoder;
using mint=atcoder::modint998244353;

template<long long mod,long long MAX_N>
struct factional_prime{
	long long inv_[MAX_N+1];
	long long fac_[MAX_N+1];
	long long fac_inv_[MAX_N+1];

	factional_prime(){
		inv_[0]=0;inv_[1]=fac_[0]=fac_[1]=fac_inv_[0]=fac_inv_[1]=1;
		for(long long i=2;i<=MAX_N;i++){
			inv_[i]=((mod-mod/i)*inv_[mod%i])%mod;
			fac_[i]=(fac_[i-1]*i)%mod;
			fac_inv_[i]=(fac_inv_[i-1]*inv_[i])%mod;
		}
	}
	long long inv(long long n){
		if(n<0)return 0;
		return inv_[n];
	}
	long long fac(long long n){
		if(n<0)return 0;
		return fac_[n];
	}
	long long finv(long long n){
		if(n<0)return 0;
		return fac_inv_[n];
	}
	long long nCr(long long n,long long r){
		if(n<r||n<0||r<0)return 0;
		return ((fac_[n]*fac_inv_[n-r])%mod*fac_inv_[r])%mod;
	}
	long long nPr(long long n,long long r){
		if(n<r||n<0||r<0)return 0;
		return (fac_[n]*fac_inv_[n-r])%mod;
	}
};

factional_prime<998244353,1000> fp;

mint dp[89][89][89];

int main(){
	int N,M,K;cin>>N>>M>>K;

	dp[0][0][0]=1;

	int T=max(M,K)+2;

	for(int i=1;i<=T;i++){
		for(int j=0;j<=min(i-1,N);j++){
			for(int k=0;k<=i-1;k++){
				//0
				dp[i][j][i]+=dp[i-1][j][k]*mint(i).pow(max(0,min(K,i)-k));
				
				if(i<=M&&j<N){
					//1
					dp[i][j+1][k]+=dp[i-1][j][k];
				}
			}
		}
	}

	mint ans=0;

	for(int i=1;i<=N;i++){
		mint sum=0;
		for(int j=0;j<=i;j++){
			sum+=mint(j).pow(N)*mint(-1).pow(i-j)*fp.nCr(i,j);
		}
		mint sum2=0;
		for(int j=0;j<=T;j++){
			sum2+=dp[T][i][j];
		}
		ans+=sum*sum2;
	}

	cout<<ans.val()<<endl;
}
0