给出一个长度为n的序列A(A1,A2…AN)。如果序列A不是非降的,你必须从中删去一个数,
这一操作,直到A非降为止。求有多少种不同的操作方案,答案模10^9+7。
对于这个题目,我们用 f [ i , j ]表示选到 i ,i 必须选,已经删了j个,构成一个不降的子序列的方案数
于是由于每次都是从前for,然后统计个数,我们考虑用树状数组优化
因为取出来有一定的顺序,所以要乘以 ( n - i ) !
最后因为可能在去i+1的时候就已经不降了,所以我们需要容斥一下
#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
#define re register
#define gc getchar()
#define ll long long
inline int read()
{
re int x(0); re char ch(gc);
while(ch>'9'||ch<'0') ch=gc;
while(ch>='0'&&ch<='9') x=(x*10)+(ch^48),ch=gc;
return x;
}
const int N=2200,mod=1e9+7;
int a[N],n,b[N];
ll tr[N],f[N][N],sum[N],jc[N],ans;
#define lowbit(x) (x&(-x));
ll query(int x)
{
ll ans=0;
while(x)
{
ans+=tr[x];
ans%=mod;
x-=lowbit(x);
}
return ans;
}
void add(int x,int v)
{
if(!x) return;
while(x<=n)
{
tr[x]+=v;
tr[x]%=mod;
x+=lowbit(x);
}
}
int main()
{
n=read();
for(int i=1;i<=n;++i)
a[i]=b[i]=read();
sort(b+1,b+1+n);
int t=unique(b+1,b+1+n)-b-1;
for(int i=1;i<=n;++i)
a[i]=lower_bound(b+1,b+1+t,a[i])-b;
f[0][0]=1;
for(int i=1;i<=n;++i)
f[i][1]=1;
for(int i=2;i<=n;++i)
{
memset(tr,0,sizeof(tr));
for(int j=1;j<=n;++j)
{
add(a[j-1],f[j-1][i-1]);
f[j][i]+=query(a[j]);
f[j][i]%=mod;
}
}
for(int i=0;i<=n;++i)
for(int j=0;j<=n;++j)
sum[i]+=f[j][i],sum[i]%=mod;
jc[0]=1;
for(int i=1;i<=n;++i)
jc[i]=jc[i-1]*i%mod;
for(int i=1;i<=n;++i)
sum[i]=sum[i]*jc[n-i]%mod;
ans=0;
for(int i=1;i<=n;++i)
ans=sum[i]+ans-sum[i+1]*(i+1)%mod,ans%=mod;
cout<<(ans+mod)%mod<<endl;
return 0;
}