apio2009 convention

话说这题非常恶心，我整整写了上午4节自习+三节晚自习
题目大意：
  给你一些线段，让你选出最多的不相交的线段，并且让这些线段的字典序最小。
题解：
  首先，如果是不要求字典序的话，应该很简单吧，按线段的结束时间从大到小排序， 然后贪心地求解。
  好，这个思想很重要，我们先用这个方法来求出最多的线段数num，然后按字典序一个一个往里放，如果这条线段放进去不与已经放进去的线段重叠并且能够构成最优解，那么就把它放进去，这个题的思想就是这样的。
  那么如果要知道它是不是与原来放进去的线段重叠很简单，只要用一棵平衡树或者线段树就可以了，我写的平衡树，那么如何判断它是否能构成最优解呢？这就是这道题的关键。
  先定义闭区间[l，r]内最多放的线段数量为find（l，r）
  每次判断一条线段能否能构成最优解时，先找到这条线段往左往右能到达的最大距离，记为l，r，然后，插入这条线段后，这段区间的最多放的线段数量=find（l，a[x]-1)+find（b[x]+1，r）+1。这个很显然吧，然后就是find的求法了。
  我们先把所有的线段按右端点升序排序，然后把所有包含小线段的大线段都删去，你会发现，连左端点也变成了升序的。这有什么好处呢，就是你以后再用这些线段搜的时候就是最优的了。
  定义f[i，j]表示从第i条线段开始取2^j条不相交的线段能达到的最近的端点。这个是不是跟解决rmq问题的st表很相似？对，就是用的st表的思想。同样，这个东西能够在nlogn的时间内解决。先预处理一个get[i]表示i这个点右边的第一条线段。
  然后f[i，j]=f[get[d[f[i,j-1]]+1]，j-1]，这样就能在nlogn内预处理出f数组了。
  这次find就应该知道怎样求了吧？枚举一个左端点i，不断减小j，看是否达到右端点，达到了就加到find上面，然后将i移动到f[i，j]后边，直到j=0为止。
  额，最后输出很简单。
  sui哥用的线段树但貌似没我快。。sbt果然神奇。
附程序，250+行：

