bzoj3270 博物馆（期望+高斯消元）

Time Limit: 30 Sec  Memory Limit: 128 MB

有一天Petya和他的朋友Vasya在进行他们众多旅行中的一次旅行，他们决定去参观一座城堡博物馆。这座博物馆有着特别的样式。它包含由m条走廊连接的n间房间，并且满足可以从任何一间房间到任何一间别的房间。

两个人在博物馆里逛了一会儿后两人决定分头行动，去看各自感兴趣的艺术品。他们约定在下午六点到一间房间会合。然而他们忘记了一件重要的事：他们并没有选好在哪儿碰面。等时间到六点，他们开始在博物馆里到处乱跑来找到对方（他们没法给对方打电话因为电话漫游费是很贵的）

不过，尽管他们到处乱跑，但他们还没有看完足够的艺术品，因此他们每个人采取如下的行动方法：每一分钟做决定往哪里走，有Pi 的概率在这分钟内不去其他地方（即呆在房间不动），有1-Pi 的概率他会在相邻的房间中等可能的选择一间并沿着走廊过去。这里的i指的是当期所在房间的序号。在古代建造是一件花费非常大的事，因此每条走廊会连接两个不同的房间，并且任意两个房间至多被一条走廊连接。

两个男孩同时行动。由于走廊很暗，两人不可能在走廊碰面，不过他们可以从走廊的两个方向通行。（此外，两个男孩可以同时地穿过同一条走廊却不会相遇）两个男孩按照上述方法行动直到他们碰面为止。更进一步地说，当两个人在某个时刻选择前往同一间房间，那么他们就会在那个房间相遇。

两个男孩现在分别处在a，b两个房间，求两人在每间房间相遇的概率。

 

第一行包含四个整数，n表示房间的个数;m表示走廊的数目;a,b (1 ≤ a, b ≤ n),表示两个男孩的初始位置。

之后m行每行包含两个整数，表示走廊所连接的两个房间。

之后n行每行一个至多精确到小数点后四位的实数 表示待在每间房间的概率。

题目保证每个房间都可以由其他任何房间通过走廊走到。

 

输出一行包含n个由空格分隔的数字，注意最后一个数字后也有空格，第i个数字代表两个人在第i间房间碰面的概率（输出保留6位小数）

注意最后一个数字后面也有一个空格

Sample Input

2 1 1 2
1 2
0.5
0.5

Sample Output

0.500000 0.500000

对于100%的数据有 n <= 20，n-1 <= m <= n(n-1)/2

---

（可以列个概率方程并暴力矩乘，然鹅并不知道什么时候收敛，精度会出现误差）

设两个人在$(x,y)$的期望经过次数为$f(x,y)$，$(x,y)$转移到$(r,w)$的概率为$a[(r,w)][(x,y)]$，

点$i$连$in[i]$条边，在点$i$停留的概率为$p[i]$

$a[(x,y)][(x,y)]=p[x]*p[y]$

$a[(x,y+1)][(x,y)]=p[x]*(1-p[y])/in[y]$

$a[(x+1,y)][(x,y)]=p[y]*(1-p[x])/in[x]$

$a[(x+1,y+1)][(x,y)]=(1-p[y])/in[y]*(1-p[x])/in[x]$

那么对于每个$f(x,y)$，我们都能得出一个方程：

$f(x,y)= \sum f(r,w)*a[(x,y)][(r,w)]$（$r=x$或存在边$(x,r)$，$w$同理）

酱紫我们就得出了$n^2$个方程。

发现方程组存在环的关系，而对于$n$元1次方程组，常规做法是套个高斯消元上去

把上面的方程移项一下

$ \sum f(r,w)*a[(x,y)][(r,w)]-f(x,y)=0$

蓝后跑高斯消元就好辣

最后对于$f(a,b)$，因为是初始点，所以要期望次数+1

终点只走一次，所以期望经过次数就等于概率

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
using namespace std;
typedef double db;
const db eps=1e-;
inline db Fabs(db x){return x<?-x:x;}
#define N 420
int n,m,w,A,B,id[N][N];
db p[N],a[N][N],sol[N];
vector <int> g[N];
void draw(int x,int y){
    int e=id[x][y],lx=g[x].size(),ly=g[y].size();//lx,ly:x,y的连边数
    for(int i=;i<lx;++i){
        int r=g[x][i];
        a[id[r][y]][e]+=(-p[x])/(1.0*lx)*p[y];
    }
    for(int i=;i<ly;++i){
        int r=g[y][i];
        a[id[x][r]][e]+=(-p[y])/(1.0*ly)*p[x];
    }
    for(int i=;i<lx;++i)
        for(int j=;j<ly;++j){
            int r1=g[x][i],r2=g[y][j];
            a[id[r1][r2]][e]+=(-p[x])/(1.0*lx)*(-p[y])/(1.0*ly);
        }
}
void gauss(){
    for(int i=,x=;i<=w;x=++i){
        for(int j=i+;j<=w;++j) if(Fabs(a[j][i])>Fabs(a[x][i])) x=j;
        if(x!=i) swap(a[x],a[i]);
        if(Fabs(a[i][i])<eps) continue;
        for(int j=;j<=w;++j){
            if(i==j) continue;
            db div=a[j][i]/a[i][i];
            for(int u=i;u<=w+;++u) a[j][u]-=a[i][u]*div;
        }
    }
    for(int i=w;i;--i){
        sol[i]=a[i][w+];
        for(int j=w;j>i;--j) sol[i]-=sol[j]*a[i][j];
        sol[i]/=a[i][i];
    }
    for(int i=;i<=n;++i) printf("%.6f ",sol[id[i][i]]+eps);//防止输出-0.00
}
int main(){
    scanf("%d%d%d%d",&n,&m,&A,&B);
    for(int i=,u,v;i<=m;++i){
        scanf("%d%d",&u,&v);
        g[u].push_back(v);
        g[v].push_back(u);
    }
    for(int i=;i<=n;++i) scanf("%lf",&p[i]);
    for(int i=;i<=n;++i)
        for(int j=;j<=n;++j){
            id[i][j]=++w;
            if(i!=j) a[w][w]=p[i]*p[j];
        }
    for(int i=;i<=n;++i)
        for(int j=;j<=n;++j)
            if(i!=j) draw(i,j);
    for(int i=;i<=w;++i) a[i][i]-=1.0;//移项系数为-1
    a[id[A][B]][w+]=-1.0;//起点+1
    gauss(); return ;
}