hdu2262 高斯消元

题目：有一个地图，一个人从某个点出发，问走到花园的期望步数为多少

设某点的期望步数为Ei。

那么目标的Ei=0。

Ei=(Enext1+Enext2……Enextk)/k+1。

为什么是这个公式 因为 如果 Ei-1的期望应该等于所有从Ei出发的点的期望总和

bfs完 套一个高斯消元就ok了

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <string.h>
#include <cmath>
#include <queue>
using namespace std;
#define eps 1e-9
const int MAXN=;
double a[MAXN][MAXN],x[MAXN];//方程的左边的矩阵和等式右边的值，求解之后x存的就是结果
int equ,var;//方程数和未知数个数
int Gauss(int cequ, int cvar)
{
    equ=cequ; var=cvar;
    int i,j,k,col,max_r;
    for(k=,col=;k<equ&&col<var;k++,col++)
    {
        max_r=k;
        for(i=k+;i<equ;i++)
          if(fabs(a[i][col])>fabs(a[max_r][col]))
            max_r=i;
        if(fabs(a[max_r][col])<eps)continue;
        if(k!=max_r)
        {
            for(j=col;j<var;j++)
              swap(a[k][j],a[max_r][j]);
            swap(x[k],x[max_r]);
        }
        x[k]/=a[k][col];
        for(j=col+;j<var;j++)a[k][j]/=a[k][col];
        a[k][col]=;
        for(i=;i<equ;i++)
          if(i!=k)
          {
              x[i]-=x[k]*a[i][col];
              for(j=col+;j<var;j++)a[i][j]-=a[k][j]*a[i][col];
              a[i][col]=;
          }
    }
    return ;
}
char str[][];
bool falg[][];
int dx[]={,,,-};
int dy[]={,-,,};
void bfs(int n,int m)
{
    queue<int>Qy,Qx;
    memset(falg,false,sizeof(falg));
    while(!Qy.empty())Qy.pop();
    while(!Qx.empty())Qx.pop();
    for(int i=; i<n; i++)
    for(int j=; j<m; j++)
        if(str[i][j]=='$'){
            falg[i][j]=true;
            Qx.push(i);
            Qy.push(j);
        }

    while(!Qy.empty()){
        int ux=Qx.front();
        int uy=Qy.front();
        Qx.pop();
        Qy.pop();
        for(int i=; i<; i++)
        {
        int tx=ux+dx[i];
        int ty=uy+dy[i];
        if(tx>=&&ty>=&&tx<n&&ty<m&&str[tx][ty]!='#'&&falg[tx][ty]==false){
                falg[tx][ty]=true;
                Qy.push(ty);
                Qx.push(tx);
          }
        }
    }
}
int main()
{
    int n,m;
    while(scanf("%d%d",&n,&m)==){
        for(int i=; i<n; i++){
            scanf("%s",str[i]);
        }
        memset(a,,sizeof(a));
        memset(x,,sizeof(x));
        bfs(n,m);
        int sx=-,sy=-;
        for(int i=; i<n; i++)
            for(int j=; j<m; j++)
            {

                if(str[i][j]=='#')continue;
                 int cnt=;
                if(str[i][j]=='$'){
                    a[i*m+j][i*m+j]=;
                    x[i*m+j]=;
                    continue;
                }
                if(str[i][j]=='@'){
                    sx=i; sy=j;
                }
                for(int k=; k<; k++){
                    int tx=i+dx[k];
                    int ty=j+dy[k];
                    if(tx>=&&ty>=&&tx<n&&ty<m&&str[tx][ty]!='#'&&falg[tx][ty]){
                        a[i*m+j][tx*m+ty]=-;
                        cnt++;
                    }
                }
                a[i*m+j][i*m+j]=cnt;
                x[i*m+j]=cnt;
            }
        if(sx!=-&&sy!=-&&falg[sx][sy]&&Gauss(n*m,n*m)){
            printf("%.6lf\n",x[sx*m+sy]);
        }else{
          puts("-1");
        }
    }
    return ;
}