function fval=genetic(number);
X= [ ] ; Y= [ ] ; n = 1 ;lchrom=53;
distmatrix=[0	152	191	400	235	332	405	557	555	619	396	389	244	188	184	373	305	352	342	397	318	238	228	139	359	233	155	234	268	394	246	327	173	218	300	103	59	112	273	477	511	330	287	270	342	339	450	258	271	161	317	189	60
152	0	48	257	83	180	253	405	403	576	411	404	259	203	199	403	457	367	357	353	395	390	380	291	317	298	307	345	300	449	313	394	325	370	452	255	211	127	130	325	468	345	302	285	357	436	298	197	339	193	372	221	92
191	48	0	209	131	228	233	385	383	556	363	356	211	155	151	355	423	319	309	401	443	368	358	269	365	346	285	364	348	497	361	442	364	365	447	294	189	79	82	305	448	297	254	237	309	411	346	245	387	241	420	269	140
400	257	209	0	240	337	278	430	428	601	492	565	420	364	360	564	632	528	518	478	552	577	567	478	474	555	494	573	557	706	570	651	573	574	656	503	398	288	127	350	493	506	463	446	518	593	423	354	496	450	629	478	349
235	83	131	240	0	97	170	322	320	493	384	487	342	286	282	486	451	445	363	270	312	398	388	374	234	361	390	428	383	532	396	477	408	453	535	338	294	210	113	242	385	428	385	368	440	353	215	114	256	276	455	304	175
332	180	228	337	97	0	73	225	223	396	287	410	439	383	379	456	391	385	266	173	252	352	342	369	215	342	420	499	480	629	493	574	505	550	632	435	391	307	210	145	288	353	482	338	347	293	118	95	237	373	552	401	272
405	253	233	278	170	73	0	152	150	323	214	337	366	388	384	483	418	412	293	200	279	379	369	442	265	415	493	572	553	702	566	647	578	598	680	508	422	312	151	72	215	280	439	365	346	320	145	168	310	446	625	474	345
557	405	385	430	322	225	152	0	158	331	222	345	374	430	536	625	570	436	435	352	431	531	521	594	417	567	645	724	705	854	718	799	730	750	832	660	574	464	303	80	223	288	447	507	354	472	297	320	462	598	777	626	497
555	403	383	428	320	223	150	158	0	173	220	187	351	388	534	623	568	434	433	350	429	529	519	592	415	565	643	722	703	852	716	797	728	748	830	658	572	462	301	78	65	265	289	505	352	470	295	318	460	596	775	624	495
619	576	556	601	493	396	323	331	173	0	374	230	394	431	579	676	741	574	573	523	602	702	692	697	588	738	713	792	830	952	843	924	792	793	875	722	617	507	474	251	108	308	332	558	492	643	468	491	633	723	875	751	622
396	411	363	492	384	287	214	222	220	374	0	144	152	208	356	403	404	214	213	306	347	407	417	474	371	564	490	569	607	729	620	701	569	570	652	499	394	284	365	142	266	66	246	285	132	306	359	382	459	500	652	528	399
389	404	356	565	487	410	337	345	187	230	144	0	164	201	349	446	514	344	343	436	477	537	547	467	501	561	483	562	600	722	613	694	562	563	645	492	387	277	438	265	122	78	102	328	262	436	482	505	589	493	645	521	392
244	259	211	420	342	439	366	374	351	394	152	164	0	56	204	282	350	180	179	272	313	373	383	322	337	416	338	417	455	577	468	549	417	418	500	347	242	132	293	294	286	86	155	164	98	272	327	445	425	348	500	376	247
188	203	155	364	286	383	388	430	388	431	208	201	56	0	148	338	406	236	235	328	369	365	355	266	393	360	282	361	399	521	412	493	361	362	444	291	186	76	237	350	323	142	99	220	154	328	383	389	398	292	444	320	191
184	199	151	360	282	379	384	536	534	579	356	349	204	148	0	204	272	168	158	251	301	339	349	262	316	356	278	357	395	517	408	489	357	358	440	287	182	72	233	456	471	290	247	86	239	260	306	385	394	288	440	316	187
373	403	355	564	486	456	483	625	623	676	403	446	282	338	204	0	68	200	190	283	207	135	145	234	285	382	410	489	599	649	579	660	506	490	572	436	314	276	437	545	568	368	437	118	271	166	338	405	277	487	572	520	391
305	457	423	632	451	391	418	570	568	741	404	514	350	406	272	68	0	190	258	218	139	67	77	166	217	314	342	421	535	581	511	592	438	422	504	368	246	344	505	490	633	436	505	186	272	98	273	337	209	419	504	456	365
352	367	319	528	445	385	412	436	434	574	214	344	180	236	168	200	190	0	154	212	133	193	203	292	211	404	446	525	563	685	576	657	525	526	608	455	350	240	401	356	466	266	335	82	82	92	267	331	299	456	608	484	355
342	357	309	518	363	266	293	435	433	573	213	343	179	235	158	190	258	154	0	93	172	272	262	351	158	351	429	508	553	668	566	647	515	516	598	445	340	230	391	355	465	265	334	72	81	213	148	278	246	446	591	474	345
397	353	401	478	270	173	200	352	350	523	306	436	272	328	251	283	218	212	93	0	79	179	169	258	65	258	336	415	515	575	603	678	530	514	596	460	338	323	351	272	415	358	427	165	174	120	55	185	153	363	498	512	383
318	395	443	552	312	252	279	431	429	602	347	477	313	369	301	207	139	133	172	79	0	100	90	179	78	271	349	428	528	588	524	605	451	435	517	381	259	369	425	351	494	399	468	215	215	41	134	198	166	376	511	469	378
238	390	368	577	398	352	379	531	529	702	407	537	373	365	339	135	67	193	272	179	100	0	10	99	178	247	275	354	468	514	444	525	371	355	437	301	179	289	450	451	594	459	464	253	275	101	234	284	142	352	437	389	298
228	380	358	567	388	342	369	521	519	692	417	547	383	355	349	145	77	203	262	169	90	10	0	89	168	237	265	344	458	504	434	515	361	345	427	291	169	279	440	441	584	469	454	263	285	111	224	274	132	342	427	379	288
139	291	269	478	374	369	442	594	592	697	474	467	322	266	262	234	166	292	351	258	179	99	89	0	220	237	176	255	369	415	345	426	272	256	338	202	80	190	351	514	589	408	365	348	374	200	313	274	132	300	338	290	199
359	317	365	474	234	215	265	417	415	588	371	501	337	393	316	285	217	211	158	65	78	178	168	220	0	193	271	350	450	510	539	613	492	476	558	422	300	388	347	337	480	423	492	230	239	119	120	120	88	298	433	447	318
233	298	346	555	361	342	415	567	565	738	564	561	416	360	356	382	314	404	351	258	271	247	237	237	193	0	78	157	257	317	427	420	366	350	432	296	174	284	428	487	630	502	459	423	432	312	313	247	105	105	240	335	206
155	307	285	494	390	420	493	645	643	713	490	483	338	282	278	410	342	446	429	336	349	275	265	176	271	78	0	79	234	239	361	342	288	272	354	218	96	206	367	565	605	424	381	364	436	376	391	325	183	183	162	306	215
234	345	364	573	428	499	572	724	722	792	569	562	417	361	357	489	421	525	508	415	428	354	344	255	350	157	79	0	155	160	326	263	367	351	433	297	175	285	446	644	684	503	460	443	515	455	470	404	262	152	83	234	253
268	300	348	557	383	480	553	705	703	830	607	600	455	399	395	599	535	563	553	515	528	468	458	369	450	257	234	155	0	149	171	252	237	282	364	167	289	323	430	625	722	541	498	481	553	569	570	406	362	152	72	79	208
394	449	497	706	532	629	702	854	852	952	729	722	577	521	517	649	581	685	668	575	588	514	504	415	510	317	239	160	149	0	184	103	257	431	513	316	335	445	579	774	844	663	620	603	675	615	630	555	422	301	77	228	357
246	313	361	570	396	493	566	718	716	843	620	613	468	412	408	579	511	576	566	603	524	444	434	345	539	427	361	326	171	184	0	81	73	258	340	143	265	336	443	638	735	554	511	494	566	545	611	419	477	322	243	92	221
327	394	442	651	477	574	647	799	797	924	701	694	549	493	489	660	592	657	647	678	605	525	515	426	613	420	342	263	252	103	81	0	154	339	421	224	346	417	524	719	816	635	592	575	647	626	692	500	525	403	180	173	302
173	325	364	573	408	505	578	730	728	792	569	562	417	361	357	506	438	525	515	530	451	371	361	272	492	366	288	367	237	257	73	154	0	185	267	70	192	285	446	650	684	503	460	443	515	472	585	431	404	334	309	158	233
218	370	365	574	453	550	598	750	748	793	570	563	418	362	358	490	422	526	516	514	435	355	345	256	476	350	272	351	282	431	258	339	185	0	82	115	176	286	447	670	685	504	461	444	516	456	569	476	388	379	354	203	278
300	452	447	656	535	632	680	832	830	875	652	645	500	444	440	572	504	608	598	596	517	437	427	338	558	432	354	433	364	513	340	421	267	82	0	197	258	368	529	752	767	586	543	526	598	538	651	558	470	461	436	285	360
103	255	294	503	338	435	508	660	658	722	499	492	347	291	287	436	368	455	445	460	381	301	291	202	422	296	218	297	167	316	143	224	70	115	197	0	122	215	376	580	614	433	390	373	445	402	515	361	334	264	239	88	163
59	211	189	398	294	391	422	574	572	617	394	387	242	186	182	314	246	350	340	338	259	179	169	80	300	174	96	175	289	335	265	346	192	176	258	122	0	110	271	494	509	328	285	268	340	280	393	317	212	220	258	210	119
112	127	79	288	210	307	312	464	462	507	284	277	132	76	72	276	344	240	230	323	369	289	279	190	388	284	206	285	323	445	336	417	285	286	368	215	110	0	161	384	399	218	175	158	230	332	378	313	322	216	368	244	115
273	130	82	127	113	210	151	303	301	474	365	438	293	237	233	437	505	401	391	351	425	450	440	351	347	428	367	446	430	579	443	524	446	447	529	376	271	161	0	223	366	379	336	319	391	466	296	227	369	323	502	351	222
477	325	305	350	242	145	72	80	78	251	142	265	294	350	456	545	490	356	355	272	351	451	441	514	337	487	565	644	625	774	638	719	650	670	752	580	494	384	223	0	143	208	367	427	274	392	217	240	382	518	697	546	417
511	468	448	493	385	288	215	223	65	108	266	122	286	323	471	568	633	466	465	415	494	594	584	589	480	630	605	684	722	844	735	816	684	685	767	614	509	399	366	143	0	200	224	450	384	535	360	383	525	615	767	643	514
330	345	297	506	428	353	280	288	265	308	66	78	86	142	290	368	436	266	265	358	399	459	469	408	423	502	424	503	541	663	554	635	503	504	586	433	328	218	379	208	200	0	180	250	184	358	413	448	511	434	586	462	333
287	302	254	463	385	482	439	447	289	332	246	102	155	99	247	437	505	335	334	427	468	464	454	365	492	459	381	460	498	620	511	592	460	461	543	390	285	175	336	367	224	180	0	319	253	427	482	488	497	391	543	419	290
270	285	237	446	368	338	365	507	505	558	285	328	164	220	86	118	186	82	72	165	215	253	263	348	230	423	364	443	481	603	494	575	443	444	526	373	268	158	319	427	450	250	319	0	153	174	220	350	318	374	526	402	273
342	357	309	518	440	347	346	354	352	492	132	262	98	154	239	271	272	82	81	174	215	275	285	374	239	432	436	515	553	675	566	647	515	516	598	445	340	230	391	274	384	184	253	153	0	174	229	359	327	446	598	474	345
339	436	411	593	353	293	320	472	470	643	306	436	272	328	260	166	98	92	213	120	41	101	111	200	119	312	376	455	569	615	545	626	472	456	538	402	280	332	466	392	535	358	427	174	174	0	175	239	207	417	538	490	399
450	298	346	423	215	118	145	297	295	468	359	482	327	383	306	338	273	267	148	55	134	234	224	313	120	313	391	470	570	630	611	692	585	569	651	515	393	378	296	217	360	413	482	220	229	175	0	213	208	418	553	519	390
258	197	245	354	114	95	168	320	318	491	382	505	445	389	385	405	337	331	278	185	198	284	274	274	120	247	325	404	406	555	419	500	431	476	558	361	317	313	227	240	383	448	488	350	359	239	213	0	142	299	478	327	198
271	339	387	496	256	237	310	462	460	633	459	589	425	398	394	277	209	299	246	153	166	142	132	132	88	105	183	262	362	422	477	525	404	388	470	334	212	322	369	382	525	511	497	318	327	207	208	142	0	210	345	422	311
161	193	241	450	276	373	446	598	596	723	500	493	348	292	288	487	419	456	446	363	376	352	342	300	298	105	183	152	152	301	322	403	334	379	461	264	220	216	323	518	615	434	391	374	446	417	418	299	210	0	224	230	101
317	372	420	629	455	552	625	777	775	875	652	645	500	444	440	572	504	608	591	498	511	437	427	338	433	240	162	83	72	77	243	180	309	354	436	239	258	368	502	697	767	586	543	526	598	538	553	478	345	224	0	151	280
189	221	269	478	304	401	474	626	624	751	528	521	376	320	316	520	456	484	474	512	469	389	379	290	447	335	306	234	79	228	92	173	158	203	285	88	210	244	351	546	643	462	419	402	474	490	519	327	422	230	151	0	129
60	92	140	349	175	272	345	497	495	622	399	392	247	191	187	391	365	355	345	383	378	298	288	199	318	206	215	253	208	357	221	302	233	278	360	163	119	115	222	417	514	333	290	273	345	399	390	198	311	101	280	129	120


    ];
