免费注册 查看新帖 |

Chinaunix

  平台 论坛 博客 文库
12345
最近访问板块 发新帖
楼主: 发仔很忙
打印 上一主题 下一主题

[算法] 如何实现返回0和1的概率为50% [复制链接]

论坛徽章:
1
2015年迎新春徽章
日期:2015-03-04 09:49:45
41 [报告]
发表于 2012-10-30 00:14 |只看该作者

  1. int func() {
  2.     return foo() * 2 + 1 - foo() <= 1
  3. }
复制代码
目测俺的方法最简单{:3_189:}

论坛徽章:
0
42 [报告]
发表于 2012-10-30 06:07 |只看该作者
对于独立事件,P(AB)=P(A)P(B),所以只需调用两次foo,P(0,1)与P(1,0)的概率都是0.24,由此就可得到一个等概率。这样当(0,1)出现时返回0,当(1,0)出现时,返回1。
  1. /**
  2. * File:MyFoo.java
  3. * -----------------
  4. * P(0,1)=P(0)*P(1)=0.6*0.4=0.24
  5. * P(1,0)=P(1)*P(0)=0.4*0.6=0.24
  6. * P(0,1)=P(1,0)
  7. */
  8. package org.cudemo;

  9. /**
  10. * @author isaacxu
  11. *
  12. */
  13. public class MyFoo {

  14.        
  15.         public static void main(String[] args) {
  16.                 int val = myfoo();
  17.         System.out.println(val);
  18.         }
  19.        
  20.         // foo that returns 0 by 0.6 and 1 by 0.4
  21.         static int foo()
  22.         {
  23.                 int i=0;
  24.             // some code here
  25.                 return i;
  26.         }
  27.        
  28.         // returns  0 and 1 by 0.5 probability
  29.         static int myfoo()
  30.         {
  31.                 int returnVal= 0;
  32.             int firstCall = foo();
  33.             int secondCall = foo();
  34.             //get pair (0,1) or (1,0)
  35.             while(firstCall==secondCall){
  36.                     firstCall = foo();
  37.                     secondCall = foo();
  38.             }
  39.           
  40.             if (firstCall == 0 && secondCall == 1)
  41.                 returnVal= 0;   // get 0.24 probability
  42.             if (firstCall == 1 && secondCall == 0)
  43.                     returnVal = 1;   // get 0.24 probability
  44.             return returnVal;  
  45.         }

  46. }
复制代码

论坛徽章:
0
43 [报告]
发表于 2012-10-30 14:10 |只看该作者
本身此题目就提供了基本概率函数,没必要重新编写其他函数,不然笔试肯定不合格,这个题目更多的是考察如何构造算法:
因为 0:60% 1:40%
所以 00:36% 01:24% 10:24% 11:16%
可以看出 01 和10出现的概率是一样的
所以函数中若出现11和00则重新产生随机数
出现10或01则可以

论坛徽章:
0
44 [报告]
发表于 2012-10-30 15:09 |只看该作者
A1:A2= x:y
A1+A2 : A2+A1 =1:1

论坛徽章:
0
45 [报告]
发表于 2012-10-30 15:28 |只看该作者
本帖最后由 FaintKnowledge 于 2012-10-30 15:34 编辑

回复 1# 发仔很忙


论坛徽章:
0
46 [报告]
发表于 2012-10-30 17:23 |只看该作者
本帖最后由 fallening_cu 于 2012-10-30 17:56 编辑

其实可以这样:
  1. #include <cstdint>
  2. #include <cstddef>

  3. #include <vg.hpp>//在这里 https://github.com/fengwang/random_variate_generator

  4. int foo()
  5. {
  6.     static vg::vg<double> gen(0.0, 1.0);
  7.     if ( gen() > 0.4 ) return 1;
  8.     return 0;
  9. }

  10. int gen_50_50()
  11. {
  12.     std::uint64_t const a = 6364136223846793005ULL;
  13.     std::uint64_t const b = 1442695040888963407ULL;
  14.     std::uint64_t const mask = 256;
  15.     std::uint64_t s = 1;

  16.     for ( std::size_t i = 0; i != 64; ++i )
  17.     {
  18.         s <<= 1;
  19.         s += foo();
  20.     }

  21.     static std::uint64_t ans = s;
  22.     ans *= a;
  23.     ans += b;

  24.     return ans & mask ? 1 : 0;
  25. }
  26. //test
  27. int main()
  28. {
  29.     const std::size_t n = 8888888;

  30.     int acc = 0;

  31.     for ( std::size_t i = 0; i != n; ++i )
  32.         acc += gen_50_50();

  33.     std::cout << "accumulation of " << n << " 0-1 random number is " << acc << "\n";

  34.     return 0;
  35. }

复制代码
运行了3次,分别输出 4444442、4444443 和 4444445,应该可以接受作为一个 50-50 的随机数生成器
您需要登录后才可以回帖 登录 | 注册

本版积分规则 发表回复

  

北京盛拓优讯信息技术有限公司. 版权所有 京ICP备16024965号-6 北京市公安局海淀分局网监中心备案编号:11010802020122 niuxiaotong@pcpop.com 17352615567
未成年举报专区
中国互联网协会会员  联系我们:huangweiwei@itpub.net
感谢所有关心和支持过ChinaUnix的朋友们 转载本站内容请注明原作者名及出处

清除 Cookies - ChinaUnix - Archiver - WAP - TOP