【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);
}