popsize=53;
n = 1;pop = zeros(popsize) ;
while n<=popsize ,
  pop(n,:) = randperm(lchrom);
  n = n + 1 ;
end ;
countnum=0;
xg=[1:popsize];
for i=1:number
     maxn=0;
     while maxn<1000,
      num=randperm(popsize);
      for j=1:popsize-8
          [pnew1,pnew2]=crossover(pop(num(j),:),pop(num(j+1),:));
          if Fitness(pnew1,distmatrix)<Fitness(pop(num(j),:),distmatrix),
              pop(num(j),:)=pnew1;
          end;
          if Fitness(pnew2,distmatrix)<Fitness(pop(num(j+1),:),distmatrix),
              pop(num(j+1),:)=pnew1;
          end;
         %[pop(num(j),:),pop(num(j+1),:)]=
      end;
     for j=1:8
        pnew1=mutation(pop(num(j),:),num(j),num(j+1));
        if Fitness(pnew1,distmatrix)<Fitness(pop(num(j),:),distmatrix),
          pop(num(j),:)=pnew1;
        end;
         pnew2=inverse(pop(num(j),:),num(j+2),num(j+1),popsize);
         if Fitness(pnew2,distmatrix)<Fitness(pop(num(j),:),distmatrix),
          pop(num(j),:)=pnew1;
        end;
     end;
      maxn=maxn+1;
   end;
    for j=1:popsize
        if Fitness(pop(j,:),distmatrix)<Fitness(xg,distmatrix),
            xg=pop(j,:);
        end;
    end;
    %for j=1:popsize
       % if Fitness(pop(j,:),distmatrix)>Fitness(xg,distmatrix),
           % pop(j,:)=xg;
        %end;
    %end;
   fval = Fitness(xg,distmatrix);
end; 

fval