【HDU】4742 Pinball Game 3D cdq分治+DP

传送门：【HDU】4742 Pinball Game 3D

题目分析：找rabbit要cdq分治的入门题，然后他就给了这题。。。TUT哭瞎啊，调了好久呢。。巨巨的入门题太难了TUT，当然我太渣了也是原因。。。

求3维的最长上升子序列的长度以及对应长度里的方案数。

一共是三维，排序x维（x维就相当于编号，插入顺序），cdq分治解决剩下的两维，用cdq_solve ( l , m )解决了子问题后更新对 ( m + 1 , r )内的影响，然后递归解决cdq_solve ( m + 1 , r )。向右区间递归比用两个子区间更新自己的方式更方便，左区间解决完以后的答案是可以直接用来更新右区间的左区间的。

令dp为一个结构体，用元素v表示最长上升子序列，cnt表示对应的方案数。每次用树状数组维护dp的值就行了。

今天简直没治。。。把z轴当编号扔到dp里面去更新解。。。然后一直wa了十数发以后才发现。。。。都怪自己太懒不造数据TUT。。

代码如下：

#include <cstdio>
#include <vector>
#include <cstring>
#include <algorithm>
using namespace std ;
 
#define REP( i , a , b ) for ( int i = ( a ) ; i <  ( b ) ; ++ i )
#define FOR( i , a , b ) for ( int i = ( a ) ; i <= ( b ) ; ++ i )
#define REV( i , a , b ) for ( int i = ( a ) ; i >= ( b ) ; -- i )
#define travel( e , H , u ) for ( Edge* e = H[u] ; e ; e = e -> next )
#define CLR( a , x ) memset ( a , x , sizeof a )
#define CPY( a , x ) memcpy ( a , x , sizeof a )
#define ls ( o << 1 )
#define rs ( o << 1 | 1 )
#define lson ls , l , m
#define rson rs , m + 1 , r
#define root 1 , 1 , cnt
#define mid ( ( l + r ) >> 1 )

const int MAXN = 100005 ;
const int mod = 1 << 30 ;

struct Node {
	int v , cnt ;
	Node () {}
	Node ( int v , int cnt ) : v ( v ) , cnt ( cnt ) {}
} dp[MAXN] , T[MAXN] ;

struct Point {
	int x , y , z ;
} p[MAXN] ;

int hash[MAXN] , Len ;
int n ;

int cmp1 ( const Point& a , const Point& 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 ;
}

int cmp2 ( const Point& a , const Point& b ) {
	return a.y != b.y ? a.y < b.y : a.z < b.z ;
}

int unique ( int a[] , int n ) {
	int cnt = 1 ;
	sort ( a + 1 , a + n + 1 ) ;
	FOR ( i , 2 , n ) if ( a[i] != a[cnt] ) a[++ cnt] = a[i] ;
	return cnt ;
}

int search ( int x , int l , int r ) {
	while ( l < r ) {
		int m = ( l + r ) >> 1 ;
		if ( hash[m] >= x ) r = m ;
		else l = m + 1 ;
	}
	return l ;
}

void fun ( Node& a , Node b ) {
	if ( a.v < b.v ) a = b ;
	else if ( a.v == b.v ) {
		a.cnt += b.cnt ;
		if ( a.cnt >= mod ) a.cnt -= mod ;
	}
}

void modify ( int x , Node a ) {
	for ( ; x <= Len ; x += x & -x ) fun ( T[x] , a ) ;
}

Node get_max ( int x , Node res = Node ( 0 , 0 ) ) {
	for ( ; x ; x -= x & -x ) fun ( res , T[x] ) ;
	return res ;
}

void clear ( int x ) {
	for ( ; x <= Len ; x += x & -x ) {
		if ( !T[x].v ) return ;
		T[x] = Node ( 0 , 0 ) ;
	}
}

void cdq_solve ( int l , int r ) {
	if ( l == r ) return ;
	int m = ( l + r ) >> 1 , j = l ;
	cdq_solve ( l , m ) ;
	sort ( p + l , p + m + 1 , cmp2 ) ;
	sort ( p + m + 1 , p + r + 1 , cmp2 ) ;
	FOR ( i , m + 1 , r ) {
		while ( j <= m && p[j].y <= p[i].y ) {
			modify ( p[j].z , dp[p[j].x] ) ;
			++ j ;
		}
		Node tmp = get_max ( p[i].z ) ;
		tmp.v ++ ;
		fun ( dp[p[i].x] , tmp ) ;
	}
	sort ( p + m + 1 , p + r + 1 , cmp1 ) ;
	FOR ( i , l , m ) clear ( p[i].z ) ;
	cdq_solve ( m + 1 , r ) ;
}

void solve () {
	scanf ( "%d" , &n ) ;
	FOR ( i , 1 , n ) {
		scanf ( "%d%d%d" , &p[i].x , &p[i].y , &p[i].z ) ;
		hash[i] = p[i].z ;
	}
	sort ( p + 1 , p + n + 1 , cmp1 ) ;
	Len = unique ( hash , n ) ;
	FOR ( i , 1 , n ) {
		p[i].x = i ;
		p[i].z = search ( p[i].z , 1 , Len ) ;
	}
	FOR ( i , 1 , Len ) T[i] = Node ( 0 , 0 ) ;
	FOR ( i , 1 , n ) dp[i] = Node ( 1 , 1 ) ;
	cdq_solve ( 1 , n ) ;
	Node ans = Node ( 0 , 0 ) ;
	FOR ( i , 1 , n ) fun ( ans , dp[p[i].x] ) ;
	printf ( "%d %d\n" , ans.v , ans.cnt ) ;
}

int main () {
	int T ;
	scanf ( "%d" , &T ) ;
	while ( T -- ) solve () ;
	return 0 ;
}