一、问题简介
这是一个类似于麦穗问题的小游戏,具体规则就是一共有100个礼物被杯子覆盖,纸杯下的奖品价值是均匀分布的,小明要从中选择价值最高礼物,这时候小明就有两种可能:
(1)拿走任意一个纸杯下的奖品;
(2)从左往右依次打开纸杯,小明可以选择要或者不要纸杯下的奖品,如果小明选择不要,那小明可以继续打开后续纸杯,但是不能再要之前放弃的物品,如果小明选择要,则游戏结束。
二、问题思路
对于第一种情况显然可能性为百分之一选到最贵重的物品;所以我们来讨论第二种情况,模拟第二种情况的过程:
设定一个值K,小明选择在前K个礼物中价值最高的一件礼物为标准A,在打开第K个礼物之后,只要选择到价值比A高的礼物就选择该礼物。然后我们通过模拟来考量当K取多少时,小明选到价值最高的礼物的概率是最大的。
三、代码实现
#include<iostream>
#include<cmath>
#include<cstdlib>
#include<ctime>
#include "matplotlibcpp.h"
#define N 100
#define MAX 100000
namespace plt = matplotlibcpp;
using namespace std;
//随机打乱权重数组
void random(int w[])
{
int index, tmp, i;
srand(time(NULL));//随机数种子为time
for (i = 0; i < N; i++)
{
index = rand() % (N - i) + i;
if (index != i)
{
tmp = w[i];
w[i] = w[index];
w[index] = tmp;
}
}
}
int main()
{
//开辟两个动态数组分别记录权重与k次取得最大值的次数
int* w = new int[N];
int* k_times = new int[N + 1];
for (int i = 0; i < N; i++) {
w[i] = i;
k_times[i] = 0;
}
k_times[N] = 0;
for (int s = 0; s < MAX; s++)
{
//打乱w权重的次序
random(w);
//进行抽奖模拟
for (int k = 1; k < 101; k++)
{
int flag = 0;
//设置最大值
int max = 0;
//遍历K个礼物,从中选出价值最大的礼物
for (int stop = 0; stop < k; stop++)
{
if (w[stop] > max) max = w[stop];
}
//遍历后100-K个礼物,选出第一个价值大于max的礼物
for (int stop = k; stop < N; stop++)
{
if (w[stop] > max)
{
flag = 1;
max = w[stop];
break;
}
}
//当flag == 0,没有找到大于max的数,只能选择最后那个
if (flag == 0)max = w[N - 1];
//当k==100时,没得选择,只能选最后那个
if (k == 100) max = w[N - 1];
//如果选出的价值最大,则记录K的值+1
if (max == 99)
{
k_times[k]++;
}
}
}
for (int i = 0; i < N + 1; i++)
{
cout << k_times[i] << " ";
}
//对绘制K等于1到99的概率曲线
std::vector<int> x(N), y(N);
//给横坐标与纵坐标赋值
for (int i = 0; i < N; i++)
{
x.at(i) = i + 1;
y.at(i) = k_times[i + 1];
}
//开始绘图
plt::figure_size(1200, 780);
plt::plot(x, y);
plt::xlim(0, 100);
plt::title("K_times");
plt::show();
delete[]w;
delete[]k_times;
return 0;
}
四、运行结果
在这里我们运用到了C++的一个matplotlib-cpp库,这个需要单独配置运行环境。通过每次如果取到了最大值时候K的次数的分布,可以为小明抽取礼物的时设置的K值作为一个参考。
五、注意
这份代码因为没有结果对照,所以无法确认代码正确性,如果代码中有错误的地方还希望读者指出,这也是我发这篇博客的初心所在,我也会在第一时间更正代码。
github仓库地址:GitHub - beihaidaoxingwan/AI_project