#include <stdio.h>
#include <stdint.h>
#include <ctime>
#include <vector>
#include <map>

template<typename T>
T *RandomSelect(const std::vector<T*> &itemPtrs);
uint32_t Random(uint32_t range);
uint32_t Random2(uint32_t range);

template<typename T>
T * RandomSelect(const std::vector<T*> &itemPtrs)
{
    if (itemPtrs.empty()){ return NULL; }
    uint32_t factor = ;
    std::vector<uint32_t> totalProbs;
    totalProbs.resize(itemPtrs.size());
    totalProbs[] = itemPtrs[]->Prob() * factor;
    ; i < itemPtrs.size(); ++i)
    {
        totalProbs[i] = totalProbs[i - ] + itemPtrs[i]->Prob() * factor;
    }
    uint32_t randInt = Random(totalProbs[totalProbs.size() - ]);
    ; i < totalProbs.size(); ++i)
    {
        if (randInt < totalProbs[i])
        {
            return itemPtrs[i];
        }
    }
    return NULL;
}

uint32_t Random(uint32_t range)
{
    ){ ; }
    return uint32_t((double)rand() / RAND_MAX * range) % range;
}

uint32_t Random2(uint32_t range)
{
    ){ ; }
    uint32_t rand1 = Random();
    uint32_t rand2 = Random();
    uint32_t rand3 = Random();
    uint32_t rand4 = Random();
    uint32_t num = (rand4 << ) | (rand3 << ) | (rand2 << ) | rand1;
    return uint32_t(double(num) / 0xFFFFFFFF * range) % range;
}

struct Item
{
    int32_t prob;
    Item(int32_t p) :prob(p){}
    int32_t Prob(){ return prob; }
};

int32_t main()
{
    std::map<Item*, int32_t> statistic;
    ;
    Item arrItem[] = { , , , , , , , , , ,  };
    std::vector<Item*> vecProbs;
    float sum = 0.0;
    ; j < ]); ++j)
    {
        sum += arrItem[j].Prob();
        vecProbs.push_back(&arrItem[j]);
    }
    ; i < loop; ++i)
    {
        Item *ret = RandomSelect(vecProbs);
        if (ret)
        {
            //printf("%d\n", ret->Prob());
            statistic[ret] += ;
        }
    }
    printf("statistic size = %u\n", statistic.size());
    double error = 0.0;
    for (const std::pair<Item*, int32_t> &p : statistic)
    {
        printf("%-8d %-8d %-12f %-12f\n", p.first->Prob(), p.second, (double)p.second / loop * 100.0, p.first->Prob() / sum * 100.0);
        error += fabs((double)p.second / loop * 100.0 - p.first->Prob() / sum * 100.0);
    }
    printf("error = %f\n", error);
    getchar();
    ;
}

事实证明,Random2要比Random的随机误差略小一点点,多随机几次还是有好处的.

