求助：单纯形法的程序(c或其它语言都可)

lauqt1980612

急需啊，谁有啊。感激不尽啊!!!

lauqt1980612

最好是C语言的!

julysea

以前遍了一个小程序，可惜只能解某个特定问题。。。。

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
" ]/ @: r/ W% w) }#include&lt;math.h&gt;
0 Q* F" J2 r& Y#define X 5
#define Y 7
void xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
4 a8 Z/ P: }/ u& _{
1 l1 ~  ~( X; Q  int j;
& \2 Z" h, b1 t; R2 K  *c=0;
  for(j=1;j&lt;=*is;j++)
3 E4 V" U. X) O0 \% ~    {
! u' P( T" r$ V! I: ]. q* V4 z  x/ i% p      if((*a)[*ir][j]-*c&gt;0)
5 A" F- T. Q: ?! P- X# R    {
      *c=(*a)[*ir][j];
      *j0=j;
9 _! F% [) G* y& o( |    }
9 j: M) i! r$ v( z# X# J* C9 T9 V' n/ o     }
% Y& G- U; p. q( s1 q}</P><>      /***************** 参数说明 **********************/
      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
5 z+ B" B; w" q& E8 I! ^9 M- f: }; u      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
1 r: }$ H! V! h) N) H2 v3 U2 s      /********* <a href="http://happyyangxu.home.sunbo.net*******/" target="_blank" >http://happyyangxu.home.sunbo.net*******/</A>
int xi_sm(int m,int n,int m2,int mn1,int l1,float (*a)[X][Y],
- B  Y, N! I- P% O+ c% E6 b8 W       int(*k)[],float(*x)[])
# U* W6 O1 x' y; H2 Y3 \2 z{
$ R2 J: @* ~7 L/ d# d  int i,m1,mn,j0,i0,j;
  float c,g;
  m1=m+1;
  mn=m+n;
! `" F( q7 V7 V' ~/ m1 `) L' `) K  for(i=1;i&lt;=m;i++)
    (*k)=n+i;
