【NOIP2016提高A组模拟7.17】寻找

题目

Bob和Alice出去度蜜月，但Alice不慎走失，Bob在伤心过后，决定前去寻找Alice。
他们度蜜月的地方是一棵树，共有N个节点，Bob会使用下列DFS算法对该树进行遍历。
starting_time是一个容量为n的数组
current_time = 0
dfs(v):
current_time = current_time + 1
starting_time[v] = current_time
将children[v]的顺序随机排列 (每个排列的概率相同)
// children[v]v的直接儿子组成的数组
for u in children[v]:
dfs(u)
1是这棵树的根，Bob会从1出发，即运行dfs(1)，现在他想知道每个点starting_time的期望值。

分析

发现，当遍历到某个节点是，以该节点为根的子树一定比它的兄弟早遍历。
如果要求x节点，它的兄弟一定会对它有贡献。
对于一个父亲节点F，它的儿子分别是S1、S2、S3······
假设现在求S1，
因为children[v]的顺序随机排列，所以，它的兄弟在它前面的几率是(dfrac{1}{2})。
总贡献
starting_time[S1]=starting_time[F]+(size[S2]+size[S3]+···)/2+1=starting_time[F]+(size[F]-1-size[S1])/2+1

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const int maxlongint=2147483647;
const int mo=1000000007;
const int N=100005;
using namespace std;
double f[N];
int size[N],fa[N],next[N*2],last[N*2],to[N*2],n,m,tot;
int bj(int x,int y)
{
	next[++tot]=last[x];
	last[x]=tot;
	to[tot]=y;
}
int dg(int x)
{
	size[x]=1;
	for(int i=last[x];i;i=next[i])
	{
		int j=to[i];
		if(j!=fa[x])
		{
			dg(j);
			size[x]+=size[j];
		}
	}
}
int dg1(int x)
{
	for(int i=last[x];i;i=next[i])
	{
		int j=to[i];
		if(j!=fa[x])
		{
			f[j]=f[x]+double(size[x]-1-size[j])*0.5+1;
			dg1(j);
		}
	}
}
int main()
{
	scanf("%d",&n);
	for(int i=2;i<=n;i++)
	{
		scanf("%d",&fa[i]);
		bj(fa[i],i);
		bj(i,fa[i]);
	}
	dg(1);
	f[1]=1;
	dg1(1);
	for(int i=1;i<=n;i++)
	{
		printf("%.1lf ",f[i]);
	}
}