HDU-6348 序列计数(dp+树状数组)

题目:度度熊了解到,1,2,…,n 的排列一共有 n!=n×(n−1)×⋯×1 个。现在度度熊从所有排列中等概率随机选出一个排列 p1,p2,…,pn,你需要对 k=1,2,3,…,n 分别求出长度为 k 的上升子序列个数,也就是计算满足 1≤a1 < a2 < … < ak ≤n 且 pa1 <pa2< … < pak 的 k 元组 (a1,a2,…,ak) 的个数。    由于结果可能很大,同时也是为了 ruin the legend, 你只需要输出结果对 1000000007(=109+7) 取模后的值。

思路:因为是随机排列,循环长度用树状数组维护长度为len-1时的答案,滚动数组更新dp就行了。n*nlogn,随机序列能过的。。

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
const int maxn=1e4+10;
const ll mod=1e9+7;
int t,n,a[maxn];
ll dp[2][maxn],c[maxn],ans[maxn];
int lowbit(int x){return x&-x;}
void update(int pos,int val)
{
    while(pos<=n)
    {
        (c[pos]+=val)%=mod;
        pos+=lowbit(pos);
    }
}
ll query(int pos)
{
    ll ans=0;
    while(pos>0)
    {
        (ans+=c[pos])%=mod;
        pos-=lowbit(pos);
    }
    return ans;
}
int main()
{
    scanf("%d",&t);
    int cas=0;
    while(t--)
    {
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
            dp[1][i]=1;
        }
        memset(ans,0,sizeof ans);
        ans[1]=n;
        for(int len=2;len<=n;len++)
        {
            if(ans[len-1]==0) break;
            for(int i=1;i<=n;i++)
            {
                ll tmp=query(a[i]-1);
                (ans[len]+=tmp)%=mod;
                dp[len&1][i]=tmp;
                update(a[i],dp[(len-1)&1][i]);
            }
            memset(c,0,sizeof c);
        }
        printf("Case #%d:",++cas);
        for(int i=1;i<=n;i++)
            printf(" %lld",ans[i]);
        puts("");
    }
    return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值