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 //提示:自动阅卷结束唯一标识,请勿删除或增加。


阅读更多
个人分类: hero.pongo
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

不良信息举报

pongo(英雄会)题解之理想随机数发生器

最多只允许输入30个字

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