BZOJ 1103: [POI2007]大都市meg( 树链剖分 )

早上数学考挂了...欲哭无泪啊下午去写半个小时政治然后就又可以来刷题了..

树链剖分 , 为什么跑得这么慢...

---------------------------------------------------------------------

#include<cstdio>#include<algorithm>#include<iostream>#include<cstring> #define rep( i , n ) for( int i = 0 ; i < n ; ++i )#define clr( x , c ) memset( x , c , sizeof( x ) )#define REP( x ) for( edge* e = head[ x ] ; e ; e = e -> next ) using namespace std; const int maxn = 250000 + 5;const int maxnode = 525000; int n; struct edge { int to; edge* next;}; edge* pt , E[ maxn << 1 ];edge* head[ maxn ]; void init() { pt = E; clr( head , 0 );} inline void Add( int u , int v ) { pt -> to = v; pt -> next = head[ u ]; head[ u ] = pt++;} #define add_edge( u , v ) Add( u , v ) , Add( v , u ) int top[ maxn ] , fa[ maxn ] , size[ maxn ] , dep[ maxn ] , son[ maxn ];int id[ maxn ] , id_cnt = 0 , TOP; int dfs( int x ) { size[ x ] = 1; son[ x ] = -1; REP( x ) { int to = e -> to; if( to == fa[ x ] ) continue; fa[ to ] = x; dep[ to ] = dep[ x ] + 1; dfs( to ); size[ x ] += size[ to ]; if( son[ x ] == -1 || size[ son[ x ] ] < size[ to ] ) son[ x ] = to; }} void DFS( int x ) { id[ x ] = ++id_cnt; top[ x ] = TOP; if( son[ x ] != -1 ) DFS( son[ x ] ); REP( x ) { int to = e -> to; if( to == fa[ x ] || to == son[ x ] ) continue; DFS( TOP = to ); }} void DFS_init() { dfs( dep[ 0 ] = 0 ); DFS( TOP = 0 );} int sum[ maxnode ] , L , R;bool tag[ maxnode ]; #define L( x ) ( ( x ) << 1 )#define R( x ) ( L( x ) ^ 1 )#define M( l , r ) ( ( l + r ) >> 1 ) void ADD( int x , int l , int r ) { if( tag[ x ] ) return; if( L <= l && r <= R ) { tag[ x ] = true; sum[ x ] = 0; } else { int m = M( l , r ); if( L <= m ) ADD( L( x ) , l , m ); if( m < R ) ADD( R( x ) , m + 1 , r ); sum[ x ] = sum[ L( x ) ] + sum[ R( x ) ]; }} int query( int x , int l , int r ) { if( L <= l && r <= R ) return sum[ x ]; int m = M( l , r ); return ( L <= m ? query( L( x ) , l , m ) : 0 ) + ( m < R ? query( R( x ) , m + 1 , r ) : 0 );} void modify( int x , int y ) { while( top[ x ] != top[ y ] ) { if( dep[ top[ x ] ] < dep[ top[ y ] ] ) swap( x , y ); L = id[ top[ x ] ] , R = id[ x ]; ADD( 1 , 1 , n ); x = fa[ top[ x ] ]; } if( x == y ) return; if( dep[ x ] < dep[ y ] ) swap( x , y ); L = id[ y ] + 1 , R = id[ x ]; ADD( 1 , 1 , n );} int Q( int x , int y ) { int ans = 0; while( top[ x ] != top[ y ] ) { if( dep[ top[ x ] ] < dep[ top[ y ] ] ) swap( x , y ); L = id[ top[ x ] ] , R = id[ x ]; ans += query( 1 , 1 , n ); x = fa[ top[ x ] ]; } if( x == y ) return ans; if( dep[ x ] < dep[ y ] ) swap( x , y ); L = id[ y ] + 1 , R = id[ x ]; return ans + query( 1 , 1 , n ); } void build( int x , int l , int r ) { sum[ x ] = r - l + 1; if( r > l ) { int m = M( l , r ); build( L( x ) , l , m ); build( R( x ) , m + 1 , r ); }} int main() {// freopen( "test.in" , "r" , stdin ); init(); cin >> n; rep( i , n - 1 ) { int u , v; scanf( "%d%d" , &u , &v ); --u , --v; add_edge( u , v ); } DFS_init(); clr( tag , 0 ); build( 1 , 1 , n ); int m; cin >> m; m += n; char c; while( --m ) { int u , v; scanf( " %c" , &c ); scanf( "%d" , &u ); u--; if( c == 'W' ) { printf( "%d\n" , Q( 0 , u ) ); } else { scanf( "%d" , &v ); --v; modify( u , v ); } } return 0;}

--------------------------------------------------------------------

1103: [POI2007]大都市meg

Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 1334 Solved: 690
[Submit][Status][Discuss]

Description

在经济全球化浪潮的影响下,习惯于漫步在清晨的乡间小路的邮递员Blue Mary也开始骑着摩托车传递邮件了。不过，她经常回忆起以前在乡间漫步的情景。昔日，乡下有依次编号为1..n的n个小村庄，某些村庄之间有一些双向的土路。从每个村庄都恰好有一条路径到达村庄1（即比特堡）。并且，对于每个村庄，它到比特堡的路径恰好只经过编号比它的编号小的村庄。另外，对于所有道路而言，它们都不在除村庄以外的其他地点相遇。在这个未开化的地方，从来没有过高架桥和地下铁道。随着时间的推移，越来越多的土路被改造成了公路。至今，Blue Mary还清晰地记得最后一条土路被改造为公路的情景。现在，这里已经没有土路了——所有的路都成为了公路，而昔日的村庄已经变成了一个大都市。 Blue Mary想起了在改造期间她送信的经历。她从比特堡出发，需要去某个村庄，并且在两次送信经历的间隔期间,有某些土路被改造成了公路.现在Blue Mary需要你的帮助：计算出每次送信她需要走过的土路数目。（对于公路，她可以骑摩托车；而对于土路，她就只好推车了。）

Input

第一行是一个数n(1 < = n < = 2 50000).
以下n-1行，每行两个整数a，b（1 < = a以下一行包含一个整数m（1 < = m < = 2 50000），表示Blue Mary曾经在改造期间送过m次信。
以下n+m-1行，每行有两种格式的若干信息，表示按时间先后发生过的n+m-1次事件:
若这行为 A a b(a若这行为 W a, 则表示Blue Mary曾经从比特堡送信到村庄a。

Output

有m行，每行包含一个整数，表示对应的某次送信时经过的土路数目。

Sample Input

5
1 2
1 3
1 4
4 5
4
W 5
A 1 4
W 5
A 4 5
W 5
W 2
A 1 2
A 1 3

Sample Output

2
1
0
1

HINT

Source