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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
2 o3 A, U: ~+ ^; c) J0 z9 I% r) C; z) L#include&lt;math.h&gt;
#define X 5
#define Y 7
5 e( [. l& Z+ S, X$ ?1 Dvoid xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
) a; L3 H+ R! i{
1 h  C2 i( G  d# R5 k  int j;
  *c=0;
  for(j=1;j&lt;=*is;j++)
, A+ B8 U' J6 A    {
      if((*a)[*ir][j]-*c&gt;0)
5 z$ R% b$ }0 F) r0 ~5 z3 z' Q    {
      *c=(*a)[*ir][j];
      *j0=j;
    }
     }
}</P><>      /***************** 参数说明 **********************/
8 a; U5 A7 u. x      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
" F. u& H  t4 {  B      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
) i4 _+ p5 x2 }% U5 e% O% o/ s/ t      /********* <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],
' w- `/ g/ {3 P* n1 c1 D       int(*k)[],float(*x)[])
; l. \% r! z3 h, \* C; y3 b{
  int i,m1,mn,j0,i0,j;
+ Y- c9 a4 z. t3 k  float c,g;
5 t8 t+ T; C6 a% I& N9 I  m1=m+1;
' E0 `: F2 r6 w/ @, F% `5 q  mn=m+n;
  for(i=1;i&lt;=m;i++)
2 F. ]7 o3 P: K# `5 A8 S    (*k)=n+i;
- D3 g" f5 b& ?1 K8 E6 n' u5 \( Xleap1:
$ a' a; Z6 E3 }2 D1 z( m' j/ c  if(l1-1==0)
# [. S9 ?( y1 a  ^' c! Z" o2 t     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  else
     {
leap2:   if(l1-50==0)
  @: p/ R3 T# N) r. T8 k2 y      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
       else
      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
     }
: {& X/ ~% j7 V$ h; [   c-=1e-8;
   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
      {
    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
    return(-1);
" F' w( y5 L; P0 y) f      }
! Y  |/ X8 I; G   if((c&lt;0)&amp;&amp;(l1==1))
       {
8 H* Q4 r& ~) P& m, i7 }9 R- V  \     l1=50;
     goto leap2;
       }
+ _& |3 v  H4 J8 R. ]    if(c&lt;=0)
       {
      for(i=1;i&lt;=m;i++)
$ r- @5 y% ]+ \$ b# f         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
      for(i=1;i&lt;=m;i++)
+ n! ^+ L& q; M0 @8 I, ]/ l  b7 u4 C         printf("\ni=%d,x=%f\n",(*k),(*x));
5 K& ]6 `& M  j      printf("\nF=%f\n",(*a)[m2][mn1]);
      return 0;
    }
     c=1e8;
" X1 p$ p! T; S+ P  I     for(i=1;i&lt;=m;i++)
       {
0 V( \; T& o& p  V- E2 T     if((*a)[j0]&gt;1e-8)
, l- `2 w5 `6 d( [( M. x        {
3 Y5 J: C/ w. {, Z          g=(*a)[mn1]/(*a)[j0];
          if(g-c&lt;0);
1 c$ M0 ?+ J7 b8 G- [9 }0 L         {
           c=g;
           i0=i;
8 e: w% A+ |( e% H          }
: C) {; w5 |) j* X% E" [$ k         }
    }
$ Q& F1 B$ E, X! `0 {) z      if(c==1e+8)
) y* O* R2 B+ U2 ], u    {
& y: r" g( K6 O  U9 O      printf("\n\t*******Lp no solution********\n");
      return -3;
& E/ h& B+ G; J0 D     }
       (*k)[i0]=j0;
# g1 y. O0 a6 g) Y9 ?       for(j=1;j&lt;=mn1;j++)
      {
# h7 ~4 b! U% d; E2 b' U        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
           continue;
" C  @" i& e& g- T' R$ p         g=(*a)[i0][j]/(*a)[i0][j0];
0 h; e+ M0 a* F: U8 D           (*a)[i0][j]=g;
# ^# S& a3 B+ |% E9 a- [, i- ?' z. f         for(i=1;i&lt;=m2;i++)
           {
         if((i==i0)||(!(l1==1)&amp;(i==m1)))
            continue;
- u8 s4 E) J: o         (*a)[j]=(*a)[j]-(*a)[j0]*g;
           }
       }
       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
       (*a)[i0][j0]=1;
       goto leap1;
}  
- l9 P5 b% H* t8 j- J5 ^* X7 Amain()
{
' P9 @5 J7 x3 H2 {% E) L  float a[5][7]={{0,0,0, 0,  0,0,0},
                 {0,3,-4,3,  1,0,12},
+ N; @: L2 d1 V                 {0,3,0, 6,  0,1,12},
                 {0,0,0, 0,  0,0,0},
                 {0,69,0,144,0,0,300}};
  int k[3];
  float x[3];
  clrscr();
" \" H8 ]) C9 @" m% b$ P  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
8 U2 {$ ?" F" t* c( {- L}
; {/ T  |6 l* D( s3 j" ^9 H8 ]$ v</P>

plgatc

<>/*************************************************************************
8 B7 L8 @& a) S                    单纯型法解线性规划问题（两阶段法）        
1 |6 G; T, k* P4 a
编程环境：VC++6.0     
方程组输入说明：
变量非负，按提示输入相关参数。
*************************************************************************/
#include &lt;stdio.h&gt;
#include &lt;stdlib.h&gt;
#define MAX 100
6 ?% U7 G: [" J# l: R0 b1 p# [( O; X#define STP 100</P>
<>int stop=1; //迭代记数变量
$ Q" L# J* j: R1 Z/ p0 J( jint status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
: J1 ]3 ^; A7 G8 e. u  S$ nint step=1; //目前阶段</P>
<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
int num_x; //变量个数
int num_st; //约束方程数
2 H; v9 `4 o( _. `# @int num_ar=0; //人工变量个数
int arti[MAX]; //人工变量下标
: U+ e* k3 l+ k: ^, q$ p+ Pint base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
<>void create(); //建立方程组
void iterative(); //单纯型法迭代
9 C* i- g, B+ S  B1 e' B. zvoid output(); //输出结果
) B% W+ Q8 u' U* \void banner(); //打印程序标题
void exchange(); //交换两阶段价值系数
void show(); //输出方程组</P>
<>void main() {
3 ?- A; b3 {2 r: E! |, Z int i,j,k;
% [! ^) m/ w) m; V' O$ M banner();
  s$ N* s. A  @0 P create();
+ p* l$ n( q* @9 i //保存原价值系数，转换为第一阶段价值系数
7 Y7 }. S3 z( _8 Y) |0 v% L for(i=1;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
  if(k==1) temp_c=-1;
  else temp_c=0;
}
exchange(c,temp_c);</P>
7 o$ T5 ~0 F6 V3 f' P! W  F& T<> printf("\n\n第一阶段问题为：\n\n");
# r" c3 X0 Y/ g  A" E+ ` show();
step++;
printf("\n\n按回车开始第一阶段迭代");
getchar();
getchar();
iterative();
. b  o% B4 O& }( j; h if(status==-2) {
; ^! W! q. s' P  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
  getchar();
  R0 R) E, O7 }3 D2 h: u  u' _  exit(0);
}
output();</P>
  p& a: r2 Z3 {8 n+ U7 x# ]7 ~4 W<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
4 J1 X; X" Y% L0 D3 }1 @  puts("\n按回车结束");
- ^- ]+ E; _& W$ P  getchar();
  exit(0);
