和Leo一起做爱数学的好孩子之容斥【BZOJ3622】已经没有什么好害怕的了

已经使 Modoka 有签订契约, 和自己一起战斗的想法后 , Mami 忽然感到自己不再是孤单一人了呢.
于是, 之前的谨慎的战斗作风也消失了 , 在对 Charlotte 的傀儡使用终曲——Tiro Finale后 , Mami 面临着即将被 Charlotte 的本体吃掉的局面.
这时, 已经多次面对过 Charlotte 的 Homura 告诉了学 OI 的你这样一个性质——Charlotte的结界中有两种具有能量的元素——一种是“糖果” , 另一种是“药片” , 每种各有 n 个. 在Charlotte 发动进攻前, “糖果”和“药片”会两两配对，若恰好“糖果”比“药片”能量大的组数比“药片”比“糖果”能量大的组数多 k 组, 则在这种局面下, Charlotte 的攻击会丢失,从而 Mami 仍有消灭 Charlotte 的可能. 你必须根据 Homura 告诉你的“糖果”和“药片”的能量的信息迅速告诉 Homura 这种情况的个数.

额实际上比集合选数难一点吧

你需要DP一下容斥系数乘上的全值

和集合选数一样思考

为至少有i个的情况数

这里需要一个DP

首先稍有常识的排一下序

定义

为i号糖果干涉成功的最右端点值

易知这个值是单调不减的

定义为前i个中选j个干涉的状态数

所以前面的所有的都可以选但是有些被选了

F[i][j]=(F[i][j]+F[i-1][j-1]*(Lim[i]-j+1)%mod)%mod;

但是这只是确定了干涉点

于是乘上一个非干涉点的全排列

容斥系数和集合选数一样为

#include<bits/stdc++.h>
using namespace std;
typedef int INT;
#define int long long
const int mod=1e9+9;
inline void read(int &x){
	x=0;
	int f=1;
	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();
	}
	x*=f;
}
int C[2001][2001];
int F[2001][2001];
int A[2001];
int B[2001];
int Lim[2001];
int fac[2001];
int n,k;
INT main(){
	read(n);
	read(k);
	k=(n+k)/2;
	for(int i=1;i<=n;++i)read(A[i]);
	for(int i=1;i<=n;++i)read(B[i]);
	fac[0]=1;
	for(int i=1;i<=2000;++i)fac[i]=fac[i-1]*i%mod;
	sort(A+1,A+1+n);
	sort(B+1,B+1+n);
	for(int i=1;i<=n;++i){
		Lim[i]=Lim[i-1];
		while(Lim[i]+1<=n&&B[Lim[i]+1]<A[i])Lim[i]++;
	}
	F[0][0]=1;
	for(int i=1;i<=n;++i){
		F[i][0]=1;
		for(int j=1;j<=i;++j){
			F[i][j]=F[i-1][j];
			if(j&&Lim[i]>j-1){
				F[i][j]=(F[i][j]+F[i-1][j-1]*(Lim[i]-j+1)%mod)%mod;
			}
		}
	}
	C[0][0]=1;
	for(int i=1;i<=2000;++i){
		C[i][0]=1;
		for(int j=1;j<=i;++j){
			C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
		}
	}
	int ans=0;
	for(int i=k;i<=n;++i)F[n][i]=(F[n][i]*fac[n-i])%mod;
	for(int i=k;i<=n;++i){
		int cur=(k-i)%2?-1:1;
		ans=((ans+(cur*F[n][i]*C[i][k])%mod)%mod+mod)%mod;
	}
	cout<<ans;
}

 

转载于:https://www.cnblogs.com/Leo-JAM/p/10079083.html