题目
题解
突然get到这样路径期望的题目八成是高斯消元
因为路径上的dp往往具有后效性,这就形成了一个方程组
对于本题来说,直接对权值dp很难找到突破口
但是由于异或是位独立的,我们考虑求出每一位的期望
设\\(f[i]\\)为从节点\\(i\\)出发到达N的期望值
有\\(f[i] = \\frac{f[j]}{degree[i]} + \\frac{1 - f[k]}{degree[i]} [edge(i,j) = 0,edge(i,k) = 1]\\)
因为如果出边权值为0,异或之后值不变,等于\\(f[j]\\)的值,
如果权值为1,异或后取反,等于\\(1-f[k]\\)
同时\\(f[n] = 0\\)
列出式子后就是一个n元方程组
最后要注意自环只算该点的一条边
#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define eps 1e-9
#define LL long long int
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define Redge(u) for (int k = h[u]; k; k = ed[k].nxt)
using namespace std;
const int maxn = 105,maxm = 100005,INF = 1000000000;
inline int read(){
int out = 0,flag = 1; char c = getchar();
while (c < 48 || c > 57) {if (c == \'-\') flag = -1; c = getchar();}
while (c >= 48 && c <= 57) {out = (out << 1) + (out << 3) + c - 48; c = getchar();}
return out * flag;
}
int h[maxn],ne = 2;
double de[maxn];
struct EDGE{int to,nxt,w;}ed[maxm];
inline void build(int u,int v,int w){
ed[ne] = (EDGE){v,h[u],w}; h[u] = ne++;
if (u != v){
ed[ne] = (EDGE){u,h[v],w}; h[v] = ne++;
de[u] += 1,de[v] += 1;
}
else de[u] += 1;
}
int n,m,p;
double A[maxn][maxn],ans;
void gause(){
for (int i = 1; i <= n; i++){
int j = i;
for (int k = i + 1; k <= n; k++)
if (fabs(A[k][i]) > fabs(A[j][i])) j = k;
if (fabs(A[j][i]) < eps) exit(0);
double t = A[j][i];
for (int k = i; k <= n + 1; k++) swap(A[i][k],A[j][k]),A[i][k] /= t;
for (j = i + 1; j <= n; j++){
if (fabs(A[j][i]) > eps){
t = A[j][i];
for (int k = i; k <= n + 1; k++)
A[j][k] -= A[i][k] * t;
}
}
}
for (int i = n; i; i--){
for (int j = n; j > i; j--)
A[i][n + 1] -= A[i][j] * A[j][n + 1];
A[i][n + 1] /= A[i][i];
}
}
void solve(){
for (int i = 1; i <= n; i++)
for (int j = 1; j <= n + 1; j++)
A[i][j] = 0;
for (int i = 1; i < n; i++){
A[i][i] = de[i];
Redge(i){
if ((ed[k].w >> p) & 1){
A[i][n + 1] += 1.0;
A[i][ed[k].to] += 1.0;
}else A[i][ed[k].to] -= 1.0;
}
}
A[n][n] = 1;
gause();
ans += (1 << p) * A[1][n + 1];
}
int main(){
n = read(); m = read();
int a,b,w;
for (int i = 1; i <= m; i++){
a = read(); b = read(); w = read();
build(a,b,w);
}
for (p = 0; (1 << p) <= INF; p++) solve();
printf("%.3lf\\n",ans);
return 0;
}