免费注册 查看新帖 |

Chinaunix

  平台 论坛 博客 文库
12345下一页
最近访问板块 发新帖
查看: 13114 | 回复: 45
打印 上一主题 下一主题

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

论坛徽章:
0
跳转到指定楼层
1 [收藏(0)] [报告]
发表于 2012-10-20 22:23 |只看该作者 |正序浏览
    前几天去某公司笔试,有道题说的是有个函数foo,返回0的概率是60%,返回1的概率是40%,让你自己写一个函数,使返回0和1的概率是50%,

用前面说的foo函数实现,不能用像c++里面的rand这类的函数,怎么解?

论坛徽章:
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 的随机数生成器

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

回复 1# 发仔很忙


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

论坛徽章:
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
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. }
复制代码

论坛徽章:
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
40 [报告]
发表于 2012-10-29 20:29 来自手机 |只看该作者
感觉求(6/10)∧x+(4/10)∧y=1/2

论坛徽章:
0
39 [报告]
发表于 2012-10-29 13:14 |只看该作者
我有一个想法,不知道对不对。我是这样想的:假设我们要生成N个随机数,先用foo()生成并保存在a[N]中,那么数组a中0占60%,1占40%(N足够大),然后把60%中的10%由0变成1,这样0和1就各占一半。

论坛徽章:
0
38 [报告]
发表于 2012-10-25 17:57 |只看该作者
回复 35# bruceteen


    啊 .. 这里的意思是 2n 次调用(遇01/10则停)以内获得 01 的概率,有限次内确实是不到 0.5 的
    但如果不限次数,则是等于 0.5 的。
    而获得非 01/10 的概率随调用次数是呈指数下降的
    有个形象的比喻,这种指数下降的概率两位数次方以后,是小于宇宙射线改变改变内存存储状态的概率的,因而可以忽略。
    所谓“永远不是” 也只是理论上的可能而已。
  

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

清除 Cookies - ChinaUnix - Archiver - WAP - TOP