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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
) k0 f7 q% \! ~) H$ Y#include&lt;math.h&gt;
# b7 g- m' P& ^3 `+ _& b#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])
; B5 K; ^2 T* c) d- m{
5 U# R! \! r8 `# w- r0 _  int j;
  *c=0;
  for(j=1;j&lt;=*is;j++)
( h2 R9 W8 O1 j: d  S    {
      if((*a)[*ir][j]-*c&gt;0)
3 p& {& S: h) J/ i; a5 {( C2 T    {
      *c=(*a)[*ir][j];
# u3 v2 U$ @$ I) G& v% E+ r, x, [) M      *j0=j;
    }
' k0 C  D6 T0 R     }
}</P><>      /***************** 参数说明 **********************/
      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
6 N; y3 _2 u4 q# E2 j+ @      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
      /********* <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],
       int(*k)[],float(*x)[])
{
  int i,m1,mn,j0,i0,j;
5 f, v9 U7 h: i/ H2 f. ~9 j  u  float c,g;
  m1=m+1;
- v% T1 d, e- |; ~) ~& y  mn=m+n;
1 q2 O" a% ]5 r2 O7 V8 O7 W  @  for(i=1;i&lt;=m;i++)
    (*k)=n+i;
3 c! `* A& M5 v7 [( ]# \leap1:
" ]6 Y- c6 f3 P9 B/ M8 w2 [4 I5 K% e  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
( l2 s0 G% L7 G8 @9 Q4 T9 u  else
     {
5 v: s7 T# g, x+ v; {2 R) Rleap2:   if(l1-50==0)
( R# P8 P# u% X      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
! w) [4 \# W  Y# n       else
* h& v$ k; Z+ @; U' F+ X5 @( Z      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
     }
% d9 y3 M7 H$ u0 W, I  _- k, H   c-=1e-8;
5 [! m) G+ h. z6 [' V' f# L! p* m4 v   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);
      }
5 a( j& x& N9 Q* \   if((c&lt;0)&amp;&amp;(l1==1))
       {
     l1=50;
; a) E+ Q7 S& P% E1 W     goto leap2;
       }
    if(c&lt;=0)
       {
      for(i=1;i&lt;=m;i++)
' v; M0 h) M2 U! M4 I         (*x)=(*a)[mn1];
4 w8 P# k! G, b" ], r6 ]. ]         printf("\n\t********Optimal solution********\n");
      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
      printf("\nF=%f\n",(*a)[m2][mn1]);
6 G, T6 F: h1 g0 H      return 0;
    }
     c=1e8;
     for(i=1;i&lt;=m;i++)
       {
     if((*a)[j0]&gt;1e-8)
        {
          g=(*a)[mn1]/(*a)[j0];
# u+ Z* X" O: X. l          if(g-c&lt;0);
         {
2 E# p! V$ ~3 Z& ]           c=g;
           i0=i;
3 N# s6 f3 D' w4 {% A& n& C          }
$ K' Z8 I# y, s4 _4 e         }
% u; F) L" g! P) Y* d    }
      if(c==1e+8)
* c5 ?- g6 d9 S+ j- e    {
      printf("\n\t*******Lp no solution********\n");
      return -3;
: `# c  _2 _& |/ D     }
1 \. w! H' Q9 z: f5 b       (*k)[i0]=j0;
/ q: g" J% @) I4 ?       for(j=1;j&lt;=mn1;j++)
      {
8 t2 ]8 _# c9 [6 {1 l+ k/ i6 f        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
. Y3 y4 ^& n, \           continue;
3 g8 h9 m5 ]; m6 j& Z% D8 ?         g=(*a)[i0][j]/(*a)[i0][j0];
* d2 j4 K, X" l3 c0 e8 Q8 h           (*a)[i0][j]=g;
% T% G' a* a3 Y3 U1 w; G3 K         for(i=1;i&lt;=m2;i++)
           {
" `! {2 T/ ^. n* i4 g2 Z5 q! i, H         if((i==i0)||(!(l1==1)&amp;(i==m1)))
# f5 O4 s$ n/ d, r5 i( z            continue;
$ v6 N  e+ e+ L8 H& L+ p: I         (*a)[j]=(*a)[j]-(*a)[j0]*g;
0 e* {' n; m3 R8 f0 e, e* d           }
       }
' w4 Q8 E2 v+ c       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
       (*a)[i0][j0]=1;
       goto leap1;
}  
main()
0 v/ O& I# v% r3 {& r' d# g* Z# X{
5 Z. w! y6 b& W* n  float a[5][7]={{0,0,0, 0,  0,0,0},
                 {0,3,-4,3,  1,0,12},
                 {0,3,0, 6,  0,1,12},
  [) H2 d" V# i" M                 {0,0,0, 0,  0,0,0},
2 n* P3 g+ {, x6 D- L; z7 m( s8 P                 {0,69,0,144,0,0,300}};
: f- D: _3 I+ ?+ X! S; y  int k[3];
# R4 y1 O% N* W: {) C" Y; m9 ^" @  float x[3];
5 [$ |" I8 W8 ^9 u) F  Y) P. D  clrscr();
  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
8 b* A6 p$ U/ O}
</P>

plgatc

<>/*************************************************************************
$ @, J. u& z5 c6 z2 f                    单纯型法解线性规划问题（两阶段法）        

编程环境：VC++6.0     
; r6 S+ P! ^7 o! n- {' a8 o* v 方程组输入说明：
变量非负，按提示输入相关参数。
  a5 s! @% V9 @0 z *************************************************************************/
1 ]2 q9 g( s  e3 a1 D0 g6 u* N#include &lt;stdio.h&gt;
7 J9 K7 m. L) m, Q#include &lt;stdlib.h&gt;
#define MAX 100
#define STP 100</P>
9 R; X3 i0 e, ^5 g6 p9 z8 |<>int stop=1; //迭代记数变量
! Z, b) S2 h8 d7 s$ nint status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
+ |) u: V( T6 X& X% W" f* Fint step=1; //目前阶段</P>
<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
int num_x; //变量个数
  t) v; P3 I$ J6 K; ]; uint num_st; //约束方程数
) L6 S% z4 I9 E% f" j6 xint num_ar=0; //人工变量个数
; e  k. `" d4 r' m$ y8 eint arti[MAX]; //人工变量下标
. Z5 h; }2 N6 x! vint base[MAX]; //基变量下标
" M! _9 ~4 ^: U5 o8 Z* \int ma_mi; //1为求最大值，2为求最小值</P>
<>void create(); //建立方程组
void iterative(); //单纯型法迭代
void output(); //输出结果
void banner(); //打印程序标题
, _  w; k! z8 Z) z+ ~  mvoid exchange(); //交换两阶段价值系数
- ?0 N9 D( T% t. S( w0 ^5 Evoid show(); //输出方程组</P>
<>void main() {
int i,j,k;
banner();
create();
//保存原价值系数，转换为第一阶段价值系数
( E8 k: E, |3 V+ \7 A/ |3 `$ ^ for(i=1;i&lt;=num_x;i++) {
7 @8 q9 T, f; U- l  k=0;
7 L, r- T! d! Q( Q  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
( E% |9 B7 |  m  if(k==1) temp_c=-1;
  else temp_c=0;
}
exchange(c,temp_c);</P>
: n- e! g& |1 h% n  K6 |  x<> printf("\n\n第一阶段问题为：\n\n");
1 ^9 S( V+ U% I" \" j show();
step++;
( |/ D# m/ Q5 | printf("\n\n按回车开始第一阶段迭代");
getchar();
getchar();
' Y9 L: C" _* J% g iterative();
if(status==-2) {
/ H" E- e; Z. v  S  puts("迭代超过限制次数强行终止！\n");
- K! u# {7 Y  M2 {& W* M. ^9 V- _! A  puts("\n按回车结束");
3 Y6 p( V! t# D# e  getchar();
  exit(0);
}
2 a' L. j8 \& v& k0 \/ W# G: X output();</P>
<> if(max!=0) {
; l. H- L4 F$ e0 }* z  puts("\n\n原问题无可行解。\n");
  puts("\n按回车结束");
! D% j. x0 a6 ~3 R* b  getchar();
  exit(0);
}</P>
<> //转换为第二阶段价值系数
/ J  k3 n  v% G( o exchange(c,temp_c);
5 U6 m. D: s8 Y) S //把人工变量列全设为0
* U& A5 y: T- ?: h% R! ^ for(i=1;i&lt;=num_ar;i++) {
! O& W' W. g. R, m+ R, X0 G  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
1 Q0 S1 H7 G/ J' K/ w! v2 s: k }</P>
- \- I; o& F* A1 L" K<> puts("\n\n第二阶段问题为：\n\n");
show();
8 f* ^) J3 ], {/ d' P* l9 C% c puts("\n\n按回车开始第二阶段迭代");
; N- r6 U' L( d5 M# Y% E getchar();
iterative();
! s9 m1 r" r2 I  h( T switch(status) {
$ S9 q( ^) s% g1 S. g case 1:
  output();
  puts("\n\n原问题有唯一最优解。\n");
* W8 O- y+ [4 _2 ], @! J. m  puts("\n按回车结束");
  g* C/ [2 B7 i6 h  getchar();
. `( ?! l# n( u0 ]/ O  exit(0);
% I3 i% F- P" w4 A' J+ L& o case 0:
* c( N, E  [; a' _3 p  puts("\n\n原问题为无界解。\n");
* [, ]' k2 K- J* x* F& z; J  puts("\n按回车结束");
  getchar();
  exit(0);
6 _  K% s* U* A+ [4 U8 o case -1:
  output();
  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
: ?, N+ P+ ?7 U+ e7 l! X# u  getchar();
& v# ^1 N! r* Q  exit(0);
case -2:
( H" G2 Q9 z  v  M5 W+ x! O  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
  getchar();
$ `7 m# T+ {, {  exit(0);
5 ?8 [' u5 ]0 o- | }//switch
}</P>
<>void banner() {
' a# @% d6 L) {* E printf("\t\t****************************************\n");
printf("\t\t         单纯型法解线性规划问题\n");
8 k7 o. x" D1 q$ k7 i printf("\t\t                         作者：Thunder\n");
printf("\t\t****************************************\n");
printf("\n");
}</P>
: `% k: _3 {: ?0 R: G6 G<>void show() {
* ^) c6 [4 Z' G. r- B//对方程组以自然的格式输出，系数为零的x不显示
: O) T- k' H& v3 U- r2 }//为1的不显示系数1，-1系数只显示负号
4 `, I- u: @1 w3 n int i,j,k;
switch(step) {
case 1:
  printf("min z= ");
$ U1 F- I* z' [5 w" ~: Z$ b  printf("x[%d]",arti[1]);
& f6 x6 x7 n" T8 Q* [. g! p  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
/ _) a+ c5 {$ j& N" Y0 B case 2:
printf("max z= ");
printf("%lg x[%d]",c[1],1);
+ H  z6 [" e! q for(i=2;i&lt;=num_x;i++) {
  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);
- N7 O; d. h* D# Z4 }) U9 y$ R7 E  else printf(" %lg x[%d]",c,i);
, t- t: w5 k: \/ y }
5 C, z) D7 j$ c# E8 a break;
, x& r& M2 ]2 ~" a) _ }</P>
9 T. p) U! D/ n; ~; t! n<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
4 o& Q! L6 F! A" G3 u* m( Z4 z. M  k=0;
* Y  S. s3 [  }  for(j=1;j&lt;=num_x;j++) {
1 o- K6 W% H0 x6 i4 n   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
, J$ B8 j1 t8 {, W5 s7 s& A    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
    else if(a[j]==-1) printf(" - x[%d]",j);
! K$ j8 D, e7 B$ P& W( t    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);
    else printf(" %lg x[%d]",a[j],j);
5 p- q. ?  H8 \' C+ m    k=1;
$ C: N) P. |3 G6 F  ^6 n! N6 l   }
; _8 H: L  ]3 I, v" y  }
  printf(" == %lg\n",b);
}
) i# s0 d( e) u9 g( _* Y printf(" x[1]~x[%d]&gt;=0",num_x);
' i/ a* `3 T& k0 g! `6 O}</P>
<>void exchange() {
1 M! J; O. Y0 f- h  b int i;
4 I- }1 I: m: `. F9 {2 w8 W double temp[MAX];
% b. l! X, q; }" o  p for(i=1;i&lt;=num_x;i++) {
  temp=temp_c;
  temp_c=c;
& G# f0 @+ |% a6 s  c=temp;
0 f5 D  k( P3 u* y6 T9 r% H }
5 p6 A+ L# ^0 m  Y8 ^}</P>
! Y' ], W/ G1 _$ Y4 B5 E<>void create() {
//输入方程组系数，每个方程输完后回显确认
int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
char confirm;

4 a& Z" }  d6 O* r( Q) S while(1) {
# S/ @, T. ^$ @6 Z) K* k5 \0 c. \  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  scanf("%d",&amp;ma_mi);
  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
# r4 X/ m' M6 Q  else break;
}
, n/ ^) j9 k$ j6 M9 l) x
while(1) {
, i) d$ }6 W2 e) r: `4 D  printf("指定变量个数：");
. W  G0 a4 q  Y4 Q4 M/ m- g( e  scanf("%d",&amp;num_x);
3 ~! n% m' Q5 o# R! N( T# N& y  printf("输入价值系数c1-c%d:\n",num_x);
  for(i=1;i&lt;=num_x;i++) {   
! q% D2 ]% _2 D% ~3 s( U   printf("c%d=",i);
   scanf("%lf",&amp;c);
  }
  if(ma_mi==1) printf("max z= ");
  else printf("min z= ");
  printf("%lg x[%d]",c[1],1);
  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);
  }
