C语言Folyd算法和模拟退火算法解决新增巡警平台以及优化全市平台设置程序.doc

#include stdio.h#include math.h#include malloc.h#define N 582#define T 929#define M 582#define Max 10000double *lshujuN+1;double *zdjN+1;void juzhen()int i,j;float zshujuN+14;double xx=0.0;double lucheng;int a,b;FILE *fp1,*fp2,*fp3,*fp4;printf(gajj);for(i=1;i=N;i+)lshujui=(double *)malloc(N+1)*sizeof(double);for(i=1;i=N;i+)for(j=1;j=N;j+)lshujuij=Max;for(i=1;i=N;i+)lshujuii=0;printf(ijdfjj);fp1=fopen(d1.txt,r);for(i=1;i=N;i+)fscanf(fp1,%f%f%f,&zshujui1,&zshujui2,&zshujui3);fclose(fp1);fp2=fopen(d2.txt,r);for(i=1;i %d , %d , - %fn,i,a,b,lucheng);/if(i % 50 = 0)/getchar();lshujuab=lucheng;lshujuba=lucheng;/printf(%f %f ,lshujuab,lshujuba);fclose(fp2);fp3=fopen(lcy.txt,w);for(i=1;i=N;i+)for(j=1;j=N;j+)if(lshujuij0)/printf(%f,lshujuij);fprintf(fp3,%d %d: %f,i,j,lshujuij);fprintf(fp3,n);fclose(fp3); fp4=fopen(jj.txt,w);for(i=1;i=M;i+)for(j=1;j=N;j+)xx = lshujuij;if(xx=Max)xx = 0; fprintf(fp4,%d %d: %f n,i,j,xx);fclose(fp4); main() int i,j,k,t=0,mark=0; FILE *fp,*fp1; double *bN+1,*cN+1; juzhen();for(i=1;i=N;i+) bi=(double *)malloc(N+1)*sizeof(double); ci=(double *)malloc(N+1)*sizeof(double); / printf(di 1 ci de ju zheng:n); /输出表示带权图的矩阵a for(i=1;i=N;i+) for(j=1;j=N;j+) cij=lshujuij; / printf( %f ,cij); / printf(n); for(t=1;t=N;t+) / t为矩阵a相乘的次数 mark=0; /mark为矩阵a进行t次相乘后看是否第t次的矩阵和t+1次矩阵相等，mark=0表示相等停止矩阵a做乘法，否则继续做矩阵a的乘法 /printf(di %d ci de ju zheng:n,t+1); for(i=1;i=N;i+)for(k=1;k=N;k+) bik=cik; /b矩阵用来记住上一次相乘后的矩阵 for(i=1;i=N;i+) for(k=1;k=N;k+)for(j=1;jbij+lshujujk) /找到顶点vi到vk的最短距离，放在cik中 printf( %f ,bij+lshujujk); cik=bij+lshujujk; / printf( %f,cik); if(cik!=bik) mark=1; / printf(n); if(mark=0) t=N+1; printf(dddsadasfn); /* for(i=1;i=M;i+) zdji=(double *)malloc(M+1)*sizeof(double);*/fp=fopen(zdl.txt,w);fp1=fopen(zdlj.txt,w); for(k=1;k=M;k+) for(j=1;j=M;j+) for(i=1;i 点%d 的最短路径 :点%d 的+点 %d-点%d ,k,j,i,i,j); fprintf(fp, quan he:%fn,ckj); fprintf(fp,nnn); printf(dsjfjasjgjka); for(i=1;iN+