传送门
思路:
直接dp求答案很麻烦,因为限制条件较严格
不妨放宽限制条件
1.首先想到
(n−k)mod2=0
,不然就是0
2.升序排两遍,然后对于a[i],我们记录b数组中所有比它小的元素的个数,记为pos[i]
3.f[i][j]指前i个元素中a数组有j个元素比b数组中元素大,其他i-j个元素位置不确定的方案数,显然我们有
f[i][j]=f[i−1][j]+f[i−1][j−1]×(pos[i]−j+1)
初始状态为
f[i][0]=1(i=1,2,...,n)
4.
g[i]
表示到第n个元素,糖果>药片的组数为i,且剩下的为糖果<药片的方案数,显然
g[i]
是可以通过f[i][n]和之前算出的g减出来的
g[i]=f[n][i]×(n−i)!−∑j=i+1ng[j]×Cij
#include<cstdio>
#include<algorithm>
#define N 2005
#define mo 1000000009
#define LL long long
using namespace std;
int n,k;
int fac[N],inv[N],a[N],b[N],pos[N],f[N][N],g[N];
int C(int m,int n){return (LL)fac[n]*inv[m]%mo*inv[n-m]%mo;}
main()
{
scanf("%d%d",&n,&k);
if (n-k&1) return puts("0"),0;
inv[0]=fac[0]=1;
for (int i=1;i<=n;++i) fac[i]=(LL)fac[i-1]*i%mo;
inv[1]=1;
for (int i=2;i<=n;++i) inv[i]=(LL)(mo-mo/i)*inv[mo%i]%mo;
for (int i=2;i<=n;++i) inv[i]=(LL)inv[i-1]*inv[i]%mo;
for (int i=1;i<=n;++i) scanf("%d",a+i);
for (int i=1;i<=n;++i) scanf("%d",b+i);
sort(a+1,a+n+1);
sort(b+1,b+n+1);
int now=1;
for (int i=1;i<=n;++i)
{
while (now<=n&&b[now]<a[i]) ++now;
--now;
pos[i]=now;
}
for (int i=0;i<=n;++i) f[i][0]=1;
for (int i=1;i<=n;++i)
for (int j=1;j<=i;++j)
f[i][j]=(f[i-1][j]+(LL)f[i-1][j-1]*(pos[i]-j+1)%mo)%mo;
for (int i=n;i>=(n+k)/2;--i)
{
g[i]=(LL)f[n][i]*fac[n-i]%mo;
for (int j=i+1;j<=n;++j)
g[i]=(g[i]-(LL)g[j]*C(i,j))%mo;
}
printf("%d\n",(g[n+k>>1]+mo)%mo);
}