hdu 6348 序列计数(dp+bit)

像2015ccpc的一道题..赛时直接贴代码上去，果断mle到死，因为数据范围加到了1e4，所以没法开出这么大的二维数组。一个优化思路就是，当长度大于序列的最长上升子序列时，上升序列的数量一定是0，所以dp[i][j]，定义长度为前i个数里面长度为j的上升序列个数，j这一维只需要开到最长上升子序列的个数，用vector动态开就行，超过最长上升子序列长度的直接输出0,。然后这题序列是随机的，根据结论随机序列的最长上升子序列的期望长度是sqrtn，所以这样是可做的。

太久没写，代码好挫。然后lis还求错了..555

 

 

#include <bits/stdc++.h>
using namespace std;
const int maxn=1e4+100;
const int mod=1e9+7;
int n, m;
vector<int>a, b;
//unsigned int a[maxn];
//unsigned int bit[maxn][105];
vector<int>bit[maxn];
int dp[maxn];
//int dp[maxn][maxn];
int lowbit(int x)
{
    return x&-x;
}
void add(int x, int pos, int y)
{
    while(x<n+5)
    {
        bit[x][pos]=(bit[x][pos]+y);
//        cout<<x<<endl;

        bit[x][pos]=bit[x][pos]>=mod?bit[x][pos]-mod:bit[x][pos];
        x+=lowbit(x);
    }
}

int  query(int x, int pos)
{
    int res=0;
    while(x>0)
    {
        res=(res+bit[x][pos])%mod;
        res=res>=mod?res-mod:res;
        x-=lowbit(x);
    }
    return res%mod;
}
int main()
{
    int e=1, t, i, j;
    cin>>t;
    while(t--)
    {

         scanf("%d", &n);
         a.clear();
         b.clear();
         int x;
         dp[1]=1;
         int m=1;
         for(i=0; i<n; i++)
         {
            scanf("%d", &x);
            a.push_back(x);
//            b[i]=a[i];
            b.push_back(x);
            dp[i]=1;
            for(int j=0; j<i; j++)
            {
                if(a[j]<a[i])
                {
                    dp[i]=max(dp[i], dp[j]+1);
                }
            }
            m=max(dp[i], m);
            bit[i].clear();
         }
         for(int i=n; i<=n+10; i++)bit[i].clear();
//         cout<<m<<endl;
         for(int i=0; i<=n+6; i++)
         {
             for(int j=1; j<=m+6; j++)bit[i].push_back(0);
         }
         sort(b.begin(), b.end());
         int pos, sum;
//         b.erase(unique(b.begin(), b.end()),b.end());
         for(i=1; i<=n; i++)
        {
            pos=lower_bound(b.begin(), b.end(), a[i-1])-b.begin()+1;
//            cout<<pos<<endl;
            add(pos, 1, 1);
//            cout<<2333<<endl;
            for(j=1; j<=m; j++)
            {
                sum=query(pos-1, j);
                add(pos, j+1, sum);
            }
//            printf("%d %d\n", pos, query(pos, m));
        }
        printf("Case #%d:",e++);
        for(i=1; i<=m; i++)printf(" %d", query(n, i));
        for(; i<=n; i++)printf(" 0");
        printf("\n");
        for(int i=0; i<n+10; i++)bit[i].clear();

    }
}