[BZOJ4361]isn(DP+容斥+树状数组)

先用树状数组求出g[i]表示长度为i的不降子序列个数。$O(n^2\log n)$

容斥，最终序列长度为i的方案数为g[i]*(n-i)!，但这里多计算了在删得只剩i个数之前就已经构成不降子序列的情况。

这时考虑最后是从哪个不降子序列删掉一个数的，g[i+1]*(n-i-1)!*(i+1)，这里显然不会重复斥掉方案。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define rep(i,l,r) for (int i=(l); i<=(r); i++)
typedef long long ll;
using namespace std;

const int N=2010,mod=1e9+7;
int n,ans,a[N],b[N],fac[N],c[N][N],f[N][N],g[N];
void add(int c[],int x,int k){ for (; x<=n; x+=x&-x) c[x]=(c[x]+k)%mod; }
int que(int c[],int x){ int res=0; for (; x; x-=x&-x) res=(res+c[x])%mod; return res; }

int main(){
    freopen("bzoj4361.in","r",stdin);
    freopen("bzoj4361.out","w",stdout);
    scanf("%d",&n);
    rep(i,1,n) scanf("%d",&a[i]),b[i]=a[i];
    sort(b+1,b+n+1);
    rep(i,1,n) a[i]=lower_bound(b+1,b+n+1,a[i])-b;
    fac[1]=1; rep(i,2,n) fac[i]=1ll*fac[i-1]*i%mod;
    add(c[0],1,1);
    rep(i,1,n) for (int j=i; j; j--)
        f[i][j]=(f[i][j]+que(c[j-1],a[i]))%mod,add(c[j],a[i],f[i][j]);
    rep(i,1,n) rep(j,i,n) g[i]=(g[i]+f[j][i])%mod;
    rep(i,1,n) ans=(ans+((1ll*fac[n-i]*g[i]%mod-1ll*fac[n-i-1]*g[i+1]%mod*(i+1)%mod)%mod+mod)%mod)%mod;
    printf("%d\n",ans);
    return 0;
}