[Codeforces1153F][数论]Serval and Bonus Problem

翻译

你有一个长度为$L$的数列，在上面随机选择$2\ast n$条线段
求至少被覆盖了$K$次的区间的长度和的期望

题解

数学菜逼来啦
先考虑一个子问题，在一个$[0,1]$的数轴上随机撒$i$个点，求某个间隔长度的期望
抛出结论的话就是期望是$\frac{1}{i+1}$
感性认知的话我们可以想象一下每个区间的期望其实是相同的…
理性认知的话我们考虑这$i$个点的最小值的位置，这个可以直接积分出来
所以我们现在仅需要考虑每个区间被选中的贡献
枚举左边有多少个端点，再枚举左边选了几个端点当左端点
预处理一个$msk[i]$表示$i$个端点互相匹配的方案数
然后随意组合数一下就能出答案了
其实对于这种数轴上随机撒点的题…一般如果能知道随机的期望就可以很好做了的…
重点是自己要能积分出来啊

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<ctime>
#include<map>
#include<bitset>
#include<set>
#define LL long long
#define mp(x,y) make_pair(x,y)
#define pll pair<long long,long long>
#define pii pair<int,int>
using namespace std;
inline int read()
{
	int f=1,x=0;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
int stack[20];
inline void write(int x)
{
	if(x<0){putchar('-');x=-x;}
    if(!x){putchar('0');return;}
    int top=0;
    while(x)stack[++top]=x%10,x/=10;
    while(top)putchar(stack[top--]+'0');
}
inline void pr1(int x){write(x);putchar(' ');}
inline void pr2(int x){write(x);putchar('\n');}
const int MAXN=4005;
const int mod=998244353;
int pow_mod(int a,int b)
{
	int ret=1;
	while(b)
	{
		if(b&1)ret=1LL*ret*a%mod;
		a=1LL*a*a%mod;b>>=1;
	}
	return ret;
}
int pre[MAXN],inv[MAXN],msk[MAXN];
int C(int n,int m){return 1LL*pre[n]*inv[m]%mod*inv[n-m]%mod;}
int n,K,L;
void ad(int &x,int y){x+=y;if(x>=mod)x-=mod;}
int main()
{
	msk[0]=1;for(int i=2;i<MAXN;i++)msk[i]=1LL*msk[i-2]*(i-1)%mod;
	pre[0]=1;for(int i=1;i<MAXN;i++)pre[i]=1LL*pre[i-1]*i%mod;
	inv[MAXN-1]=pow_mod(pre[MAXN-1],mod-2);
	for(int i=MAXN-2;i>=0;i--)inv[i]=1LL*inv[i+1]*(i+1)%mod;
	n=read();K=read();L=read();int sum=pow_mod(2*n+1,mod-2);
	int ans=0;
	for(int i=1;i<=2*n;i++)
	{
		if(i>=K&&2*n-i>=K)
		{
			for(int j=K;j<=min(i,2*n-i);j++)if(!((i-j)&1)&&!((2*n-i-j)&1))
				ad(ans,1LL*sum*C(i,j)%mod*C(2*n-i,j)%mod*pre[j]%mod*msk[i-j]%mod*msk[2*n-i-j]%mod);
		}
	}
	ans=1LL*ans*pow_mod(msk[2*n],mod-2)%mod*L%mod;
	pr2(ans);
	return 0;
}