// IGA.cpp: implementation of the CIGA class.
//
//////////////////////////////////////////////////////////////////////

#include "stdafx.h"
#include "ImmueGA.h"
#include "IGA.h"
#include "math.h"

#ifdef _DEBUG
#undef THIS_FILE
static char THIS_FILE[]=__FILE__;
#define new DEBUG_NEW
#endif

//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////

CIGA::CIGA()
{	srand((unsigned)time(NULL)); //產生種子
	for(int i=0;i<SUM;i++)
		for(int j=0;j<VAR;j++)
			AntiBase[i].AntiGene[j]= AverageRandom(TL,TH);	 //在[a,b]之間
}

CIGA::~CIGA()
{
}

double CIGA::DoIGA(CAntiCell *ac)
{	double oldfit=0;
	srand((unsigned)time(NULL)); //產生種子
	for(int i=0;i<LOOPNUM;i++)
	{	ResortAntiBase();
		if(i==0)	InitMemory();	//i>0后才更新記憶庫
		if(i>0)	//第1次不用添加
		{	AddMemory();
			GenerateNew();
			AddNew();
		}
		if(i>0)		UpdateMemory();
		CalAntiC();
/*		if(fabs(oldfit-GetMemoryFit())<delt)//親和度趨于穩定時結束
		{	ResortMB();
			*ac=MB[0];
			return (10/ac->Fitness-LB);
		}
		oldfit=GetMemoryFit();
*/		Updatesp();	//抗體激勵和抑制
		Genetic();	//進行遺傳操作
	}
	ResortMB();	//選擇記憶庫中的最優抗體
	*ac=MB[0];
	return (10/ac->Fitness-LB);//
}

/*待函數是y=min fun(x0,x1,x2,x3)=
(0.7*x0*(2*x1-1)*x3+1.2*x2)/(1+x0*x0+x1*x1)+sin(0.5*x0*x1)
親和度*/
void CIGA::fun(CAntiCell *ac)//保守的下界估計
{	ac->Fitness=10/(LB+(0.7*ac->AntiGene[0]*(2*ac->AntiGene[1]-1)*ac->AntiGene[3]+
				 1.2*ac->AntiGene[2])/(1+ac->AntiGene[0]*ac->AntiGene[0]+
				 ac->AntiGene[1]*ac->AntiGene[1])
				 +sin(0.5*ac->AntiGene[0]*ac->AntiGene[1]));
}

void CIGA::GenerateNew()//
{	srand((unsigned)time(NULL)); //產生種子
	for(int i=0;i<NEWSUM;i++)
		for(int j=0;j<VAR;j++)
			NB[i].AntiGene[j]= AverageRandom(TL,TH);	 //在[a,b]之間
}

double CIGA::AverageRandom(double min, double max)
{	int minInteger = (int)(min*10000);
	int maxInteger = (int)(max*10000);
	int randInteger = rand()*rand();
	int diffInteger = maxInteger - minInteger;
	int resultInteger = randInteger % diffInteger + minInteger;
	return resultInteger/10000.0;
}
void CIGA::AddMemory()//添加到總細胞
{	int i,j,k;
	ResortAntiBase();
	for(i=0;i<MSUM;i++)//MB[i]
		for(j=0;j<2*MSUM;j++)//AntiBase[j]
		{	if(fabs(MB[i].Fitness-AntiBase[j].Fitness)<0.1*delt)
				break;
			if(MB[i].Fitness>AntiBase[j].Fitness)	//MB[i]加入到總免疫細胞庫
			{	for(k=SUM-1;k>j;k--)				//右移
					AntiBase[k]=AntiBase[k-1];
				AntiBase[j]=MB[i];
				break;
			}
		}
}

void CIGA::AddNew()	//添加到總細胞末尾
{	ResortAntiBase();
	for(int i=0;i<NEWSUM;i++)	
		AntiBase[SUM-NEWSUM+i]=NB[i];
}
void CIGA::ResortAntiBase()
{	int i,j;
	CAntiCell tmp;
	for(i=0;i<SUM;i++)
		fun(AntiBase+i);
	for(i=0;i<SUM;i++)
		for(j=i;j<SUM;j++)//按親和度從大到小排序
			if(AntiBase[i].Fitness<AntiBase[j].Fitness)//swap
			{	tmp=AntiBase[i];
				AntiBase[i]=AntiBase[j];
				AntiBase[j]=tmp;
			}
}

void CIGA::ResortMB()
{	int i,j;
	CAntiCell tmp;
	for(i=0;i<MSUM;i++)
		fun(MB+i);
	for(i=0;i<MSUM;i++)
		for(j=i;j<MSUM;j++)
			if(MB[i].Fitness<MB[j].Fitness)
			{	tmp=MB[i];
				MB[i]=MB[j];
				MB[j]=tmp;
			}
}

void CIGA::UpdateMemory()//更新記憶庫
{	int i,j,k;
	ResortAntiBase();
	ResortMB();
	for(i=0;i<MSUM;i++)//AntiBase[i]
		for(j=0;j<MSUM;j++)//MB[j]
		{	if(fabs(AntiBase[i].Fitness-MB[j].Fitness)<0.1*delt)
				break;
			if(AntiBase[i].Fitness>MB[j].Fitness)//AntiBase[i]加入記憶庫
			{	for(k=MSUM-1;k>j;k--)//右移
					MB[k]=MB[k-1];
				MB[j]=AntiBase[i];
				break;
			}
		}
}

