BZOJ2337 [HNOI2011] XOR和路径
2018 年 03 月 14 日发布.
Description
$n$个点的无向图,从$1$号点出发每次从所有相连的边随机选一个走过去,直到$n$号点结束。求路径上所有边权异或和的期望。$n\leqslant100, m\leqslant10000$,可能有重边自环。
Solution
显然每一位是互不影响的;那每一位分别计算(这一位为1的边共经过奇数次)的概率即可。高消。
Code
#include <algorithm>
#include <cstdio>
const int N = 105;
const int M = 10005;
double A[N][N];
int n;
inline double abs(double x) { return x < 0 ? -x : x; }
double solve() {
for (int i = 1; i <= n; ++i) {
int j = i;
for (int k = i + 1; k <= n; ++k) if (abs(A[k][i]) > abs(A[j][i])) j = k;
for (int k = i; k <= n + 1; ++k) std::swap(A[i][k], A[j][k]);
for (int k = n + 1; k >= i; --k)
[i][k] /= A[i][i];
Afor (int j = i + 1; j <= n; ++j)
for (int k = n + 1; k >= i; --k)
[j][k] -= A[i][k] * A[j][i];
A}
for (int i = n; i; --i)
for (int j = i - 1; j; --j)
[j][n + 1] -= A[j][i] * A[i][n + 1];
Areturn A[1][n + 1];
}
int u[M], v[M], w[M], deg[N];
int main() {
int m;
("%d%d", &n, &m);
scanffor (int i = 0; i < m; ++i) {
("%d%d%d", &u[i], &v[i], &w[i]);
scanf++deg[u[i]];
if (u[i] != v[i]) ++deg[v[i]];
}
double ans = .0;
for (int i = 0; i < 32; ++i) {
std::fill(A[0], A[n], .0);
for (int j = 1; j <= n; ++j) A[j][j] = 1.0;
for (int j = 0; j < m; ++j) {
if ((w[j] >> i) & 1) {
// f[u] += (1 - f[v]) / d[u]
[u[j]][v[j]] += 1.0 / deg[u[j]];
A[u[j]][n + 1] += 1.0 / deg[u[j]];
Aif (u[j] != v[j]) {
[v[j]][u[j]] += 1.0 / deg[v[j]];
A[v[j]][n + 1] += 1.0 / deg[v[j]];
A}
} else {
// f[u] += f[v] / d[u]
[u[j]][v[j]] -= 1.0 / deg[u[j]];
Aif (u[j] != v[j]) A[v[j]][u[j]] -= 1.0 / deg[v[j]];
}
}
std::fill(A[n], A[n + 1], .0);
[n][n] = 1.0;
A+= solve() * (1 << i);
ans }
("%.3lf\n", ans);
printfreturn 0;
}