快捷導(dǎo)航

概率的問題:使用遞歸與多次試驗(yàn)?zāi)M的分析

更新時(shí)間：2013年05月08日 15:49:55 作者：

以下對(duì)概率的問題:使用了遞歸和多次試驗(yàn)?zāi)M。需要的朋友參考下

多次枚舉：

實(shí)例1

口袋中有5只紅球，4只白球。隨機(jī)從口袋中取出3個(gè)球，取出1個(gè)紅球2個(gè)白球的概率

<SPAN style="FONT-SIZE: 18px"> srand( (unsigned)time( NULL ) );
int n = 0;
for(int i=0; i<100000; i++)
{
  char x[] = {1, 1, 1, 1, 1, 2, 2, 2, 2};//5個(gè)紅球用5個(gè)1表示 4個(gè)白球用4個(gè)2表示
  int a = 0; // 取到的紅球的數(shù)目
  int b = 0; // 取到的白球的數(shù)目

  for(int j=0; j<3; j++) //取3個(gè)球進(jìn)行3次循環(huán)
  {
   int k = rand() % (9-j); //下標(biāo)的確定確定范圍 9-j 是重點(diǎn)
   if(x[k]==1)
    a++;
   else
    b++;

   x[k] = x[9-j-1]; //將取出數(shù)向后移動(dòng)
  }
  if(a==1 && b==2) n++;//取出1個(gè)紅球2個(gè)白球時(shí)進(jìn)行計(jì)數(shù)
}
printf("概率=%f\n", n/100000.0*100);</SPAN>

實(shí)例2

復(fù)制代碼代碼如下:

<SPAN style="FONT-SIZE: 18px">#define N 30
......
 int a[N];
 srand( time( NULL ) );
 int n = 0;
 for(int k=0; k<10000; k++)
 {
  for(int i=0; i<N; i++)
   a[i] = rand() % 365;
  bool tag = false; // 假設(shè)沒有相同
  for(i=1; i<N; i++)
  {
   for(int j=0; j<i; j++)
   {
    if(a[i]==a[j]) 
    {
     tag = true;
     break;
    }
   }
   if(tag) break;
  }
  if(tag) n++;
 }
 printf("%f\n", 1.0 * n / 10000 * 100);
</SPAN>

遞歸：

某個(gè)袋子中有紅球m個(gè)，白球n個(gè)?，F(xiàn)在要從中取出x個(gè)球。紅球數(shù)目多于白球的概率

下面的代碼解決了這個(gè)問題。其中的y表示紅球至少出現(xiàn)的次數(shù)。

這與前文的問題是等價(jià)的。因?yàn)槿绻?0個(gè)球，要求紅球數(shù)大于白球數(shù)，則等價(jià)于至少取出16個(gè)紅球。

復(fù)制代碼代碼如下:

<SPAN style="FONT-SIZE: 18px">/*
   m: 袋中紅球的數(shù)目
   n: 袋中白球的數(shù)目
   x: 需要取出的數(shù)目
   y: 紅球至少出現(xiàn)的次數(shù)
*/

double pro(int m, int n, int x, int y)
{
if(y>x) return 0;
if(y==0) return 1; //對(duì)y沒有要求
if(y>m) return 0;
if(x-n>y) return 1; //把白球全部取出，剩下就是紅球紅球比至少取出還多，概率為1
double p1 = pro(m-1,n,x-1,y-1) ;
double p2 = pro(m,n-1,x-1,y);
return (double)m/(m+n) * p1 + (double)n/(m+n) * p2;
}</SPAN>