【高斯消元求期望】LightOJ - 1151 Snakes and Ladders

Step1 Problem:

有 100 个格子,从 1 开始走,每次抛骰子走 1~6,如果所走位置大于 100,则重新抛骰子。i 会被传送至 nex[i]。
问你从 1 走到 100 所抛骰子的期望次数。

Step2 Ideas:

由期望公式可得: d p [ i ] = ( d p [ i + 1 ] + d p [ i + 2 ] + d p [ i + 3 ] + d p [ i + 4 ] + d p [ i + 5 ] + d p [ i + 6 ] ) 6 + 1 dp[i] = \frac{(dp[i+1] + dp[i+2] + dp[i+3] + dp[i+4] + dp[i+5] + dp[i+6])}{6} + 1 dp[i]=6(dp[i+1]+dp[i+2]+dp[i+3]+dp[i+4]+dp[i+5]+dp[i+6])+1.
由上述式子可以看出,得先求出后续状态才能推出 dp[i],所以当后续状态由前面状态推出来的时候也就是 i 会被传送至 nex[i],就很明显的需要列方程了。
第二个方程: d p [ n e x t [ i ] ] = d p [ i ] dp[next[i]] = dp[i] dp[next[i]]=dp[i].
然后高斯消元解方程即可。

Step3 Code:

#include<bits/stdc++.h>
using namespace std;
const int N = 110;
double data[N][N];
double ans[N];
double eps = 1e-8;
int lib[N], nex[N], n;
bool eq(double x, double y)
{
    if(fabs(x-y) <= eps) return 1;
    else return 0;
}
int Gauss()
{
    int rst = 0, r = 0;
    memset(lib, 0, sizeof(lib));
    for(int i = 0; i < n; i++)
    {
        for(int j = r; j < n; j++)
        {
            if(eq(data[j][i], 0)) continue;
            for(int k = 0; k <= n; k++)
                swap(data[j][k], data[r][k]);
            break;
        }
        if(eq(data[r][i], 0)) continue;
        for(int j = 0; j < n; j++)//利用r消去其他行,和01矩阵不同的地方
        {
            if(j == r || eq(data[j][i], 0)) continue;
            double t = data[r][i]/data[j][i];
            for(int k = 0; k <= n; k++)
                data[j][k] = data[j][k]*t - data[r][k];
        }
        lib[i] = 1; r++; rst++;
    }
    for(int i = 0; i < n; i++)
    {
        if(eq(data[i][n], 0) == 0) {
            int flag = 0;
            for(int j = 0; j < n; j++)
            {
                if(eq(data[i][j], 0) == 0) {
                    flag = 1;
                    break;
                }
            }
            if(!flag) {
                return -1;
            }
        }
        if(!lib[i]) continue;
        if(eq(data[i][i], 0)) continue;
        ans[i] = data[i][n] / data[i][i];
    }
    return rst;
}
int main()
{
    int T, m, u, v;
    scanf("%d", &T);
    for(int Case = 1; Case <= T; Case++)
    {
        scanf("%d", &m);
        n = 100;
        memset(nex, -1, sizeof(nex));
        for(int i = 1; i <= m; i++) {
            scanf("%d %d", &u, &v);
            nex[u-1] = v-1;
        }
        memset(data, 0, sizeof(data));
        // 100个解,需要100个方程组
        data[n-1][n-1] = 1, data[n-1][n] = 0;//你在 100 这个格子,期望步数就是 0 了.
        for(int i = 0; i < n-1; i++) {
            if(nex[i] != -1) {
                data[i][i] = 1;
                data[i][nex[i]] = -1;
                data[i][n] = 0;
            }
            else {
                int cnt = 0;
                for(int j = 1; j <= 6; j++) {
                    if(i+j < n) {
                        data[i][i+j] = -1;
                        cnt++;
                    }
                }
                data[i][i] = cnt;
                data[i][n] = 6;
            }
        }
        Gauss();
        printf("Case %d: %lf\n", Case, ans[0]);
    }
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值