bzoj2337 [HNOI2011]XOR和路径（概率期望dp+gauss）

分析：
之前写过图上的线性基，所以一下子就想到了线性基

我在想是不是要把所有的路径都求出来*概率？
这显然不是很现实。。。

可以尝试着把二进制的每一位分开考虑（以下的讨论都是针对二进制的一位）

所以如果路径上经过了奇数个1，那么 “$2^i$ * 路径上经过奇数个1的概率” 就是这一位的贡献

设计状态：$f[x]$为从$x$点到$n$点，XOR和为1的概率
$f[x]->{f[y] \over deg(x)}$ （$w(x,y)=0$：这条边是0，我们可以放心大胆的走这条边）
$f[x]->{1-f[z] \over deg(x)}$ （$w(x,z)=1$：这条边是1，所以我们希望下一步有偶数个1，概率即为$1-f[z]$）

$f[x]=\sum {f[y] \over deg(x)}+\sum {1-f[z] \over deg(x)}$
$f[x]-\sum {f[y] \over deg(x)}-\sum {1-f[z] \over deg(x)}=0$

因为图上有自环和重边，不能像DAG一样直接推，所以我们要用高斯消元求解

tip

自环只要连一条
方程组中，n号点不算：$f[n]=1$

位运算的问题我们的一般处理办法就是拆位，按位处理

#include<bits/stdc++.h>

using namespace std;

const int N=102;
const double eps=1e-8;
struct node{
    int y,v,nxt;
};
node way[N*N*2];
int n,m,st[N],tot=0;
double a[N][N],du[N];

void add(int u,int w,int z) {
    du[u]++;
    tot++;way[tot].y=w;way[tot].v=z;way[tot].nxt=st[u];st[u]=tot;
}

void gauss() {
    int now=1,to;
    for (int i=1;i<=n;i++) {
        for (to=now;to<=n;to++)
            if (fabs(a[to][i])>eps) break;
        if (to>n) continue;
        if (to!=now)
            for (int j=1;j<=n+1;j++)
                swap(a[to][j],a[now][j]);
        for (int j=1;j<=n;j++)
            if (j!=now) {
                double t=a[j][i]/a[now][i];
                for (int k=1;k<=n+1;k++)
                    a[j][k]-=t*a[now][k];
            }
        now++;
    }
    for (int i=1;i<=n;i++) a[i][n+1]/=a[i][i];
}

int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=m;i++) {
        int u,w,z;
        scanf("%d%d%d",&u,&w,&z);
        add(u,w,z);
        if (u!=w) add(w,u,z);
    }

    double ans=0;
    for (int _=0;_<=30;_++) {
        memset(a,0,sizeof(a));
        for (int i=1;i<n;i++) {     //n结点不算 
            a[i][i]++;
            for (int j=st[i];j;j=way[j].nxt) {
                if (way[j].v>>_&1) {
                    a[i][way[j].y]+=(1.0/du[i]);
                    a[i][n+1]+=(1.0/du[i]);
                }
                else {
                    a[i][way[j].y]-=(1.0/du[i]);
                }
            }
        }
        gauss();
        ans+=(double)(1<<_)*a[1][n+1];
    }
    printf("%0.3lf\n",ans);
    return 0;
}