整数规划隐枚举法实验报告.docVIP

实验报告 名 称： 运筹学上机实习 一、实验的目的与要求 1、实现整数规划隐枚举法的计算机实现实验（本实验是通过c++编程实现的）。 2、实验成果要求 基本要求：能够完成整数规划隐枚举法的计算（给出的是标准形式的整数规划问题）。 提高要求：求解分成两种情况：（1）求最大值max。（2）求最小值min。 二、实验分析 由实验要求知，需要在运行界面里输入一个整数规划的标准型（包括系数矩阵，目标函数矩阵，b列值，变量x的取值矩阵），然后运行界面会返回输出结果。 实验实现过程： #includeiostream.h #includemath.h /*本程序为整数规划隐枚举法*/ void main() { float a[30],c[30][30],b[30],c_b[30]; int x[30][100],m,n,i,j,p=0,shuru,e,f; //a[]为目标系数，c[][]系数矩阵，c_b[]判是否符合要求的数列，b[]为b列值。 //x[][]为变量的取值，m行，n列。其中m行，n列是系数矩阵c[][]的行列值。 //shuru为选择变量(目标函数求最大max值输入1,目标函数求最小min值输入2)。 //本方法要将x1,x2,x3...的取值(0或1)以纵向形式(二进制形式)输入。 cout/*endl本程序为整数规划隐枚举法,endl; cout本设计方法要将x1,x2,x3,x4...的取值(0或1)以纵向形式(且为二进制形式)输入。endl*/endlendl; cout请输入系数矩阵c[][]的行数m为:\n; cinm; cout请输入系数矩阵c[][]的列数n为:\n; cinn; cout请输入目标函数a[]的系数:\n; for(i=0;in;i++) cina[i]; cout请输入系数矩阵c[][]:\n; for(i=0;im;i++) for(j=0;jn;j++) cinc[i][j]; coutb列值:\n; for(i=0;im;i++) cinb[i]; cout请输入变量的取值(变量x1,x2,x3......只能取整数0和1,且以二进制的形式)\n; cout以纵向形式写出x[][]矩阵的取值\n; cout即:(x1,x2,x3...)T (x1,x2,x3...)T (x1,x2,x3...)T (x1,x2,x3...)T...:\n; cout(其中T表示矩阵的转置)endl; for(i=0;in;i++) for(j=0;jpow(2,n);j++) cinx[i][j]; for(i=0;im;i++) { c_b[i]=0; } cout若目标函数求最大max值请输入1\n; cout若目标函数求最小min值请输入2\n; cinshuru; float z=0.0; //目标函数求最大max值的情况 if(shuru==1) { //用试探法找初始可行解。 cout目标函数求最大max值endl; for(i=0;im;i++) for(j=0;jn;j++) c_b[i]+=c[i][j]*x[j][p]; int t,s; for(i=0;im;i++) { if(c_b[i]b[i]p=pow(2,n)-1) {p++; for(s=0;sm;s++) { c_b[s]=0; } for(t=0;tm;t++) for(j=0;jn;j++) c_b[t]+=c[t][j]*x[j][p]; i=-1; } if(ppow(2,n)-1) { cout无可行解\n; break; } } //产生过滤条件。 for(j=0;jn;j++) z+=a[j]*x[j][p]; //改进过滤条件，求最优解的过程。 float temp; int l,w; e=p; f=p; while(p=pow(2,n)-1) { temp=0.0; for(j=0;jn;j++) { temp+=a[j]*x[j][p]; } if(tempz) { p++; for(w=0;wm;w++) { c_b[w]=0; } for(t=0;tm;t++) for(s=0;sn;s++) c_b[t]+=c[t][s]*x[s][p