【SHTSC2014】概率充电器

Description

著名的电子产品品牌SHOI刚刚发布了引领世界潮流的下一代电子产品——概率充电器：

“采用全新纳米级加工技术，实现元件与导线能否通电完全由真随机数决定！SHOI概率充电器，您生活不可或缺的必需品！能充上电吗？现在就试试看吧！”

SHOI概率充电器由n-1条导线连通了n个充电元件。进行充电时，每条导线是否可以导电以概率决定，每一个充电元件自身是否直接进行充电也由概率决定。随后电能可以从直接充电的元件经过通电的导线使得其他充电元件进行间接充电。

作为SHOI公司的忠实客户，你无法抑制自己购买SHOI产品的冲动。在排了一个星期的长队之后终于入手了最新型号的SHOI概率充电器。你迫不及待地将SHOI概率充电器插入电源——这时你突然想知道，进入充电状态的元件个数的期望是多少呢？

Input

第一行一个整数：n。概率充电器的充电元件个数。充电元件由1-n编号。

之后的n-1行每行三个整数a, b, p，描述了一根导线连接了编号为a和b的充电元件，通电概率为p%。

第n+2行n个整数：qi。表示i号元件直接充电的概率为qi%。

Output

输出一行一个实数，为能进入充电状态的元件个数的期望，四舍五入到小数点后6位小数。

Sample Input

输入1：

3

1 2 50

1 3 50

50 0 0

输入2：

5

1 2 90

1 3 80

1 4 70

1 5 60

100 10 20 30 40

Sample Output

输出1：

1.000000

输出2：

4.300000

Data Constraint

对于30%的数据，n≤5000。

对于100%的数据，n≤500000，0≤p,qi≤100。

Analysis

令人懵逼的概率题。
考虑正难则反。
f[i]表示不考虑由父亲导电的情况，i号节点不能导电的概率
显然f[i]=(1-w[i])（自身不导电）(f[j]+(1-f[j])(1-v[i]（导线不能导电）) (j为i的儿子）
g[i]则表示只考虑由父亲导电的情况
可以类似的dp转移，不过要考虑重复的情况。
那么答案累加一下(1-f[i]*g[i])即可
要小心dfs爆栈哦，用bfs来打吧

Code

#include <cstring>
#include <algorithm>
#include <cstdio>
#include <cmath>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define db double
using namespace std;
typedef long long ll;
const int maxn=50005;
const double pi=acos(-1);
struct edge{
    db x,y;
}  a[maxn];
int n,i,x,y,d,p;
db r;
int read(){
    int sum=0;
    char c=getchar(); int p=1;
    while (c<'0'||c>'9'){
        if (c=='-') p=-1;
        c=getchar();
    }
    while (c>='0'&&c<='9'){
        sum=sum*10+c-'0';
        c=getchar();
    }
    return sum*p;
}
db dis(edge a,edge b){a.x=a.x-b.x; a.y=a.y-b.y; return sqrt(a.x*a.x+a.y*a.y);}
edge rot(edge t,int p){
    db tri=1.0*p/180*pi; edge h;
    h.x=t.x*cos(tri)-t.y*sin(tri);
    h.y=t.x*sin(tri)+t.y*cos(tri);
    return h;
}
edge centre(edge A,edge B,edge C){
     db a1=B.x-A.x,b1=B.y-A.y,c1=(a1*a1+b1*b1)/2;
     db a2=C.x-A.x,b2=C.y-A.y,c2=(a2*a2+b2*b2)/2;
     db d=a1*b2-b1*a2;edge ret;
     ret.x=A.x+(c1*b2-c2*b1)/d;
     ret.y=A.y+(c2*a1-c1*a2)/d;
     return ret;
}
void get(){
    int i,j,k;
    random_shuffle(a+1,a+n+1);
    edge c=a[1];r=0;
    fo(i,1,n){
        if (dis(a[i],c)<=r) continue;
        c=a[i],r=0;
        fo(j,1,i-1) {
            if (dis(a[j],c)<=r) continue;
            c.x=(a[i].x+a[j].x)/2,c.y=(a[i].y+a[j].y)/2;
            r=dis(a[j],c);
            fo(k,1,j-1){
                if (dis(a[k],c)<=r) continue;
                c=centre(a[i],a[j],a[k]);
                r=dis(a[i],c);
            }
        }
    }
}
int main(){
    freopen("amplifier.in","r",stdin);
    freopen("amplifier.out","w",stdout);
    scanf("%d",&n);
    fo(i,1,n) {
        x=read(),y=read();
        a[i].x=x*1.0,a[i].y=y*1.0;
    }
    scanf("%d%d",&d,&p);
    fo(i,1,n) a[i]=rot(a[i],-d),a[i].x/=p;
    get(); 
    printf("%.3lf\n",r);
}