BZOJ3566 SHOI2014概率充电器(动态规划+概率期望)

  设f[i]为i在子树内不与充电点连通的概率。则f[i]=(1-pi)·∏(1-qk+qk·f[k])。

  然后从父亲更新答案。则f[i]=f[i]·(1-qfa+qfa*f[fa]/(1-qfa+qfa*f[i]))。

  比较好想的dp。注意第二个式子可能会除0。

#include<iostream> 
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
    int x=0,f=1;char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
    while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
    return x*f;
}
#define N 500010
int n,P[N],p[N],t=0;
double f[N],ans=0;
struct data{int to,nxt,len;
}edge[N<<1];
void addedge(int x,int y,int z){t++;edge[t].to=y,edge[t].nxt=p[x],edge[t].len=z,p[x]=t;}
void dfs(int k,int from)
{
    f[k]=1-P[k]/100.0;
    for (int i=p[k];i;i=edge[i].nxt)
    if (edge[i].to!=from)
    {
        dfs(edge[i].to,k);
        f[k]*=(f[edge[i].to]-1)*edge[i].len/100+1;
    }
}
void getans(int k,int from)
{
    for (int i=p[k];i;i=edge[i].nxt)
    if (edge[i].to!=from)
    {
        double q=edge[i].len/100.0;
        if (f[edge[i].to]>0) f[edge[i].to]*=1-q+q*f[k]/(1-q+q*f[edge[i].to]);
        getans(edge[i].to,k);
    }
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj3566.in","r",stdin);
    freopen("bzoj3566.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read();
    for (int i=1;i<n;i++)
    {
        int x=read(),y=read(),z=read();
        addedge(x,y,z),addedge(y,x,z);
    }
    for (int i=1;i<=n;i++) P[i]=read();
    dfs(1,1);
    getans(1,1);
    for (int i=1;i<=n;i++) ans+=1-f[i];
    printf("%.6lf",ans);
    return 0;
}

 

posted @ 2018-09-26 19:51  Gloid  阅读(152)  评论(0编辑  收藏  举报