题解:
1.因为只需要保存小数点后6位,所以并不需要把1e5这么大的数全部枚举出来
2.double的表示长度也只有2 ^ 1000,所以不能直接把2 ^ 100000算出来,要提前约分
3.double的精度为15–16位,作为常识可以顺便记下来
4.最多需要枚举最大的100个数,因为10e5 * 10e5 * 10e5 / 2 ^ 100 约= 1 / 1e15,几乎再小的数枚举也对结果造成不了影响
总结:
1.虽然解题的时候讲究弱化或者加强条件,但是如果实在想不出来的话,可以在重新考量某些忽略的条件,也许就是解题关键!
2.如果连WA几次,就必须不能再考虑这个题目了,放空自己,不纠结与细节,等过一段时间,回头再看这个题目的时候,就可以从宏观的思考这个问题了
#include<iostream>
#include<cstring>
#include<cstdio>
#include<vector>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
#define MAXN 100005
vector<int>vec[MAXN];
queue<int>q;
int n,d[MAXN];
inline double f(double a,double b)
{
return a * b / (a + b);
}
double solve()
{
double ans = 0;
for(int i = n - 1;i >= max(n - 100,0);i--)
for(int j = n - 1;j > i;j--)
ans += f(d[i],d[j]) / (pow(2.0,n - i) - (n + 1) / pow(2.0,i));
return ans;
}
int main()
{
while(scanf("%d",&n) != EOF)
{
for(int i = 0;i <= n;i++)vec[i].clear();
for(int i = 1;i < n;i++)
{
int u;
scanf("%d",&u);
vec[u - 1].push_back(i);
}
d[0] = 1;
q.push(0);
while(!q.empty())
{
int u = q.front();
q.pop();
for(int i = 0;i < vec[u].size();i++)
{
int v = vec[u][i];
d[v] = d[u] + 1;
q.push(v);
}
}
sort(d,d + n);
double ans = solve();
printf("%.6f\n",ans);
}
}