题目
题意: 每个球有x,y,z三个坐标,球i得三个坐标同事小于等于球j的坐标时,可以用i去撞击j,撞击后i消失,j继续去撞下一个球。问最多可以撞多少个球和有多少种不同的方法。
题目: 三维LIS 球的顺序不是按题目输入的顺序。
先对x排序,然后a[i].x=i; 其实a[i].x就是a[i].id
第45行对[l,r]之间按y排序。注意并不是排序a数组。。 我们46-52行就是单纯的为了计算左边对右边的影响。。。不像别的cdq(还有归并 每次cdq之后还可以保证某一维有序。) 我们每次cdq之后左右x都有序 左x<右x 因为根本没有对a数组排过序!!
t[i].x<=mid时 是在树状数组t[i].z的这个位置插入dp[t[i].x]。否则是查询。就保证了查询的x>=更新的x 而且按y排过序,也保证了查询的y>=更新的y 。直接在树状数组中查询t[i].z即可.
#include<bits/stdc++.h>
using namespace std;
typedef pair<int,int>pii;
const int N=1e5+5,mod=(1<<30);
struct node{int x,y,z;}a[N],t[N];
inline int cmpx(node A,node B){
if(A.x!=B.x) return A.x<B.x;
if(A.y!=B.y) return A.y<B.y;
return A.z<B.z;
}
inline int cmpy(node A,node B){
if(A.y!=B.y) return A.y<B.y;
return A.z<B.z;
}
pii dp[N];
int refl[N],tot;
inline void modify(pii &A,pii B){
if(A.first<B.first) A=B;
else if(A.first==B.first) A.second=(A.second+B.second)%mod;
}
struct BIT{
#define low(i) ((i)&(-i))
pii A[N];
void update(int o,pii T){
while(o<=tot) modify(A[o],T),o+=low(o);
}
pii query(int o){
pii res={0,0};
while(o) modify(res,A[o]),o-=low(o);
return res;
}
void clear(int o){
while(o<=tot){
if(!A[o].first&&!A[o].second) return;
A[o]={0,0},o+=low(o);
}
}
}B;
void cdq(int l,int r){
if(l==r) return;
int mid=(l+r)>>1;
cdq(l,mid);
int num=0;
for(int i=l;i<=r;++i) t[++num]=a[i];
sort(t+1,t+num+1,cmpy);
for(int i=1;i<=num;++i){
if(t[i].x<=mid) B.update(t[i].z,dp[t[i].x]);
else{
pii tmp=B.query(t[i].z);
++tmp.first,modify(dp[t[i].x],tmp);
}
}
for(int i=1;i<=num;++i)
if(t[i].x<=mid) B.clear(t[i].z);
cdq(mid+1,r);
}
int main(){
int T;scanf("%d",&T);
while(T--){
int n;scanf("%d",&n);
for(int i=1;i<=n;++i) scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].z),refl[i]=a[i].z;
sort(a+1,a+n+1,cmpx);
sort(refl+1,refl+n+1);tot=unique(refl+1,refl+n+1)-(refl+1);
for(int i=1;i<=n;++i) a[i].x=i,a[i].z=lower_bound(refl+1,refl+n+1,a[i].z)-refl,dp[i]={1,1};
cdq(1,n);
pii ans={0,0};
for(int i=1;i<=n;++i) modify(ans,dp[i]);
printf("%d %d\n",ans.first,ans.second);
}
}