HDU 3458

http://acm.hdu.edu.cn/showproblem.php?pid=3458

可以转换为LIS做，一开始按这个题的上一题做，想在原来的n^2的复杂度上面改成以前学过求LIS的n*log(n)的算法，一开始想以每个矩阵为一个点，按其左下角的点从前到后的排序，然后在逐个log(n)的复杂度查找更新，但是王牛写了一遍，说不行，网上面的人说也不行，一开始没明白怎么不行，刚刚想了很久其实真的不行，因为在后面找的时候是根据y单调递增性来二分查找的，如何按照前面的说话一矩形为一点，那么在保证x递增的情况下并不能保证y也是递增的，比如：

 

像这个图，按矩形为点的方法排序后结果插入F之后就会破坏长度为2 和 3 的y的单调性，就不能二分查找了。。

所以我们这里通过另外的一种方法实现，牛人说的，把一个矩形分成两个点，左下角和右上角；然后按x，y排序。这样就会保证X是不下降的。那么我们可以根据长度查找y，然后利用右上角哪一点更新y值，这样总是可以保证右上角y的最小，具有单调性。。。

#include<iostream> #include<algorithm> using namespace std; #define N 200020 struct Node { int x,y,lt,id ; }point[N]; int n,len,maxlp[N],by[N]; int cmp(Node a,Node b) { if(a.x != b.x) return a.x < b.x ; if(a.y == b.y) return a.lt > b.lt ; return a.y < b.y ; } int find(int len,int _y) { int mid, left = 1, right = len ; while(left <= right) { mid = (left+right) >> 1 ; if(by[mid] >= _y) right = mid-1 ; else left = mid+1 ; } return left ; } int main() { while(scanf("%d",&n) && n) { for(int i=0;i<n;i++) { scanf("%d%d",&point[2*i].x,&point[2*i].y); point[2*i].lt = 1 ; scanf("%d%d",&point[2*i+1].x,&point[2*i+1].y); point[2*i+1].lt = 0 ; point[2*i].id = point[2*i+1].id = i ; } n *= 2 ; sort(point,point+n,cmp) ; len = 0 ; for(int i=0;i<n;i++) { if(point[i].lt == 1) maxlp[point[i].id] = find(len,point[i].y) ; else { int tmp = maxlp[point[i].id] ; if(tmp > len) by[++len] = point[i].y ; else by[tmp] = min(by[tmp],point[i].y) ; } } printf("%d/n",len); } return 0; }