题目:给出一个数轴,有一个起点和一个终点,某个人可以
走1,2,3……m步,每一种情况有一个概率,初始有一个方向,
走到头则返回,问到达终点的期望步数为多少。
比较明显的高斯求期望问题,好坑啊,m可能非常大,以至于同一个点走了多次,WA了我好多发
#include <set>
#include <map>
#include <stack>
#include <queue>
#include <deque>
#include <cmath>
#include <vector>
#include <string>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
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;k<equ&&col<var;k++,col++)
{
max_r=k;
for(i=k+1;i<equ;i++)
if(fabs(a[i][col])>fabs(a[max_r][col]))
max_r=i;
if(fabs(a[max_r][col])<eps)return 0;
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+1;j<var;j++)a[k][j]/=a[k][col];
a[k][col]=1;
for(i=0;i<equ;i++)
if(i!=k)
{
x[i]-=x[k]*a[i][k];
for(j=col+1;j<var;j++)a[i][j]-=a[k][j]*a[i][col];
a[i][col]=0;
}
}
return 1;
}
void bfs()
{
queue<int> 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;
}