#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. int main( int argc, char **argv)

    1.参数 (有时参数是void) argc是程序运行时参数个数 argv是存储参数的数组,可以用char* argv[],也可以用char **argv. 例如编译一个hello.c的程序 1 #in ...

  2. C Primer Plus_第8章_字符输入输出和输入确认_编程练习

    1.题略 #include <stdio.h> int main(void) { ; printf("Please enter text here(end with Ctrl + ...

  3. .NET Lambda

    Lambda概述      lambda 表达式是一个可用于创建委托或表达式树类型的匿名函数. 通过使用 lambda 表达式,可以可作为参数或返回编写本地函数,该函数调用的值. Lambda 表达式 ...

  4. 51nod1079(中国剩余定理)

    题目链接: http://www.51nod.com/onlineJudge/user.html#!userId=21687 题意: 中文题诶~ 思路: 本题就是个中国剩余定理模板题,不过模拟也可以过 ...

  5. 第二章 编写与设置Servlet

    2.1 第一个Servlet package cc.openhome; import javax.servlet.ServletException; import javax.servlet.http ...

  6. GitHub使用心得

    PHP今天算是学习结束了;版本控制也了解很久了,比较熟悉的是svn,Git一直都不太熟,也可能是它是英文的缘故,感觉很费劲,所以用的不多,但不得不说,它功能真的很强大;今天试着和朋友用Git开发一个简 ...

  7. thinkphp中volist标签

    volist标签 volist标签主要用于在模板中循环输出数据集或者多维数组 volist(name,id,offset,length,key,mod,empty) name(必须):要输出的数据模型 ...

  8. 将Linux下的Android签名对pk8和pem转换为Eclipse下的签名(keystore)

    一 在github上下载工具 https://github.com/getfatday/keytool-importkeypair 二 将工具在Linux环境下解压或者解压后Copy到Linux下,运 ...

  9. WPF中如何将ListViewItem双击事件绑定到Command

    今天的博客将介绍如何实现ListViewItem双击事件绑定到ViewModel中的Command.实现方法借助了Style中的EventSetter,请看下面的详细代码: <ListView ...

  10. python安装库

    首先确保安装了pip,并且pip也加入了系统path路径: pip下载:https://pypi.python.org/pypi/pip#downloads 下载Python对应的包:(http:// ...