c++遗传算法，用bitset实现

#include 
#include 
#include 
#include 
#include 
#include 
using namespace std;

#define SIZE 20//染色体长度，根据实际情况改写

typedef bitset Chromosome;
typedef vector Population;
typedef vector MatingPool;


void cross（MatingPool& pool double crossRate）
{
	default_random_engine e（time（0））;
	uniform_real_distribution rand_r（0 1）;
	uniform_int_distribution rand_p（1 SIZE - 1）;
	for （int i = 0; i < pool.size（）; i += 2）
	{
		if （rand_r（e） <= crossRate）
		{
			int pos = rand_p（e）;
			Chromosome temp = pool[i] >> pos << pos | pool[i + 1] << （SIZE - pos） >> （SIZE - pos）;
			pool[i + 1] = pool[i + 1] >> pos << pos | pool[i] << （SIZE - pos） >> （SIZE - pos）;
			pool[i] = temp;
		}
	}
}


void mutate（Population& population double mutationRate）
{
	default_random_engine e（time（0））;
	uniform_real_distribution rand_r（0 1）;
	uniform_int_distribution rand_p（1 SIZE - 1）;
	for （int i = 0; i < population.size（）; i++）
	{
		if （rand_r（e） <= mutationRate）
		{
			population[i].flip（rand_p（e））;
		}
	}
}

double parse（Chromosome chromo）
{
	double value = chromo.to_ullong（）;
	value /= （1 << SIZE） - 1;
	value *= 3.141592654;
	return value;
}

//价值函数，根据实际情况改写
double value_of（Chromosome chromo）
{
	double value = cos（parse（chromo））;
	return 1 / value / value;
}


void Scale（vector& values）
{
	vector temp;
	for （int i = 0; i < values.size（）; i++）
	{
		int rank = values.size（）;
		for （int j = 0; j < values.size（）; j++）
		{
			if （i == j） continue;
			if （values[i] < values[j]） rank--;
		}
		temp.push_back（rank）;
	}
	values.swap（temp）;
}


void Init（Population& population int popSize）
{
	popSize = popSize >> 1 << 1;
	default_random_engine e（time（0））;
	uniform_int_distribution b（0 1）;
	for （int i = 0; i < popSize; i++）
	{
		Chromosome chromo;
		for （int i = 0; i < SIZE; i++）
		{
			chromo.set（i b（e））;
		}
		population.push_back（chromo）;
	}
}


void selection（Population population vector values MatingPool& pool）
{
	Scale（values）;
	int size = population.size（）;
	for （int i = 1; i < size; i++） values[i] += values[i - 1];
	default_random_engine e（time（0））;
	uniform_real_distribution rand（0 values[size - 1]）;
	for （int i = 0; i < siz