C++ 从一组数据中按概率选择一个的更多相关文章

  1. C#使用拉依达准则(3σ准则)剔除异常数据(.Net剔除一组数据中的奇异值)

    原文:C#使用拉依达准则(3σ准则)剔除异常数据(.Net剔除一组数据中的奇异值) 1.问题的提出: 电池生产中,遇到一批电池的测量结果数据: 电压值 电池个数 电压值 电池个数 电压值 电池个数 电 ...

  2. 工作 巧遇 sql 查询 一组数据中 最新的一条

    SELECT * FROM rsl a, (SELECT CODE, max(time_key) time_key FROM rsl GROUP BY CODE ) b WHERE a. CODE = ...

  3. 获取mysql一组数据中的第N大的值

    create table tb(name varchar(10),val int,memo varchar(20)) insert into tb values('a', 2, 'a2') inser ...

  4. [SQL]一组数据中Name列相同值的最大Je与最小je的差

    declare @t table(name varchar(),qy varchar(),je int) insert into @t union all union all union all un ...

  5. c语言:从一组数据中选出可以组成三角形并且周长最长的三个数(简单)

    题目如下: 思路分析: 写出完整的程序: /* 问题描述: 有n根棍子,棍子i的长度为ai.想要从中选出3根棍子组成周长尽可能长的三角形.请输 出最大的周长,若无法组成三角形则输出0. */ #inc ...

  6. 【剑指offer】求一组数据中最小的K个数

    题目:输入n个整数,找出其中最小的K个数.例如输入4,5,1,6,2,7,3,8这8个数字,则最小的4个数字是1,2,3,4,. *知识点:Java PriorityQueue 调整新插入元素 转自h ...

  7. oracle数据中记录被另一个用户锁住

    原因:PL/SQL里面执行语句执行了很久都没有结果,于是中断执行,于是就直接在上面改字段,在点打钩(记入改变)的时候提示,记录被另一个用户锁住. 解决方法: 第一步:(只是用于查看哪些表被锁住,真正有 ...

  8. 在MySQL中,如何计算一组数据的中位数?

    要得到一组数据的中位数(例如某个地区或某家公司的收入中位数),我们首先要将这一任务细分为3个小任务: 将数据排序,并给每一行数据给出其在所有数据中的排名. 找出中位数的排名数字. 找出中间排名对应的值 ...

  9. 在django中如何从零开始搭建一个mock服务

    mock概念 mock 就是模拟接口返回的一系列数据,用自定义的数据替换接口实际需要返回的数据,通过自定义的数据来实现对下级接口模块的测试.这里分为两类测试:一类是前端对接口的mock,一类是后端单元 ...

随机推荐

  1. .Net中使用无闪刷新控件时提示框不显示

    今天做提示框的时候一直不显示,让我郁闷好久,晚上吃饭的时候问了同事一下,他给了一个思路, 他说可能是因为由于页面中的无闪刷新导致的结果:百度了一下真找到了解决方法 在页面中存在无闪刷新控件的时候提示框 ...

  2. @Html.Partial和@Html.Action区别

    1.首先看一下它们的对等关系 @Html.Partial 对应 @{Html.RenderPartial();}@Html.Action 对应 @{Html.RenderAction();} 以上相互 ...

  3. Java基础学习(一)

    常见的dos命令 盘符: 进入指定的盘符下. dir : 列出当前目录下的文件以及文件夹 md : 创建目录 rd : 删除目录    注意:rd不能删除非空的文件夹,而且只能用于删除文件夹. cd ...

  4. Power BI for Office 365(四)Power View第一部分

    Power View是一种可以创建动态的交互式报表的工具,并且支持丰富多样的图表类型,在Power View中创建的报表可以很容易地进行分享,以及供查看报表的用户对Power View报表进行交互.从 ...

  5. flume到kafka和hbase配置

    # Flume test file# Listens via Avro RPC on port 41414 and dumps data received to the logagent.channe ...

  6. bin/sh^M: bad interpreter: No such file or directory解决

    问题:bin/sh^M: bad interpreter: No such file or directory 原因:.sh脚本在windows系统下用记事本文件编写的.不同系统的编码格式引起的. 解 ...

  7. 【CentOS】又是一篇Shell

    一.Shell 1.Shell脚本的格式 #!/bin/bash 或者是 #!/bin/sh 开头 ,这是用于解析脚本的工具 2.执行脚本的方法 (1)bash filename 或者是sh file ...

  8. 【前端开发】优化代码之减少引入,css预编译语言的优点,stylus的使用

    前言:我必须得承认在最最最开始的时候,我对于css的预编译是非常不以为然的,这是错误的.一般在页面编写过程中,我会将需要reset的css放在reset.css中,讲会需要重复用到的放置到public ...

  9. iOS @@创建NSURL的字面量

    @@ 是创建 NSURL 的字面量的绝佳方法(例如:@@"http://example.com")

  10. Node.js-视图引擎【1】-Swig集成express的安装与配置

    node.js视图引擎,选来选去发现Swig最符合我的胃口哈哈. 一.安装Swig视图引擎 npm install -g swig 二.在node.js代码中配置如下 var app = requir ...