蟻群算法c語(yǔ)言

蟻群算法（C語(yǔ)言實(shí)現(xiàn)）.en/s/blog6bblb4b001016pt0.html'這是我2012年第一個(gè)學(xué)期的期末大作業(yè)。不得不說(shuō)這段代碼參考了網(wǎng)絡(luò)上的一些代碼，以及書(shū)本上的一些相關(guān)數(shù)據(jù)。但是無(wú)論是網(wǎng)絡(luò)上還是書(shū)上難免總是有很多C++的影子，這是我自己理解之后用C編寫(xiě)的蟻群算法的代碼。這幾個(gè)城市的相關(guān)測(cè)試數(shù)據(jù)在CSDN上可以找到。//蟻群算法關(guān)于簡(jiǎn)單的TSP問(wèn)題求解//#include〈stdio.h〉#include〈stdlib.h〉#include〈string.h〉#include〈math.h〉#include〈time.h>#defineM13//螞蟻的數(shù)量#defineN144//城市的數(shù)量#defineR1000//迭代次數(shù)#defineIN1 //初始化的信息素的量#defineMAX0x7fffffff//定義最大值structcoordinate{charcity[15]; //城市名intx; //城市相對(duì)橫坐標(biāo)inty; //城市相對(duì)縱坐標(biāo)}coords[N];doublegraph[N][N]; //儲(chǔ)存城市之間的距離的鄰接矩陣自己到自己記作MAXdoublephe[N][N];//每條路徑上的信息素的量doubleadd[N][N];//代表相應(yīng)路徑上的信息素的增量doubleyita[N][N];//啟發(fā)函數(shù),yita[i][j]=l/graph[i][j]intvis[M][N];//標(biāo)記已經(jīng)走過(guò)的城市intmap[M][N];//map[K][N]記錄第K只螞蟻?zhàn)叩穆肪€doublesolution[M];//記錄某次循環(huán)中每只螞蟻?zhàn)叩穆肪€的距離intbestway[N];//記錄最近的那條路線doublebestsolution二MAX;intNcMax;//代表迭代次數(shù)，理論上迭代次數(shù)越多所求的解更接近最優(yōu)解,最具有說(shuō)服力doublealpha,betra,rou,Q;voidInitialize();//信息初始化voidInputcoords(FILE*fp);//將文件中的坐標(biāo)信息讀入voidGreateGraph();//根據(jù)坐標(biāo)信息建圖doubleDistance(int*p);//計(jì)算螞蟻所走的路線的總長(zhǎng)度voidResult();//將結(jié)果保存到out.txt中voidInitialize(){

{d+二graph[*(p+i)][*(p+i+l)];}d+=graph[*(p+i)][*(p)];returnd;}voidResult(){FILE*fl;inti;fl=fopen("out.txt","a"); //將結(jié)果保存在out.txt這個(gè)文件里面fprintf(fl,"%s\n","本次算法中的各參數(shù)如下:")；fprintf(fl,"alpha=%.3lf,betra=%.3lf,rou=%.3lf,Q=%.3lf\n",alpha,betra,rou,Q);fprintf(fl,"%s%d\n","本次算法迭代次數(shù)為:",NcMax);fprintf(fl,"%s%.4lf\n","本算法得出的最短路徑長(zhǎng)度為:",bestsolution);fprintf(fl,"%s\n","本算法求得的最短路徑為:");for(i=0;i<N;++i)fprintf(fl,"%sf",coords[bestway[i]].city);fprintf(fl,"%s",coords[bestway[0]].city);fprintf(fl,"\n\n\n");fclose(fl);return;

memset(map,T,sizeof(map));//把螞蟻?zhàn)叩穆肪€置空memset(vis,0,sizeof(vis));//O表示未訪問(wèn)，1表示已訪問(wèn)srand(time(NULL));while(NC++〈二NcMax){for(k=0;k〈M;++k){map[k][O]=(k+NC)%N;//給每只螞蟻分配一個(gè)起點(diǎn)，并且保證起點(diǎn)在N個(gè)城市里vis[k][map[k][O]]=1;//將起點(diǎn)標(biāo)記為已經(jīng)訪問(wèn)}s=1;while(s<N){for(k=0;k〈M;++k){psum=O;for(j=0;j〈N;++j){if(vis[k][j]==O){psum+二pow(phe[map[k][sT]][j],alpha)*pow(yita[map[k][sT]][j],betra);}}drand=(double)(rand()POO);drand/=5000.0; //生成一個(gè)小于1的隨機(jī)數(shù)pro=0;for(j=0;j〈N;++j){if(vis[k][j]==O)pro+二pow(phe[map[k][sT]][j],alpha)*pow(yita[map[k][sT]][j],betra)/psum;if(pro〉drand)break;}vis[k][j]=1; //將走過(guò)的城市標(biāo)記起來(lái)map[k][s]=j;//記錄城市的順序}s++;}memset(add,0,sizeof(add));for(k=0;k<M;++k) //計(jì)算本次中的最短路徑//{solution[k]二Distance(map[k]);//螞蟻k所走的路線的總長(zhǎng)度if(solution[k]〈bestsolution){bestsolution=sol