灾情巡视C语言代码
/*
"水灾巡视问题"模拟退火算法。
这是一个推销员问题,本题有53个点,所有可能性大约为exp(53),目前没有好方法求出精确解,既然求不出精确解,我们使用模拟退火法求出一个较优解,将所有结点编号为1到53,1到53的排列就是系统的结构,结构的变化规则是:从1到53的排列中随机选取一个子排列,将其反转或将其移至另一处,能量E自然是路径总长度。具体算法描述如下:
步1:设定初始温度T,给定一个初始的巡视路线。
步2:步3 --8循环K次
步3:步 4--7循环M次
步4:随机选择路线的一段
步5:随机确定将选定的路线反转或移动,即两种调整方式:反转、移动。
步6:计算代价D,即调整前后的总路程的长度之差
步7:按照如下规则确定是否做调整:
如果D<0,则调整
如果D>0,则按照EXP(-D/T)的概率进行调整
步8:T*0.9-->T,降温
*/
#include <stdio.h>
#include <math.h>
#define TFACTR 0.9
#define MAX 10000.0
#define MBIG 1000000000
#define MSEED 161803398
#define MZ 0
#define FAC (1.0/MBIG)
float weight[54][54];
int iorder[54];
void anneal (int ncity)
{ int irbit(unsigned long *iseed);
int metrop(float de,float t);
float ran(long *idum);
float reverse_len(int ncity,int n[]);
void reverse(int ncity,int n[]);
float transe_len(int ncity,int n[]);
void transe(int ncity,int n[]);
int ans,nover,nlimit,i1,i2;
int i,j,nsucc,nn,idec;
unsigned long k;
static int n[7];
long idum;
unsigned long iseed;
float path,de,t,T;
nlimit=100*ncity;/*max succ try numbers */
path =0.0;
for(i=1;i<ncity;i++){ /* primitive path length */
i1=iorder[i];
i2=iorder[i+1];
path+=weight[i1][i2];
}
i1=iorder[ncity];
i2=iorder[1];
path+=weight[i1][i2];
idum=-1;
iseed=111;
T=2.8 ;
for(j=1;j<=30;j++){
nsucc=0;
for(k=1;k<=20000;k++){
do { n[1]=1+(int)(ncity*ran(&idum));
n[2]=1+(int)((ncity-1)*ran(&idum));
if(n[2]>=n[1]) ++n[2];
nn=1+((n[1]-n[2]+ncity-1)%ncity);
} while (nn<3);
idec=irbit(&iseed);
if(idec==0){
n[3]=n[2]+(int) (abs(nn-2)*ran(&idum))+1;
n[3]=1+((n[3]-1)%ncity);
de=transe_len(ncity,n);
ans=metrop(de,T);
if (ans){ ++nsucc;
path+=de;
transe(ncity,n);
}
}
else {
de=reverse_len(ncity,n);
ans=metrop(de,T);
if (ans) {++nsucc;
path+=de;
reverse(ncity,n);
}
}
if(nsucc>=nlimit) break;
}
printf("\n %s %10.6f %s %10.6f j=%d","T=",T,"path len=",path,j);
printf("successful moves: %6d\n",nsucc);
for(i=1;i<=53;i++) printf("\%5d",iorder[i]);
T*=TFACTR;
if(nsucc==0) return;
}
}
float reverse_len(int ncity,int n[])
{
float de;
int xx[5],yy[5];
int j,ii;
n[3]=1+((n[1]+ncity-2)%ncity);
n[4]=1+(n[2]%ncity);
for(j=1;j<=4;j++) {
xx[j]=iorder[n[j]];
}
de=-weight[xx[1]][xx[3]];
de-=weight[xx[2]][xx[4]];
de+=weight[xx[1]][xx[4]];
de+=weight[xx[2]][xx[3]];
return de;
}
void reverse (int ncity,int n[])
{ int nn,j,k,l,itmp;
nn=(1+((n[2]-n[1]+ncity)%ncity))/2;
for(j=1;j<=nn;j++){ k=1+((n[1]+j-2)%ncity);
l=1+((n[2]-j+ncity)%ncity);
itmp=iorder[k];
iorder[k]=iorder[l];
iorder[l]=itmp;
}
}
float transe_len(int ncity,int n[])
{
float /*xx[7],yy[7],*/ de;
int j,ii,xx[7];
n[4]=1+(n[3]%ncity);
n[5]=1+((n[1]+ncity-2)%ncity);
n[6]=1+(n[2]%ncity);
for(j=1;j<=6;j++){
/*ii=iorder[n[j]];
xx[j]=x[ii];
yy[j]=y[ii];*/
xx[j]=iorder[n[j]];
}
de=-weight[xx[2]][xx[6]];
de-=weight[xx[1]][xx[5]];
de-=weight[xx[3]][xx[4]];
de+=weight[xx[1]][xx[3]];
de+=weight[xx[2]][xx[4]];
de+=weight[xx[5]][xx[6]];
return de;
}
void transe(int ncity,int n[])
{
int m1,m2,m3,nn,j,jj,jorder[100];
m1=1+((n[2]-n[1]+ncity) %ncity); /* n[1] to n[2] city numbers*/
m2=1+((n[5]-n[4]+ncity) %ncity); /*n[4] to n[5] */
m3=1+((n[3]-n[6]+ncity) %ncity); /*n[3].. n[6] .....*/
nn=1;
for(j=1;j<=m1;j++)
{jj=1+((j+n[1]-2)%ncity);
jorder[nn++]=iorder[jj];
}
if(m2>0) for(j=1;j<=m2;j++)
{ jj=1+((j+n[4]-2)%ncity);
jorder[nn++]=iorder[jj];
}
if(m3>0) for(j=1;j<=m3;j++){ jj=1+((j+n[6]-2)%ncity);
jorder[nn++]=iorder[jj];
}
for(j=1;j<=ncity;j++)
iorder[j]=jorder[j];
}
int metrop(float de,float t)
{
float ran(long *idum);
static long gljdum=1;
return de<0.0 ||ran(&gljdum)<exp(-de/t);
}
float ran(long *idum)
{static int inext,inextp;
static long ma[56];
static int iff=0;
long mj,mk;
int i,ii,k;
if(*idum<0||iff==0) {
iff=1;
mj=MSEED-(*idum<0?-*idum:*idum);
mj%=MBIG;
ma[55]=mj;
mk=1;
for(i=1;i<54;i++){
ii=(21*i)%55;
ma[ii]=mk;
mk=mj-mk;
if(mk<MZ)mk+=MBIG;
mj=ma[ii];
}
for(k=1;k<=4;k++)
for(i=1;i<=55;i++){
ma[i]-=ma[1+(i+30)%55];
if(ma[i]<MZ) ma[i]+=MBIG;
}
inext=0;
inextp=31;
*idum=1;
}
if(++inext==56)inext=1;
if(++inextp==56) inextp=1;
mj=ma[inext]-ma[inextp];
if(mj<MZ) mj+=MBIG;
ma[inext]=mj;
return mj*FAC;
}
int irbit(unsigned long *iseed)
{
unsigned long newbit;
newbit=(*iseed &131072)>>17
^(*iseed&16)>>4
^(*iseed &2)>>1
^(*iseed &1);
*iseed=(*iseed<<1)|newbit;
return newbit;
}
main()
{
int ncity=53;
int i,j,k,l;
for(i=0;i<=53;i++) iorder[i]=i;
for(i=1;i<=53;i++)
for(j=1;j<=53;j++)
weight[i][j]=0;
weight[53][1]=6;
weight[53][50]=10.1;
weight[1][38]=11.2;
weight[3][39]=8.2;
weight[5][39]=11.3;
weight[6][48]=9.5;
weight[8][40]=8;
weight[11][40]=14.2;
weight[12][42]=7.8;
weight[13][44]=16.4;
weight[15][44]=8.8;
weight[17][46]=9.8;
weight[19][20]=9.3;
weight[20][25]=6.5;
weight[21][46]=4.1;
weight[23][49]=7.9;
weight[25][49]=8.8;
weight[27][28]=7.9;
weight[29][51]=7.2;
weight[31][32]=8.1;
weight[33][35]=20.3;
weight[34][37]=17.6;
weight[44][45]=15.8;
weight[53][2]=9.2;
weight[53][52]=12.9;
weight[2][3]=4.8;
weight[4][8]=20.4;
weight[5][48]=11.4;
weight[7][39]=15.1;
weight[9][40]=7.8;
weight[11][42]=6.8;
weight[12][43]=10.2;
weight[13][45]=9.8;
weight[16][17]=6.8;
weight[18][44]=8.2;
weight[19][45]=8.1;
weight[20][47]=5.5;
weight[22][23]=10;
weight[24][27]=18.8;
weight[26][27]=7.8;
weight[28][50]=12.1;
weight[29][52]=7.9;
weight[31][33]=7.3;
weight[33][36]=7.4;
weight[36][37]=12.2;
weight[48][49]=14.2;
weight[53][38]=11.5;
weight[1][36]=10.3;
weight[2][5]=8.3;
weight[4][39]=12.7;
weight[6][7]=7.3;
weight[7][40]=7.2;
weight[9][41]=5.6;
weight[11][45]=13.2;
weight[13][14]=8.6;
weight[14][15]=15;
weight[16][44]=11.8;
weight[18][45]=8.2;
weight[19][47]=7.2;
weight[21][23]=9.1;
weight[22][46]=10.1;
weight[24][49]=13.2;
weight[26][49]=10.5;
weight[28][51]=8.3;
weight[30][32]=10.3;
weight[31][52]=9.2;
weight[34][35]=8.2;
weight[36][52]=8.8;
weight[53][48]=19.8;
weight[1][37]=5.9;
weight[3][38]=7.9;
weight[5][6]=9.7;
weight[6][47]=11.8;
weight[7][47]=14.5;
weight[10][41]=10.8;
weight[12][41]=12.2;
weight[13][42]=8.6;
weight[14][43]=9.9;
weight[17][22]=6.7;
weight[18][46]=9.2;
weight[20][21]=7.9;
weight[21][25]=7.8;
weight[23][24]=8.9;
weight[25][48]=12;
weight[26][50]=10.5;
weight[29][50]=15.2;
weight[30][51]=7.7;
weight[32][35]=14.9;
weight[34][36]=11.5;
weight[37][38]=11;
for(i=1;i<=53;i++)
for(j=1;j<=53;j++)
if(weight[i][j]!=0)weight[j][i]=weight[i][j];
for(i=1;i<=53;i++)
for(j=1;j<=53;j++)
if(weight[i][j]==0) weight[i][j]=MAX;
for(k=1;k<=53;k++)
for(i=1;i<=53;i++)
for(j=1;j<=53;j++)
if(!(i==j)&&!(i==k)&&!(k==j))
if( weight[i][k]+weight[k][j]<weight[i][j])
weight[i][j]=weight[i][k]+weight[k][j] ;
printf("\n");
anneal(ncity);
for(i=1;i<=53;i++) printf("\%5d",iorder[i]);
}
灾情巡视C语言代码的更多相关文章
- Latex中插入C语言代码
Latex是一个文本排版的语言,能排版出各种我们想要的效果.而且用代码排版的优点是易于修改板式,因此在文本内容的排版时,Latex应用十分广泛. 当我们需要在Latex中插入代码时,就需要用到 \us ...
- loadrnner添加C语言代码的几种方式
今天有人在群里问,想直接把自己的C语言代码让lr调用,该怎么搞. 这东西说来简单,只是对Loadrunner这工具不熟悉可能才会有这种问题吧.个人理解,一般有三种方法吧,废话不多,直接干货. 1.直接 ...
- HTML之一语言代码
HTML的lang属性可用于网页或部分网页的语言.这对搜索引擎和浏览器是有帮助的. 同时也可以是指HTTP Header中的Accept-Language/Content-Language. ISO ...
- 让你的Windows不断重启的C语言代码
原文:让你的Windows不断重启的C语言代码 没有写Linux的原因是因为搞不定Linux下的权限问题,而Windows下基本上使用电脑的用户都是管理员,所以钻个空了,不多说下面是代码#includ ...
- 国家语言,语言代码,locale id对应表
国家语言,语言代码,locale id对应表.比如 en_US对应的id为1033, 中文的locale=zh_CN,id=2052. Locale Languagecode LCIDstring L ...
- 如何提高单片机C语言代码效率
代码效率包括两个方面内容:代码的大小和代码执行速度.如果代码精简和执行速度快,我们就说这个代码效率高.一般情况下,代码精简了速度也相应提上来了.单片机的ROM和RAM的空间都很有限,当您编程时遇到单片 ...
- 使用highlight.js高亮静态页面的语言代码
显示静态的代码其实html的pre标签基本可以满足需求了,至少不会将换行的文本显示成一堆字符串. 不过能使静态的文本能高亮显示,倒更炫酷一点.其实很简单的,引入highlight.js包,可以使用cd ...
- HTML5 脚本 语言代码 URL 符号实体 ASCII码 颜色
1.HTML<noscript> 标签 <noscript> 标签提供无法使用脚本时的替代内容,比方在浏览器禁用脚本时,或浏览器不支持客户端脚本时. <noscript& ...
- 02 How to Write Go Code 如何编写go语言代码
How to Write Go Code 如何编写go语言代码 Introduction 介绍 Code organization 组织代码 Overview 概述 Workspaces ...
随机推荐
- 什么是LeapMotion
LeapMotion预览——什么是LeapMotion LeapMotion预览 这个就是LeapMotion: 原文转自: LeapMotion预览 LeapMotion 官网:http://l ...
- Machine Learning/Introducing Logistic Function
Machine Learning/Introducing Logistic Function 打算写点关于Machine Learning的东西, 正好也在cnBlogs上新开了这个博客, 也就更新在 ...
- Windows Store 手势编程小结
Windows Store 手势编程小结 最近完成了一个Windows Store上面的手势操作的页面.在这里总结了一下经验和心得,希望能和大家一起分享和讨论一下. 首先,要纠正一个误区,在Windo ...
- Organic Solar Cells - Theory and Practice
renewable engergy: Wind solar Tidal Wave Ocean Bio Etc. How much energy we gonna need as following y ...
- 压缩文件 compress files 以7z 格式及解压 或者别的格式
主要是为了能大量的减少文件使用空间,为了能节约带宽. 那么就用了7z的压缩方式. 这里,使用了7z的压缩方式,硬生生的将一个10k多的图片压缩成了3k左右的包.图片是不好压缩的,这个压缩比比zip g ...
- linux下php-5.4.8.tar.gz编译安装全攻略
首先安装基础依赖组建,注:这些依赖组建也是LINUX+PHP+MYSQL+APACHE+NGINX+MEMCACHED时必要的系统组件 LANG=C yum -y install gcc gcc-c ...
- qt的安装及连接sql使用注意
qt安装与使用 今天写下qt使用的注意事项,qt-project.org上已经更新了qt5的相关库,但是建议不是很熟系qt的朋友,或者说想使用传统qt api而不是qml的朋友, 继续下载老版本4.8 ...
- [置顶] JAVA从零单排4-----继承、封装和多态详解
继承 Java三大特征之一:继承.Java的继承具有单继承的特点,每个子类只能有一个直接父类. 继承的特点 Java的继承用extends关键字来实现,被继承的类成为父类,实现继承的类被称为子类.子类 ...
- 黑马程序员:Java基础总结----GUI图形化界面
黑马程序员:Java基础总结 GUI图形化界面 ASP.Net+Android+IO开发 . .Net培训 .期待与您交流! GUI(Graphical User Interface)图形化界 ...
- OC之类与对象
1.面向过程与面向对象. 1). 完成需求1 将大象放进冰箱. a. 把冰箱门打开. b. 把大象放进去. c. 把冰箱门关上. 这是面向过程的思路. 找1个冰箱,要求这个冰箱可以自己开门,自己把大象 ...