题意:
给出一个长度为n的序列A(A1,A2…AN)。如果序列A不是非降的,你必须从中删去一个数,
这一操作,直到A非降为止。求有多少种不同的操作方案,答案模10^9+7。
题解:
还挺好做的吧。
首先不考虑强制停止的情况,
f[i][j]
表示前i个数,选出j个不降的数的方案数(i必须选)。
那么用个树状数组就可以
n2log n
转移了。
那么对于选出k个不降的数的方案数,就是
f′[k]=∑nif[i][k]
那么删数的方案数就是
F[i]=f′[i]∗(n−i)!
然后用容斥去掉不合法的情况。
显然就是在剩i+1或之前的时候就停下来的情况,
F′[i]=F[i]−F[i+1]∗(i+1)
然后
ans=∑niF′[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);
}