【文章內容簡介】 數組,用于記錄該個體所有小方格能達到的最大高度,數組一長度為預建建筑區(qū)域的小方格個數n*/double m_fitn。/*適應度*/double m_Ps。/*選擇概率*/double m_Qs。/*累計選擇概率*/ }std::vectorindividualpopulation。//populat={individual[1],individual[2],..,individualpp_s。z一,}//產生初始種群后調用void KeepTheBest(){ (population[pop_size1])。 int cur_best=pop_size_1。 for(int i=0。ipop_size。++i) { if(population[i].m_fitnesspopulation[pop_size].m_fitness) { cur_best=i。 population[pop_size].m_fitness=population[i].m_fitness。 } } population[pop_size]=population[cur_best]。}//進行完遺傳操作后調用void Elitist(){ double bestFit=population[0].m_fitness,worseFit=population[0].m_fitness。 int bestIdx=0。 std::vectorintworseIdxs。 for(int i=0。ipop_size1。++i) { if(population[i].m_fitnesspopulation[i+1].m_fitness) { if(population[i].m_fitness=bestFit) { bestFit=population[i].m_fitness。 bestIdx=i。 } if(population[i+1].m_fitnessworseFit) { worseFit=population[i].m_fitness。 ()。 (i+1)。 } else if(population[i+1].m_fitness==worseFit) (i+1)。 else continue。 } else { if(population[i+1].m_fitness=bestFit) { bestFit=population[i].n_fitness。 bestIdx=i+1。 } if(population[i].m_fitnessworseFit) { worseFit=population[i].m_fitness。 ()。 (i)。 } else if(population[i].m_fitness==worseFit) (i)。 else continue。 } } if(bestFit=population[pop_size].m_fitness) population[pop_size]=population[bestIdx]。 else { for(i=0。i()。++i) population[worseIdxs[i]]=population[pop_size]。 }}//適應度函數/*void Evaluate(double sumhmax){ for(i=O。i(pop_size。++i) { for(j=O。jn。++j) { population[i].m_fitness+=population[i].m_ChrH[j]。 } population[i].m_fitness/=Sumhmax。 }}*///選擇函數Selection(){ int i,j。 double sumFitness=O。 for(i=O。ipop_size。++i) sumFitness+=populationt[i].m_fitness。 double qs=0。//累計選擇概率 for(i=O。ipop_size。++i) { population[i].mee_Ps=population[i].m_fitness/sumFitness。 qs+=population[i].m_Ps。 population[i].m_Qs=qs。 } //根據Qs選擇生存者 std::vectorindividual newpopulation。 double r。 for(i=O。ipop_size。++i) { r=rand()%1000/1000。 for(j=O。jpop_size。++j) { if(r=population[j].m_Qs) break。 } (population[j])。 } (newpopulation)。//用新的一代替換上一代}double GetVariance(individual *indiv1,individual *indiv2){ double variance=0。 for(int i=0。im。++i) { for(int j=0。jsegment。++j) variance+=([i][j][i][j])^2。 } return sqrt(variance)。}float Pc(double fitness,double maxfitness,double avgfitness){ float pc。 if(fitnessavgfitness) pc=float()。 else pc=()*(fitnessavgfitness)/float((maxfitnessavgfitness))。 return pc。}float Pm(double fitness,double maxfitness,double avgfitness){ float pm。 if(fitnessavgfitness) pm=Pm。 else pm=Pm()*(maxfitnessfitness)/float(maxfitnessavgfitness)。 return pm。}//評價函數double Evaluate(double Sumhmax){ double sunfitness=0。 for(i=O。i(pop_size。++i) { for(j=O。jn。++j) { population[i].m_fitness+=population[i].m_ChrH[j]。 } population[i].m_fitness/=Sumhmax。 sumfitness+=population[i].m_fitness。 } return sumfitness/pop_size。}//交叉函數void Crossover(double avgfitness){ int parentId1,parentId2,first=0。 float lastfitness=0。 for(int