hdu4742

题意：给定3维的n（<=100000）个点，求最长不下降子序列长度(对于1和2两个点，2可以放在1后面，x1<=x2,y1<=y2,z1<=z2 )，并求出有多少种方案。

思路：裸的cdq分治。。

首先可以先对x排序，就降成二维了。。

定义solve（l,r）为处理[l,r]的区间

那么solve（l，r）为如下：{

solve（l, mid）

对于l,r的点按照y关键字排序

然后从左到右（y增大方向）扫描，对于（l, mid）的点插入将值插入z离散化的数据结构里维护，

对于（mid+1，r）的点直接查询数据结果，更改结果

清空数据结构。

solve（mid+1， r）

}

大致就这样，可以做到nlog^2n的复杂度。。数据结构可以选用bit或者线段树。。不过大家都选用bit吧。。

cdq果然强大，还需慢慢体会呀。。

 /*
  * Author:  Yzcstc
  * Created Time:  2014年10月03日 星期五 22时57分13秒
  * File Name: hdu4742.cpp
  */
 #include <iostream>
 #include <cstdio>
 #include <cstdlib>
 #include <cstring>
 #include <cmath>
 #include <algorithm>
 #include <string>
 #include <vector>
 #include <stack>
 #include <queue>
 #include <set>
 #include <utility>
 #define M0(x) memset(x, 0, sizeof(x))
 #define MP make_pair
 #define PB push_back
 #define repf(i, a, b) for (int i = (a); i <= (b); ++i)
 #define Inf 0x3fffffff
 #define eps 1e-8
 #define pi acos(-1)
 #define maxn 120000
 #define X first
 #define Y second
 using namespace std;
 typedef pair<int, int> pii;
 struct point{
     int x, y, z,id;
     void input(){
         scanf("%d%d%d", &x, &y, &z);
     }
     bool operator<(const point& p) const{
         if (x < p.x) return true;
         if (x == p.x && y < p.y) return true;
         if (x == p.x && y == p.y && z < p.z) return true;
         return false;
     }
 } p[maxn], pt[maxn];
 int t[maxn], n, m;
 pii dp[maxn], v[maxn], zero(, );

 void init(){
      scanf("%d", &n);
      for (int i = ; i <= n; ++i)
           p[i].input(), t[i-] = p[i].z;
      sort(p + , p +  + n);
      sort(t, t + n);
      m = unique(t, t + n) - t;
 //     cout << m << endl;
      for (int i = ; i <= n; ++i)
           p[i].z = lower_bound(t, t + m, p[i].z) - t + ;
 //     repf(i, 1, n) printf("%d\n" ,p[i].z);
      for (int i = ; i <= n; ++i)
           p[i].id = i;
 }

 /** BIT **/
 inline int lowbit(const int& x){
      return x & (-x);
 }

 inline void update(pii &a, const pii &b){
      if (a.X < b.X) a = b;
      else if (a.X == b.X) a.Y += b.Y;
 }

 void modify(int x, const pii& val){
      for ( ; x <= m; x += lowbit(x))
          update(v[x], val);
 }

 pii query(int x){
      pii ret(, );
      for ( ; x >  ; x -= lowbit(x))
           update(ret, v[x]);
      return ret;
 }

 void recover(int x){
      for (; x <= m; x += lowbit(x))
           v[x] = zero;
 }
 /** BIT-end **/
 /** Divide and conquer**/
 pii tmp;
 void solve(const int& l,const int& r){
      if (l == r) return;
      int mid = (l + r) >> ;
      solve(l, mid);
      int sz = ;
      for (int i = l; i <= r; ++i)
          pt[sz] = p[i], pt[sz++].x = ;
      sort(pt, pt+sz);
      for (int i = ; i < sz; ++i)
          if (pt[i].id <= mid)
                modify(pt[i].z, dp[pt[i].id]);
          else
                tmp = query(pt[i].z), ++tmp.X, update(dp[pt[i].id], tmp);
      for (int i = ; i < sz; ++i)
          if (pt[i].id <= mid)
                recover(pt[i].z);
      solve(mid + , r);
 }

 void solve(){
      for (int i = ; i <= n; ++i)
           dp[i].X = dp[i].Y = ;
      solve(, n);
      pii ans(, );
      for (int i = ; i <= n; ++i)
          update(ans, dp[i]); // printf("i = %d: %d %d\n", i, dp[i].first, dp[i].second);
      printf("%d %d\n", ans.X, ans.Y);
 }

 int main(){
  //   freopen("a.in", "r", stdin);
 //    freopen("a.out", "w", stdout);
     int cas;
     scanf("%d", &cas);
     while (cas--){
          init();
          solve();
     }
     return ;
 }