bzoj2152-[国家集训队]聪聪可可

Description

聪聪和可可是兄弟俩，他们俩经常为了一些琐事打起来，例如家中只剩下最后一根冰棍而两人都想吃、两个人都想玩儿电脑（可是他们家只有一台电脑）……遇到这种问题，一般情况下石头剪刀布就好了，可是他们已经玩儿腻了这种低智商的游戏。他们的爸爸快被他们的争吵烦死了，所以他发明了一个新游戏：由爸爸在纸上画n个“点”，并用n-1条“边”把这n个“点”恰好连通（其实这就是一棵树）。并且每条“边”上都有一个数。接下来由聪聪和可可分别随即选一个点（当然他们选点时是看不到这棵树的），如果两个点之间所有边上数的和加起来恰好是3的倍数，则判聪聪赢，否则可可赢。聪聪非常爱思考问题，在每次游戏后都会仔细研究这棵树，希望知道对于这张图自己的获胜概率是多少。现请你帮忙求出这个值以验证聪聪的答案是否正确。

Input

输入的第1行包含1个正整数n。后面n-1行，每行3个整数x、y、w，表示x号点和y号点之间有一条边，上面的数是w。

Output

以即约分数形式输出这个概率（即“a/b”的形式，其中a和b必须互质。如果概率为1，输出“1/1”）。

Sample Input

5

1 2 1

1 3 2

1 4 1

2 5 3

Sample Output

13/25

【样例说明】

13组点对分别是(1,1) (2,2) (2,3) (2,5) (3,2) (3,3) (3,4) (3,5) (4,3) (4,4) (5,2) (5,3) (5,5)。

【数据规模】

对于100%的数据，n<=20000。

Solution

[模板] 树的重心/点分治/动态点分治

点分治模板题. 似乎还可以dp...

Code

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<set>
#include<map>
using namespace std;
#define rep(i,l,r) for(register int i=(l);i<=(r);++i)
#define repdo(i,l,r) for(register int i=(l);i>=(r);--i)
#define il inline
typedef double db;
typedef long long ll;

//---------------------------------------
const int nsz=20050;

ll n,ans=0;

struct te{int t,pr,v;}edge[nsz*2];
int hd[nsz],pe=1;
void adde(int f,int t,int v){edge[++pe]=(te){t,hd[f],v};hd[f]=pe;}
void adddb(int f,int t,int v){adde(f,t,v);adde(t,f,v);}
#define forg(p,i,v) for(int i=hd[p],v=edge[i].t;i;i=edge[i].pr,v=edge[i].t)

int vi[nsz];
int szp[nsz],sum,maxp[nsz]{2e4+5},rt;

void getrt(int p,int fa){
	szp[p]=1,maxp[p]=0;
	forg(p,i,v){
		if(vi[v]||v==fa)continue;
		getrt(v,p);
		szp[p]+=szp[v];
		maxp[p]=max(maxp[p],szp[v]);
	}
	maxp[p]=max(maxp[p],sum-szp[p]);
	if(maxp[p]<maxp[rt])rt=p;
}

ll cnt[4];
void getdis(int p,int fa,int d){
	++cnt[d];//dis[++pd]=d;
	forg(p,i,v){
		if(v==fa||vi[v])continue;
		getdis(v,p,(d+edge[i].v)%3);
	}
}

void sol1(int p,int v0,int fl){
	memset(cnt,0,sizeof(cnt));
	getdis(p,0,v0%3);
	ans+=cnt[1]*cnt[2]*2*fl;
	ans+=cnt[0]*cnt[0]*fl;
}

void divide(int p){
	int sum0=sum;
	vi[p]=1;
	sol1(p,0,1);
	forg(p,i,v){
		if(vi[v])continue;
		sol1(v,edge[i].v,-1);
		rt=0,sum=sum0-maxp[v],getrt(v,0);
		divide(rt);
	}
}

ll sol(){
	rt=0,sum=n,getrt(1,0);
	divide(rt);
	return ans;
}

ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}

int main(){
	ios::sync_with_stdio(0),cin.tie(0);
	cin>>n;
	int a,b,c;
	rep(i,2,n){
		cin>>a>>b>>c;
		adddb(a,b,c);
	}
	ll ans=sol(),g=gcd(ans,n*n);
	cout<<ans/g<<'/'<<n*n/g;
	return 0;
}