模拟退火算法解决工作指派问题c++源程序.doc

/*********************************************************************** 模拟退火算法解决工作指派问题 输入格式(工作时间)：矩阵 行i；工人 列j；工作 ***********************************************************************/ #include iostream #include cmath #include time.h #include stdio.h using namespace std; const int MAXN = 501; //最大约束数n const double INIT_T =2000; //初始温度 const double RATE = 0.95; //温度衰减率 const double FINAL_T = 1E-10; //终止温度 const int IN_LOOP = 13000; //内层循环次数 const int OUT_LOOP = 2000; //外层循环次数 struct work { //定义工作结构类型 int p[MAXN]; //工人序号 double time; //时间 }; int N; //工人数量 double D[MAXN][MAXN]; //工作时间 int x[MAXN][ MAXN]; //选择工人与工作 work bestwork; //最优的遍历路径 void swap(int a,int b) { int t; t=a;a=b;b=t; } inline double totaldist(work p) //计算总时间 { int i,j; double cost = 0,sumi=0, sumj=0; for (j=1; j=N; j++) { sumj+= x[i][j]; } for (i=1; i=N; i++) { sumi+= x[i][j]; } if(sumj==1 sumi==1) { for (j=1; j=N; j++) { for (i=1; i=N; i++) { cost += D[i][j]*x[i][j]; } } return cost; } } void init() //读入数据，并初始化 { cin N; int i,j; for(i=0;i=N;i++) { for(j=0;j=N;j++) cin D[i][j]; } for (i=j=1; i=j=N; i++,j++) //最优解的初始状态 { bestwork.p[i] = i; bestwork.time = totaldist(bestwork); } } /*********************************************************************** 产生指派方法 ***********************************************************************/ work getnext(work p) //新解产生函数 { int x, y,i,j; work ret; ret = p; do { x = rand() % N + 1; y = rand() % N + 1; }while(x == y); swap(ret.p[x], ret.p[y]); //交换两工人位置顺序 ret.time = totaldist(bestwork); return ret; } /*****************************************************************