program convention; var ctree:array[0..200000] of record size,l,r,a,b:longint; end; a,b,c,d,id,pos:array[0..200000] of longint; nb,get:array[0..400000] of longint; f:array[0..200000,0..30] of longint; v:array[0..200000] of boolean; n,len,l,num,root,tot,fin:longint; {-----------------------------------------} procedure qsort(x,y:longint); var g,t,s,l,r:longint; begin t:=b[(x+y) div 2];l:=x; s:=a[(x+y) div 2];r:=y; repeat while (b[l]<t) or ((b[l]=t) and (a[l]>s)) do inc(l); while (b[r]>t) or ((b[r]=t) and (a[r]<s)) do dec(r); if l<=r then begin g:=a[l];a[l]:=a[r];a[r]:=g; g:=b[l];b[l]:=b[r];b[r]:=g; g:=id[l];id[l]:=id[r];id[r]:=g; pos[id[l]]:=l;pos[id[r]]:=r; inc(l);dec(r); end; until l>r; if l<y then qsort(l,y); if x<r then qsort(x,r); end; {-----------------------------------------} procedure qusort(x,y:longint); var g,t,l,r:longint; begin t:=nb[(x+y) div 2]; l:=x;r:=y; repeat while nb[l]<t do inc(l); while nb[r]>t do dec(r); if l<=r then begin g:=nb[l];nb[l]:=nb[r];nb[r]:=g; inc(l);dec(r); end; until l>r; if l<y then qusort(l,y); if x<r then qusort(x,r); end; {-----------------------------------------} procedure change(var x:longint); var l,r,mid:longint; begin l:=1;r:=len; repeat mid:=(l+r) div 2; if x>nb[mid] then l:=mid+1 else r:=mid-1; until l>r; x:=l; end; {-----------------------------------------} procedure tl(var u:longint); var t:longint; begin t:=ctree[u].r; ctree[u].r:=ctree[t].l; ctree[t].l:=u; ctree[t].size:=ctree[u].size; ctree[u].size:=ctree[ctree[u].l].size+ctree[ctree[u].r].size+1; u:=t; end; {-----------------------------------------} procedure tr(var u:longint); var t:longint; begin t:=ctree[u].l; ctree[u].l:=ctree[t].r; ctree[t].r:=u; ctree[t].size:=ctree[u].size; ctree[u].size:=ctree[ctree[u].l].size+ctree[ctree[u].r].size+1; u:=t; end; {-----------------------------------------} procedure maintain(var u:longint;flag:boolean); begin if not flag then if ctree[ctree[ctree[u].l].l].size>ctree[ctree[u].r].size then tr(u) else if ctree[ctree[ctree[u].l].r].size>ctree[ctree[u].r].size then begin tl(ctree[u].l);tr(u); end else exit else if ctree[ctree[ctree[u].r].r].size>ctree[ctree[u].l].size then tl(u) else if ctree[ctree[ctree[u].r].l].size>ctree[ctree[u].l].size then begin tr(ctree[u].r);tl(u); end else exit; maintain(ctree[u].l,false); maintain(ctree[u].r,true); maintain(u,true); maintain(u,false); end; {-----------------------------------------} procedure makeget; var i:longint; begin for i:=1 to l do get[c[i]]:=i; for i:=c[l] downto 1 do if get[i]=0 then get[i]:=get[i+1]; end; {-----------------------------------------} function find(x,y:longint):longint; var i,j:longint; begin if x>=y then exit(0); find:=0; i:=get[x]; j:=fin; while j>=0 do begin while ((f[i,j]=0) or (d[f[i,j]]>y)) and (j>=0) do dec(j); if j>=0 then inc(find,1 shl j); i:=f[f[i,j],1]; end; end; {-----------------------------------------} function cmp(x,st,ed:longint):boolean; begin if find(st,a[x]-1)+find(b[x]+1,ed)+1=find(st,ed) then exit(true) else exit(false); end; {-----------------------------------------} function check(x,st,ed,u:longint):boolean; begin if u=0 then begin if cmp(x,st,ed) then exit(true) else exit(false); end; if a[x]>ctree[u].b then begin check:=check(x,ctree[u].b+1,ed,ctree[u].r); exit; end else if b[x]<ctree[u].a then begin check:=check(x,st,ctree[u].a-1,ctree[u].l); exit; end else exit(false); end; {-----------------------------------------} procedure insert(x:longint; var u:longint); begin if u=0 then begin inc(tot); ctree[tot].size:=1; ctree[tot].a:=a[x]; ctree[tot].b:=b[x]; u:=tot; exit; end; if a[x]>ctree[u].a then begin insert(x,ctree[u].r); inc(ctree[u].size); maintain(u,true); end else begin insert(x,ctree[u].l); inc(ctree[u].size); maintain(u,false); end; end; {-----------------------------------------} procedure init; var i:longint; begin readln(n); for i:=1 to n do begin readln(a[i],b[i]); pos[i]:=i;id[i]:=i; inc(len);nb[len]:=a[i]; inc(len);nb[len]:=b[i]; end; qusort(1,len); len:=1; for i:=1 to n*2 do if nb[len]<>nb[i] then begin inc(len);nb[len]:=nb[i]; end; for i:=1 to n do begin change(a[i]);change(b[i]); end; end; {-----------------------------------------} procedure makef; var i,j:longint; begin for i:=1 to l do f[i,0]:=i; for i:=1 to l do f[i,1]:=f[get[d[f[i,0]]+1],0]; for j:=2 to fin do for i:=1 to l-1 shl j+1 do f[i,j]:=f[f[f[i,j-1],1],j-1]; end; {-----------------------------------------} procedure prepare; var i,tmp:longint; begin root:=0; qsort(1,n); c[1]:=a[1];d[1]:=b[1]; l:=1;i:=1; for i:=1 to n do if (a[i]>c[l]) and (b[i]>d[l]) then begin inc(l);c[l]:=a[i];d[l]:=b[i]; end; fin:=trunc(ln(l)/ln(2)); makeget; makef; tmp:=0;num:=0; for i:=1 to l do if c[i]>tmp then begin inc(num);tmp:=d[i]; end; end; {-----------------------------------------} procedure main; var i:longint; begin prepare; for i:=1 to n do if check(pos[i],1,len,root) then begin v[i]:=true; insert(pos[i],root); end; writeln(num); for i:=1 to n do if v[i] then write(i,' ');writeln; end; {-----------------------------------------} begin assign(input,'convention.in');reset(input); assign(output,'convention.out');rewrite(output); init; main; close(output); end.