蒙提霍尔问题（三门问题）程序

http://baike.baidu.com/view/657835.htm

 

蒙提霍尔问题：参赛者会看见三扇关闭了的门，其中一扇的后面有一辆汽车，选中后面有车的那扇门就可以赢得该汽车，而另外两扇门后面则各藏有一只山羊。当参赛者选定了一扇门，但未去开启它的时候，节目主持人会开启剩下两扇门的其中一扇，露出其中一只山羊。主持人其后会问参赛者要不要换另一扇仍然关上的门。问题是：换另一扇门是否会增加参赛者赢得汽车的机率？

 

程序：

 

static int n[3];//统计产生的随即数0、1、2的次数 //产生随机数0、1或2 //方法：先生成一个随机数，把这个数的各个10进制位加起来，把这个和除以3的余数返回。 int GetRandom() { int nRandom = rand(); int nDigitalSum = 0; int nInt = nRandom; while( nInt > 0 ) { int nDec = nInt % 10; nInt /= 10; nDigitalSum += nDec; } int nRan = nDigitalSum % 3; ++ n[nRan]; return nRan; } //一次试验（一次猜测），bChange表示主持人打开一个门后，参与者是否改变选择 //返回值表示参与者是否选中了奖品 bool guess( bool bChange ) { int nSet = GetRandom();//奖品在第nSet个门里面 int nGuess = GetRandom();//参与者猜测奖品在第nGuess个门里面 int nPublic = nGuess;//主持人会打开第nPublic个门（nPublic与nSet不同，与nGuess也不同） while( nPublic == nSet || nPublic == nGuess ) nPublic = GetRandom(); if( bChange ) { //参与者改变选择，他只能选择nGuess和nPublic以外的那个门 vector<int> vecPos( 3 ); vecPos[0] = 0; vecPos[1] = 1; vecPos[2] = 2; int nMax = max( nGuess, nPublic ); int nMin = min( nGuess, nPublic ); vecPos.erase( vecPos.begin() + nMax ); vecPos.erase( vecPos.begin() + nMin ); nGuess = vecPos[0]; } return nSet == nGuess; } //做nTimes次试验（nTimes次猜测） int guess( int nTimes, bool bChange ) { int nRight = 0; for( int i = 0; i < nTimes; ++ i ) { bool bRight = guess( bChange ); if( bRight ) ++ nRight; } return nRight; } void main() { int nTimes = 1000;//试验次数（猜测次数） srand( (unsigned)time( NULL ) ); bool bChange = false; n[0] = n[1] = n[2] = 0; int nRight = guess( nTimes, bChange ); TRACE( "随机数: 0:%d次, 1:%d次, 2:%d次/n", n[0], n[1], n[2] ); TRACE( "不改变选择/n共%d次，猜中%d次/n/n", nTimes, nRight ); bChange = true; n[0] = n[1] = n[2] = 0; nRight = guess( nTimes, bChange ); TRACE( "随机数: 0:%d次, 1:%d次, 2:%d次/n", n[0], n[1], n[2] ); TRACE( "改变选择/n共%d次，猜中%d次/n", nTimes, nRight ); }

 

执行输出：

随机数: 0:1519次, 1:1458次, 2:1469次
不改变选择
共1000次，猜中331次

随机数: 0:1553次, 1:1508次, 2:1461次
改变选择
共1000次，猜中669次