hdu4742 Pinball Game 3D

        杭州网赛的1005，题意挺简单，每个球有x,y,z三个坐标，球i得三个坐标同事小于等于球j的坐标时，可以用i去撞击j，撞击后i消失，j继续去撞下一个球。问最多可以撞多少个球和有多少种不同的方法。比赛结束后也一直没什么好的思路，刚开始想着用拓扑排序去搞结果发现极限情况的边太多记录不下..后来看了 http://troysky.com/hdu-4742-pinball-game-3d/  这篇题解，又看了看一些大神的代码，大致是明白怎么做了。

        总体思路其实还是个Dp,dp[k]是一个pair<int,int>用来表示当前激活的球是k时，撞击的球数和撞击的方法数.具体做法用的是归并排序的思想，假定现在要搞区间[l,r]，（l,r）表示当前要处理的球的编号是l到r  （这里的球编号可以看做x坐标的哈希值，编号大的球x坐标一定不小于编号小的球）），区间[l,mid]已经处理完毕，那么对区间l,r按照y坐标排个序，这样y也变得有序了，之后还剩下z坐标，这里可以用线段树或者树状数组去维护一个最优值，对于排序以后的序列，从前到后依次遍历，若dt[i].id小于等于mid,那么dp[dt[a[i].id]]已经处理过了，直接按dt[i].z的值更新到树状数组里即可，对于dt[i].id大于mid的球，根据dt[i].z，需要在0--dt[i].z中中找到dp值最大的，直接用树状数组查询即可。之后递归去搞区间[mid+1,r]就ok了。因为用到了树状数组，所以读取数据的时候要对z坐标做一个哈希。之后再按x为第一关键字排个序，记录一下每个球的id。

        代码如下：

        

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <memory.h>
#include <cmath>
#include <algorithm>
#include <queue>
#include <stack>
#include <map>
using namespace std;
typedef long long ll;
const ll mod=(1<<30)-1;
const int maxn=100000+100;
inline int getint() 
{
    char c=getchar();
    int t=0;
    while (c<'0' || c>'9') c=getchar();
    while (c >='0' && c<='9')
    {
        t=t*10+c-'0';
        c=getchar();
    }
    return t;
}
int tt,n,m,p,q;

struct node
{
    int x,y,z,id;
    node(){}
    node(int a,int b,int c)
    {
        x=a;
        y=b;
        z=c;
    }
    bool operator<(const node& p)const
    {
//        return x<p.x||(x==p.x && y<p.y)||(x&&p.x&&y==p.y && z<p.z);
        if (x!=p.x) return x<p.x;
        if (y!=p.y) return y<p.y;
        return z<p.z;
    }
}dt[maxn],tp[maxn];
int valz[maxn];
typedef pair<int,int> pii;
pii dp[maxn];
void merge(pii& p,pii q)
{
    if (p.first<q.first) p=q;
    else if (p.first==q.first)
    {
        p.second+=q.second;
        p.second=p.second&mod;
    }
}
struct BIT
{
    pii c[maxn];
    void init()
    {
        memset(c,0,sizeof c);
    }
    inline int lowbit(int x)
    {
        return x&-x;
    }
    void add(int pos,pii t)
    {
        for(;pos<=n; pos+=lowbit(pos))
        {
            merge(c[pos],t);
        }
    }
    void clear(int pos)
    {
        for(;pos<=n;pos+=lowbit(pos))
        {
            c[pos]=pii(0,0);
        }
    }
    pii query(int pos)
    {
        pii ret(0,0);
        for(;pos; pos-=lowbit(pos))
        {
            merge(ret,c[pos]);
        }
        return ret;
    }

}bit;

void slove(int l,int r)
{
    if (l==r) return;
    int mid=(l+r)>>1;
    slove(l,mid);

    int sp=0;
    for (int i=l; i<=r;i++)
    {
        tp[sp]=dt[i];
        tp[sp].x=0;
        sp++;
    }
    sort(tp,tp+sp);
    for (int i=0; i<sp; i++)
    {
        if(tp[i].id<=mid)
        bit.add(tp[i].z,dp[tp[i].id]);
        else
        {
            pii res=bit.query(tp[i].z);
            res.first++;
            merge(dp[tp[i].id],res);
        }
    }
    for (int i=0; i<sp; i++)
    {
        if (tp[i].id<=mid) bit.clear(tp[i].z);
    }
    slove(mid+1,r);
}

int main()
{
//   freopen("in.txt","r",stdin);
    tt=getint();
    while(tt--)
    {
        bit.init();
        n=getint();
        int tolz=0;
        for (int i=0; i<n; i++)
        {
            dt[i].x=getint();
            dt[i].y=getint();
            dt[i].z=getint();
            valz[tolz++]=dt[i].z;
        }
        sort(valz,valz+tolz);
        tolz=unique(valz,valz+tolz)-valz;
        for (int i=0; i<n;i++)
        {
            dt[i].z=lower_bound(valz,valz+tolz,dt[i].z)-valz+1;
        }
        sort(dt,dt+n);
        for (int i=0;i<n;i++)
        {
            dt[i].id=i;
        }
        for (int i=0;i<=n;i++)
        dp[i]=pii(1,1);
        slove(0,n-1);
        pii ans(0,0);
        for (int i=0; i<n; i++)
        merge(ans,dp[i]);
        printf("%d %d\n",ans.first,ans.second);
    }
    return 0;
}