8 t" x+ d8 }: j6 \( `leap1:
+ }. V" L9 S9 F# v4 ?: O  if(l1-1==0)
' x, \/ k% x: y: S* p- u     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  else
     {
leap2:   if(l1-50==0)
, W% Y7 O! J3 w) u      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
       else
) y/ M( ~. \& p- t0 M: _( o      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
8 m: j# K" F! H0 y  r     }
$ Q/ P& B# `- C! z   c-=1e-8;
   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
6 }/ O9 R' q1 e' e5 V9 J      {
9 x9 n, w0 o* q# e9 g) }2 i    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
$ \: q, g% j7 F9 i8 D4 y$ D    return(-1);
* T7 R" h, o- V( j, c      }
   if((c&lt;0)&amp;&amp;(l1==1))
( Y* Y8 q7 u5 {6 ?, d" Q       {
     l1=50;
2 U2 z. o# z8 j     goto leap2;
       }
    if(c&lt;=0)
       {
) P2 Z) S- P) F/ `      for(i=1;i&lt;=m;i++)
         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
* c$ Q4 h7 }2 W      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
" C, j9 m+ o/ j4 N% N3 @- }      printf("\nF=%f\n",(*a)[m2][mn1]);
% e  l4 j+ M0 C" k      return 0;
; L" ^2 |$ W3 Y6 a) Q    }
0 ]6 c# m. g7 i3 }& O- u6 u     c=1e8;
     for(i=1;i&lt;=m;i++)
       {
     if((*a)[j0]&gt;1e-8)
        {
          g=(*a)[mn1]/(*a)[j0];
! ]$ l: M# B) S. D: J          if(g-c&lt;0);
- J0 h3 }8 h' m( H% T- ^         {
           c=g;
5 z/ q# V* @( _; y4 ~           i0=i;
          }
         }
    }
      if(c==1e+8)
) T1 }, F5 f# U    {
6 T! R; Y" j; v* U  u      printf("\n\t*******Lp no solution********\n");
      return -3;
3 ?6 o8 e& ]8 ?/ A/ B     }
, u" d; [. x# k1 e7 D+ S' O       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
0 I& x0 b/ E, T0 r: n6 ]/ f      {
        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
- Z$ S7 I( I6 E8 N! D2 F7 I1 V           continue;
* p$ B+ o4 f; [  v         g=(*a)[i0][j]/(*a)[i0][j0];
           (*a)[i0][j]=g;
+ E# g; d+ X$ w  U         for(i=1;i&lt;=m2;i++)
           {
         if((i==i0)||(!(l1==1)&amp;(i==m1)))
7 S  `3 Z5 F% a: @; d2 l- h            continue;
         (*a)[j]=(*a)[j]-(*a)[j0]*g;
% o7 f4 [5 T6 @9 M# R( c7 A           }
       }
       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
       (*a)[i0][j0]=1;
+ c3 w& _" c( ?. P       goto leap1;
, V6 ]$ n8 [7 b- Z& _9 G# e}  
main()
+ E3 S: s  ^) [; l. m{
  float a[5][7]={{0,0,0, 0,  0,0,0},
7 K# ~/ @$ E8 ]3 r1 T& ^0 l; b+ F                 {0,3,-4,3,  1,0,12},
                 {0,3,0, 6,  0,1,12},
                 {0,0,0, 0,  0,0,0},
: T( l* I# y9 x9 O( s- R# h                 {0,69,0,144,0,0,300}};
  int k[3];
+ i$ r; p7 d3 Y0 E  float x[3];
7 L0 I; b6 z1 X- H3 G* [  clrscr();
  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
( J7 ]- G- S# g' u  a9 E. w}
) e2 h/ @( W0 N* `# z1 g4 N) Z</P>

plgatc

<>/*************************************************************************
                    单纯型法解线性规划问题（两阶段法）        

编程环境：VC++6.0     
方程组输入说明：
变量非负，按提示输入相关参数。
! g/ a5 ?2 Z8 V. I" D% W3 ] *************************************************************************/
3 m3 W4 G# {% I6 G  b  P#include &lt;stdio.h&gt;
#include &lt;stdlib.h&gt;
#define MAX 100
8 V% h2 t+ L( f, i#define STP 100</P>
% F0 b; ^; t7 b& |& g% F<>int stop=1; //迭代记数变量
; Z/ I$ x" `* z% [8 j7 l1 g3 G2 iint status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
int step=1; //目前阶段</P>
<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
" {2 P( ^/ S/ O# e* Uint num_x; //变量个数
) V% P$ K) i& Z* t: R9 e' Dint num_st; //约束方程数
int num_ar=0; //人工变量个数
! j; S+ _& N% h- w. O/ O0 m0 a* iint arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
" @: R* ?# @! o# o/ @<>void create(); //建立方程组
void iterative(); //单纯型法迭代
void output(); //输出结果
void banner(); //打印程序标题
void exchange(); //交换两阶段价值系数
void show(); //输出方程组</P>
<>void main() {
int i,j,k;
# z# T( l* ]/ ^( N banner();
create();
//保存原价值系数，转换为第一阶段价值系数
' h0 x0 R! L# z for(i=1;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
: K' h, h, J2 r5 M$ }  if(k==1) temp_c=-1;
4 d3 i$ z: L& S; @0 g  else temp_c=0;
( `3 e( @4 v8 `& G8 a0 ^; G }
+ ~( Q! I) @& k" X exchange(c,temp_c);</P>
<> printf("\n\n第一阶段问题为：\n\n");
show();
6 Z' ?; K) J2 E% M1 ? step++;
+ V8 Z2 {' D1 c) P9 g# p printf("\n\n按回车开始第一阶段迭代");
getchar();
9 R8 T' g- {; a( Q- t6 C7 k; b1 w getchar();
* s0 p, |+ t( A; `1 }- z iterative();
. w8 x4 Z  v: {& ]+ ?7 y' b if(status==-2) {
  puts("迭代超过限制次数强行终止！\n");
$ l* W, F3 G3 b& C) v  puts("\n按回车结束");
1 y& j1 Q& I- [' _  ]  G  getchar();
/ J3 ~9 k" v- e  exit(0);
# O! ~( n) G! f: m$ P8 F: [ }
7 l' o% h; W* S5 } output();</P>
* v+ R7 ?- ?. P( K- h' T<> if(max!=0) {
8 u* A  A$ e! |( S7 o' E+ }  puts("\n\n原问题无可行解。\n");
9 n, x9 H7 Q7 ?! b2 B0 e  puts("\n按回车结束");
6 P8 }: [4 i' {$ }  getchar();
  exit(0);
# \+ Z+ c2 c3 ~4 @+ q! s0 S }</P>
<> //转换为第二阶段价值系数
! c; H+ R% n3 h! ^ exchange(c,temp_c);
//把人工变量列全设为0
, o0 V# G8 X, E# ~/ B for(i=1;i&lt;=num_ar;i++) {
* a, u  w4 v. u  J) b4 U  c[arti]=0;
" _0 t3 U* S% n% r" f  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
}</P>
0 {( O& u. n8 I; m9 m* b5 J0 G<> puts("\n\n第二阶段问题为：\n\n");
. D# r( h5 b! [! o/ F$ \ show();
$ d3 z! Q- V7 e; b9 r puts("\n\n按回车开始第二阶段迭代");
getchar();
iterative();
switch(status) {
  G: |, q! M6 C' e! |1 G case 1:
  output();
  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
, U- L( b+ \2 q3 t# x, {  getchar();
: n2 Q; V; o: X' ?; K  exit(0);
case 0:
  puts("\n\n原问题为无界解。\n");
! p$ b9 Q# \0 C3 }1 d# L  puts("\n按回车结束");
  getchar();
8 F) J6 Q( s' A  @9 f  exit(0);
case -1:
% ?) I" S2 o5 Z. T  output();
  puts("\n\n原问题有无穷多最优解。\n");
- s: ]8 r* \4 E/ |% F) q  puts("\n按回车结束");
  getchar();
  exit(0);
case -2:
, W' e* g' e% w, Y1 i4 h) J  puts("迭代超过限制次数强行终止！\n");
  M& e  r+ s/ J4 u( N* n# u  puts("\n按回车结束");
  getchar();
  exit(0);
}//switch
}</P>
<>void banner() {
5 T1 ]! \9 i' B. F8 q, J* L: M printf("\t\t****************************************\n");
! K. B( W4 P. d; }! @, t! P printf("\t\t         单纯型法解线性规划问题\n");
printf("\t\t                         作者：Thunder\n");
; `; c+ x3 O- Z. A% A" M& _! A7 R printf("\t\t****************************************\n");
printf("\n");
4 l7 A$ ?6 L! h9 _/ Q' |$ c}</P>
<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
! I8 \. W$ _4 r( s2 K//为1的不显示系数1，-1系数只显示负号
0 {( G; J. ]8 g2 h2 K) W int i,j,k;
' F* X6 W8 k( w switch(step) {
7 v. @! r/ A+ w/ `# y- ^9 d/ { case 1:
  printf("min z= ");
  printf("x[%d]",arti[1]);
9 \: @  J: k6 C0 x3 G  ?  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
7 _1 Y! u4 K. [4 M4 W  break;
" ?8 h7 h6 _, x# B! x& D case 2:
+ P: y5 U, l) F9 V printf("max z= ");
printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
5 ?7 o4 r- B4 ~- o6 m( a  if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
4 [7 y, p, s. i. y- Q  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
  else printf(" %lg x[%d]",c,i);
* L( e7 x: l; p, X3 L8 j( ~ }
break;
. M% D7 v( l. k! a9 i; d4 p* Q3 ` }</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
2 d) @: ~# K0 y5 j   if(a[j]!=0) {
8 P& ~0 m4 }( g! i0 }0 m9 L    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
2 ~; E5 G- B  X# `$ ?8 B' {( s    else if(a[j]==-1) printf(" - x[%d]",j);
    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
! ]# X* P' T2 C! N( t' k    else printf(" %lg x[%d]",a[j],j);
$ |& Z% S0 B& n+ p" ~    k=1;
   }
) f7 G. z# Q" i" z; h  }
  printf(" == %lg\n",b);
% C+ T- J- m8 ~# z) I. e- y" J( _. R }
( P8 L/ C9 f# F7 N+ f& k printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
<>void exchange() {
int i;
2 V' @$ y. h8 b- T9 g' E% m4 h double temp[MAX];
for(i=1;i&lt;=num_x;i++) {
1 {6 G- T+ A! @3 P  temp=temp_c;
  temp_c=c;
  c=temp;
: O( L  H& H  K4 L }
}</P>
<>void create() {
% I# ^6 N' {: g+ Y( F//输入方程组系数，每个方程输完后回显确认
int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
9 |$ }' d4 v# p  S6 N* F char confirm;
. s; z0 H" y$ c: Z; F) }: s
! u8 A' U* v. z& Q9 X8 Z9 J0 G7 t8 B while(1) {
  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  scanf("%d",&amp;ma_mi);
  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
  else break;
0 P0 j3 t$ Q: w4 M2 b }
$ U3 z  m  i1 i9 h( b9 ~
; U% p3 j9 [) `, R while(1) {
7 T. \) Z, ^( N: e) G1 B  printf("指定变量个数：");
  scanf("%d",&amp;num_x);
) ~# J" C$ k, l# g  printf("输入价值系数c1-c%d:\n",num_x);
  for(i=1;i&lt;=num_x;i++) {   
/ R; E$ z2 c, h7 a% D$ p# W   printf("c%d=",i);
& }  g) ~$ D9 D2 s* w   scanf("%lf",&amp;c);
  }
2 l/ e+ \) Q8 ]% Q4 m3 |  if(ma_mi==1) printf("max z= ");
9 f3 Q. O2 m0 n! N2 O: Q2 G  else printf("min z= ");
  printf("%lg x[%d]",c[1],1);
  Y* z) ^( G6 i0 m% @; F  a  for(i=2;i&lt;=num_x;i++) {
   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
   else printf("  %lg x[%d]",c,i);
0 Q& `; S& q- r% F  }
printf("\n正确吗？：(y/n)");
getchar();
confirm=getchar();
: I; d0 A2 l& S5 w if (confirm=='y') break;
; u& Q4 @8 ~- p else if(confirm=='n') continue;
}</P>
<> printf("输入约束方程组个数：");
* r7 p5 [6 a  C3 o) u scanf("%d",&amp;num_st);
! V" I2 v! N$ Z( U for(i=1;i&lt;=num_st;i++) {
4 i( ^: c5 ~3 ^! m0 o! E. t  printf("st.%d:\n",i);
1 q5 y( K6 Y$ W; g  while(1) {
   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
   scanf("%d",&amp;re_st);
% J7 E$ Y% P9 ~   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
4 O. y- p) A# I# m   else break;
  }
9 }. [  w7 U, `7 I$ p7 g  E( F  printf("输入技术系数：\n");
  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
- e( U! X2 e- {( }3 y9 A   scanf("%lf",&amp;a[j]);
  }
4 |' }' J$ Q5 u) g- ^; r  printf("输入资源拥有量：\nb%d=",i);
  scanf("%lf",&amp;b);
  
) l* H; G- l$ v# X  printf("st.%i:\n",i);
: O1 Z+ v7 s1 z) x- r% B  printf("%lg x[%d]",a[1],1);
! W  L' T+ n0 {- d( N  for(j=2;j&lt;=num_x;j++) {
   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
3 d8 a) U7 a! T0 p( m! v   else printf(" %lg x[%d]",a[j],j);
& J6 k6 g8 s* k( q! d  }
  switch(re_st) {
   case 1: printf(" == %lg",b); break;
* U: C; D* u/ V, V) j% v& I   case 2: printf(" &gt;= %lg",b); break;
5 v4 O0 P6 e( E: ^" _% O# H   case 3: printf(" &lt;= %lg",b); break;
  }</P>
% U% N, w* A% Y9 Q% @' ]<>  while(1) {
) y7 O' e6 Y9 O2 x. o   printf("\n正确吗？(y/n)");
  b! w- u/ {  }- ?! X   getchar();
   confirm=getchar();
6 u& F! A5 N$ z# j   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
! O- V1 q2 i% ]: @& k  }
}</P>
( ]3 z; E/ Q( w<>//显示输入的方程组
4 B+ r% X1 ?6 o+ W printf("\n原问题为：\n\n");
if(ma_mi==1) printf("max z= ");
else printf("min z= ");
printf("%lg x[%d]",c[1],1);
. v: Y* U3 ]! _% Y3 ` for(i=2;i&lt;=num_x;i++) {
3 N7 }1 H+ w) B6 b' W  if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
( v  b% p8 o2 q* J  else printf(" %lg x[%d]",c,i);
}</P>
<> printf("\nst:\n");
+ Y. L& {: D5 d: S for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
9 a6 E0 a7 m. p& c: [   if(a[j]!=0) {
: H! u# g4 @/ {  K8 v! ?    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
' J$ m" ?! j, _9 h, H6 X    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
    else if(a[j]==-1) printf(" - x[%d]",j);
' v  J- l  g& |3 _4 w# n5 h+ S    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
& V7 M) @) G* D    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
    else printf(" %lg x[%d]",a[j],j);
3 n% u; P& }5 K  V+ K& L. g    k=1;
* \+ c; u1 y3 o; X4 Q/ d) D7 `) f) R   }
  }
% d9 @" V- A2 C- W: ?6 s+ k- R7 r) b  switch(re_st) {
   case 1:
    printf(" == %lg\n",b);
    break;
   case 2:
    printf(" &gt;= %lg\n",b);
    break;
5 p1 W0 T0 N0 I# ~   case 3:
1 e0 T  i4 ?! P) b1 h9 |9 a' _7 F    printf(" &lt;= %lg\n",b);
% O! Y8 u3 ], N1 Y4 \+ E    break;
  }
}
/ n; d6 h5 Z/ t  H  @& w+ B1 r printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
<> tnum_x=num_x;
* t6 ]( x6 w, l" ^ for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
  case 1:
( [( [# ]9 f3 `6 T$ ~8 r  S  case 3:
7 Z1 o& S2 |. p, B5 l% o7 }   num_x+=1;
! r9 u/ H+ ~+ Y' Y+ A   break;
  case 2:
   num_x+=2;
; w5 E, f& Q; H0 H) S1 t$ V; l   break;
' |' u3 f  q; V  H- p  }
}</P>
<>//化为标准形式
6 W: K  N1 `6 s5 g# ~ if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
4 N1 z4 ?, [1 y, G6 i1 H: z for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
   case 1:
) O% x' t3 A) k    num_addv++;
+ X8 D% C1 C" v1 A1 A    num_ba++;
    num_ar++;
    c[tnum_x+num_addv]=0;
    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
* Q: y! X3 ?) `$ l0 P, h4 v    for(j=tnum_x+1;j&lt;=num_x;j++)
     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
