已经使 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;
}