提出 #69906067

---

ソースコード 拡げる

/*
*            /$$           /$$
*           |__/          |__/
*  /$$$$$$$$ /$$ /$$$$$$$$ /$$  /$$$$$$
* |____ /$$/| $$|____ /$$/| $$ /$$__  $$
*    /$$$$/ | $$   /$$$$/ | $$| $$  \ $$
*   /$$__/  | $$  /$$__/  | $$| $$  | $$
*  /$$$$$$$$| $$ /$$$$$$$$| $$|  $$$$$$$
* |________/|__/|________/|__/ \____  $$
*                                   | $$
*                                   | $$
*                                   |__/
*/
//hj23308保佑我
//Missile保佑我
/*
* 醒了在梦里挣扎，不觉黯淡了朝霞
*/
/*
* 我很高兴你没有忘了我，但是我现在更希望你已经忘了我了。
* 希望在你的记忆中，我只是尘土一撮，从你的全世界路过，然后四散飞扬不留下一点痕迹，而你要不回头的往前走。
* 我更希望我只是从你的全世界路过，只是路过
*/
/*
* 只是我在十字路口守了太久，守到黄沙如雨掩埋一切痕迹，才发现自己等的人已经离开了。
*/
/*
* 听我的 别回头 回头就可能会泪流满面，会被黄沙掩埋，所以即使痛苦也要向前走
*/
/*
* 我听到了「天行健」的回响，这是一个伟大斗士的不息自强；
* 我听到了「破万法」的回响，这是一个黑道打手的守护欲望；
* 我看见了「生生不息」的激荡，这是一个骗子的伟大乐章！
*/
/*
* 我用虚假的面具照顾着细腻的感情；
* 我以华丽的衣物下藏着腐烂的血肉；
* 当我摘下面具，褪去衣物，即便是我最亲近的人，也无法直视我
*/
#include<bits/stdc++.h>
using namespace std;
mt19937 engine(chrono::steady_clock().now().time_since_epoch().count());
const int MOD=998244353;
const int MAXN=105;
void AddC(int &x,int y)
{
	if((x+=y)>=MOD) x-=MOD;
}
int Mul(int x,int y)
{
	return 1ll*x*y%MOD;
}
int n,x;
long long res;
int pre,a[MAXN],sum[MAXN];
long long fac[MAXN],ifac[MAXN];
int f[MAXN][MAXN*MAXN],g[MAXN][MAXN*MAXN];
long long qpow(long long u,long long k)
{
	long long Ans=1;
	while(k) {
		if(k&1) Ans=Ans*u%MOD;
		k/=2;
		u=u*u%MOD;
	}
	return Ans;
}
long long inv(long long u)
{
	return qpow(u,MOD-2);
}
long long C(int N,int M)
{
	return fac[N]*ifac[M]%MOD*ifac[N-M]%MOD;
}
void Init()
{
	fac[0]=ifac[0]=1;
	for(int i=1;i<=n;i++) {
		fac[i]=fac[i-1]*i%MOD;
		ifac[i]=ifac[i-1]*inv(i)%MOD;
	}
}
void update(int pos,int tmp)
{
	memcpy(g,f,sizeof(g)),memset(f,0,sizeof(f));
	for(int i=0;i<=n;i++) {
		for(int j=0;j<=n*n;j++) {
			int val=g[i][j];
			if(!val) continue;
			int use=pos-1-(pre-tmp-i);
			for(int k=0;k<=min(tmp,use);k++) {
				int Val=Mul(val,Mul(Mul(C(use,k),C(tmp,k)),fac[k]));
				if(i) AddC(f[i+tmp-k-1][j+i],Mul(i,Val)); 
				if(tmp-k) AddC(f[i+tmp-k-1][j+i],Mul(tmp-k,Val));
				AddC(f[i+tmp-k][j+i],Val);
			}
		}
	}
//	for(int i=0;i<=n;i++) {
//		for(int j=0;j<=n*n;j++) {
//			if(f[i][j]) cout<<i<<" , "<<j<<" : "<<f[i][j]<<"\n";
//		}
//	}
//	cout<<"-------------\n";
}
int solve()
{
	f[0][0]=1;
	for(int i=1;i<=n;i++) {
		pre+=sum[i];
		update(i,sum[i]);
	}
	res-=1ll*n*(n-1)/2;
	int Ans=0;
	for(int j=0;j<=n*n;j++) {
		if(res+j<=x) AddC(Ans,f[0][j]);
	}
	return Ans;
}
int main()
{
//	freopen(".in","r",stdin);
//	freopen(".out","w",stdout);
	std::ios::sync_with_stdio(false);
	cin.tie(0),cout.tie(0);
	cin>>n>>x;
	for(int i=1;i<=n;i++) {
		cin>>a[i];
		res+=a[i];
		sum[min(a[i]+1,n)]++;
	}
	Init();
	cout<<solve()<<"\n";
	return 0;
}
/*
input:
3 1
0 1 2
output:
3
*/

提出情報

ジャッジ結果