[bzoj2337]XOR和路径

[bzoj2337]XOR和路径

---

位运算不会进位所以就直接一位一位算。

然后列方程，移项然后高斯消元

• 代码

#include<bits/stdc++.h>
using namespace std;
const int N=110,M=20010;
int cnt,hed[N],nxt[M],to[M],val[M],d[N];
inline void adde(int u,int v,int w){
    cnt++;nxt[cnt]=hed[u],to[cnt]=v,val[cnt]=w;hed[u]=cnt;
}
typedef double db;
int n,m;
db f[N][N];

inline void guass(){
    for(int i=1;i<=n;i++){
        int pos=i;
        for(int j=i+1;j<=n;j++)if(fabs(f[j][i])>fabs(f[pos][i]))pos=j;
        if(pos!=i)for(int k=1;k<=n+1;k++)swap(f[i][k],f[pos][k]);
        for(int j=i+1;j<=n;j++){
            db c=f[j][i]/f[i][i];
            for(int k=i;k<=n+1;k++){
                f[j][k]-=c*f[i][k];
            }
        }
    }
    for(int i=n;i;i--){
        db res=0;
        for(int j=i+1;j<=n;j++){
            res+=f[j][n+1]*f[i][j];
        }
        f[i][n+1]=(f[i][n+1]-res)/f[i][i];
    }
    return;
}
db ans=0;
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++){
        int u,v,w;scanf("%d%d%d",&u,&v,&w);
        d[u]++;adde(u,v,w);
        if(u!=v){d[v]++;adde(v,u,w);}
    }
    for(int p=0;p<=31;p++){
        memset(f,0,sizeof(f));
        for(int i=1;i<=n;i++)f[i][i]=1;
        for(int u=1;u<n;u++){
            db POS=1.0/(db)d[u];
            for(int i=hed[u];i;i=nxt[i]){
                int v=to[i];
                if((val[i]>>p)&1){
                    f[u][v]+=POS;
                    f[u][n+1]+=POS;
                }else{
                    f[u][v]-=POS;
                }
            }
        }
        guass();
        ans+=f[1][n+1]*(1<<p);
    }
    printf("%.3lf\n",ans);
}