实验题目:写一Sherwood算法C,与算法A, B, D比较,给出实验结果
算法的思想很简单,因为经过计算,算法B是从val前个数中找到一个不大于x的数y,然后从y开始寻找,直到找到x返回x的下标。那么作为一个概率算法,算法C采用与B相似的思维,只不过算法C不是直接从val前个数中找,而是随机产生个数,接下来做与算法B相同工作,就是找这数中不大于x的最大的数,然后再从这个数开始search。因为val作为一个完全随机的数组,随机选取个数和选取前个数效果是一样的。
这里给出算法A、B、C、D代码。
// HomeWork_P67.cpp : 此文件包含 "main" 函数。程序执行将在此处开始并结束。
//
#include "pch.h"
#include <iostream>
#include <random>
#include <vector>
#include <ctime>
#include <algorithm>
#include<fstream>
#define MAXN 10000
using namespace std;
//产生一个随机数组val
vector<int> getVector() {
vector<int> val(1000);
for (auto i = 0; i < val.size(); ++i)
val[i] = i;
uniform_int_distribution<unsigned> u(0, 999);
default_random_engine e(time(0));
for (auto i = 0; i < val.size(); ++i)
{
int tmpa = u(e);
int tmpb = u(e);
swap(val[tmpa], val[tmpb]);
}
return val;
}
//打印数组
void print(vector<int> v) {
for (auto i : v)
cout << i << " ";
cout << endl;
}
//生成数组ptr[]
vector<int> getPtr(vector<int> val) {
//算法思想:每次选出val的最小值,将其标记为MAXN,那么下一个最小值就是当前值的next。当前的ptr值就应该是下一个的下标
vector<int> ptr(val.size(),-1);
auto min_it = min_element(val.begin(), val.end());
int pos = min_it-val.begin();
int head = pos;
while (*min_it != MAXN) {
*min_it = MAXN;
pos = min_it - val.begin();
min_it = min_element(val.begin(), val.end());
ptr[pos] = min_it - val.begin();
}
ptr[pos] = head;
return ptr;
}
//查找x在数组val的下标.返回pair(下标,比较次数)
pair<int,int> Search(int x, int i,vector<int> val,vector<int> ptr) {
int count = 1;//统计比较次数
while (x > val[i]) {
i = ptr[i];
++count;
}
return { i,count };
}
//确定性算法A(x),在表val中查找x的下标
pair<int, int> A(int x, vector<int> val, vector<int> ptr) {
int head = min_element(val.begin(), val.end()) - val.begin();
return Search(x, head, val, ptr);
}
//概率算法D(x),在表val中查找x的下标
pair<int, int> D(int x, vector<int> val, vector<int> ptr) {
default_random_engine e(time(0));
uniform_int_distribution<int> u(0, val.size());
int head = min_element(val.begin(), val.end()) - val.begin();
int i = u(e);
int y = val[i];
if (x < y)
return Search(x, head, val, ptr);
else if (x > y)
return Search(x, ptr[i], val, ptr);
else
return { i ,0};
}
//确定性算法B(x),在表val中查找x的下标
pair<int, int> B(int x, vector<int> val, vector<int> ptr) {
int head = min_element(val.begin(), val.end()) - val.begin();
int max = val[head];
int i = head;
int y = -1;
for (int j = 0; j < sqrt(val.size()); ++j) {
y = val[j];
if (max < y&&y <= x) {
i = j;
max = y;
}
}
return Search(x, i,val,ptr);
}
//概率算法C(x),先随机取L个下标,找到一个不大于x的数,从这个数开始查找。
pair<int, int> C(int x, vector<int> val, vector<int> ptr) {
int L = sqrt(val.size());
uniform_int_distribution<int> u(0, 999);
default_random_engine e(time(0));
vector<int> xb;
for (auto i = 0; i < L; ++i)
xb.push_back(u(e));
int i = xb[0];
int y = -1;
int max = i;
for (int j = 0; j < L; ++j) {
y = val[xb[j]];
if (max < y&&y < x) {
i = xb[j];
max = y;
}
}
return Search(x, i, val, ptr);
}
//测试向文件输入数据 用于测试数组
void fwrite(vector<int> v) {
ofstream os("a.txt");
for (auto i = 0; i < v.size(); ++i) {
os << v[i] << " ";
if ((i+1) % 10 == 0)
os << endl;
}
}
//测试随机数组 从文件中读取val
vector<int> getVal() {
vector<int> val;
ifstream is("a.txt");
int tmp;
while (is >> tmp)
val.push_back(tmp);
return val;
}
//随机抽10个数字,
vector<int> randVal() {
uniform_int_distribution<int> u(0, 999);
default_random_engine e(time(0));
vector<int> rV;
for (auto i = 0; i < 100; ++i)
rV.push_back(u(e));
return rV;
}
//各算法查找10个随机数,的查找次数
int searchTime(pair<int, int>(*pf)(int , vector<int> , vector<int> ), vector<int> val,vector<int> ptr,vector<int> RV) {
int sum = 0;//统计总共比较次数
for (auto si : RV) {
pair<int, int> p = pf(si, val, ptr);
sum += p.second;
}
return 1.0*sum / RV.size();
}
int main() {
vector<int> val = getVector();
vector<int> ptr = getPtr(val);
//随机抽100个数字,测试他们用四种方法所需要查询的次数
vector<int> RV = randVal();
cout << "随机抽100个数字如下:\n";
for (auto i : RV)
cout << i << " ";
cout << endl;
//用算法
cout <<"算法A 查找这100个数字的平均查找次数 " << searchTime(A,val,ptr,RV) <<endl;
cout << "算法B 查找这100个数字的平均查找次数 " << searchTime(B, val, ptr, RV) << endl;
cout << "算法C 查找这100个数字的平均查找次数 " << searchTime(C, val, ptr, RV) << endl;
cout << "算法D 查找这100个数字的平均查找次数 " << searchTime(D, val, ptr, RV) << endl;
}