% 蠟傈磊 舅絆府硫闌 撈儈茄 鑒雀 魄概盔 鞏力(TSP)狼 秦搬
% N: 檔矯狼 薦. 檔矯甸籃 1俊輯 N欄肺 扼駭傅. 
clear allN=input('鑒雀且 檔矯狼 薦: ');%======================================================================
% N俺狼 檔矯俊 措茄 TSP 鞏力 積己 pos=rand(N,2);  % 竄困 瀝葷阿屈郴俊 檔矯狼 困摹甫 烙狼肺 汲瀝 
distance=dist(pos,pos); %檔矯埃狼 芭府 青紡  figure(1),clfplot(pos(:,1),pos(:,2),'o') %檔矯 困摹 釬矯axis([0 1 0 1]), axis square, hold on%弊覆惑俊 檔矯 鍋齲 扼駭傅
for i=1:N   text(pos(i,1)+.03,pos(i,2),int2str(i));endhold off%========================================================================
% 蠟傈磊 漚禍 disp('********** 蠟傈磊 舅絆府硫*****************');ngen=input('柳拳矯懦 技措薦: ');ngpool=input('# 蠟傈磊 錢(gene pool)郴狼 俺眉(chromosoms)薦: '); % 蠟傈磊 錢狼 農扁gpool=zeros(ngpool,N); % 蠟傈磊 錢(gene pool)% 蠟傈磊 錢 檬扁拳 
for i=1:ngpool, 
   gpool(i,:)=[1 randomize([2:N]')'];end
oldstep=0;
oldcostmin=99999;costmin=N, tourmin=zeros(1,N); cost=zeros(1,ngpool); oldtourmin=gpool(2,:);
%========================================================================
for step=1:ngen,    % 1竄拌. 泅犁 俺眉狼 利欽檔 乞啊, 鑒雀 芭府啊 利欽檔啊 登哥 陋闌薦廢 亮促.     for i=1:ngpool,      cost(i)=sum(diag(distance(gpool(i,:)',rshift(gpool(i,:))')));   end   % 泅犁 彌惑狼 菩膠甫 扁廢    [costmin,idx]=min(cost);   tourmin=gpool(idx,:);   % 2竄拌. 背瞞(crossover)客 倒楷函撈(mutation) 楷魂 薦青    [csort,ridx]=sort(cost); % 厚儈撈 利籃巴俊輯 奴巴欄肺 瀝紡茄促.    
   for i=2:ngpool/2 
       sameidx=[gpool(2*i-1,:)==gpool(2*i,:)]; 
       diffidx=find(sameidx==0); 
       if length(diffidx)<=2 
           % 倒楷函撈積己
           gpool(2*i,:)=[1 randomize([2:N]')']; 
       else
           if i%3 == 0 % 何盒背眉背硅
             [gpool(2*i-1,:),gpool(2*i,:)]=PMX(gpool(ridx(2*i-1),:),gpool(ridx(2*i),:));         
           end
           if i%3 == 1  % 鑒輯背硅
             [gpool(2*i-1,:),gpool(2*i,:)]=OX(gpool(ridx(2*i-1),:),gpool(ridx(2*i),:));
           end
           if i%3 == 2 % 林扁背硅
             [gpool(2*i-1,:),gpool(2*i,:)]=CX(gpool(ridx(2*i-1),:),gpool(ridx(2*i),:));
           end
       end 
   end 
   
   %pause;
   %彌利鑒雀 風飄啊 函竅綽 何盒父 免仿秦林扁 困茄 貿府 
   sameidx=[oldtourmin==tourmin]; % 鞍籃 鎊俊綽 #1肺 付膠農甫 競快絆                            
   diffidx=find(sameidx==0); % 促弗 牢鄲膠甫 茫酒輯 
   if length(diffidx) ~= 0 & costmin < oldcostmin     
       figure;
       plot(pos(:,1),pos(:,2),'o'),axis([0 1 0 1]), axis square, hold on 
       trip=pos(tourmin',:); trip=[trip; trip(1,:)]; 
       plot(trip(:,1),trip(:,2),'-');
       title(['技措:' num2str(step) ',  蠟傈磊 八禍, 厚儈:' num2str(costmin)])
       hold off;
       drawnow
       disp(['技措' num2str(step) '::彌利鑒雀:'  num2str(tourmin)  '::厚儈:'  num2str(costmin) ]);
      oldcostmin = costmin;
      oldtourmin = tourmin;
      oldstep=step;
   end
enddisp(['厚儈 竊薦 乞啊 : ' int2str(ngen*2+ngpool) ' times!'])disp(['彌家 鑒雀 芭府 : ' num2str(oldcostmin)])disp('彌利 鑒雀 : ') 
disp(num2str(oldtourmin))figure
plot(pos(:,1),pos(:,2),'o'),axis([0 1 0 1]), axis square, hold ontrip=pos(oldtourmin',:); trip=[trip; trip(1,:)];plot(trip(:,1),trip(:,2),'-')title(['技措:' num2str(oldstep) ',  蠟傈磊 八禍, 厚儈:' num2str(oldcostmin)])
hold off