codeforces1153F Serval and Bonus Problem

题面

题意

在一条线段上随机放n个线段,则被至少k条线段覆盖的区间总长度期望为多少.

做法

首先可以发现这n个线段的$2\ast n$个端点期望把线段等分为$2\ast n+1$段,那么只要计算每一段的贡献即可.这个可以用dp求出每段区间被大于等于k条线段覆盖的方案数,再除以总方案数,即可得到每段区间对答案的贡献,然后就能统计答案了.

代码

#include<bits/stdc++.h>
#define ll long long
#define N 2010
#define MN 4000
#define M 998244353
using namespace std;

ll n,m,L,tot,ans,dp[N<<1][N],jc[N<<1],nj[N<<1];

inline void Add(ll &u,ll v){u=(u+v)%M;}
inline ll po(ll u,ll v)
{
	ll res=1;
	for(;v;)
	{
		if(v&1) res=res*u%M;
		u=u*u%M;
		v>>=1;
	}
	return res;
}

int main()
{
	ll i,j,len;
	jc[0]=1;for(i=1;i<=MN;i++) jc[i]=jc[i-1]*i%M;
	nj[MN]=po(jc[MN],M-2);for(i=MN-1;i>=0;i--) nj[i]=nj[i+1]*(i+1)%M;
	cin>>n>>m>>L;
	len=po(n*2+1,M-2);
	dp[1][1]=1;
	for(i=1;i<n*2;i++)
	{
		for(j=0;j<=min(i,n);j++)
		{
			if(!dp[i][j]) continue;
			if(j+1<=n*2-i-1) Add(dp[i+1][j+1],dp[i][j]);
			if(j) Add(dp[i+1][j-1],dp[i][j]*j%M);
		}
	}
	tot=dp[n*2][0];
	for(i=1;i<n*2;i++)
	{
		for(j=m;j<=n;j++)
		{
			Add(ans,dp[i][j]*dp[n*2-i][j]%M*jc[j]%M);
		}
	}
	ans=ans*po(tot,M-2)%M*len%M*L%M;
	cout<<ans;
}