pongo（英雄会）题解之理想随机数发生器

点击打开转载地址

不得感慨下本科没好好学习概率论，这道题用概率论很简单。

首先还是题目：

题目详情:

一个理想随机数发生器R，给定参数正整数x，它可以均匀随机产生一个闭区间[0,x]之间的实数（注意是实数,每个实数出现的概率相同）R(x)。现在给定3个正整数a,b,c，我们使用a,b产生两个随机实数R(a), R(b)，问R(a) + R(b) <= c的概率有多大？

输入: 正整数a,b,c。正整数均不超过10000

输出: R(a) + R(b) <= c的概率。注意输出是一个有理数（分数），请化为最简分数。（分子分母的最大公约数是1）。

例如 输入a = 1, b = 1, c = 4，输出"1/1"

这道题转换为概率论题目就是：x与y分别在[0,a]和[0,b]上服从均匀分布，并且x与y相互独立，求x+y<=z ，当z=c时的概率。

x与y分别在[0,a]和[0,b]上的概率密度为f(x)=1/a,f(y)=1/b；x与y相互独立，那么联合概率密度f(x,y)=1/(a*b)  x属于[0,a]、y属于[0,b]，

....

....

用图示如下（假如a>=b)：

就是求阴影部分的面积与矩形面积的比例，分c<=b、b<c<=a、a<c<a+b和c>=a+b四种情况求，记得最后需要约去公约数，即结果中的分子与分母的最大公约数为1，代码如下：

#include <stdio.h>  
#include <iostream>  
#include <string>  
#include<sstream>  
using namespace std;  
class Test {  
public:  
    static int maxG(int a,int b){  
        while(b){  
            int mid=b;  
            b=a%b;  
            a=mid;  
        }  
        return a;  
    }  
    static string calculate (int   a,int   b,int   c){  
        if(c>=a+b)return "1/1";  
        if(a<b){//for a>=b  
                int mid =a;  
                a=b;  
                b=mid;  
        }     
        int res=0;  
        int sqr=2*a*b;  
        if(c<=b){  
            res=c*c;  
        }else if(c>b&&c<=a){  
            //res=c*c-(c-b)*(c-b);  
            res=b*(2*c-b);  
        }else {//a<c<a+b  
            //res=2*a*b-(a+b-c)*(a+b-c);  
            res=2*c*(a+b)-a*a-b*b-c*c;  
        }  
        int maxGYS=1;  
        while((maxGYS=maxG(sqr,res))>1){  
            res/=maxGYS;  
            sqr/=maxGYS;  
        }  
        std::stringstream ss;  
        ss<<res<<"/"<<sqr;  
        std::string ret=ss.str();     
        return ret;  
    }  
};  
//start 提示：自动阅卷起始唯一标识，请勿删除或增加。  
int main(){     
    cout<<Test::calculate(1,3,2)<<endl;     
}   
//end //提示：自动阅卷结束唯一标识，请勿删除或增加。