bzoj 4361: isn

27 篇文章 0 订阅
21 篇文章 0 订阅

题意:

给出一个长度为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](ni)!
然后用容斥去掉不合法的情况。
显然就是在剩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);
}
  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值