jzoj 1382. 区间

Description

　　Alice收到一些很特别的生日礼物：区间。即使很无聊，Alice还是能想出关于区间的很多游戏，其中一个是，Alice从中选出最长的不同区间的序列，其中满足每个区间必须在礼物中，另序列中每个区间必须包含下一个区间。
　　编程计算最长序列的长度。

Input

　　输入文件第一行包含一个整数N(1<=N<=100000)，表示区间的个数。
　　接下来N行，每行两个整数A和B描述一个区间(1<=A<=B<=1000000)。

Output

　　输出满足条件的序列的最大长度。

Sample Input

输入1：
3
3 4
2 5
1 6

输入2：
5
10 30
20 40
30 50
10 60
30 40

输入3：
6
1 4
1 5
1 6
1 7
2 5
3 5

Sample Output

输出1：
3

输出2：
3

输出3：
5

Data Constraint

Hint

【样例解释】
　　例3中可以找到长度为5的区间序列是：[1,7]、[1,6]、[1,5]、[2,5]、[3,5]

正解

一眼最长不下降子序列，不过N=100000,过不了N^2,只能用n log n的方法（除非RP爆表）

对于处理二维信息，我们可以在保证前者单调的情况下，讨论后者的单调性。我们令单调不增，那么当且仅当所选序列的值单调不减才存在覆盖关系，于是我们以为第一关键字，以为第二关键字排序，再以为关键字求一遍最长不下降子序列即可。

CODE

uses math;
var
        a,b,f:array[0..100010] of longint;
        len,i,j,k,n,ans,l,r,mid:longint;
procedure sort(l,r:longint);
var
        i,j,t,mid,mid2:longint;
begin
        mid:=a[(l+r) div 2];
        mid2:=b[(l+r) div 2];
        i:=l;
        j:=r;
        repeat
                while (a[i]<mid)or((a[i]=mid)and(b[i]>mid2)) do inc(i);
                while (a[j]>mid)or((a[j]=mid)and(b[j]<mid2)) do dec(j);
                if i<=j then begin
                        t:=a[i];
                        a[i]:=a[j];
                        a[j]:=t;
                        t:=b[i];
                        b[i]:=b[j];
                        b[j]:=t;
                        inc(i);
                        dec(j);
                end;
        until i>j;
        if l<=j then sort(l,j);
        if i<=r then sort(i,r);
end;
begin
        readln(n);
        for i:=1 to n do readln(a[i],b[i]);
        sort(1,n);
        f[0]:=maxlongint; 
        f[1]:=b[1];
        len:=1;
        for i:=2 to n do begin
                if b[i]<=f[len] then begin
                        inc(len);
                        f[len]:=b[i];
                end
                else begin
                        l:=1;
                        r:=len;
                        while l<r do begin
                                mid:=(l+r) div 2;
                                if f[mid]>=b[i] then
                                        l:=mid+1
                                else
                                        r:=mid;
                        end;
                        while b[i]>f[l] do dec(l);
                        f[l+1]:=b[i];
                end;
        end;
        writeln(len);
end.