#include?
#include?
#include?
#include?
#include?
#include?
#include?

#define?ulong?unsigned?long?	//?For?KISS?pseudo?random?numbers?generator
#define?RAND_MAX_KISS?4294967295ul?//?idem
					//?You?don‘t?need?KISS?if?you?have?a?good?rand（）?function


double?x[100][2]?v[100][2]?fitness_lbest[100]?X_lbest[100][2]?X_gbest[2]?fitness_X[100];
double?fitness_gbest?c1?c2r1r2w;
double?fitness_gbest_old?gamma;

double?f_rand（）;
ulong	rand_kiss（）;
void	seed_rand_kiss（ulong?seed）;
double	alea（double?adouble?b）;
double	phi?s?rmax;		//?Control?exploration/exploitation

void?main（）
{
	FILE???*chkdsk;
	
	int?n?nitit?i?jmax_fitness_indexrt;
	double?c_min?c_max?g_min?g_maxmax_fitnessvl?=?5;
	
	char?pr[]?=?“svmtrain?“;
	//char?fin[]?=?“?heart_scale“;
	
	//char?fin[]?=?“?train1.scale“;
	//char?fin[]?=?“?train2.scale“;
	//char?fin[]?=?“?train3.scale“;
	char?fin[]?=?“?sptest.scale“;
	//char?fin[]?=?“?spect_train.txt“;
	c_min?=?1?/?pow（25）;?c_max?=?pow（28）;
	g_min?=?1?/?pow（215）;?g_max?=?pow（23）;
	c1?=?1.5;?c2?=?1.5;w=0.5;gamma?=?0.1;
	//c1?=?2;?c2?=?2;w=0.8;gamma?=?0.1;
	srand（?（unsigned）time（?NULL?）?）;
	//seed_rand_kiss（1）;?//?Initialize?KISS
	
	fflush（stdin）;
	cout<<“Enter?the?number?of?particles?:?“;
	cin>>n;
	cout<<“Enter?the?number?of?iteration?is?:?“;
	cin>>it;
	cout<<“Enter?the?number?of?exactly?is?is:?“;
	cin>>gamma;
	//?Initialize?the?particle?positions?and?their?velocities
	for?（i=0;i	{
		//r1?=?alea（c_minc_max）;
		//r2?=?f_rand（）;
		//x[i][0]?=?alea（c_minc_max）;
		//x[i][1]?=?alea（g_ming_max）;
		r1?=?f_rand（）;
		r2?=?f_rand（）;
		x[i][0]?=?c_min?+?（c_max?-?c_min）?*?r1;
		x[i][1]?=?g_min?+?（g_max?-?g_min）?*?r2;
		v[i][0]?=?0;?
		v[i][1]?=?0;

	
	}
	int?d1d2pq;
	d1?=?0;	d2?=?1;
	
	s?=?sqrt（pow（x[d1][0]?-?x[d2][0]2）?+?pow（x[d1][1]?-?x[d2][1]2））;

	for?（p=0;?p	{
		for?（q=p+1;q		{
			if?（sqrt（pow（x[p][0]?-?x[q][0]2）?+?pow（x[p][1]?-?x[q][1]2））?>?s）
				s?=?sqrt（pow（x[p][0]?-?x[q][0]2）?+?pow（x[p][1]?-?x[q][1]2））;
		}
	}


	//?Initialize?the?global?and?local?fitness?to?the?worst?possible
	fitness_gbest?=?0;
	for?（i?=?0;?i?	{
		fitness_lbest[i]?=?0;
	}

	//?Loop?until?convergence?in?this?example?a?finite?number?of?iterations?chosen
	nit?=?0;
	do//for?（p?=?0;?p?	{
		//fitness_gbest_old?=?fitness_gbest;
		//?evaluate?the?fitness?of?each?particle
		nit?+=1;
		for?（i?=?0;?i?		{
			char???psBuffer[128];
			char?cmd[200];
			int?k;
			double?cx;
			if?（x[i][0]?>=?c_min?&&?x[i][0]?<=?c_max?&&?x[i][1]?>=?g_min?&&?x[i][1]?<=?g_max）
			{
			k?=?sprintf（cmd?“%s“?pr）;
			k?+=?sprintf（cmd?+?k?“?-c?%f“?x[i][0]）;
			k?+=?sprintf（cmd?+?k?“?-g?%f“?x[i][1]）;
			k?+=?sprintf（cmd?+?k?“?-v?%f“?vl）;
			k?+=?sprintf（cmd?+?k?“%s“?fin）;
???
			if（?（chkdsk?=?_popen（cmd?“rt“?））?==?NULL?）
			{//	exit（?1?）;