题目地址：CDOJ...但是我进不去(;´༎ຶД༎ຶ`)这道题目前还没有交给评测鸭测试。。不过搜一下网上的代码差不多也都是这个亚子

题目：

众所周知，柱爷的数学非常好，尤其擅长概率论！某日柱爷在喵蛤蛤村散步，无意间踏入了远古法阵！法阵很奇怪，是一个长度为N的走廊，初始时柱爷在最左边，现在柱爷要到最右边去!

柱爷的行动方式如下:

每个回合柱爷会投一次骰子，根据骰子上的点数X，柱爷会相应的往右边移动X步.

骰子的数值是 1到6，取到每面的概率相同

在某些位置可能有传送门，一旦柱爷在该回合结束后在这个位置上，会被强制传送到传送门的另外一边

传送门是单向的,同时每个位置不会有超过1个传送门,同时不会存在a→b,b→c这种情况

在任意时刻柱爷都必须保证在法阵内，也就说如果在这一回合结束后柱爷的位置在法阵外，那么这回合柱爷将什么都不做

那么请问柱爷到达最右边的期望回合数是多少呢？或者是永远都无法到达？ 
Input

第一行两个整数N ,M，分别表示法阵的长度和传送门的数量

接下来M 行，每行两个整数u,v，表示从u到v有一扇传送门

output:

输出柱爷到法阵最右边的期望回合数，如果无法到达输出-1

数据保证:

1≤N≤300
10≤M≤[N−22]

解题思路： 

听卿学姐讲高斯消元

这道题很卡精度（虽然我没做(._.)），eps都已经到了1e-14，还用到了long double ,但是最后用%.12lf输出的时候一定要强制转成double，否则结果都是0.0000000

感觉思路我有点想不到：

dp[i]表示从i到n期望的回合数

  （1）

移项有：

（2）

（3）

这样的话i从1到n可以列n个等式，最后一个dp[n]=0;(n到n不需要经过任何回合)，高斯消元之后再回代，用a[i][j]表示第i个方程第j个位置的系数

因为有传送门i->j，所以（1）式变成了dp[i]=dp[j],dp[i]-dp[j]=0,因为我们令每一个dp[i]前面的系数为6，所以dp[j]前面的系数为-6，且右边的常系数为0

如果i+x（掷出的骰子数）>n,那么就没有dp[i+x]这一项了，dp[i]前面的系数要-1，同时右侧常数不变。因为本来（2）式就是从1式导出的，且对于越界的i+x,也需要一个回合，所以右侧的常数是不变的，但是因为dp[i+x]不存在了，所以dp[i]也要少一个，这点从（3）式可以推出。

这题用读入挂读的话第二个样例本地都跑不对，还是scanf吧(´･Д･)」很迷

ac代码： 

感觉网上的很多题解博客中的代码也是用的卿学姐的代码，但是他代码中的那个p变量我觉得可以没有吧，感觉作用不是很大。因为目前还无法评测，所以代码的正确性不完全保证，但是大概思路就是酱紫。

#include<bits/stdc++.h>
using namespace std;
const int maxn=303;
const long double eps=1e-14;
long double a[maxn][maxn];//a[i][j]表示第i个方程组的第j个系数的多少
int n,m,u,v,f[maxn];//i->f[i]：从i可以传送到f[i]
int main()
{
    scanf("%d %d",&n,&m);
    for(int i=1;i<=n;i++)
        f[i]=i;
    for(int i=1;i<=m;i++)
    {
        scanf("%d %d",&u,&v);
        f[u] = v;//有传送
    }
    //建立增广矩阵
    for(int i=1;i<n;i++)
    {
        a[i][i]=6.0;//6dp[i]-dp[i+1]-dp[i+2]+...+dp[i+6]=6
        if(f[i]!=i)
        {
            a[i][f[i]]=-6.0;//6dp[u]=6dp[v];
        }
        else
        {
            a[i][n+1]=6.0;//方程右边的常数
            for(int j=1;j<=6;j++)
            {
                if(i+j<=n)
                    a[i][i+j]=-1.0;
                else
                    a[i][i]-=1.0;//i+j越界
            }
        }
    }
    a[n][n]=1.0;
    a[n][n+1]=0.0;
    //高斯消元的过程
    for(int i=1;i<=n;i++)
    {
             if(fabs(a[i][i])>eps)
             {
                 for (int j = i + 1; j <= n; j++)//把下面的所有第i位的系数都调0
                 {
                     if (fabs(a[j][i]) > eps) {
                         long double k = a[j][i] / a[i][i];//消元
                         for (int t = i; t <= n + 1; t++)
                             a[j][t] -= a[i][t] * k;
                     }
                 }
             }
    }
    //回代过程
    for(int i=n;i>=1;i--)
    {
        for(int j=i+1;j<=n;j++)//i+1开始
            if(fabs(a[i][j])>eps)//有系数
                a[i][n+1]-=a[i][j]*a[j][n+1];//先减去代入的值
        if(fabs(a[i][i])<=eps && fabs(a[i][n+1])>eps)//系数为0，但乘积不为0，无法求解
        {
            printf("-1\n");
            return 0;
        }
        a[i][n+1]/=a[i][i];
    }
    printf("%.12lf\n",(double)a[1][n+1]);//最后long double 要转成double 否则结果都是0
    return 0;
}