三维偏序



太神奇的，二分的魅力

给定N个没有顺序三元组(xi,yi,zi),求其中能构成最长上升序列的长度，大于的定义为（xi>xj && yi>yj && zi>zi）

1.在处理区间[L,R]的时候，先二分区间[L, (L+R)/ 2]，递归求这个左区间（二分的原因是我在维护y维的时候难免破坏x维的性质，但是二分之后我还是可以保证左区间的x全都大于右区间的x）。

(注意对x排序的时候，如果x相等，则根据y降序来排，以避免之后x相同，yi>yj的情况被计算在内)

(同理对y排序相等时，也要对z进行降序)

2.这个时候左区间y维已经有序，我们只好把右区间[(L+R)/2，R]，快排一遍，进行归并式的转移（左区间一个指针，右区间一个指针，保证左区间指针所扫过的y一定小于右区间指针所扫过的y），然后在用常规的树状数组维护z。

3.当然，最后还应该把右区间[(L+R)/2，R]按x快排回去，并递归处理这一段。

4.最后归并或者快排维护整个区间的y有序。

详见http://blog.csdn.net/cjoilmd/article/details/6884432

#include<iostream>
#include<stdlib.h>
#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
int N;
struct node
{
	int x,y,z,ori;
}a[300011];
int f[300011];
int ans=0;
int Max=-2147483647,Min=2147483647;
int d[1000002];
bool cmpx(node e1,node e2)
{
		 if(e1.x<e2.x)return 1;
	else if(e1.x>e2.x)return 0;
	else if(e1.y>e2.y)return 1;
	else return 0;
}
bool cmpy(const node e1,const node e2)
{
		 if(e1.y<e2.y)return 1;
	else if(e1.y>e2.y)return 0;
	else if(e1.z>e2.z)return 1;
	else return 0;
}
int lowbit(int x)
{
	return x&(-x);
}
int getmax(int x)
{
	int maxn=-2147483647;
	while(x>=Min)
	   {
	   	    maxn=max(d[x],maxn);
	   	    x-=lowbit(x);
	   }
	return maxn;
}
void updata(int x,int num)
{
	while(x<=Max)
	    {
	    	d[x]=max(num,d[x]);
	    	x+=lowbit(x);
	    }
}
void two(int L,int R)
{
	int mid=(L+R)>>1;
	if(L==R)
	return;
	two(L,mid);
	
	sort(a+mid+1,a+R+1,cmpy);
	int p1=L,p2=mid+1,temp=0;
	memset(d,0,sizeof(d));
	int s=L;
	for(;p2<=R;p2++)
	  {
	  	 while(a[p1].y<a[p2].y && p1<=mid)
	  	   {
	  	   	  updata(a[p1].z,f[a[p1].ori]);
	  	   	  p1++;
	  	   }
	  	   f[a[p2].ori]=max(f[a[p2].ori],getmax(a[p2].z-1)+1);
	  	   
	  }
	ans+=temp;
	sort(a+mid+1,a+R+1,cmpx);
	two(mid+1,R);
	sort(a+L,a+R+1,cmpy);
}
int main()
{
	freopen("A.in","r",stdin);
	freopen("A.out","w",stdout);
	while(scanf("%d",&N)!=EOF)
	  {
	  	 for(int i=0;i<=10000;i++)
	  	  f[i]=1;
	  	 memset(d,0,sizeof(d));
	  	 for(int i=0;i<N;i++)
	  	   {
	  	   	  scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].z);
	  	   	  a[i].ori=i;
			  a[i].z++;
	  	   	  Max=max(a[i].z,Max);
	  	   	  Min=min(a[i].z,Min);
	  	   }
	  	 sort(a,a+N,cmpx);
	  	 two(0,N-1);
	  	 for(int i=0;i<N;i++)
	  	   {
	  	   	  ans=max(f[i],ans);
	  	   }
	  	 cout<<ans<<endl;
	  }
	
}