bzoj 4361: isn

题意：

给出一个长度为n的序列A(A1,A2…AN)。如果序列A不是非降的，你必须从中删去一个数，
这一操作，直到A非降为止。求有多少种不同的操作方案，答案模10^9+7。

题解：

还挺好做的吧。
首先不考虑强制停止的情况，$f[i][j]$表示前i个数，选出j个不降的数的方案数（i必须选）。
那么用个树状数组就可以$n^2log\ n$转移了。
那么对于选出k个不降的数的方案数，就是$f'[k]=\sum_{i}^n f[i][k]$
那么删数的方案数就是$F[i]=f'[i]*(n-i)!$
然后用容斥去掉不合法的情况。
显然就是在剩i+1或之前的时候就停下来的情况，$F'[i]=F[i]-F[i+1]*(i+1)$
然后$ans=\sum_{i}^n F'[i]$
code：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define LL long long
const LL mod=1000000007;
using namespace std;
LL n,tr[2010];
LL f[2010][2010],g[2010],F[2010];
struct node{
    LL a,b,num;
}a[2010];
bool cmp(node a,node b) {return a.b<b.b;}
void change(LL k,LL c) {for(LL i=k;i<=n;i+=(i&-i)) (tr[i]+=c)%=mod;}
LL get(LL k) {LL ans=0;for(LL i=k;i>=1;i-=(i&-i)) (ans+=tr[i])%=mod;return ans;}
int main()
{
    scanf("%lld",&n);
    for(LL i=1;i<=n;i++) scanf("%lld",&a[i].b),a[i].num=i;
    g[0]=1;for(LL i=1;i<=n;i++) g[i]=g[i-1]*i%mod;
    sort(a+1,a+n+1,cmp);
    LL tot=0;
    for(LL i=1;i<=n;i++)
    {
        if(i==1||a[i].b!=a[i-1].b) tot++;
        a[a[i].num].a=tot;
    }
    for (int i=1;i<=n;++i)
    {
        memset(tr,0,sizeof(tr));
        if (i==1) change(1,1);
        for (int j=1;j<i;++j) change(a[j].a,f[j][i-1]);
        for (int j=i;j<=n;++j)
        {
            f[j][i]=get(a[j].a);
            change(a[j].a,f[j][i-1]);
        }
    }
    LL ans=0;
    for(LL i=n;i>=1;i--)
    {
        F[i]=0;
        for(LL j=1;j<=n;j++) (F[i]+=f[j][i])%=mod;
        F[i]=F[i]*g[n-i]%mod;ans=(ans+F[i]-F[i+1]*(i+1))%mod;
        (ans+=mod)%=mod;
    }
    printf("%lld",ans);
}