[CTS2019]氪金手游（容斥+树形背包DP）

降智好题。本蒟蒻VP时没想到怎么做被题面迷惑了，只会20分的“好”成绩。简直自闭了。

首先显然度为0的点是白给的，根据等比数列求和公式即可求得。然后考虑这个树如果是一颗外向树，就是每个点先父亲再自己。然后直接DP，令f[i][j]表示子树i内Σw=j的概率，转移时直接用背包转移一发即可。边是正向的直接转移，反向的加上去掉该限制的答案，并减去反向的答案。复杂度显然是O(n²)

#include<bits/stdc++.h>
using namespace std;
const int N=,mod=;
int n,cnt,ans,p[N][],f[N][N*],g[N*],sz[N],inv[N*],hd[N],v[N<<],nxt[N<<],w[N<<];
void add(int&x,int y){x=x+y>=mod?x+y-mod:x+y;}
int qpow(int a,int b)
{
    int ret=;
    while(b)
    {
        if(b&)ret=1ll*ret*a%mod;
        a=1ll*a*a%mod,b>>=;
    }
    return ret;
}
void adde(int x,int y)
{
    v[++cnt]=y,w[cnt]=,nxt[cnt]=hd[x],hd[x]=cnt;
    v[++cnt]=x,w[cnt]=,nxt[cnt]=hd[y],hd[y]=cnt;
}
void dfs(int u,int fa)
{
    f[u][]=;
    for(int i=hd[u];i;i=nxt[i])
    if(v[i]!=fa)
    {
        dfs(v[i],u);
        for(int j=;j<=sz[u]+sz[v[i]];j++)g[j]=;
        for(int j=;j<=sz[u];j++)
        for(int k=;k<=sz[v[i]];k++)
        {
            int val=1ll*f[u][j]*f[v[i]][k]%mod;
            if(w[i])add(g[j+k],val);else add(g[j+k],mod-val),add(g[j],val);
        }
        sz[u]+=sz[v[i]];
        for(int j=;j<=sz[u];j++)f[u][j]=g[j];
    }
    memset(g,,sizeof g);
    for(int i=;i<=sz[u];i++)
    for(int j=;j<=;j++)
    add(g[i+j],1ll*f[u][i]*p[u][j]%mod*j%mod*inv[i+j]%mod);
    sz[u]+=;
    for(int i=;i<=sz[u];i++)f[u][i]=g[i];
}
int main()
{
    scanf("%d",&n);
    for(int i=,x,y,z,sum;i<=n;i++)
    {
        scanf("%d%d%d",&x,&y,&z);
        sum=qpow(x+y+z,mod-);
        p[i][]=1ll*x*sum%mod,p[i][]=1ll*y*sum%mod,p[i][]=1ll*z*sum%mod;
    }
    inv[]=inv[]=;for(int i=;i<=n*;i++)inv[i]=1ll*(mod-mod/i)*inv[mod%i]%mod;
    for(int i=,x,y;i<n;i++)scanf("%d%d",&x,&y),adde(x,y);
    dfs(,);
    for(int i=;i<=*n;i++)add(ans,f[][i]);
    printf("%d",ans);
}