题目链接：

here

题意：

和hdu4742类似。区别就是一部分三元组是直接给出的。还有一部分是用他给的那个函数生成的。还有就是这里的大于是严格的大于a>b必须ax>bx,ay>by,az>bz。

思路:

思路也和hdu4742here类似。只是有几个比较棘手的问题。现在变成严格大于了。对于y还是很好办的。我们在排序y的时候可以使的标号大的排在前面这样就可以防止y和它一样的更新它了。感觉比较麻烦的是x一样怎么办。这个真没想出什么好办法。就只有x和mid+1的x不一样的建一个树状数组bv。然后所有的不管一不一样再建个树状数组av。这样。遇到x等于mid+1的就用bv更新。其它的用av更新就好了。

详细见代码:

#include
  
   
#include
   
     #include
    
      #include
     
       using namespace std; const int INF=0x3f3f3f3f; const int maxn=300030; typedef long long ll; struct node { int x,y,z; inline bool operator <(const node& tt) const { if(x!=tt.x) return x
      
       b; } void update(int mv[],int x,int v,int mn) { for(int i=x;i<=mn;i+=i&-i) mv[i]=max(mv[i],v); } int qu(int mv[],int x) { int ret=0; for(int i=x;i>0;i-=i&-i) ret=max(ret,mv[i]); return ret; } void solve(int l,int r) { if(l==r) { ans=max(ans,dp[l]); sz[l]=po[l].z; return; } int mid=(l+r)>>1,le=l,ri=mid+1,len=r-l+1,ml=mid-l+1,lim,h,i; memmove(tmp+ptr,sy+l,len*sizeof(int)); for(i=0;i
       
        =ml||sz[l+h]>po[tmp[i]].z-1) h--; if(h>=0) { int *mv=po[tmp[i]].x==po[mid+1].x?bv:av; dp[tmp[i]]=max(dp[tmp[i]],qu(mv,h+1)+1); } } solve(mid+1,r); merge(sz+l,sz+l+ml,sz+l+ml,sz+l+len,ty); memmove(sz+l,ty,len*sizeof(int)); } int A,B,C = ~(1<<31),M = (1<<16)-1; int rd() { A = 36969 * (A & M) + (A >> 16); B = 18000 * (B & M) + (B >> 16); return (C & ((A << 16) + B)) % 1000000; } int main() { int m,n,i; while(scanf("%d%d%d%d",&n,&m,&A,&B),m||n||A||B) { for(i=1;i<=n;i++) scanf("%d%d%d",&po[i].x,&po[i].y,&po[i].z); for(i=n+1,n+=m;i<=n;i++) po[i].x=rd(),po[i].y=rd(),po[i].z=rd(); ptr=ans=0; sort(po+1,po+n+1); for(i=1;i<=n;i++) sy[i]=i,dp[i]=1; sort(sy+1,sy+n+1,cmp); solve(1,n); printf("%d\n",ans); } return 0; }