printf("\n正确吗？：(y/n)");
( D9 D/ x; d% _4 Q5 f getchar();
confirm=getchar();
  X& N9 Y" K2 }7 @ if (confirm=='y') break;
else if(confirm=='n') continue;
}</P>
<> printf("输入约束方程组个数：");
scanf("%d",&amp;num_st);
# \4 B1 V- b: a/ I( e% v for(i=1;i&lt;=num_st;i++) {
  printf("st.%d:\n",i);
  while(1) {
2 }$ ^, `* @/ n" w' F; S   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
6 A# q1 x/ |% p: O- F0 }   scanf("%d",&amp;re_st);
' {% ]) B. i6 {7 T   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
   else break;
  }
  printf("输入技术系数：\n");
6 j$ c# d0 W0 _  for(j=1;j&lt;=num_x;j++) {  
* p* S: O9 t6 F+ B/ L, I   printf("a%d=",j);
0 Z2 u+ w; l, d7 d1 s   scanf("%lf",&amp;a[j]);
& p/ e, ^! C. ^/ o$ i) H) W# L  }
2 f. c- G+ _3 @* f: s1 d  printf("输入资源拥有量：\nb%d=",i);
; k. D; q7 Z/ j) R, V! q  scanf("%lf",&amp;b);
. _  `' h; \; a- \; v) c3 F  
# v$ h- J& _; l8 R( J  printf("st.%i:\n",i);
  printf("%lg x[%d]",a[1],1);
  d- o' r6 H0 q  e( g  for(j=2;j&lt;=num_x;j++) {
  j0 g( o' k- s' u( }+ R& Q5 N   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
' m: w3 L( N( a+ Q% W9 n' x   else printf(" %lg x[%d]",a[j],j);
  }
