結果

問題 No.2243 Coaching Schedule
ユーザー 沙耶花沙耶花
提出日時 2023-03-11 21:21:58
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 569 ms / 4,000 ms
コード長 1,070 bytes
コンパイル時間 4,783 ms
コンパイル使用メモリ 269,356 KB
実行使用メモリ 8,584 KB
最終ジャッジ日時 2024-09-18 06:48:36
合計ジャッジ時間 11,647 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 42 ms
5,248 KB
testcase_01 AC 43 ms
5,248 KB
testcase_02 AC 42 ms
5,376 KB
testcase_03 AC 43 ms
5,376 KB
testcase_04 AC 42 ms
5,376 KB
testcase_05 AC 126 ms
8,456 KB
testcase_06 AC 153 ms
8,456 KB
testcase_07 AC 152 ms
8,580 KB
testcase_08 AC 152 ms
8,456 KB
testcase_09 AC 173 ms
8,456 KB
testcase_10 AC 174 ms
8,452 KB
testcase_11 AC 176 ms
8,456 KB
testcase_12 AC 569 ms
8,452 KB
testcase_13 AC 455 ms
8,584 KB
testcase_14 AC 462 ms
8,456 KB
testcase_15 AC 458 ms
8,580 KB
testcase_16 AC 124 ms
8,584 KB
testcase_17 AC 74 ms
5,760 KB
testcase_18 AC 122 ms
8,120 KB
testcase_19 AC 159 ms
8,016 KB
testcase_20 AC 131 ms
6,592 KB
testcase_21 AC 90 ms
6,528 KB
testcase_22 AC 86 ms
6,272 KB
testcase_23 AC 49 ms
5,376 KB
testcase_24 AC 106 ms
6,528 KB
testcase_25 AC 52 ms
5,376 KB
testcase_26 AC 66 ms
5,760 KB
testcase_27 AC 113 ms
6,596 KB
testcase_28 AC 53 ms
5,376 KB
testcase_29 AC 133 ms
7,864 KB
testcase_30 AC 139 ms
8,160 KB
testcase_31 AC 146 ms
8,460 KB
testcase_32 AC 82 ms
6,136 KB
testcase_33 AC 67 ms
5,888 KB
testcase_34 AC 122 ms
7,800 KB
testcase_35 AC 138 ms
7,952 KB
testcase_36 AC 135 ms
8,312 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>
#include <atcoder/all>
#include <bits/stdc++.h>
using namespace std;
using namespace atcoder;
using mint = modint998244353;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf32 1000000001
#define Inf64 4000000000000000001

mint iv[300000];
int main(){
	for(int i=1;i<300000;i++)iv[i] = mint(i).inv();
	
	int M,N;
	cin>>M>>N;
	
	vector<int> c(M);
	rep(i,N){
		int a;
		cin>>a;
		c[a-1]++;
	}
	
	vector<int> d(N+1);
	rep(i,c.size()){
		d[c[i]]++;
	}
	
	vector<mint> num(N+1,1);
	for(int i=1;i<=N;i++){
		if(d[i]==0)continue;
		mint t = 1;
		rep(j,i)t *= j+1;
		for(int j=i;j<=N;j++){
			num[j] *= t.pow(d[i]);
			t *= j+1;
			t *= iv[j-i+1];
		}
		rep(j,i)num[j] = 0;
		
	}
	vector<mint> x(N+1,0);
	x[0] = 1;
	{
		mint c = 1;
		for(int i=1;i<=N;i++){
			c /= i;
			if(i%2==0)x[i] = c;
			else x[i] = -c;
		}
	}
	{
		mint c = 1;
		for(int i=1;i<=N;i++){
			c /= i;
			num[i] *= c;
		}
	}
	x = convolution(x,num);
	mint ans= 0;
	{
		mint c = 1;
		for(int i=1;i<=N;i++){
			c *= i;
			ans += x[i]*c;
		}
	}
	
	cout<<ans.val()<<endl;
	
	return 0;
}
0