void CIGA::CalAntiC()//計算濃度
{	int i,j,k;
	int genenum,num;
	for(i=0;i<SUM;i++)	//計算總免疫細胞庫的細胞濃度
	{	num=0;
		for(j=0;j<SUM;j++)
		{	genenum=0;
			for(k=0;k<VAR;k++)
				if(fabs(AntiBase[i].AntiGene[k]-AntiBase[j].AntiGene[k])<10*delt)
					genenum++;
			if(genenum>=VAR-1)	num++;//相似
		}
		AntiBase[i].Concentration=num*10.0/SUM;
	}
}

void CIGA::CalMemoryC()
{	int i,j,k;
	int genenum,num;
	for(i=0;i<MSUM;i++)	//計算記憶免疫細胞庫的細胞濃度
	{	num=0;
		for(j=0;j<MSUM;j++)
		{	genenum=0;
			for(k=0;k<VAR;k++)
				if(fabs(MB[i].AntiGene[k]-MB[j].AntiGene[k])<5*delt)
					genenum++;
			if(genenum>=VAR-1)	num++;//相似
		}
		MB[i].Concentration=num*1.0/SUM;
	}
}

double CIGA::GetMemoryFit()
{	double sumfit=0;
	for(int j=0;j<MSUM;j++)
		sumfit+=MB[j].Fitness;
	return sumfit;
}

void CIGA::Updatesp()//抗體的激勵和抑制
{	int i,j;
	double fave,cfave,sumc=0,sumfit=0;
	CalAntiC();
	ResortAntiBase();
	for(j=0;j<SUM;j++)
	{	sumfit+=AntiBase[j].Fitness;
		sumc+=AntiBase[j].Concentration;
	}
	fave=sumfit/SUM;	//計算總免疫細胞庫的平均親和度
	cfave=sumc/SUM;		//計算總免疫細胞庫的平均細胞濃度
	for(i=0;i<SUM;i++)	//用濃度調整選擇概率
		if(AntiBase[i].Fitness>fave)
			AntiBase[i].sp=(1+beita*(cfave-AntiBase[i].Concentration)
			/(cfave+AntiBase[i].Concentration))*AntiBase[i].Fitness;
		else 
			AntiBase[i].sp=AntiBase[i].Fitness;
}

void CIGA::Genetic()//交叉和變異
{	int i,j;
	int i1,i2;
	CAntiCell tmpAnti[SUM];
	double psum=0;
	for(j=0;j<SUM;j++)	psum+=AntiBase[j].sp;
	for(i=0;i<SUM;i++,i++)
	{	SelectAnti(psum,i1,i2);
		CrossOver(i1,i2,tmpAnti+i,tmpAnti+i+1);
		Mutation(tmpAnti+i,tmpAnti+i+1);
	}
	for(i=0;i<SUM;i++,i++)
		AntiBase[i]=tmpAnti[i];
	ResortAntiBase();//just for dubug
}

void CIGA::SelectAnti(double psum,int&i1,int&i2)
{	double r1,r2,tmpsum=0;
	r1=AverageRandom(0,1);
	r2=AverageRandom(0,1);
	for(i1=0;i1<SUM && tmpsum<r1;i1++)//輪賭法
		tmpsum+=AntiBase[i1].sp/psum;
	tmpsum=0;
	for(i2=0;i2<SUM && tmpsum<r2;i2++)//輪賭法
		tmpsum+=AntiBase[i2].sp/psum;
}

void CIGA::CrossOver(int i1, int i2, CAntiCell *c1, CAntiCell *c2)
{	int i,k;
	double r,tmpsum=0;
//	srand((unsigned)time(NULL)); //產生種子
	r=AverageRandom(0,1);
	if(r<1-PC)	//直接克隆
	{	*c1=AntiBase[i1];
		*c2=AntiBase[i2];
	}
	else
	{	k=rand()%4;	//產生交叉位置
		for(i=0;i<k;i++)//
		{	c1->AntiGene[i]=AntiBase[i1].AntiGene[i];
			c2->AntiGene[i]=AntiBase[i2].AntiGene[i];
		}
		for(i=k;i<VAR;i++)//
		{	c1->AntiGene[i]=AntiBase[i2].AntiGene[i];
			c2->AntiGene[i]=AntiBase[i1].AntiGene[i];
		}
	}
}

void CIGA::Mutation(CAntiCell *ac1,CAntiCell *ac2)
{	int k;
	double r1,r2,anti;
//	srand((unsigned)time(NULL)); //產生種子
	r1=AverageRandom(0,1);
	r2=AverageRandom(0,1);
	if(r1>1-PM)	//判斷是否產生變異
	{	anti=AverageRandom(TL,TH);
		k=rand()%VAR;	//產生交叉位置
		ac1->AntiGene[k]=anti;
	}
	if(r2>1-PM)	//判斷是否產生變異
	{	anti=AverageRandom(TL,TH);
		k=rand()%VAR;	//產生交叉位置
		ac2->AntiGene[k]=anti;
	}
}

void CIGA::InitMemory()
{	for(int i=0;i<MSUM;i++)
		MB[i]=AntiBase[i];
}