+ j; D% y+ F# M" \6 M    break;
   case 2:
    num_addv++;
    c[tnum_x+num_addv]=0;
    num_addv++;
% q" ~6 M9 y; y    num_ba++;
' S% L# W* G% u0 d) b1 M    num_ar++;
* j  N* N5 r' M* K- v& B    c[tnum_x+num_addv]=0;
    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
    for(j=tnum_x+1;j&lt;=num_x;j++)
6 T* V; A) P$ }& i) p     if(j==tnum_x+num_addv-1) a[tnum_x+num_addv-1]=-1;
/ V) Z8 _& n- M) o( J8 H     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
7 A5 q5 }0 R6 Y4 J     else a[j]=0;
    break;
   case 3:
, d- W6 r$ D+ f5 g+ {/ t% J    num_addv++;
    num_ba++;
    c[tnum_x+num_addv]=0;
    base[num_ba]=tnum_x+num_addv;
    for(j=tnum_x+1;j&lt;=num_x;j++)
' l7 u1 @8 b* e$ u, Y' m3 Z     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
, j: K" l- y2 C/ f4 M+ I% V     else a[j]=0;
+ k* Y2 U0 a" X' ^" q5 |5 q    break;
) `$ Y. C( ]# Z6 E& D0 @) R  }//switch
& u5 W# O: ?6 } }//增加松弛变量、剩余变量、人工变量、确定基变量</P>
<>//显示标准化后的方程组
& O0 d6 j: t6 n; y% I( o printf("\n化为标准形式后：\n\n");
if(ma_mi==1) printf("max z= ");
3 C$ P1 g& A0 e% C1 K/ `/ A0 m* r else printf("max z'= ");
printf("%lg x[%d]",c[1],1);
' v" D& U, S  f0 g# w3 ]5 Q6 W2 k for(i=2;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_ar;j++)
8 ~5 P+ Z, ^; r9 `9 U  if(i==arti[j]) k=1;
+ c  b* ~# H5 u, @1 B+ B6 W- w  if(k==1) printf(" -M x[%d]",i);
  else if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
  else printf(" %lg x[%d]",c,i);
}</P>
& `6 C, a. r2 k0 S' P' f6 N& K, E1 S<> printf("\nst:\n");
! V0 \* v; t) s  Z) K for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
7 {2 F& |% }+ A# R9 [5 s   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
# [2 p" W2 V  @& z0 U6 @' m    else if(a[j]==-1) printf(" - x[%d]",j);
    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
