A - A Dangerous Maze LightOJ - 1027(概率/期望)

A - A Dangerous Maze LightOJ - 1027(概率/期望)

题意

1. 一个人被困在一个迷宫中，在他面前有 n 个门，对于第 i 个门有 ai,
2. 如果 ai>0 表示这个人消耗 ai 分钟之后，传送出迷宫，
3. 如果 ai<0 表示这个人消耗 ai 分钟之后，又被传送到这个 n 个门前，继续随机选择一个门，
4. 问逃出迷宫，所需要的时间期望是多少？

思路

1. 首先选择每个门是等概率的，
2. s1 = 所有整数传送门所用的时间之和，
3. s2= 所有负数，取绝对值之后的和，
4. 设期望为 E, 设时间大于 0 的传送门有 n1 个，小于 0 的有 n1 个（n1+n2==n）,
5. 对于每个传送门有两种情况：
   • 将人传送出去的传送门：t/n，
   • 将人重新传送到迷宫中：（t+E）/n
     传送消耗时间为 t
6. 那么有下面的等式 $$E=s1/n+(s2+n2\ast E)/n$$
7. 化简一下得到：$E=(s1+s2)/(n-n2)$

这题有个明显的特点是，让求的 E，出现的在的了在等式的右边中，让求解的过程类似是一个 递归意味的过程。

代码

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <math.h>
#include <algorithm>
#include <string>
#include <queue>
#include <map>
#include <bitset>
#include <vector>
using namespace std;
void fre() { system("clear"), freopen("A.txt", "r", stdin); freopen("Ans.txt","w",stdout); }
void Fre() { system("clear"), freopen("A.txt", "r", stdin);}
void Run(int x = 0) {     
#ifdef ACM  //宏定义免注释 freopen
    if(! x) fre(); else Fre();
#endif
}
#define ios ios::sync_with_stdio(false)
#define Pi acos(-1)
#define pb push_back
#define fi first
#define se second
#define db double
#define ll long long
#define ull unsigned long long
#define Pir pair<ll, ll>
#define m_p make_pair
#define for_(i, s, e) for(ll i = (ll)(s); i <= (ll)(e); i ++)
#define rep_(i, e, s) for(ll i = (ll)(e); i >= (ll)(s); i --)
#define memset(a, b, c) memset(a, (int)b, c);
#define size() size() * 1LL
#define sc scanf
#define pr printf
#define sd(a) sc("%lld", &a)
#define ss(a) sc("%s", a)
#define __  pr( "------------------------\n" );
#define ___ pr("\n------------------------\n");
#define inf 0x3f3f3f3f
#define INF 0x3f3f3f3f3f3f3f3f
#define esp 1e-7
#define mod (ll)(1e9 + 7)
/*=========================ACMer===========================*/
ll gcd(ll a, ll b) { return b ? gcd(b, a % b) : a; }

int main()
{
    Run();
    ll T, cas = 1; sd(T);
    while(T --)
    {
        ll n; sd(n);
        ll s1 = 0, s2 = 0, t;
        ll c1 = 0, c2 = 0;
        for_(i, 1, n)
        {
            sd(t);
            if(t > 0)
                s1 += t, c1 ++;
            else
                s2 += -t, c2 ++;
        }

        ll x = (s1 + s2);
        ll y = c1;
        ll d = gcd(x, y);

        if(! c1)
        {
            pr("Case %lld: inf\n", cas ++);
        }
        else
            pr("Case %lld: %lld/%lld\n", cas ++, x / d, y / d);
    }


    return 0;
}