CRB and Apple HDU - 5406（双线LIS+多颗数组数组优化）

思路

首先，考虑第1个LIS的最后一项是i,第二个LIS的最后一项是j。 首先，考虑第1个LIS的最后一项是i,第二个LIS的最后一项是j。

d[i][j]=max(dp[i]k+1

可以这样考虑状态转移也可以考虑i变化其实都一样，因为谁的第一个LIS根本不重要。这样的状态方程，要求每次转移j时，都要去找前j个的最大值，由此我们想到，可以用数状数组来优化。
可是问题还有一个，不同的i值，对应着不同的值，（意思是：当i=1时，我们查dp[1][j]的最大值时，要找dp[1][k],k<j的最大值，而当i=2时，则找的是dp[2][k]的最大值），为了维护这样的性质，所以我们需要多颗树状数组来维护，每颗分别维护dp[1],dp[2]。。。
最后一个问题：可能大家都纠结这个问题，这样选择到底能不能保证不选择重复的（1,3，2,4，5，这个序列，第一个LIS选了1,3，4,5，第二个会不会选了2,4，5？？）。如果有这样的疑问，说明对离散化树状数组理解不深，我们是一边更新一边查询，对于每个dp[i]都不会有重复选择，所以对整个问题都不会选择到重复元素。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e3+7;
const int inf=0x3f3f3f3f;
int n;
int dp[maxn][maxn];
int tmp[maxn];
struct Point
{
    int h,d;
}p[maxn];
bool cmp(Point p0,Point p1)
{
    return p0.h==p1.h?p0.d<p1.d:p0.h>p1.h;
}
int b[maxn];
int lowbit(int x){ return x&-x;}
void updata(int *bit,int v,int x)
{
    for(int i=x;i<maxn;i+=lowbit(i))
    {
        bit[i]=max(bit[i],v);
    }
}
int query(int *bit,int x)
{
    int ans=0;
    for(int i=x;i;i-=lowbit(i))
    {
        ans=max(bit[i],ans);
    }
    return ans;
}
int main()
{
    #ifndef ONLINE_JUDGE
        freopen("in.txt","r",stdin);
        freopen("out.txt","w",stdout);
    #endif
    int T;
    cin>>T;
    while(T--)
    {
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
        {
            scanf("%d%d",&p[i].h,&p[i].d);
            b[i]=p[i].d;
        }
        sort(p+1,p+1+n,cmp);
        sort(b+1,b+1+n);
        int m=unique(b+1,b+1+n)-&b[1];
        int ans=0;
        memset(dp,0,sizeof(dp));
        memset(tmp,0,sizeof(tmp));
        for(int i=1;i<=n;i++)
        {
            int l=lower_bound(b+1,b+1+m,p[i].d)-b;
            for(int j=1;j<=n;j++)
            {
                tmp[j]=query(dp[j],l)+1;
                ans=max(ans,tmp[j]);
            }
            for(int j=1;j<=n;j++)
            {
                updata(dp[j],tmp[j],l);//更新dp[j][i]
                updata(dp[l],tmp[j],j);//更新dp[i][j]
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}