正題

題目鏈接:https://www.luogu.com.cn/problem/P3211

---

題目大意

一個(gè)$n$個(gè)點(diǎn)$m$條邊的無(wú)向圖，從$1$到$n$隨機(jī)游走。求期望路徑異或和。

$2\le n\le 100,1\le m\le 10^4$

---

解題思路

因?yàn)槭钱惢虻钠谕?#xff0c;很難直接處理，所以考慮按位考慮每一位是$1$的概率。

然后$n$很小就是一個(gè)很顯然的高斯消元了。設(shè)$f_i$表示$i～n$是$1$的概率。

$$f_x=\frac{1}{de{g}_x}(\sum _{x?>y,w=1}(1?f_y)+\sum _{x?>y,w=0}{f}_y)$$

時(shí)間復(fù)雜度$O(n^3\log w_i)$

---

code

#include<cstdio> #include<cstring> #include<algorithm> using namespace std; const int N=110; struct node{int to,next,w; }a[N*N*2]; int n,m,tot,deg[N],ls[N]; double f[N],ans; void addl(int x,int y,int w){a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;a[tot].w=w;return; } namespace G{double a[N][N],b[N];void init(){for(int i=1;i<=n;i++){for(int j=1;j<=n;j++)a[i][j]=0;b[i]=0;}return;}void solve(double *f){for(int i=1;i<=n;i++){int z=i;for(int j=i+1;j<=n;j++)if(a[j][i]>a[z][i])z=i;swap(a[i],a[z]);swap(b[i],b[z]);double inv=a[i][i];for(int j=i;j<=n;j++)a[i][j]=a[i][j]/inv;b[i]=b[i]/inv;for(int j=i+1;j<=n;j++){double rate=-a[j][i];for(int k=i;k<=n;k++)a[j][k]+=a[i][k]*rate;b[j]+=b[i]*rate;}}for(int i=n-1;i>=1;i--){for(int j=i+1;j<=n;j++)b[i]-=b[j]*a[i][j]/a[j][j];f[i]=b[i];}return;} }; void solve(int w){G::init();G::a[n][n]=1;for(int x=1;x<n;x++){for(int i=ls[x];i;i=a[i].next){int y=a[i].to;if(a[i].w&w)G::a[x][y]++,G::b[x]++;else G::a[x][y]--;}G::a[x][x]+=deg[x];}G::solve(f);ans+=(double)w*f[1];return; } int main() {scanf("%d%d",&n,&m);for(int i=1;i<=m;i++){int x,y,w;scanf("%d%d%d",&x,&y,&w);deg[x]++;addl(x,y,w);if(x!=y)deg[y]++,addl(y,x,w);}for(int i=0;i<=30;i++)solve(1<<i);printf("%.3lf\n",ans);return 0; } 創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎(jiǎng)勵(lì)來(lái)咯，堅(jiān)持創(chuàng)作打卡瓜分現(xiàn)金大獎(jiǎng)

總結(jié)

以上是生活随笔為你收集整理的P3211-[HNOI2011]XOR和路径【高斯消元】的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。