hdu5125 树状数组+dp

 hdu5125 他说的是n个人每个人都有两个气球a,b，气球各自都有相应的体积，现在让他们按照序号排列好来，对他们的a气球体积值计算最长上升子序列，对于这整个排列来说有m次机会让你将a气球替换成b气球(允许不使用完)，问最后的最长上升子序列 的长度是多少，哈哈，当然用dp的思想我们很容易就能知道状态的转移 dp[1000][1000][2],但是苦于状态转移的复杂度太大了达到了 n*n*m肯定受不了，那好我们可以列出这个方程的转移方法

(0表示a气球1为b气球)

dp[i][j][0]=max(  dp[k][j][0]     +1(a[k]<a[i]) , dp[k][j][j]     +1(b[k]<a[i]) )0<=k<i 

dp[i][j][1]=max(  dp[k][j-1][0] +1(a[k]<b[i]) , dp[k][j-1][1] +1(b[k]<b[i])  )0<=k<i

想想优化方法 看来还是需要有换个想法的能力 啊 ！

通过建立m棵树状数组 第j棵树 表示 使用了 j 个 能 量 的 时 候 每个位置所能到达的最高点，说清楚一点就是讲这n*2个气球体积进行离散，得到了树状数组的每个节点从第0个人枚举到第n-1个，比如到达了第i个人 那么 到达他时使用了j个机会的 方案是不是 就是去找比ai(假设离散后在k这个位置)小的那些气球使用j个机会时达到的最长序列加1呢，好那么现就使用树状数组去计算出前k-1 项的 最大值S，然后将S这个值插入当前这个树状数组中 ，那么现在 我们考虑使用b气球可以去j-1 这棵树上去找 然后得到的值加1 插入第j棵树上，哇这样很完美啊！ 瞬间将复杂度减了2个0，这样一直不断地做下去直到结束.

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <string.h>
using namespace std;
const int maxn=;
int a[maxn],b[maxn];
int c[maxn][maxn*],h[maxn*],L;
int lowbit(int x){
 return x&(-x);
}
void add(int loc, int floor,int v){
      while(loc<=L){
         c[floor][loc]=max(v,c[floor][loc]);
         loc+=lowbit(loc);
      }
}
int sum(int loc, int floor){
     int ans=;
     while(loc>){
         ans=max(c[floor][loc],ans);
         loc-=lowbit(loc);
     }
     return ans;
}

int main()
{
    int cas;
    scanf("%d",&cas);
    while(cas--){
         int n,m;
         scanf("%d%d",&n,&m);
         L=;
         for(int i=; i<n; ++i){
            scanf("%d%d",&a[i],&b[i]);
            h[L++]=a[i]; h[L++]=b[i];
         }
         memset(c,,sizeof(c));
         sort(h,h+L);
         L=unique(h,h+L)-h;
         int ans=;
         for(int i=; i<n; ++i){
            int loca = lower_bound(h,h+L,a[i])-h+;
            int locb = lower_bound(h,h+L,b[i])-h+;
            int val;
            for(int j=min(i+,m);j>; j--){
                   val = sum(loca-,j);
                   ans=max(ans,val+);
                  add(loca,j,val+);
                    val = sum(locb-,j-);
                   ans=max(ans,val+);
                add(locb,j,val+);
            }
                 val = sum(loca-,);
                 ans=max(ans,val+);
                add(loca,,val+);
         }
         printf("%d\n",ans);
    }
    return ;
}