【题解】洛谷P3211XOR和路径[HNOI2011] 概率DP+高斯消元

题目链接
有点像以前写过的一道题，但是写着写着感觉不对劲……大佬博客写的很好，学习了。写完此题赶紧去复习一波高斯消元

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define _rep(i,a,b) for(int i=(a);i<=(b);i++)
const int N=1e4+10;
const int M=110;
int n,m;
struct Edge{
    int v,val,nx;
}edge[N<<1];
int head[M],tot,deg[M];
double dp[M][M];
void addedge(int u,int v,int w)
{
    edge[tot].v=v;
    edge[tot].val=w;
    edge[tot].nx=head[u];
    head[u]=tot++;
}
void gauss()//高斯消元 
{
    _rep(i,1,n)
    {
        int maxn=i;
        while(!dp[maxn][i])maxn++;
        if(i!=maxn)swap(dp[i],dp[maxn]);
        double mul=dp[i][i];
        _rep(j,i,n+1)dp[i][j]/=mul;
        _rep(k,1,n)if(k!=i&&dp[k][i])
        {
            mul=dp[k][i];
            _rep(j,i,n+1)dp[k][j]-=mul*dp[i][j];
        }
    }
}
int main()
{
    //freopen("in.txt","r",stdin);
    scanf("%d%d",&n,&m);
    memset(head,-1,sizeof(head));
    int u,v,w;
    _rep(i,1,m)
    {
        scanf("%d%d%d",&u,&v,&w);
        if(u!=v)addedge(u,v,w),addedge(v,u,w),deg[u]++,deg[v]++;
        else addedge(u,v,w),deg[u]++;
    }
    double ans=0.0;
    _rep(i,0,30)
    {
        memset(dp,0,sizeof(dp));
        _rep(j,1,n-1)
        {
            dp[j][j]=1.0;
            for(int k=head[j];k!=-1;k=edge[k].nx)
            {
                v=edge[k].v;
                w=edge[k].val;
                if((w>>i)&1)dp[j][v]+=1.0/deg[j],dp[j][n+1]+=1.0/deg[j];
                else dp[j][v]-=1.0/deg[j];
            }
        }
        dp[n][n]=1.0;
        gauss();
        ans+=dp[1][n+1]*1.0*(1<<i);
    }
    printf("%.3lf\n",ans);
    return 0;
}