/ R* ^% I' G' D8 B! ?    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
    else printf(" %lg x[%d]",a[j],j);
1 _% Y- p5 ]: c0 }& m    k=1;
$ F7 {# |8 O$ H. E! j& x, c   }
  }
; G* }2 u1 `  Z. ~( c  printf(" == %lg\n",b);
" M7 V2 N- U, Y- E5 d }
, k1 ~  H1 i0 S/ m printf(" x[1]~x[%d]&gt;=0",num_x);
. p6 z2 o# L5 S9 t$ d+ V( i- A% O}</P>
<>void iterative() {
! u' K7 N9 O" ^7 C4 u2 | int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
int base_elem;
7 {  D1 ~9 \9 y( }1 d5 X int base_out,base_in;
double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
<> printf("\n\n第%d次迭代：\n\n",stop);
for(i=1;i&lt;=num_st;i++) {
1 o) o  ^$ U( c% h. S, t  printf("c%d=%lg\t",base,c[base]);
3 x! q) s: V; y8 E+ N  printf("b%d=%lg\t",i,b);</P>
<>  switch(step) {
   case 1:
, [  R) {: L9 p. [4 f; A$ L& f/ j    for(j=1;j&lt;=num_x;j++)
    {
- n7 n/ }2 N- g     printf("a[%d][%d]=%lg\t",i,j,a[j]);
6 q- o1 J& y4 O. {7 {    }
" U6 U! s4 Q1 t0 V2 J    printf("\n");
    break;
' U2 Q- r/ I3 b" T( |   case 2:
    for(j=1;j&lt;=num_x;j++) {
     k_a=0;
     for(l=1;l&lt;=num_ar;l++) if(j==arti[l])k_a=1;
0 T$ z5 K6 F. d2 A! ], K     if(k_a!=1) printf("a[%d][%d]=%lg\t",i,j,a[j]);
    }
    printf("\n");
/ X1 B( A6 Z# L+ o7 l& v    break;
$ ]3 q- L% @9 M+ O' O, y( k  }
3 f$ |% q3 x# I0 X }
//求检验数sigma
: n  B" _2 t+ }1 s- J: E1 [ for(i=1;i&lt;=num_x;i++) {
7 `1 K( O2 `1 `5 D% ?  sigma=c;
  for(j=1;j&lt;=num_st;j++) sigma-=c[base[j]]*a[j];
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
) t4 u0 o7 G( k4 [, _* g  switch(step) {
   case 1:
    printf("sigma[%d]=%lg\t",i,sigma);
    break;
4 f" Z. G3 t! W# @7 h   case 2:
    k_a=0;
& \! x- E, f: I' {" J5 c    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
4 E% h" r" I$ T    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
    break;
$ N0 [5 H$ l; N& l1 T  J$ O  }
/ Z- d: T/ f' j. N" W }
putchar('\n');
3 S" ^# X0 g0 m- u5 }//检验检验数sigma是否全小于等于0
k=0;
for(i=1;i&lt;=num_x;i++) {
" n: N7 L+ k. W) ~& ^5 F9 ^9 Y  if(sigma&gt;0)
1 p- {; ^; ?: h' ?  k=1;
}
if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
# v6 j5 B$ X4 K  for(i=1;i&lt;=num_x;i++) {
   k_f=k_a=0;
/ b3 `6 k2 l$ F   for(j=1;j&lt;=num_ar;j++)
; ^& H! Q4 _) T  l: Z- T   if(i==arti[j]) k_a=1;
   if(sigma==0&amp;&amp;k_a!=1) {
/ M  d! a) J% \7 i- v    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
   if(k_f==0) {status=-1; return;}
% t8 H3 ?$ |; c   }
  }
9 K) T  P! n* |' w& k- g  status=1;
  r5 ]) i; f7 d  return;
}
//检查是否为无界解
for(i=1;i&lt;=num_x;i++) {
  k_f=0;
  if(sigma&gt;0) {
   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
   if(k_f!=1) {status=0; return;}
& W/ P9 _+ L% X2 S9 h' f* c  }
}</P>
<>//确定换入变量
for(i=1;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  |* A9 r0 w9 b4 d5 R8 }  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
}//temp赋初值
for(i=1;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  if(k==0)
   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
    base_in=i;
    temp=sigma;
   }
3 _/ @& a# W- G" a' ]! c2 o& A% M }</P>
( f* D" r+ I7 r  D<>//确定换出变量
( ^# i: N6 K3 x% e( u for(i=1;i&lt;=num_st;i++)
( M# K3 r- j+ C  if(a[base_in]&gt;0) {
1 p! J, l. _; g8 P   temp=b/a[base_in]+1;
   break;
, `4 [' k, r9 g  |* l' Q& T  }//temp赋初值
; G; F# M' V% L7 j3 C. E for(i=1;i&lt;=num_st;i++) {
  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
' Z0 P, O; B8 `9 [7 W   for(j=1;j&lt;=num_ar;j++)
    if(base==arti[j]) {
     base_out=base;
     base_elem=i;
& q, H$ p% s0 ^* ]     temp=b/a[base_in];
8 E7 |- r2 z+ R( p     break;
    }
: ^3 s* L8 ^, V( \! G& ], A4 u5 {  }//人工变量优先换出
  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
   base_out=base;
1 r) K; `# {5 j5 [* @   base_elem=i;
   temp=b/a[base_in];
  }
}</P>
<> printf(" 基变量：");
$ L, R; k) A9 i for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
4 k4 V  H+ v2 [; g8 l6 D+ t//基变量变换，进行新方程初始化后迭代
' x  I. H! f+ @- o" n for(i=1;i&lt;=num_st;i++) {
$ t3 m& }6 E1 ]  P0 F  j9 s2 |  if(base==base_out) base=base_in;
- W% a1 B: [# h0 }3 {5 K }
//初始化主元素行系数
value_be=a[base_elem][base_in];
9 d: E. {  @2 s4 X  r/ B b[base_elem]/=value_be;
# C  j/ l) L' {8 r& [5 u1 e for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
. l5 i2 h4 j( s* w3 M. N<> for(i=1;i&lt;=num_st;i++) {
2 ^( v. f2 ^+ @  if(i!=base_elem) {
0 ?  a6 a' {; D* t/ e/ C: V   b-=b[base_elem]*a[base_in];
* J& d6 Y# M( l* t   value_be=a[base_in];
   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  }
}
: W# Z* \- _# T& c stop++;
if(stop&gt;STP) {status=-2; return;}
iterative();
}</P>
% A/ |# `0 P7 N* _<>void output() {
8 ]; B3 I/ g& {) n+ U1 w% T int i,j;
( l- d0 l) A. A" r5 h, Z double X[MAX];
printf("\n结果如下：\n");
printf("\nX=(");
: C& k1 o$ J9 w! e* E  w for(i=1;i&lt;=num_x;i++) {
  for(j=1;j&lt;=num_st;j++)
. h, Q1 {- t8 E1 O  R& l2 ~  if(i==base[j]) {X=b[j];break;}
: d, F  Q3 e+ J# O  else X=0;
  V0 V/ [, [1 @) [  e* U% n  printf("%lg ",X);
}
printf(")");
# ?6 ?9 C- T) s4 b for(i=1;i&lt;=num_x;i++) max+=c*X;
! b2 e+ y8 s; u3 h0 n if(ma_mi==1) printf("\nMax z= %lf\n",max);
' i; c: @0 u8 _$ } else printf("\nMin z= %lf\n",-max);
1 x3 \) h* c  b  |, |}</P>

sg47

好详细！

yqm10507

<>very good.Well  done.</P>

tw1982

先谢谢了，我看看能不能用

tw1982

我用.net要出错，能不能提供一个.net的啊

winzipftp

是你吗？thunder？帅呆了