P4859 已经没有什么好害怕的了

传送门

见计数想容斥

首先题目可以简单转化一下， 求 糖果比药片能量大的组数比药片比糖果能量大的组数多 $k$ 组 的方案数

因为所有能量各不相同，所以就相当于求 糖果比药片能量大的组数为 $(n+k)/2$ 组的方案数，如果 $(n+k)$ 为奇数则无解

发现这个 '恰好' 很不好算，考虑先算出 '至少'，设 $F[i]$ 表示至少有 $i$ 对糖果比药片大的方案数

那么就是要强制选 $i$ 对糖果比药片大，然后再随便选，发现这个强制选 $i$ 对糖果比药片大的方案数也不好算..

考虑先把糖果和药片排序，然后 $dp$

设 $f[i][j]$ 表示从小到大前 $i$ 个糖果，和药片匹配了 $j$ 对的方案数

那么 $f[i][j]=f[i-1][j]+f[i-1][j-1]*(g[i]-j+1)$，其中 $g[i]$ 表示比糖果 $i$ 能量小的药片的数量

然后 $f[n][i]$ 就是强制选 $i$ 对糖果比药片大的方案数，因为剩下的随便选，所以剩下方案数就是 $(n-i)!$

所以 $F[i]=f[n][i]*(n-i)!$（注意这里的 $F[i]$ 中其实有些方案是重复算了）

设恰好 $i$ 对糖果比药片大的方案数为 $ans[i]$，可以（不能）发现，对于 $F[i]=f[n][i]*(n-i)!$

$F[i]=\sum_{j=i}^{n}C_{j}^{i}\cdot ans[j]$ （乘 $C_{j}^{i}$ 是因为有重复算）

所以逆推 $ans[i]$，$ans[i]=F[i]-\sum_{j=i+1}^{n}C_{j}^{i}ans[j]$

然后就可以算了

具体看代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
typedef long long ll;
inline int read()
{
    int x=,f=; char ch=getchar();
    while(ch<''||ch>'') { if(ch=='-') f=-; ch=getchar(); }
    while(ch>=''&&ch<='') { x=(x<<)+(x<<)+(ch^); ch=getchar(); }
    return x*f;
}
const int N=,mo=1e9+;
int n,K;
int A[N],B[N],fac[N],C[N][N];
int cnt[N],f[N][N],ans[N];//cnt[i]是比糖果i小的药片的数量
inline int fk(int x) { return x>=mo ? x-mo : x; }
int main()
{
    n=read(),K=read();
    for(int i=;i<=n;i++) A[i]=read();
    for(int i=;i<=n;i++) B[i]=read();
    sort(A+,A+n+); sort(B+,B+n+);
    if((n+K)&) { printf(""); return ; }//判断无解
    int p=;
    for(int i=;i<=n;i++)
    {
        cnt[i]=cnt[i-];
        while(p<=n&&A[i]>B[p]) cnt[i]++,p++;
    }
    for(int i=;i<=n;i++) f[i][]=;//初始化
    for(int i=;i<=n;i++)
        for(int j=;j<=n;j++) f[i][j]=fk(f[i-][j] + 1ll*(cnt[i]-j+)*f[i-][j-]%mo );//dp
    fac[]=; for(int i=;i<=n;i++) fac[i]=1ll*fac[i-]*i%mo;//求阶乘
    C[][]=;
    for(int i=;i<=n;i++)
        for(int j=;j<=n;j++)
        {
            C[i+][j]=fk(C[i+][j]+C[i][j]);
            C[i+][j+]=fk(C[i+][j+]+C[i][j]);//求组合数
        }
    for(int i=n;i>=(n+K)>>;i--)//逆推ans
    {
        ans[i]=1ll*f[n][i]*fac[n-i]%mo;
        for(int j=i+;j<=n;j++) ans[i]=fk(ans[i]-1ll*C[j][i]*ans[j]%mo+mo);
    }
    printf("%d",ans[(n+K)>>]);
    return ;
}