bzoj 3622 已经没有什么好害怕的了

我好害怕这种题 φ(≧ω≦*)?

两个数组排序后，处理出next[i]表示满足tang[i]>yao[j]的最大的j。

f[i][j]前i种糖果，有j个糖比药多，不考虑剩余情况的方案数//也就是剩余n-j个糖果的放法算一种，最后乘上阶乘。

f[i][j]=f[i-1][j]+f[i-1][j-1]*max(next[i]-j+1,0) 开始忘记max了。因为f[i][j]不一定都是合法状态，所以要判一下。

套容斥：

ans[j]=f[n][j]*jie[n-j]-sigma(ans[k]*c[k][j])

总结一下，容斥原理对于等于k个的可以使用上面的公式，而不是+1-1.spring我算明白了?

吐槽一下，这题样例也太强了吧，随便怎么写都能过。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>

#define md 1000000009
#define ll long long
#define inf (int) 1e9
#define eps 1e-8
#define N 2010
using namespace std;
int tang[N],yao[N],next[N];
ll f[N][N],c[N][N],jie[N];
int main()
{
int n,k; 
scanf("%d%d",&n,&k);
if ((n+k)&1) { printf("0\n"); return 0; }
k=(n+k)>>1;
for (int i=1;i<=n;i++) scanf("%d",&tang[i]);
for (int i=1;i<=n;i++) scanf("%d",&yao[i]);
sort(tang+1,tang+n+1); sort(yao+1,yao+n+1);
int w=1;
for (int i=1;i<=n;i++)
{
while (tang[i]>yao[w]&&w<=n) w++;
next[i]=w-1;
}
//for (int i=1;i<=n;i++) printf("%d ",next[i]); printf("\n%d\n",k);
f[0][0]=1;
for (int i=1;i<=n;i++)
  for (int j=0;j<=n;j++)
    f[i][j]=(f[i-1][j]+f[i-1][j-1]*max(next[i]-j+1,0))%md;
c[0][0]=1;
for (int i=1;i<=n;i++)
  for (int j=0;j<=i;j++)
    c[i][j]=(c[i-1][j]+c[i-1][j-1])%md;
jie[0]=1;
for (int i=1;i<=n;i++) jie[i]=(jie[i-1]*i)%md;
for (int i=n-1;i>=k;i--)
{
f[n][i]=(f[n][i]*jie[n-i])%md;
ll sum=0;
for (int j=i+1;j<=n;j++) sum=(sum+f[n][j]*c[j][i])%md;
f[n][i]=(f[n][i]-sum+md)%md;
}
printf("%lld\n",f[n][k]);
return 0;
}