hdu4605Magic Ball Game 树状数组

//给一棵树。树的每个节点的子节点个数是0或2
//对于每个节点都有一个权值w[i]
//一个权值为x的球在每个节点的情况有
//x=w[i] 这个球在该点不向下掉
//x<w[i] 这个球往左节点和右节点掉的概率各为1/2
//x>w[i] 这个球往左节点掉的概率为1/8 ， 往右掉的概率为7/8
//问对于每个权值为x其掉到节点为v的概率


//对于一颗树从一点到还有一点的路径是确定的，仅仅须要记录这条路径中
//往左走的路径中大于x的节点个数 l_grt，小于x的节点个数l_less
//往右走中大于x的节点个数 r_grt , 小于x的节点个数r_less
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
using namespace std ;
const int maxn = 1e6+10 ;
int map[maxn][2] ;
vector<int> vec[maxn] ;
int x[maxn] , y[maxn] ;
int root[maxn] ;
int w_x[maxn] ;
int ?w[maxn] ;
int a[maxn] ;int len ;
int tree[2][maxn] ;
void update(int x , int flag , int dx)
{
? ? while(x < maxn)
? ? {
? ? ? ? tree[flag][x] += dx ;
? ? ? ? x += x&(-x) ;
? ? }
}
int getsum(int x ,int flag)
{
? ? int sum = 0 ;
? ? while(x)
? ? {
? ? ? ? sum += tree[flag][x] ;
? ? ? ? x -= x&(-x) ;
? ? }
? ? return sum ;
}
void dfs(int u)
{
? ? //cout<<u<<endl;
? ? for(int i = 0;i < vec[u].size() ;i++)
? ? {
? ? ? ? int id = vec[u][i] ;
? ? ? ? int pos = lower_bound(a, a + len - 1, w_x[id]) - a + 1;
? ? ? ? int l_less = getsum(pos - 1 , 0) ;
? ? ? ? int r_less = getsum(pos - 1, 1) ;
? ? ? ? int l_all = getsum(len , 0) ;
? ? ? ? int r_all = getsum(len , 1) ;
? ? ? ? int l_grt = l_all - getsum(pos , 0) ;
? ? ? ? int r_grt = r_all - getsum(pos , 1) ;
? ? ? ? if(l_less + r_less + l_grt + r_grt != l_all + r_all)
? ? ? ? {
? ? ? ? ? ? x[id] = y[id] = -1 ;
? ? ? ? ? ? continue ;
? ? ? ? }
? ? ? ? x[id] = r_less ;
? ? ? ? y[id] = l_grt + r_grt + (l_less + r_less)*3 ;
? ? }
? ? for(int i = 0;i < 2;i++)
? ? {
? ? ? ? if(!map[u][i])continue ;
? ? ? ? int v = map[u][i] ;
? ? ? ? int pos = lower_bound(a , a + len - 1 , w[u]) - a + 1 ;
? ? ? ? update(pos , i , 1) ;
? ? ? ? dfs(v) ;
? ? ? ? update(pos , i , -1) ;
? ? }
}
int main()
{
? ? //freopen("in.txt" ,"r" , stdin) ;
? ? int T ;
? ? int N ;int M , Q ;
? ? scanf("%d" , &T) ;
? ? while(T--)
? ? {
? ? ? ? scanf("%d" , &N) ;
? ? ? ? len = 0 ;
? ? ? ? for(int i = 1;i <= N;i++)
? ? ? ? {
? ? ? ? ? ? scanf("%d" , &w[i]) ;
? ? ? ? ? ? a[len++] = w[i] ;
? ? ? ? }
? ? ? ? sort(a , a + N) ;
? ? ? ? len = unique(a , a + len) - a ;
? ? ? ? scanf("%d" ,&M) ;
? ? ? ? memset(map , 0 , sizeof(map)) ;
? ? ? ? memset(root , 0 , sizeof(root)) ;
? ? ? ? while(M--)
? ? ? ? {
? ? ? ? ? ? int u , l , r ;
? ? ? ? ? ? scanf("%d%d%d" , &u , &l , &r) ;
? ? ? ? ? ? map[u][0] = l ; map[u][1] = r ;
? ? ? ? ? ? root[l] = root[r] = 1;
? ? ? ? }
? ? ? ? scanf("%d" , &Q) ;
? ? ? ? for(int i = 1;i <= Q ;i++)
? ? ? ? {
? ? ? ? ? ? int v ;
? ? ? ? ? ? scanf("%d%d" , &v , &w_x[i]) ;
? ? ? ? ? ? vec[v].push_back(i) ;
? ? ? ? }
? ? ? ? int pos ;
? ? ? ? for(int i = 1;i <= N;i++)
? ? ? ? if(!root[i]){pos = i ;break;}
? ? ? ? dfs(pos) ;
? ? ? ? for(int i = 1;i <= Q;i++)
? ? ? ? {
? ? ? ? ? ? if(x[i]==-1)puts("0") ;
? ? ? ? ? ? else printf("%d %d\n" ?, x[i] , y[i]) ;
? ? ? ? }
? ? }
}