- E7 q1 F. E5 X4 F; O: P) B }</P>
<> //转换为第二阶段价值系数
; U$ T5 {( x. M+ Y9 {1 @ exchange(c,temp_c);
+ o1 o: g/ {: R6 r* b! E //把人工变量列全设为0
for(i=1;i&lt;=num_ar;i++) {
) l# v( n: H5 o  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
. b1 O! Y  }7 @, c2 t7 q4 c }</P>
4 J/ s5 N* F+ o% A0 ?<> puts("\n\n第二阶段问题为：\n\n");
show();
  L* y* v8 i4 T$ D+ R! \ puts("\n\n按回车开始第二阶段迭代");
% M5 ?$ d3 C' M, Z: A$ l+ H8 A getchar();
. g" v9 d, \7 |) u iterative();
: [- |6 W  }( Q- n switch(status) {
case 1:
  output();
  puts("\n\n原问题有唯一最优解。\n");
9 S" a6 W5 J% Q7 I" L  puts("\n按回车结束");
  getchar();
  exit(0);
case 0:
  F; Y% ?  ~- s' I8 O, q  puts("\n\n原问题为无界解。\n");
: q, M: Y: r0 B4 j  puts("\n按回车结束");
  getchar();
& K% s/ o" T( X0 w, Z$ {  exit(0);
, `0 _8 O. X. g; o case -1:
  output();
  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
+ r5 ]# L5 c1 A  getchar();
  exit(0);
  o, N& _% x0 C5 X/ O/ O* b case -2:
; O; x9 c% {2 {6 Y0 ~. U5 M  puts("迭代超过限制次数强行终止！\n");
, v) T# {, t/ Q% ^: Z5 n  puts("\n按回车结束");
6 w/ G- w! g% e+ g0 r7 R  getchar();
& x4 J7 K# y9 y2 p8 w2 o% R  exit(0);
}//switch
}</P>
6 G0 K& B' U" M<>void banner() {
; v, |+ P7 S  Q3 t; J printf("\t\t****************************************\n");
printf("\t\t         单纯型法解线性规划问题\n");
printf("\t\t                         作者：Thunder\n");
5 T7 @/ e8 c, V! `& N printf("\t\t****************************************\n");
7 R, a+ z: `! B' l+ y5 D printf("\n");
5 `3 R1 a$ X) w$ b( i$ c; \}</P>
4 S; I& P0 ]: W  O<>void show() {
) Y" Q% Z+ {( s" J( v. n5 N//对方程组以自然的格式输出，系数为零的x不显示
" l% o) U2 O( M! `( a; @4 J//为1的不显示系数1，-1系数只显示负号
int i,j,k;
switch(step) {
case 1:
. C  d$ f# J  J/ t  printf("min z= ");
3 p% O' s8 X6 G! s2 v( |5 a) c; r8 \3 g  printf("x[%d]",arti[1]);
  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
4 f# e% o& @/ F8 J! R2 [6 p) h7 Q case 2:
9 V3 |3 }, H% N4 d, ~ printf("max z= ");
9 Z6 Z9 j% W" l" M! s printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
3 c& c5 Y% D8 l& Y& Z  if(c==1) printf(" + x[%d]",i);
& ]. ]$ U& |9 H! Z. f& Y  else if(c==-1) printf(" - x[%d]",i);
  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
( ?" l9 `/ {: ~: C" ?( {5 U  else printf(" %lg x[%d]",c,i);
}
" s; E% \  _: Z7 e  F- G; \$ a break;
* Q  a9 @* s/ k/ d/ D& ~- Q }</P>
* ~. S$ V1 p' i. H4 a- c<> printf("\nst:\n");
6 w4 i1 R6 ]/ D3 O" I for(i=1;i&lt;=num_st;i++) {
# q, l, u1 M1 A9 d3 W2 `  k=0;
& |. [4 g2 K+ L/ R; R  for(j=1;j&lt;=num_x;j++) {
* ^, l5 f, b; O5 Z   if(a[j]!=0) {
, w# @, J* r/ C    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);
    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);
3 X' O1 V( K+ P5 g0 s9 n2 E    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
, B& r& t$ @" }3 h5 x% F' @; f    else printf(" %lg x[%d]",a[j],j);
+ u0 j/ a: w/ ~: s    k=1;
   }
% \4 }) H8 Y2 z% c7 ?$ J/ l1 C( z3 H  }
1 [' _8 o9 @5 M& Q! x  printf(" == %lg\n",b);
}
# m* D$ W0 A, u+ W* d3 A5 n- } printf(" x[1]~x[%d]&gt;=0",num_x);
9 {( K6 Y, n# G7 w5 C" @! z8 O}</P>
; K/ H* q# i. g9 J) Z7 {! i<>void exchange() {
int i;
* x  [5 W4 }) \+ I7 j double temp[MAX];
* C/ I" Z. c2 A6 L& _1 { for(i=1;i&lt;=num_x;i++) {
" w6 B% A. F7 M+ x7 f1 g  temp=temp_c;
" L8 L2 J& L: b  temp_c=c;
  c=temp;
8 W" u0 A/ Q1 E3 ]) h }
}</P>
. b9 c5 L) V. [$ D1 G<>void create() {
! W% E$ P+ {5 I% y  g3 L//输入方程组系数，每个方程输完后回显确认
0 f. P* ]  ?# s7 ]3 x1 D0 P int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
3 ]8 {* Y- ]8 Y9 F char confirm;
8 @5 M. q6 Z2 n' R
while(1) {
* J; ~% |) S/ m8 W' S8 P% K2 R  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  scanf("%d",&amp;ma_mi);
  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
/ Q6 V. a% F# l0 h( h+ a  else break;
}
4 k2 n  r9 S( F8 v& b
  _/ P& Y- q% @2 B8 L* B$ O+ K while(1) {
  printf("指定变量个数：");
4 s: M$ \# X2 n$ y  scanf("%d",&amp;num_x);
  printf("输入价值系数c1-c%d:\n",num_x);
  for(i=1;i&lt;=num_x;i++) {   
   printf("c%d=",i);
   scanf("%lf",&amp;c);
$ O3 S% q# {+ A- g& p, G0 H  }
  if(ma_mi==1) printf("max z= ");
  else printf("min z= ");
  printf("%lg x[%d]",c[1],1);
' Q! C8 X2 O' U- F/ \  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);
! g6 o1 ^. }" ^+ o6 s0 V+ a  }
printf("\n正确吗？：(y/n)");
5 X! z; k/ W  e) F3 F' o getchar();
confirm=getchar();
if (confirm=='y') break;
else if(confirm=='n') continue;
& m2 `. Z/ ]# d) a, E) j) N }</P>
<> printf("输入约束方程组个数：");
scanf("%d",&amp;num_st);
for(i=1;i&lt;=num_st;i++) {
  printf("st.%d:\n",i);
3 G' o3 ]1 W% z# C5 {2 T' Z  while(1) {
" O3 w& x3 f8 B   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
& x8 h- A- F0 o  D" v   scanf("%d",&amp;re_st);
   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
   else break;
9 p2 }* \: n) ~! a; ^9 x  }
- W0 i3 R2 p/ q  printf("输入技术系数：\n");
; ?; S4 u! m+ G2 ^: d+ I6 |  for(j=1;j&lt;=num_x;j++) {  
  D0 d" R* S+ g4 h1 Z2 [   printf("a%d=",j);
  ?6 q3 k& b; R' X   scanf("%lf",&amp;a[j]);
  }
0 l' [5 i" W* m1 \- ~  printf("输入资源拥有量：\nb%d=",i);
  scanf("%lf",&amp;b);
  
4 ~" f4 L* _0 h5 |3 y2 s  printf("st.%i:\n",i);
  printf("%lg x[%d]",a[1],1);
  for(j=2;j&lt;=num_x;j++) {
1 N# ]. }9 ]) }% r   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
   else printf(" %lg x[%d]",a[j],j);
) [" C$ f& n( ~  }
0 h" E% T, z1 h; \: _  switch(re_st) {
& h& f! m$ [5 e$ l0 k   case 1: printf(" == %lg",b); break;
   case 2: printf(" &gt;= %lg",b); break;
   case 3: printf(" &lt;= %lg",b); break;
  }</P>
% `* ^- `4 `7 V" {& a<>  while(1) {
, y$ Z+ S5 p4 ?8 L   printf("\n正确吗？(y/n)");
   getchar();
   confirm=getchar();
   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