. |% ?3 Y( [, C" e, Y; u* z  switch(re_st) {
   case 1: printf(" == %lg",b); break;
   case 2: printf(" &gt;= %lg",b); break;
   case 3: printf(" &lt;= %lg",b); break;
( N& W  }1 `: M  }</P>
+ M7 |' K  \- h2 B% P<>  while(1) {
   printf("\n正确吗？(y/n)");
& U, o. N4 S4 v2 _; S( V, {5 z$ }   getchar();
   confirm=getchar();
! F$ r4 X: J' p" w' X, I) w   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
  }
: k: N* B: A3 i; S, ^2 E }</P>
) z& a4 c- b. T0 x+ K7 U<>//显示输入的方程组
printf("\n原问题为：\n\n");
/ V9 x( ]% P9 ]: \+ }; g if(ma_mi==1) printf("max z= ");
7 V9 d3 a" k( U: ]$ e else printf("min z= ");
printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
! U# e/ F3 B! [. y) ?  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);
- t8 z8 D! Q0 x& |! E }</P>
/ ?4 n8 P, @6 {6 [" p+ D* h- Z7 x<> printf("\nst:\n");
7 A! D% D: F, E2 H) S5 r for(i=1;i&lt;=num_st;i++) {
' W( U1 ~$ B# b( @) A; L% ^& @  k=0;
  for(j=1;j&lt;=num_x;j++) {
& @, o7 U1 G# l4 V4 A   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
' Q+ `' h, \: y% d! }    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
    else if(a[j]==-1) printf(" - x[%d]",j);
+ {8 r+ e: z0 A+ _    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
8 q' L- ]5 ]* P; @( Y3 @9 s    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
  v& ?0 h" ]8 O2 q: I6 b+ s1 |    else printf(" %lg x[%d]",a[j],j);
    k=1;
   }
  }
! k$ c1 V0 g) [  switch(re_st) {
7 L/ F3 S) _2 h+ j' b: d   case 1:
3 i( i  F" K! ^+ _6 c    printf(" == %lg\n",b);
5 u3 ?) _; g2 w* ^  u8 O* e3 \    break;
+ t: B# e+ z' M* K! A0 ^   case 2:
    printf(" &gt;= %lg\n",b);
+ B! M8 h/ a+ D+ a2 y% d    break;
4 b/ @' y( ?8 E+ s. e; d   case 3:
! ^* a' z1 \% s9 R; p0 v* `/ w    printf(" &lt;= %lg\n",b);
    break;
  }
