NOIP2016提高A组模拟7.17寻找

Posted 无尽的蓝黄

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了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]);
    }
}

以上是关于NOIP2016提高A组模拟7.17寻找的主要内容,如果未能解决你的问题,请参考以下文章

NOIP2016提高A组模拟10.15总结

NOIP2016提高A组模拟8.14总结

NOIP2016提高A组模拟9.17序列

NOIP2016提高A组模拟10.15打膈膜

NOIP2016提高A组模拟9.14数列编辑器

NOIP2016提高A组模拟8.15Throw