比较Sherwood算法与确定性算法

实验题目：写一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;
	
}