bzoj2337 [HNOI2011]XOR和路径（期望dp+按位+高斯消元）

我们把异或拆开，一位位的考虑，这样边权就只有0/1了。
考虑期望dp，f[i]表示i节点到n的路径异或和为1的概率。
那么我们有转移：如果i指向y，且边权为1，$f[i]+=(1-f[y])/du[i]$
如果i指向y，且边权为0，$f[i]+=f[y]/du[i]$
特殊的，f[n]=0.
因为转移带环，所以高斯消元来解一下就好啦。
对答案的贡献就是$f[1]*2^k$
注意自环只加一次。
复杂度$O(logwn^3)$

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 110
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,h[N],num=0,du[N];
double ans=0,a[N][N];
struct edge{
    int to,next,val;
}data[20010];
inline void Gauss(){
    for(int i=1;i<=n;++i){
        int r=i;
        for(int j=i+1;j<=n;++j) if(fabs(a[j][i])>fabs(a[r][i])) r=j;
        if(r!=i) for(int j=1;j<=n+1;++j) swap(a[r][j],a[i][j]);
        for(int j=i+1;j<=n+1;++j){
            double t=a[j][i]/a[i][i];
            for(int k=i;k<=n+1;++k) a[j][k]-=a[i][k]*t;
        }
    }for(int i=n;i>=1;--i){
        for(int j=i+1;j<=n;++j) a[i][n+1]-=a[i][j]*a[j][n+1];
        a[i][n+1]/=a[i][i];
    }
}
inline void solve(int id){
    memset(a,0,sizeof(a));
    for(int x=1;x<n;++x){
        a[x][x]=du[x];
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;
            if(data[i].val>>id&1) a[x][y]+=1,a[x][n+1]+=1;
            else a[x][y]-=1;
        }
    }a[n][n]=1;Gauss();
    ans+=(1<<id)*a[1][n+1];
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();
    while(m--){
        int x=read(),y=read(),val=read();
        data[++num].to=y;data[num].next=h[x];h[x]=num;data[num].val=val;du[x]++;
        if(x==y) continue;
        data[++num].to=x;data[num].next=h[y];h[y]=num;data[num].val=val;du[y]++;
    }for(int i=0;i<=30;++i) solve(i);
    printf("%.3lf\n",ans);
    return 0;
}