HDU 4118 Time travel (高斯消元+概率dp)


题目:给出一个数轴,有一个起点和一个终点,某个人可以
走1,2,3……m步,每一种情况有一个概率,初始有一个方向,
走到头则返回,问到达终点的期望步数为多少。

比较明显的高斯求期望问题,好坑啊,m可能非常大,以至于同一个点走了多次,WA了我好多发

#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
using namespace std;
#define L(i) i<<1
#define R(i) i<<1|1
#define INF  0x3f3f3f3f
#define pi acos(-1.0)
#define eps 1e-9
#define maxn 2000010
#define MOD 1000000007

const int MAXN=220;


double a[MAXN][MAXN],x[MAXN];//方程的左边的矩阵和等式右边的值,求解之后x存的就是结果
int equ,var;//方程数和未知数个数
int n,m,s,t,d;
double p[MAXN];
int num[MAXN];
int vis[MAXN];
int Gauss()
{
    int i,j,k,col,max_r;
    for(k=0,col=0;kfabs(a[max_r][col]))
            max_r=i;
        if(fabs(a[max_r][col]) q;
    memset(vis,0,sizeof(vis));
    vis[s] = 1;
    q.push(s);
    while(!q.empty())
    {
        int x = q.front();
        q.pop();
        for(int i = 1; i <= m; i++)
        {
            if(fabs(p[i]) < eps)
                continue;
            int y = (x+i) % n;
            if(!vis[y])
            {
                vis[y] = 1;
                q.push(y);
            }
        }
    }
}
int main(void)
{
    int T;
    scanf("%d",&T);
    while (T--)
    {
        scanf("%d%d%d%d%d",&n,&m,&t,&s,&d);
        double cnt = 0;
        for(int i = 1; i <= m; i++)
        {
            scanf("%lf",&p[i]);
            p[i] /= 100;
            cnt += p[i] * i;
        }
        if(s == t)
        {
            printf("0.00\n");
            continue;
        }
        n = 2*n-2;
        if(d == 1)
            s = n - s;
        bfs();
        if(!vis[t] && !vis[n-t])
        {
            printf("Impossible !\n");
            continue;
        }
        memset(a,0,sizeof(a));
        memset(x,0,sizeof(x));
        for(int i = 0; i < n; i++)
        {
            a[i][i] = 1;
            if(i == t || (n-i) == t)
                continue;
            if(!vis[i])
            {
                x[i] = INF;
                continue;
            }
            for(int j = 1; j <= m; j++)
                a[i][(i+j)%n] -= p[j];
            x[i] = cnt;
        }
        equ = n;
        var = n;
        if(!Gauss())
            printf("Impossible !\n");
        else
            printf("%.2f\n",x[s]);
    }
    return 0;
}


你可能感兴趣的:(动态规划—概率dp,杂七杂八—高斯消元)