POJ 1065 Wooden Sticks

http://blog.csdn.net/acdreamers/article/details/7626671

学习一下Dilworth定理 推荐一篇写得很好博客

要求最少的覆盖，按照Dilworth定理

最少链划分 = 最长反链长度

所以最少系统 = 最长导弹高度上升序列长度。

之前写的LIS模板不对。。。。。。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
#include<stack>
#include<cmath>
#include<algorithm>
using namespace std;
#define clc(a,b) memset(a,b,sizeof(a))
#define inf 0x3f3f3f3f
const int N=10010;
#define LL long long
const double eps = 1e-5;
const double pi = acos(-1);
// inline int r(){
//     int x=0,f=1;char ch=getchar();
//     while(ch>'9'||ch<'0'){if(ch=='-') f=-1;ch=getchar();}
//     while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
//     return x*f;
// }
const int  maxn = 500010;
int a[maxn];

struct Node{
    int x,y;
};
Node A[N];
int d[N];

bool cmp(Node a,Node b){
    return a.x<b.x||(a.x==b.x&&a.y<b.y);
}

// int Search(int num,int low,int high){
//     int mid;
//     while(low<=high){
//         mid=(low+high)/2;
//         if(num>=b[mid]) low=mid+1;
//         else high=mid-1;
//     }
//     return low;
// }

// int DP(int n){
//     int len,pos;
//     b[1]=a[1];
//     len=1;
//     for(int i=2;i<=n;i++){
//         if(a[i]>=b[len]){
//             b[++len]=a[i];
//         }
//         else {
//             pos=Search(a[i],1,len);
//             b[pos]=a[i];
//         }
//     }
//     return len;
// }

int BSearch(int l,int r,int x)  
{  
    while(l < r)  
    {  
        int m = (l + r) >> 1;  
        if(x <= d[m]) r = m;  
        else l = m + 1;  
    }  
    return l;  
}  
  
int DP(int a[],int n)  
{  
    d[0] = -1;    
    int len = 1;  
    for(int i=1;i<=n;i++)  
    {  
        d[len] = inf;  
        int j = BSearch(0,len,a[i]);  
        if(j == len) len++;  
        d[j] = a[i];  
    }  
    return len - 1;  
} 
int main(){
    int n,T;
    scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
            scanf("%d%d",&A[i].x,&A[i].y);
        sort(A+1,A+n+1,cmp);
        // for(int i=1;i<=n;i++){
        //     printf("x:%d y:%d
",A[i].x,A[i].y);
        // }
        for(int i=1;i<=n;i++)
            a[n-i+1]=A[i].y;
        // for(int i=1;i<=n;i++)
        //     printf("x:%d ",a[i]);
        printf("%d
",DP(a,n));
    }
    return 0;
}