题意:
给两个长度为n,所有数互不相同的序列,然后两两配对
题解:
好题啊,几乎没什么思路。
首先假如n,k奇偶性不一样显然无解,要求有(n+k)/2对每组要求。
看到这种计数题考虑容斥,对两个序列排序后,
f[i][j]
表示前i个a列的配对,有至少j个符合要求。
预处理个
next[i]
,表示第i个数最多比b列第几个数大。
所以
f[i][j]=f[i−1][j]+f[i−1][j−1]∗(next[i]−j+1)
然后就要减去大于k对的情况,即求
f′[n][k]
于是可以推出
f′[n][i]=f[n][i]∗(n−i)!−∑nj=i+1f′[n][j]∗Cij
这个挺好理解把,就是其他的随便匹配,然后减去>i的情况。
code:
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define LL long long
using namespace std;
const LL mod=1000000009;
LL n,k,a[2010],b[2010],next[2010];
LL g[2010],C[2010][2010];
LL f[2010][2010];
void pre()
{
g[0]=1;
for(LL i=1;i<=2000;i++) g[i]=g[i-1]*i%mod;
for(LL i=0;i<=2000;i++)
{
C[i][0]=1;
for(LL j=1;j<=i;j++) C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
}
}
int main()
{
scanf("%lld %lld",&n,&k);
if((n+k)%2!=0) {printf("0");return 0;}
k=(n+k)>>1;
for(LL i=1;i<=n;i++) scanf("%lld",&a[i]);
for(LL i=1;i<=n;i++) scanf("%lld",&b[i]);
sort(a+1,a+n+1);sort(b+1,b+n+1);
for(LL i=1;i<=n;i++)
{
LL j=0;
while(j<n&&b[j+1]<a[i]) j++;
next[i]=j;
}
pre();
memset(f,0,sizeof(f));
f[0][0]=1;
for(LL i=1;i<=n;i++)
for(LL j=0;j<=i;j++)
{
f[i][j]=f[i-1][j];
if(j>=1&&next[i]>=j) f[i][j]=(f[i][j]+f[i-1][j-1]*(next[i]-j+1)%mod)%mod;
}
for(LL i=n;i>=k;i--)
{
f[n][i]=f[n][i]*g[n-i]%mod;
for(LL j=i+1;j<=n;j++) f[n][i]=(f[n][i]-f[n][j]*C[j][i]%mod)%mod;
f[n][i]=(f[n][i]+mod)%mod;
}
printf("%lld",f[n][k]);
}