8 b/ m( j5 d' E- u& u }
printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
# d/ z% I* ~4 e$ Q<> tnum_x=num_x;
for(i=1;i&lt;=num_st;i++) {
- i# D! G5 t* C9 V  switch(re_st) {
  case 1:
  case 3:
' n1 r$ z, _( m8 Z   num_x+=1;
   break;
  case 2:
   num_x+=2;
   break;
  }
' g. I! y3 |# T! n2 t+ R1 a0 B }</P>
6 J& G* {: F# V" g! X; }* Q<>//化为标准形式
+ s% R2 J8 L# f; O if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
, ?1 A7 H: W# W/ K) E  S- N* P for(i=1;i&lt;=num_st;i++) {
- Q0 i3 @: w6 u- B  switch(re_st) {
   case 1:
    num_addv++;
3 J. e  [5 m/ @& @% i    num_ba++;
; n" U  d8 T2 Y5 C. e6 A5 [    num_ar++;
    c[tnum_x+num_addv]=0;
( s) c9 m- @: w1 p0 K    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
    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;
    break;
- M% I& z, X; T' L   case 2:
    num_addv++;
    c[tnum_x+num_addv]=0;
  s0 Z+ }( Q6 l9 x4 |; L1 h2 U/ a    num_addv++;
7 g3 W# i1 u6 I# P1 [- u    num_ba++;
; J/ G+ [% e- f9 U5 N4 E5 ^    num_ar++;
    c[tnum_x+num_addv]=0;
' j( C% e1 A& D7 F    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
    for(j=tnum_x+1;j&lt;=num_x;j++)
5 C/ [" m* Z( p! D, S     if(j==tnum_x+num_addv-1) a[tnum_x+num_addv-1]=-1;
) T- V( S9 @8 x! B5 m3 ?% V     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
' Q+ z2 K1 W: d9 q% I& H* q# m: W     else a[j]=0;
    break;
) ^0 R- R; d  }7 ?; _! G   case 3:
( a! ?7 y6 i2 ^, X    num_addv++;
  x& E! B( _: |- B6 G* T  v: x    num_ba++;
    c[tnum_x+num_addv]=0;
    base[num_ba]=tnum_x+num_addv;
! S& ]! [2 W0 ]3 Z' Z8 O2 A; r    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;
    break;
  }//switch
}//增加松弛变量、剩余变量、人工变量、确定基变量</P>
<>//显示标准化后的方程组
0 r9 R: ~, {3 n4 b printf("\n化为标准形式后：\n\n");
6 s$ x% z" F3 u) C; P3 S8 S+ I9 \& b if(ma_mi==1) printf("max z= ");
else printf("max z'= ");
% g3 k6 X. o4 Z7 t$ k3 n8 ] printf("%lg x[%d]",c[1],1);
& t' i, S, W6 i: ~ for(i=2;i&lt;=num_x;i++) {
  k=0;
! N4 {8 v1 J  O2 ?" s3 q' q7 ^2 T  for(j=1;j&lt;=num_ar;j++)
, \0 G& X1 \, u  if(i==arti[j]) k=1;
  if(k==1) printf(" -M x[%d]",i);
  else if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
2 b1 D- B7 H- L/ h& Z  o3 I  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
6 `* H1 M1 l: T2 Q6 g) l5 T  else printf(" %lg x[%d]",c,i);
}</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  k=0;
# F  _4 n% O! Y" }7 y2 x  for(j=1;j&lt;=num_x;j++) {
1 X: H2 J, R- h% n   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);
, G7 l* s7 e" V4 m' 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);
4 ^& [, ?# p% H    else printf(" %lg x[%d]",a[j],j);
! O' e/ x4 G8 `    k=1;
   }
  }
( [9 e) g9 t4 x( f  printf(" == %lg\n",b);
}
" Z) r/ \7 o% A0 Z* e1 { printf(" x[1]~x[%d]&gt;=0",num_x);
9 |6 o* B$ |; G$ B+ ]}</P>
<>void iterative() {
int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
! u/ ^- I) [& P, h int base_elem;
int base_out,base_in;
3 d, |3 }: N4 x% H# O# Z double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
& q* f9 `# N2 p. C( @1 s( Z<> printf("\n\n第%d次迭代：\n\n",stop);
for(i=1;i&lt;=num_st;i++) {
: \3 B% |! ]8 V' y5 Y& j  printf("c%d=%lg\t",base,c[base]);
3 j. H0 k& u9 }0 }- i$ u  printf("b%d=%lg\t",i,b);</P>
' \; u8 \) E+ N<>  switch(step) {
; Y1 P+ Z& |3 Y   case 1:
    for(j=1;j&lt;=num_x;j++)
    {
2 O8 w, \& a. D# m- w0 G2 [1 X* a' I     printf("a[%d][%d]=%lg\t",i,j,a[j]);
    }
    printf("\n");
6 v) E! y" x6 _    break;
" r" g- N! [' v& ^+ U4 l   case 2:
' G, h/ D0 {0 o% _1 @    for(j=1;j&lt;=num_x;j++) {
4 w( c% u3 k- k3 E& p; @     k_a=0;
! \6 q; n; \) c9 Z5 r     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]);
) V4 t4 v3 |- V  a    }
) S% W! W/ B/ Y) J$ }3 Y* S    printf("\n");
; ^- H& {) D7 N' i- c( K1 v    break;
  }
}
* s( x$ S6 y; O- D//求检验数sigma
* w7 ?' y' i7 z3 B& Z( P, P+ L for(i=1;i&lt;=num_x;i++) {
( s$ I* f/ R1 |: i  sigma=c;
  for(j=1;j&lt;=num_st;j++) sigma-=c[base[j]]*a[j];
9 s5 M1 a; f8 ?8 N9 }  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
  switch(step) {
! e5 {# X+ ~0 o. G5 z, k   case 1:
) W: E+ A8 _% c. h: K    printf("sigma[%d]=%lg\t",i,sigma);
! U' Y0 ^* w7 Q' u9 l' @7 x7 y% M; g    break;
   case 2:
6 E7 X8 e9 o# _1 ?    k_a=0;
    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
3 S; I' j; @& L; y# F/ b0 m% g) W    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
    break;
6 v5 ?0 r9 T' q" s+ m. `; U  w+ O  }
}
putchar('\n');
//检验检验数sigma是否全小于等于0
, @' `8 e( N+ J& ~" q, w* ^ k=0;
for(i=1;i&lt;=num_x;i++) {
( B0 H' r5 |- q& D1 _3 b  if(sigma&gt;0)
  k=1;
/ ?% @5 s# t/ k7 {7 l7 l& Y }
; |$ o8 s6 h- f1 v& @ if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
6 T5 V* O" E, x8 t4 Y  for(i=1;i&lt;=num_x;i++) {
   k_f=k_a=0;
   for(j=1;j&lt;=num_ar;j++)
4 @( F: F9 U2 {6 j  X; p   if(i==arti[j]) k_a=1;
   if(sigma==0&amp;&amp;k_a!=1) {
1 R; ]& q8 ?1 G$ a7 P. o: I    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
   if(k_f==0) {status=-1; return;}
   }
  }
  status=1;
  return;
}
* R, @9 G  v. m7 Z1 _; u1 O5 P* U//检查是否为无界解
! {, F3 H/ x: L$ b for(i=1;i&lt;=num_x;i++) {
$ Q( z5 [( E3 r! v  k_f=0;
) |- k% n6 P" k: b4 C7 P  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;}
  }
}</P>
<>//确定换入变量
5 y9 a# J9 U1 ~2 U0 z for(i=1;i&lt;=num_x;i++) {
  k=0;
6 o( y  D3 P9 P% V' h  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赋初值
for(i=1;i&lt;=num_x;i++) {
  k=0;
# V) _6 I6 q4 J  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
9 I' p0 s; k2 N1 u3 n( Y; W  if(k==0)
   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
6 M& h. {% c( y    base_in=i;
    temp=sigma;
   }
}</P>
<>//确定换出变量
; _6 i* Z$ u' P  n1 ] for(i=1;i&lt;=num_st;i++)
; N- Z! a" q" i) e8 F  if(a[base_in]&gt;0) {
* p2 h# M# C9 k6 L+ `   temp=b/a[base_in]+1;
. L% p6 l+ U; m- o. Y. x/ i   break;
" q* F7 i# l  @( ~  }//temp赋初值
. `5 W+ }* f9 g/ }5 h+ A* ~) ] for(i=1;i&lt;=num_st;i++) {
0 j$ O  b3 W4 a# j  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
5 L# a2 K2 m5 u& V# U   for(j=1;j&lt;=num_ar;j++)
: R* F: w$ z8 v3 S    if(base==arti[j]) {
     base_out=base;
, K; c, w, Q# O5 r" j     base_elem=i;
     temp=b/a[base_in];
( i- p4 u7 \% p! K; D- {0 \3 n     break;
    }
  }//人工变量优先换出
4 l9 p* b: k7 ~* L. v6 r/ S- _  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
: v( C: r$ W! ~   base_out=base;
   base_elem=i;
   temp=b/a[base_in];
  {+ F0 [; A8 U0 n  }
}</P>
5 t# c( y4 r% a5 ]<> printf(" 基变量：");
for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
3 _/ ]8 R+ Y& k' }, C7 `2 `% B printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
* y6 g% c/ T/ {2 `4 @//基变量变换，进行新方程初始化后迭代
3 b+ m) N. i7 r% k7 f' E for(i=1;i&lt;=num_st;i++) {
  if(base==base_out) base=base_in;
4 u. ^3 G; s" L, p6 W( S. E4 | }
% i7 l5 L: A* G" J: f9 f//初始化主元素行系数
2 \1 W. m! p  ]9 `; t value_be=a[base_elem][base_in];
1 y" Q4 |# p% c4 ^, T* ~) Y b[base_elem]/=value_be;
for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
0 V4 \1 v6 p. d" Z" x  D<> for(i=1;i&lt;=num_st;i++) {
  if(i!=base_elem) {
   b-=b[base_elem]*a[base_in];
   value_be=a[base_in];
  x/ M- H; j8 X0 Q* {% J# v   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  }
}
8 }! e: P( w+ F- b2 F* B6 u stop++;
if(stop&gt;STP) {status=-2; return;}
: }" a1 G8 r6 ^" f8 X4 e1 I iterative();
}</P>
* U. [# v/ q2 ^/ ]( V1 C% v, E<>void output() {
int i,j;
" G+ E( n* S. W- w% ]# Q3 |% S double X[MAX];
; ~9 r& C/ }/ S& ? printf("\n结果如下：\n");
printf("\nX=(");
for(i=1;i&lt;=num_x;i++) {
  for(j=1;j&lt;=num_st;j++)
  if(i==base[j]) {X=b[j];break;}
! [0 H& a! f% z8 {  else X=0;
  printf("%lg ",X);
}
printf(")");
0 R$ g7 T6 h' X8 x; R; H" q for(i=1;i&lt;=num_x;i++) max+=c*X;
, D/ z0 r/ i$ P6 v8 c- P if(ma_mi==1) printf("\nMax z= %lf\n",max);
" Z/ O$ b: L5 R else printf("\nMin z= %lf\n",-max);
}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了