2 M9 y2 S% H, b/ D5 Y: z2 F  }
}</P>
<>//显示输入的方程组
printf("\n原问题为：\n\n");
. z6 O2 p1 N! x if(ma_mi==1) printf("max z= ");
else printf("min z= ");
& [: Y# M3 d3 X printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
* r% I8 b+ A- M; \7 m+ |# m& ~1 ~  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);
* w/ C6 a/ A0 [, C }</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
& J' w8 l8 g/ v% h  k=0;
  for(j=1;j&lt;=num_x;j++) {
7 y4 o  r7 E& f2 ?0 u" _   if(a[j]!=0) {
% f# F, r$ u' d1 K* x& Q4 V/ F    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
: l1 P, |- s( u, S3 E+ x7 q$ t# g# i    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
6 l! q3 i0 G7 v/ y5 d4 d    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);
# M8 |0 @9 |0 z) U2 T, ~    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
$ j$ v" d  e! @9 I1 x2 Y  U- R    else printf(" %lg x[%d]",a[j],j);
2 s" T7 O6 }2 M# [3 f8 h6 S    k=1;
   }
  }- k  G4 e9 n- x. m  }
- }4 h% _; O+ [  switch(re_st) {
( M% i+ `6 n" U' Q1 w   case 1:
% R7 i2 m2 p8 B9 ^6 m' q* S. s    printf(" == %lg\n",b);
    break;
6 ?4 }6 D3 j6 j* H. W  `! R! a   case 2:
    printf(" &gt;= %lg\n",b);
0 ?# [) y* I, |' L    break;
3 `5 Z7 @: f' O6 {( ?   case 3:
    printf(" &lt;= %lg\n",b);
* S0 Q- R2 G1 a, Z: ~" Z. K    break;
  }
5 B0 ~, `' O* b$ j, X( ^5 ~ }
printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
<> tnum_x=num_x;
for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
  case 1:
  case 3:
   num_x+=1;
7 y. a" q' \# V5 i" P7 b   break;
  case 2:
) X. Y: D6 V; d# A) Q8 R) m% @( S   num_x+=2;
1 h1 W8 y% {3 l9 o- _: O   break;
2 X4 p+ Y% K9 l9 `( z  }
# j% c( Q0 Y% q }</P>
<>//化为标准形式
if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
0 ?! h1 e' I8 j for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
   case 1:
8 v6 t( k- y( |& W    num_addv++;
) X% k- l3 C& U    num_ba++;
    num_ar++;
) u8 b& Z( C4 Q9 n' r    c[tnum_x+num_addv]=0;
    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
5 x# Y% R, j  l& |: a    for(j=tnum_x+1;j&lt;=num_x;j++)
: H( a/ _# a* p5 X     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
: l$ X+ l6 Q  O9 _     else a[j]=0;
6 K0 y+ c) b+ I5 B+ U* n2 X0 _/ ^    break;
. u) ^( n' w! ]% W2 j. _   case 2:
: k! }0 T' M9 v  e, R6 c: o% F  g( |    num_addv++;
    c[tnum_x+num_addv]=0;
$ t! n. K0 D0 k    num_addv++;
    num_ba++;
( |) R: b: L. u1 S3 C" r    num_ar++;
1 D% Q% E7 A  Z9 B, S    c[tnum_x+num_addv]=0;
( v' F. I( d+ a% |    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
) ]) n+ X4 N# V+ \8 {    for(j=tnum_x+1;j&lt;=num_x;j++)
# \9 Y. w- y) B: Q' D1 m     if(j==tnum_x+num_addv-1) a[tnum_x+num_addv-1]=-1;
     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
    break;
   case 3:
    num_addv++;
7 u( v) C( i/ T0 s" d% M    num_ba++;
    c[tnum_x+num_addv]=0;
6 g! b1 d& X( D: u    base[num_ba]=tnum_x+num_addv;
0 e* I2 s6 F' S  T' A2 g    for(j=tnum_x+1;j&lt;=num_x;j++)
$ L7 ]; U2 }; m     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
    break;
  }//switch
}//增加松弛变量、剩余变量、人工变量、确定基变量</P>
8 ]9 M( m; V/ E' O7 Z* K<>//显示标准化后的方程组
printf("\n化为标准形式后：\n\n");
/ ~8 @& y0 B9 ?4 r, o$ `' F if(ma_mi==1) printf("max z= ");
- B# [% R5 M" p1 d* j. U  V else printf("max z'= ");
printf("%lg x[%d]",c[1],1);
8 U; I% O1 \* H  e8 b9 n+ r for(i=2;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_ar;j++)
  if(i==arti[j]) k=1;
  D4 Z" U/ J+ R, Y# [  N4 @* a9 m* [0 `  if(k==1) printf(" -M x[%d]",i);
4 p5 b- C5 [6 d( R4 ]- p6 M0 O/ A  else if(c==1) printf(" + x[%d]",i);
6 l2 x$ r$ a0 C5 d1 c  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>
( R! `) m9 k& r8 i2 d$ @2 U7 F+ c1 ?<> printf("\nst:\n");
0 L: A& g) z, v/ d1 A) ]9 ]! ? for(i=1;i&lt;=num_st;i++) {
+ y* v3 g4 w/ b% A8 Y! h; t  k=0;
  for(j=1;j&lt;=num_x;j++) {
; j& A8 k5 K1 t: `5 N- u3 d: E2 @   if(a[j]!=0) {
; d6 `# M' ?+ a% A    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
( g, T* i  e6 w7 H# u, L7 D    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
; B" ?" |% z4 U% ^) ?9 g    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);
1 P' @1 W. l' m8 l% k2 P2 c    else printf(" %lg x[%d]",a[j],j);
    k=1;
! x1 ^5 P  ]: t' a: ]   }
$ A6 j2 A) J$ q! E' w' W) I, M  }
' J" y% O+ i0 m- N* t, C  printf(" == %lg\n",b);
  B1 @1 t0 p3 ^9 c) h: v- _ }
- d8 Y# Y' G0 S6 X( z/ v printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
<>void iterative() {
int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
int base_elem;
int base_out,base_in;
double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
<> printf("\n\n第%d次迭代：\n\n",stop);
* w; s0 T: c7 p9 E* Z, o$ \! T; p& v for(i=1;i&lt;=num_st;i++) {
' E2 k/ ]% o( j" ~  printf("c%d=%lg\t",base,c[base]);
  printf("b%d=%lg\t",i,b);</P>
2 f( T8 C# ^3 K1 R2 R$ ?3 J<>  switch(step) {
   case 1:
* b$ ^; [3 {( T    for(j=1;j&lt;=num_x;j++)
    {
     printf("a[%d][%d]=%lg\t",i,j,a[j]);
2 S( R7 J5 Q& y* N/ S, h0 u    }
    printf("\n");
    break;
( h  C: _6 E" O+ [4 ~6 x   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;
     if(k_a!=1) printf("a[%d][%d]=%lg\t",i,j,a[j]);
/ N  X) C( z' d1 S    }
; ?6 z! F: f( B- }  C    printf("\n");
    break;
  }
}
: C5 o& i0 T( ^8 ?6 I' Y//求检验数sigma
, W3 ]" H6 a  z4 y/ S for(i=1;i&lt;=num_x;i++) {
3 }5 L6 C( n+ B$ K& R# b  sigma=c;
" m1 e% B$ p9 {  for(j=1;j&lt;=num_st;j++) sigma-=c[base[j]]*a[j];
3 F9 A% N/ K( }4 A* F" o$ V  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
  C4 y: I+ g- h$ p* h  Q  switch(step) {
   case 1:
3 e& u! ?4 }; W6 y9 [2 A+ O    printf("sigma[%d]=%lg\t",i,sigma);
0 O. S' S# f. u' x+ v1 _    break;
. n: j+ z8 W0 J' k   case 2:
- T) k, P2 ?, N6 u    k_a=0;
& L- R& O; q2 ?! Z# N1 f    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
; B2 b! h' n& d$ K: T' o3 U( u) ^    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
    break;
3 U) T6 ~; b" G, F  }
}
# L, P; Y3 |9 @+ z; A5 t; r putchar('\n');
//检验检验数sigma是否全小于等于0
9 O/ b7 K" ?+ `- L# } k=0;
  ~. x  w8 U1 e for(i=1;i&lt;=num_x;i++) {
  if(sigma&gt;0)
  k=1;
}
if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
  K1 W7 V: l0 a0 E: L  for(i=1;i&lt;=num_x;i++) {
7 S8 [2 r3 z' G' d: O1 @, Y   k_f=k_a=0;
   for(j=1;j&lt;=num_ar;j++)
6 r) f& J2 O( X* ^   if(i==arti[j]) k_a=1;
3 [& Z5 g- b! o% Q3 v   if(sigma==0&amp;&amp;k_a!=1) {
    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
) E2 H7 `3 f2 r. R  A4 @6 W. w$ k   if(k_f==0) {status=-1; return;}
$ D6 t, d8 @1 z: z   }
  }
3 e: h' r2 C+ d  status=1;
  return;
}
9 J$ B+ T- l# q//检查是否为无界解
/ |& d% F7 |- H6 Y% a for(i=1;i&lt;=num_x;i++) {
5 u0 p4 [7 G- _  k_f=0;
; d' P: Z8 b9 C6 d  if(sigma&gt;0) {
   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
4 f2 r' l1 S$ A7 y, A+ t5 c   if(k_f!=1) {status=0; return;}
0 \! C, ~+ f) S2 I  }
}</P>
0 M, E8 v% G1 Y( ]1 Q<>//确定换入变量
for(i=1;i&lt;=num_x;i++) {
  k=0;
3 Y9 ?! X8 f) J/ G% y8 Y+ q  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
}//temp赋初值
7 m% ]1 j0 f6 A for(i=1;i&lt;=num_x;i++) {
6 \! @. p& e% M6 o4 M  k=0;
4 H1 |% R9 [' z% a! }3 T0 n! a  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
0 n# g* J3 s$ d0 r3 X1 W& w" S% D- [; W  if(k==0)
   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
    base_in=i;
3 u$ f6 P; h6 {; {$ b    temp=sigma;
   }
7 N6 c( P0 {( [% w+ c2 K }</P>
8 W7 z  J' n- \6 V8 V) |<>//确定换出变量
/ v: t# N6 W' i& e: i! T: E for(i=1;i&lt;=num_st;i++)
4 Z+ }5 J% t  j$ P% A: v% L  if(a[base_in]&gt;0) {
! z3 ~& ^8 P  \! a   temp=b/a[base_in]+1;
9 w8 I! b9 P1 a' }2 o/ h   break;
: {! P- U9 R0 o1 P  }//temp赋初值
for(i=1;i&lt;=num_st;i++) {
  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
   for(j=1;j&lt;=num_ar;j++)
    if(base==arti[j]) {
     base_out=base;
0 q' ~  v! K/ `6 j9 N4 o     base_elem=i;
     temp=b/a[base_in];
; o1 b  ~8 R+ [     break;
6 \2 x% V. Q) [7 |    }
  a; @* Q9 P& T* U# C  }//人工变量优先换出
) }9 }6 h) D0 f- b( N% Q  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
   base_out=base;
   base_elem=i;
   temp=b/a[base_in];
  }
- m6 u5 T& A2 ~7 M$ _+ m }</P>
; Z  c$ l5 E; w1 _3 X7 L2 x  G<> printf(" 基变量：");
for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
//基变量变换，进行新方程初始化后迭代
for(i=1;i&lt;=num_st;i++) {
  if(base==base_out) base=base_in;
6 J5 u+ }" P6 C) w6 n6 Y }
- `1 h& U. r6 {6 g* L9 I//初始化主元素行系数
0 }9 K/ p; l3 F value_be=a[base_elem][base_in];
; L/ N- n* }8 z6 W0 `( P9 N4 {( T+ _ b[base_elem]/=value_be;
* t3 |4 W6 E% F7 g for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
; G/ K) g' c, G' I, c9 g<> for(i=1;i&lt;=num_st;i++) {
' U# u7 u6 z5 n% y' Z; |  if(i!=base_elem) {
- E  \; @" F* y% m2 ?) b$ |   b-=b[base_elem]*a[base_in];
& z: o8 r+ T& d, s; }   value_be=a[base_in];
5 ^0 M! }' o. N) j; B8 v) l   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  }
}
) Q0 P- O6 D' u* {  v6 M, J stop++;
8 d5 d1 F' a. V, f# H  a if(stop&gt;STP) {status=-2; return;}
iterative();
}</P>
<>void output() {
+ b9 n* f7 j7 N$ D+ D6 b" O int i,j;
# E8 i8 i: [: M2 R- P double X[MAX];
printf("\n结果如下：\n");
printf("\nX=(");
for(i=1;i&lt;=num_x;i++) {
  for(j=1;j&lt;=num_st;j++)
9 Z! d) V8 ^- F0 u7 o) X: T  if(i==base[j]) {X=b[j];break;}
  else X=0;
/ k: n1 _& P6 o" g+ N  printf("%lg ",X);
}
printf(")");
4 l+ x" J: D, x# E- w+ V for(i=1;i&lt;=num_x;i++) max+=c*X;
8 g( d% A/ S1 T' g. [+ a2 T if(ma_mi==1) printf("\nMax z= %lf\n",max);
else printf("\nMin z= %lf\n",-max);
: E4 q% T" m, F; {) X* Q' ~}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了