通用对偶单纯形法C语言程序.doc

《通用对偶单纯形法C语言程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、#include#include#defineMAX_N100#defineM1000000.0;intm,n;floatA[MAX_N][MAX_N],C[MAX_N],b[MAX_N],seta[MAX_N],CZ[MAX_N];intnum[MAX_N];intGetOutCow()//求换出基{inti,k;intflag;floatmin=0;for(i=0;i=0)flag=1;else{flag=0;break;}}if(flag==1)r

2、eturn-1;for(i=0;ib[i]){min=b[i];k=i;}}returnk;}intGetInCow(intp)//求换入基{inti,j;intflag=0;floatmin;for(j=0;j=0)flag=1;else{flag=0;break;}}if(flag==1){printf("原线性规划问题无可行解！");return-1;}for(j=0;j

3、j]=CZ[j]/A[p][j];elseseta[j]=M;}min=M;for(j=0;j=seta[j]){min=seta[j];i=j;}}num[p]=i+1;returni;}voidchange(intp,intq)//计算新的单纯形表{inti,j;floattemp1,temp2,temp3;temp1=A[p][q];for(i=0;i

4、)A[i][j]=A[i][j]-A[p][j]*temp2;b[i]=b[i]-b[p]*temp2;}}}temp3=CZ[q]/temp1;for(i=0;i

5、n");printf("t");for(i=1;i<=n;i++)printf("X(%d)t",i);printf("RHS");printf("-----------------------------------------------------------------");printf("t");for(i=0;i

6、-------------");}voidprint2(){inti,j;printf("---------------------------------------------------------------");for(i=0;i

7、----------------------------------");printf("cj-zjt");for(i=0;i

8、%d",&n);printf("请输入方程组的系数矩阵A(%d行%d列):",m,n);for(i=0;i