解题思路:
当
n+k
是奇数时肯定无解。
然后
k=(n+k)/2
,直接算糖果大于药片的对数刚好有
k
个。
恰好为
首先对两个队的数组分别从小到大排序, 记 cnti 表示 ai 比 B中多少元素大.
设
f(i,j)
表示考虑了
a1,...,ai
, 在其中选出
j
个, 且这
显然有
f(i,j)=f(i−1,j)+(cnti−j+1)f(i−1,j−1)
可能有的同学觉得
f(n,i)=f(n,i)∗(n−i)!
就是对数大于
i
的方案数,然后恰好
如:
糖果:
a1,a2,a3
药片:
b1,b2,b3
a1>b1,a2>b2,a3>b3
那么
((a1,b1),(a2,b2),a3不明)((a1,b1),(a3,b3),a2不明)
就会被视为两种答案,其实是同一种。
考虑
g(j)
是恰好有
j
对的方案,那么它在
g(i)=f(n,i)−∑j=i+1ng(j)
#include<bits/stdc++.h>
#define ll long long
using namespace std;
int getint()
{
int i=0,f=1;char c;
for(c=getchar();(c!='-')&&(c<'0'||c>'9');c=getchar());
if(c=='-')c=getchar(),f=-1;
for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
return i*f;
}
const int N=2005,p=1e9+9;
int n,k,a[N],b[N],cnt[N];
int f[N][N],c[N][N],fac[N];
int main()
{
//freopen("lx.in","r",stdin);
n=getint(),k=getint();
if ((n+k)&1) {cout<<0;return 0;}
k=(n+k)>>1;
for(int i=0;i<=n;i++)c[i][0]=1;
for(int i=1;i<=n;i++)
for(int j=1;j<=i;j++)
c[i][j]=(c[i-1][j]+c[i-1][j-1])%p;
fac[0]=1;
for(int i=1;i<=n;i++)fac[i]=(ll)fac[i-1]*i%p;
for(int i=1;i<=n;i++)a[i]=getint();
for(int i=1;i<=n;i++)b[i]=getint();
sort(a+1,a+n+1);
sort(b+1,b+n+1);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
if(a[i]>b[j])cnt[i]++;
else break;
f[0][0]=1;
for(int i=1;i<=n;i++)
for(int j=0;j<=i;j++)
{
f[i][j]=f[i-1][j];
if(j&&cnt[i]-j+1>0)f[i][j]=(f[i][j]+(ll)(cnt[i]-j+1)*f[i-1][j-1]%p)%p;
}
for(int i=n;i>=k;i--)f[n][i]=(ll)f[n][i]*fac[n-i]%p;
for(int i=n;i>=k;i--)
for(int j=i+1;j<=n;j++)f[n][i]=(f[n][i]-(ll)c[j][i]*f[n][j]%p+p)%p;
cout<<f[n][k];
return 0;
}