YBROJ洛谷P3211:XOR和路径(线性基,期望dp)
生活随笔
收集整理的這篇文章主要介紹了
YBROJ洛谷P3211:XOR和路径(线性基,期望dp)
小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
解析
不難想到第一步利用期望線性性逐位考慮。
然后就變成求一個(gè)布爾變量的期望了,可以直接轉(zhuǎn)化為求概率。
我一開始的想求從1出發(fā)異或和為0/1的概率,然而這個(gè)東西在原點(diǎn)1附近的轉(zhuǎn)移特別別扭…老出現(xiàn)概率大于1的迷惑情況。
然后我就不會(huì)了
正解是反過來想,設(shè) fxf_xfx? 為從 xxx 出發(fā)到 nnn 的路徑異或和為1的概率,然后像類似游走的方法轉(zhuǎn)移就行了。
為什么要這樣?
首先,本題的圖必然聯(lián)通,此時(shí),從一個(gè)點(diǎn) xxx 必然可以走到 nnn,但從 111 不一定能走到 xxx。
比如說 1-3-2 這樣的圖中,1就是走不到2的。
那么在題解的定義中,1?fx1-f_x1?fx? 就使 x?>nx->nx?>n 的路徑異或和為 000 的概率,而在我的定義中,1?fx1-f_x1?fx? 則是 111 走不到 xxx 和 1?>x1->x1?>x 路徑異或和為0的概率和。
代碼
#include<bits/stdc++.h> using namespace std; #define ll long long #define debug(...) fprintf(stderr,__VA_ARGS__) #define OK debug("OK\n") inline ll read(){ll x(0),f(1);char c=getchar();while(!isdigit(c)){if(c=='-') f=-1;c=getchar();}while(isdigit(c)){x=(x<<1)+(x<<3)+c-'0';c=getchar();}return x*f; } const int N=100+20; const int mod=1e9+7; int n,m; struct node{int to,nxt,w; }p[N*N<<1]; int fi[N],cnt; inline void addline(int x,int y,int w){p[++cnt]=(node){y,fi[x],w};fi[x]=cnt; } int mi[40],o=30; double a[N][N],ans,f[N]; int d[N]; void gauss(){for(int i=1;i<=n;i++){int id=i;for(int j=i+1;j<=n;j++) if(abs(a[j][i])>abs(a[id][i])) id=j;if(id!=i) swap(a[id],a[i]);for(int j=i+1;j<=n;j++){for(int k=i+1;k<=n+1;k++) a[j][k]-=a[j][i]/a[i][i]*a[i][k];}}for(int i=n;i>=1;i--){f[i]=a[i][n+1]/a[i][i];for(int j=1;j<i;j++) a[j][n+1]-=a[j][i]*f[i];}return; } void work(int k){memset(a,0,sizeof(a));for(int i=1;i<=n;i++){a[i][i]=1;if(i==n) continue;for(int o=fi[i];~o;o=p[o].nxt){int j=p[o].to;int op=(p[o].w&mi[k])?1:0;if(op){a[i][n+1]+=1.0/d[i];a[i][j]+=1.0/d[i];}else{a[i][j]+=-1.0/d[i];}}}//for(int i=1;i<=n;i++){//for(int j=1;j<=n+1;j++) printf("%.2lf ",a[i][j]);//puts("");//}gauss();ans+=mi[k]*f[1]; //printf("k=%d f=%lf\n",k,f[1]); }signed main(){ #ifndef ONLINE_JUDGEfreopen("a.in","r",stdin);freopen("a.out","w",stdout); #endifmemset(fi,-1,sizeof(fi));cnt=-1;mi[0]=1;for(int i=1;i<=o;i++) mi[i]=mi[i-1]<<1;n=read();m=read();for(int i=1;i<=m;i++){int x=read(),y=read(),w=read();d[x]++;addline(x,y,w);if(x^y) d[y]++,addline(y,x,w);}for(int i=0;i<=o;i++) work(i);printf("%.3lf\n",ans); } /* */總結(jié)
以上是生活随笔為你收集整理的YBROJ洛谷P3211:XOR和路径(线性基,期望dp)的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 客厅与餐厅一体怎么隔 居然可以这样做隔断
- 下一篇: 液化天然气lng和天然气的区别 区别盘点