遗传算法解决tsp问题

1. /*baixzh SYSU 2008-12-16 */
2. #include<iostream>
3. #include<string>
4. #include<math.h>
5. #include<cstdlib>
6. #include<algorithm>
7. using namespace std;
9. int i;int j;    //循环变量
10. int t=0;//选择代数
11. struct city//城市节点
12. {
13.     double x;
14.     double y;
15.     city(); 
16. };
17. city::city() 
18. {
19.     x=0.0;
20.     y=0.0;
21. }
22. double pos[20]=//初始城市坐标数据
23.     {
24.         0.400,  0.4439,//代表x,y
25.         0.2439, 0.1463,
26.         0.1707, 0.2293,
27.         0.2293, 0.7610,
28.         0.5171, 0.9414,
29.         0.8732, 0.6536,
30.         0.6878, 0.5219,
31.         0.8488, 0.3609,
32.         0.6683, 0.2536,
33.         0.6195, 0.2634
35.     };
36. city c_t[10];
37. double dis[10][10];//城市距离矩阵
38. struct c_gen
39. {
40.     string route;//旅行商路径也代表基因编码
41.     double lenth;//路径长度
42.     double p;//基因的概率
43.     c_gen();
44. };
45. c_gen::c_gen()
46. {
47.     route="";
48.     lenth=0.0;
49.     p=0.0;
50. };
51. c_gen group[100];//种群
52. void init_population();//初始化种群
53. void cal_p();//种群中各个基因概率计算
54. void select();//选择函数 淘汰33%
55. void mating();//交配
56. void aberrance();//变异
57. void replace_population();//代替旧种群
58. bool cmp(c_gen a,c_gen b)//自定义排序比较函数
59. {
60.     return a.lenth <b.lenth ;
61. }
62. void init_population()
63. {
64.     for(i=0;i<100;i++)
65.     {
66.         int an[10]={2,5,4,9,8,7,0,1,6,3}    ;//给出一个原始祖先
67.         int ran1=rand()%10;
68.         int ran2=rand()%10;
69.         if(ran1==ran2)
70.         
71.         {
72.             ran2=(ran2+rand()%5+1)%10;
73.         }
74.         swap(an[ran1],an[ran2]);
75.         swap(an[(ran1+9)%10],an[(ran1+4)%10]);
76.         for(j=0;j<10;j++)
77.         {
78.             group[i].route+=char(an[j]+'0');
79.         }
80.         //产生100个随机的初始祖先 组成种群 
81.         for(j=0;j<9;j++)
82.         {
83.             group[i].lenth +=dis[group[i].route[j]-'0'][group[i].route[j+1]-'0'];
84.         }
85.         group[i].lenth +=dis[group[i].route[0]-'0'][group[i].route[9]-'0'];
86.     }
87.      
88.     return ;
90. }
91. void cal_p()
92. {
93.     double sum=0.0;
94.     for(i=0;i<100;i++)
95.     {
96.         sum+=group[i].lenth ;
97.     }
98.     for(i=0;i<100;i++)
99.     {
100.         group[i].p=group[i].lenth /sum;
101.     }
103. return ;
104. }
105. void select()
106. {
107.     sort(group,group+100,cmp);//种群按照评价值进行排序
108.     cout<<"第"<<t<<"代最优品种:  " ;
109.     for(int temp=0;temp<10;temp++)
110.     {
111.         cout<<char(group[0].route [temp]-'0'+'A')<<"-->";
112.     }
113.     cout<<char(group[0].route [0]-'0'+'A')<<endl;
114.     cout<<"路径长度为:   "<<group[0].lenth <<endl;
116. return ;
117. }
118. void mating()//种群规模100 0-65的优秀品种交配 66保留 67-99品种淘汰 被交配的子代代替
119. {
120.      
121. for(i=0;i<33;i++)
122. {
123.     int ran1=rand()%10;
124.     int ran2=rand()%10;
125.     if(ran1==ran2)
126.     {
127.         ran2=(ran2+2)%10;
128.     }
129.     char c1=group[i*2].route [ran1];
130.     char c2=group[i*2].route [ran2];    
131.     //偶数基因中随机选出2段染色体然后从奇数基因中除去这2段染色体 再将这2段随机插入产生子代
132.     string t=group[i*2+1].route ;
133.     int t1,t2;
134.     int temp;
135.     bool k1=false;
136.     bool k2=false;
137.     for(temp=0;temp<10;temp++)
138.     {   
139.         if(t[temp]==c1)
140.         {
141.         t.erase ( temp,1);
142.          
143.         }
144.     }
146.     for(temp=0;temp<9;temp++)
147.     {   
148.         if(t[temp]==c2)
149.         {
150.         t.erase ( temp,1);
151.          
152.         }
153.     }
155.     int ran3=rand()%8;
156.     t.insert (t.begin ()+ran3,c1);
157.     t.insert (t.begin ()+(ran3+rand())%8,c2);   
158.     group[67+i].route =t;//子代产生 加入种群
159. }
163. return ;
164. }
165. void aberrance()
166. {
167.     int ran1=rand()%100;
168.     int ran2=rand()%10;
169.     int ran3=rand()%10;
170.     if(ran2==ran3)
171.     {
172.     ran3=(ran3+rand()%5+1)%10;
173.     }
174.     swap(group[ran1].route [ran2],group[ran1].route [ran3]);
175.     //随机选出一段基因,然后随机交换2个染色体位置 产生变异
177. return ;
178. }
179. void replace_population()
180. {
181.     for(i=0;i<100;i++)
182.     {
183.         group[i].lenth =0.0;
184.     for(j=0;j<9;j++)
185.         {
186.             
187.             group[i].lenth +=dis[group[i].route[j]-'0'][group[i].route[j+1]-'0'];
188.         }
189.         group[i].lenth +=dis[group[i].route[0]-'0'][group[i].route[9]-'0'];
190.     }
191.     //新种群的评价值重新计算,代替原种群
192. }
193. int main()
194. {
195.         for(i=0;i<10;i++)
196.         {
197.             c_t[i].x=pos[i*2];
198.             c_t[i].y=pos[i*2+1];
200.         }//城市节点初始化
201.         memset(dis,0.0,sizeof(dis));
202.         for(i=0;i<10;i++)
203.             for(j=0;j<10;j++)
204.             {
205.                 dis[i][j]=sqrt(pow((c_t[i].x -c_t[j].x),2) +pow((c_t[i].y-c_t[j].y),2)) ; 
206.             }
207.         //距离矩阵初始化
209.     init_population();//种群初始化
210. //遗传算法主体,经过50代的最优解完全符合最优解的条件
211.     for(t=0;t<50;t++)
212.     {
213.             select();
214.             mating();
215.             aberrance();
216.             replace_population();
217.              
218.     }
221. return 0;
223. }