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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
7 S2 W- o% {+ `#include&lt;math.h&gt;
# y, h# A" w8 [* w$ I2 X3 n#define X 5
#define Y 7
- U" h) b* |& x, h! M8 qvoid xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
3 T/ v; N5 J' g2 l{
: x! u6 K, i' Y5 Z7 K( R, y* g  int j;
0 I- a; Y; n1 m" `  *c=0;
& ~) _3 @# D6 p/ V  for(j=1;j&lt;=*is;j++)
  f* K# }4 {" k8 }# D8 p7 p% V    {
      if((*a)[*ir][j]-*c&gt;0)
8 F. V/ }( t  k1 W; ^# b# r6 D    {
      *c=(*a)[*ir][j];
: s% `4 E  w, h! j. I8 Z      *j0=j;
    }
  l2 h. j3 s. g3 i     }
}</P><>      /***************** 参数说明 **********************/
      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
      /*   (*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],
1 C, V5 N0 y* i/ d7 F+ B6 u1 F. O" k       int(*k)[],float(*x)[])
2 E& Q  d# D5 a6 d1 J' A{
' B3 T4 L9 p, C) J8 v  s5 v3 `. ~7 h  int i,m1,mn,j0,i0,j;
! D: A, N( ~# Y0 H) y  float c,g;
  m1=m+1;
  mn=m+n;
  for(i=1;i&lt;=m;i++)
& c- a+ h6 ], w    (*k)=n+i;
leap1:
, F3 X: Y& f# s) x( b7 @  if(l1-1==0)
! C( b  W  ]( J, D, r4 {1 _     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
# ?5 m* o% f! N  else
" v/ O1 J4 n, N- B% b3 D2 O4 }4 e     {
( {/ P# [1 J9 g8 M1 W  j" jleap2:   if(l1-50==0)
      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
       else
2 i! x: u" X' v/ j! a  P4 e2 a5 K" g( J      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
" u* L# l3 k& @     }
; S* I% M" F* C5 j/ A   c-=1e-8;
( n( f0 O! _7 w& \  L, {   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
1 X2 Z# @% a$ C; H      {
2 V* Q' Z* U9 m9 f# }  |    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
    return(-1);
, d; N. B3 Z& O- Q) d      }
   if((c&lt;0)&amp;&amp;(l1==1))
  h  H" r% y3 ?7 e  r6 {5 \# `       {
     l1=50;
     goto leap2;
% t) ~8 C6 Y' K' t2 e/ g+ Q! P       }
" G, F7 d3 ~7 R7 J" b2 ~4 x! v& u    if(c&lt;=0)
( e, D% H3 t# J# b) r$ d       {
( ]* L9 o2 |8 @7 c- b& I      for(i=1;i&lt;=m;i++)
         (*x)=(*a)[mn1];
9 q" Z* C# X  M% ^         printf("\n\t********Optimal solution********\n");
& O8 t$ ?6 e3 o" k      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
      printf("\nF=%f\n",(*a)[m2][mn1]);
# b, `+ |/ O; H" V. v      return 0;
" h3 @) C& u* ^# u  j9 a    }
     c=1e8;
     for(i=1;i&lt;=m;i++)
       {
     if((*a)[j0]&gt;1e-8)
        {
          g=(*a)[mn1]/(*a)[j0];
          if(g-c&lt;0);
+ k" q3 z, w8 P- M3 h         {
           c=g;
           i0=i;
          }
         }
    }
% P0 t2 e1 V# J, u9 u5 R      if(c==1e+8)
/ |, Q. R6 b7 i4 o9 `. N# @    {
# t& Q: Z6 }2 V" f- j" \      printf("\n\t*******Lp no solution********\n");
4 x$ z( s- V" ^      return -3;
     }
       (*k)[i0]=j0;
1 E* h/ T% x3 q9 C       for(j=1;j&lt;=mn1;j++)
! Y# z3 Y& U$ k/ A9 r      {
: l0 f% H# _. P, m$ G' \: e        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
( q0 m8 n8 x* `, V" u           continue;
  g" Y' f7 y/ @         g=(*a)[i0][j]/(*a)[i0][j0];
" L# s& K; F- ~5 u; Y, _$ e2 |* D8 c           (*a)[i0][j]=g;
& j) @) U/ z; c9 t. t         for(i=1;i&lt;=m2;i++)
           {
         if((i==i0)||(!(l1==1)&amp;(i==m1)))
$ v7 T+ i: l& e7 y( _+ g6 K            continue;
         (*a)[j]=(*a)[j]-(*a)[j0]*g;
           }
       }
       for(i=1;i&lt;=m2;i++)
" E/ Y. c" @) j( K+ }; b$ w! \      (*a)[j0]=0;
       (*a)[i0][j0]=1;
       goto leap1;
}  
; U# }+ U5 p- j1 Q- R2 G- kmain()
{
  float a[5][7]={{0,0,0, 0,  0,0,0},
) Z& q: |% Y9 M/ R6 G7 V6 i                 {0,3,-4,3,  1,0,12},
% A; m  x! q3 B* u6 ^! i. t2 ^/ k3 a8 K' c                 {0,3,0, 6,  0,1,12},
                 {0,0,0, 0,  0,0,0},
3 B' B, X0 n& p/ V  R! n                 {0,69,0,144,0,0,300}};
  int k[3];
4 T9 y* A3 U5 q4 w& v3 }* H4 Y  float x[3];
  clrscr();
4 ?+ ^3 o. V9 `! W  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
9 R, d7 L2 b+ \: e}
9 h6 T5 E* E- f2 K9 ]: F, B</P>

plgatc

<>/*************************************************************************
2 A3 O) \! X7 {4 ^7 D+ l3 u                    单纯型法解线性规划问题（两阶段法）        
+ X. t- M' P. q7 q' i- Y6 T/ x
* l  {) b( V' J& P 编程环境：VC++6.0     
方程组输入说明：
变量非负，按提示输入相关参数。
*************************************************************************/
#include &lt;stdio.h&gt;
#include &lt;stdlib.h&gt;
#define MAX 100
#define STP 100</P>
<>int stop=1; //迭代记数变量
int status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
int step=1; //目前阶段</P>
<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
$ T7 ?3 m9 {, {6 zint num_x; //变量个数
1 O7 d: l2 ?' \( n2 Wint num_st; //约束方程数
3 \) A" n  w* o: Z: Jint num_ar=0; //人工变量个数
int arti[MAX]; //人工变量下标
3 o; k; R- T9 hint base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
' r! N5 H& j' ~1 ]( A* x<>void create(); //建立方程组
void iterative(); //单纯型法迭代
/ I, D% M! {2 V- v. U: {void output(); //输出结果
void banner(); //打印程序标题
void exchange(); //交换两阶段价值系数
void show(); //输出方程组</P>
% B6 i6 u6 P! A# Q" B<>void main() {
( r; U9 ~1 [; z. o+ d int i,j,k;
& M. j0 Y1 h* `7 O banner();
create();
//保存原价值系数，转换为第一阶段价值系数
for(i=1;i&lt;=num_x;i++) {
0 D8 P4 M! }1 M" ~3 w  k=0;
0 K9 M: p7 z! H: ?  S! e- @6 L2 h  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>
. d% {$ o" F8 o4 P- H+ C<> printf("\n\n第一阶段问题为：\n\n");
# j  i. R$ x* i; {  O% a' M- ? show();
" i- F+ r4 U3 W; | step++;
printf("\n\n按回车开始第一阶段迭代");
getchar();
getchar();
iterative();
: T! _4 y2 c9 _5 W) X/ i0 l$ B if(status==-2) {
& f0 x1 n, @4 J8 r9 z  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
3 j# p& D+ t& r0 q$ t6 h7 _  getchar();
/ i- b+ w( n# o+ ^" [' ?" Z  exit(0);
9 E5 Z0 D3 L0 W3 k  H9 I, U6 B }
output();</P>
. L  \% @! K/ f1 O3 i<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
) m9 w: O" W; R  J# G$ x  puts("\n按回车结束");
  getchar();
  exit(0);
: v- E: X9 ]' W9 }" q }</P>
<> //转换为第二阶段价值系数
5 q. G3 ~; l$ H5 U4 z7 _: A exchange(c,temp_c);
4 ~1 {1 E0 f" T1 N" E" ] //把人工变量列全设为0
. _) g1 t: _- j for(i=1;i&lt;=num_ar;i++) {
5 R7 _4 }6 W* Y  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
}</P>
<> puts("\n\n第二阶段问题为：\n\n");
show();
. X9 a8 C) \1 _" Z puts("\n\n按回车开始第二阶段迭代");
getchar();
/ |0 F9 I. Y9 ^$ E' D: ` iterative();
: S: ^2 ~  S$ l5 ]5 }4 S, P, d switch(status) {
case 1:
  output();
, W/ ~+ W% L4 E& E1 ~/ {7 ?  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
0 }4 O1 a( y5 q7 C) M8 H% G% @8 M; x  getchar();
  exit(0);
0 Q5 V, x" j, d8 U  D& C9 n& B case 0:
  puts("\n\n原问题为无界解。\n");
& N& i' ]! K# g# ]+ t  puts("\n按回车结束");
  X/ B  |$ I. J  getchar();
9 F# R6 l9 ~0 k, d1 u! b+ {  exit(0);
case -1:
  output();
  puts("\n\n原问题有无穷多最优解。\n");
; X% L: e" {3 W4 B4 s$ }+ w  puts("\n按回车结束");
  getchar();
# B# S- B2 I* t$ M/ X6 Y  exit(0);
case -2:
  puts("迭代超过限制次数强行终止！\n");
  \0 d4 ^4 a/ }; R0 A' d1 a) \, Q# G  puts("\n按回车结束");
; A- b# i# q4 \+ Z$ t2 i1 i% ~2 v$ Z  getchar();
! ]/ U. O4 z8 |/ W7 ]  exit(0);
6 A4 A* D9 L3 I% \6 G0 g. |2 z }//switch
/ t" F$ J' A) R& G! W}</P>
<>void banner() {
0 U) C. ^7 [7 ?/ z: t printf("\t\t****************************************\n");
8 B; ]" N8 [) U1 h+ R printf("\t\t         单纯型法解线性规划问题\n");
& E$ {1 B/ ]) e* H$ _: f printf("\t\t                         作者：Thunder\n");
printf("\t\t****************************************\n");
" P6 O  V7 h3 P  X7 f: G/ L printf("\n");
}</P>
<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
8 M# X: J+ k* \; ?; U//为1的不显示系数1，-1系数只显示负号
int i,j,k;
switch(step) {
case 1:
  printf("min z= ");
  printf("x[%d]",arti[1]);
  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
case 2:
* B- }' W9 ^5 J. a printf("max z= ");
* i2 Q9 N, L* e: h% u. g" Y printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
; r  C9 n- S8 k5 e  else if(c==-1) printf(" - x[%d]",i);
  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
3 t+ L' J8 {' `  ^. Q  else printf(" %lg x[%d]",c,i);
6 \- t5 B" d  P9 g( `$ M# J4 p5 P9 M5 y }
. ~1 `' o& y0 Q0 C9 j break;
, ]0 k7 f* X2 Y1 e6 a5 X- O }</P>
: \5 E. `+ C# ?) t6 U7 Z<> printf("\nst:\n");
4 N3 V- k2 `* D# D' V for(i=1;i&lt;=num_st;i++) {
  k=0;
' m* g8 T0 \# }. X8 v: g# w6 k4 e  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
: M& [9 b% Y% O" e$ D4 b  u: S" d    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);
: H4 t0 y2 k6 h% K; f/ x- y4 Y& w    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);
    else printf(" %lg x[%d]",a[j],j);
9 ~0 g4 P; Q; F; l& d: ~& s    k=1;
: ]3 s# I0 ~+ A' E: Y   }
' U9 Q) I/ Y8 `; r* x, Z" p  }
9 V' w5 f+ }, l' Y  printf(" == %lg\n",b);
, ~& x% E+ j/ _+ P2 b1 i }
printf(" x[1]~x[%d]&gt;=0",num_x);
: `1 f4 y, T7 T7 [}</P>
+ E7 D$ Y* ^: s8 [2 r* q2 V! Z<>void exchange() {
int i;
& X( E4 }; x/ u' Y. D2 {0 a double temp[MAX];
for(i=1;i&lt;=num_x;i++) {
  temp=temp_c;
  temp_c=c;
  c=temp;
$ ]) g- O9 `# X6 g, k0 A1 C' D  o }
}</P>
0 B6 O7 D) s; h0 Y4 Z6 i, `4 f9 A( v<>void create() {
//输入方程组系数，每个方程输完后回显确认
int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
char confirm;
2 m1 v" m. r. l, R! T; h% n
while(1) {
  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  scanf("%d",&amp;ma_mi);
) M/ K" Q' k1 y7 h: q  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
  else break;
}

& B5 w7 I+ z7 ^2 r: _$ V while(1) {
. l5 t& C4 f% X5 m  printf("指定变量个数：");
  scanf("%d",&amp;num_x);
8 s0 J, C2 y7 B* a/ g5 N) h  printf("输入价值系数c1-c%d:\n",num_x);
2 i% s# S2 z+ [  for(i=1;i&lt;=num_x;i++) {   
   printf("c%d=",i);
9 }1 U; A3 q; q' |0 }   scanf("%lf",&amp;c);
& j7 H* L9 j9 t& v' x% d! I: w  }
7 ?' O6 T/ P1 V5 M! M4 U  if(ma_mi==1) printf("max z= ");
4 ^1 Q& ^! w- T4 H4 L1 B+ R/ X  else printf("min z= ");
5 Q- ]2 B, o, x5 Y% c  P  printf("%lg x[%d]",c[1],1);
7 V( U" n7 ^" A( D+ X. N  for(i=2;i&lt;=num_x;i++) {
   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
; o$ b, ]: n* @. D; x# f   else printf("  %lg x[%d]",c,i);
5 ~7 R5 @2 ^0 V+ U* ?  }
printf("\n正确吗？：(y/n)");
. G/ ]' j; K/ O1 a- k% Y getchar();
confirm=getchar();
$ ?9 n9 o% ^0 n/ y; \9 ~ if (confirm=='y') break;
else if(confirm=='n') continue;
" X( u, @, V" P6 \6 ^6 z; U, i% J }</P>
<> printf("输入约束方程组个数：");
/ q, m5 I9 n5 k  o( e scanf("%d",&amp;num_st);
for(i=1;i&lt;=num_st;i++) {
8 U) q0 q" y4 C& s  e/ K: s" Z* H  printf("st.%d:\n",i);
4 a5 k5 d1 |7 ^/ F% {. V8 [  while(1) {
4 w8 k$ b4 a4 A2 B5 n/ r   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
   scanf("%d",&amp;re_st);
   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
0 {# }$ m  X* r( j( C7 p* e1 h! B   else break;
  }
  printf("输入技术系数：\n");
  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
  }
  printf("输入资源拥有量：\nb%d=",i);
  scanf("%lf",&amp;b);
/ \5 i8 A3 }, s& ^+ g  
# X" Q( o# o1 Z5 K: t/ D" `0 ?  printf("st.%i:\n",i);
  printf("%lg x[%d]",a[1],1);
7 c( o$ z8 N2 M& H% h3 K9 l8 s  for(j=2;j&lt;=num_x;j++) {
   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
7 q" x3 s- L2 m" [   else printf(" %lg x[%d]",a[j],j);
  }
3 L& Z/ S0 o7 a) X  switch(re_st) {
   case 1: printf(" == %lg",b); break;
  ?& Z5 w! f0 a- f! ^# {& |: ~# @( @! \   case 2: printf(" &gt;= %lg",b); break;
4 r3 r) L: c$ |: l- a   case 3: printf(" &lt;= %lg",b); break;
  }</P>
<>  while(1) {
   printf("\n正确吗？(y/n)");
   getchar();
   confirm=getchar();
   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
  }
}</P>
<>//显示输入的方程组
printf("\n原问题为：\n\n");
if(ma_mi==1) printf("max z= ");
3 R* H: D' ]; ~* b, z; g. G8 T else printf("min z= ");
printf("%lg x[%d]",c[1],1);
2 K- Q! Z% o# H- a2 {* r8 M7 p# ?1 E for(i=2;i&lt;=num_x;i++) {
7 A+ r4 c* n% F- S  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>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
+ I) a6 F( w3 ^5 P: O2 Q6 q( M# r  k=0;
  for(j=1;j&lt;=num_x;j++) {
5 s$ T9 X* Q" A% t' \   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
" K( p( f$ m- G  h. b4 V$ v% O  ]    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
% G2 r& ~$ G& j$ m" s5 Y+ @  r    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);
: z& l; s/ {! e    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
! ]- L* J6 ]( i' j+ @    else printf(" %lg x[%d]",a[j],j);
    k=1;
   }
9 X6 e2 W5 T8 N( f) G  }
8 X- [6 T( e3 B! F+ d5 J* N5 S  K  switch(re_st) {
   case 1:
    printf(" == %lg\n",b);
5 s* m5 j( B$ s4 ?$ Q3 C    break;
   case 2:
    printf(" &gt;= %lg\n",b);
    break;
" P9 q( q7 W$ i1 I+ ?2 `   case 3:
: v! d+ n, p) E& E2 P0 t    printf(" &lt;= %lg\n",b);
( R1 A; I2 I# I) M    break;
  }
}
printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
<> tnum_x=num_x;
7 Z/ A: P6 y; b/ Q9 A; F/ {/ y  G for(i=1;i&lt;=num_st;i++) {
/ K' T8 e0 f% \% A9 f  switch(re_st) {
  case 1:
& u+ o9 n6 J$ A) h9 E  Y  case 3:
* ?$ c" ~2 E  V   num_x+=1;
  H9 D$ ^8 z& [5 B- w   break;
# n, R5 E% m$ X5 ~; ~1 R9 ]  case 2:
   num_x+=2;
   break;
  }
8 {) L9 S! t# D$ H" K* [ }</P>
# v/ N6 J/ q' @8 b<>//化为标准形式
if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
% i3 h7 q) E6 F6 p1 y. y   case 1:
& T6 q6 P+ X8 e8 b    num_addv++;
1 }, h; }' T4 S    num_ba++;
9 v' }9 d& g+ X4 M* z* i+ Q    num_ar++;
% }% A# S" g# \% l" s- ^/ M+ L    c[tnum_x+num_addv]=0;
) Z$ ]( H/ G1 I7 d. Q( ~& e    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
7 Q* Q( ]$ B5 N    for(j=tnum_x+1;j&lt;=num_x;j++)
     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
& q' d; V7 G9 Z0 J+ `     else a[j]=0;
    break;
   case 2:
7 _& k6 O% ?* Z3 [2 Y. D    num_addv++;
- P7 g  b7 |# Y    c[tnum_x+num_addv]=0;
. g- x! D# d1 @    num_addv++;
    num_ba++;
! |! n3 Y1 j1 J! _4 M% k    num_ar++;
0 F% S, O' a$ I3 s/ P    c[tnum_x+num_addv]=0;
  H+ J4 @) |5 ?, h    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-1) a[tnum_x+num_addv-1]=-1;
     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
* Z( D" {, j2 N3 T# M8 u0 o     else a[j]=0;
    break;
$ V9 e2 K1 W. W' x6 [) ~   case 3:
    num_addv++;
; r" d7 y; j; a$ z& Y) V    num_ba++;
9 n: H7 D+ C9 L( i# [* T: N    c[tnum_x+num_addv]=0;
    base[num_ba]=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;
6 I* e7 ?/ G) k" \6 n5 i5 n1 A: l! h" g    break;
  }//switch
}//增加松弛变量、剩余变量、人工变量、确定基变量</P>
6 k" F( v. ?4 p! r<>//显示标准化后的方程组
printf("\n化为标准形式后：\n\n");
' w( X% M4 }1 Y$ m5 p if(ma_mi==1) printf("max z= ");
- _8 m& R& J0 |/ @3 M else printf("max z'= ");
printf("%lg x[%d]",c[1],1);
6 z3 w' P1 J  ~2 D( J% D for(i=2;i&lt;=num_x;i++) {
. y# k1 m% D4 L" f; m- e7 q  k=0;
- t: N7 I0 n/ d- Z/ c( ?% Y  for(j=1;j&lt;=num_ar;j++)
/ a5 B$ }. R: L( H  if(i==arti[j]) k=1;
: h0 R8 L5 ]4 R! ?& g; O3 i. F  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);
  z5 Z) E1 ]3 [& @$ t }</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
. ]1 f) y4 s. z, @, P   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
( Z/ O% Y- O! G: |    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);
2 q6 c# c# W, A- 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);
    k=1;
   }
/ T1 V: I# y) L6 Z" s  }
  printf(" == %lg\n",b);
% E. b3 M* ^9 X- i2 v }
printf(" x[1]~x[%d]&gt;=0",num_x);
- h# T& ^/ }/ m# w  f: L9 q" ~& J- ^}</P>
<>void iterative() {
0 p  e0 D% J, c+ H2 i8 U int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
* i6 V; r$ i  ~5 w# j" e3 \- I% O int base_elem;
int base_out,base_in;
double sigma[MAX],temp;
  x* S0 c; r# d% Y double value_be; //高斯消元里保存主元素值</P>
<> printf("\n\n第%d次迭代：\n\n",stop);
for(i=1;i&lt;=num_st;i++) {
1 }& N0 r& _" A/ c: E! W  printf("c%d=%lg\t",base,c[base]);
  printf("b%d=%lg\t",i,b);</P>
4 D+ K$ f# S7 z! h5 n<>  switch(step) {
( R3 ^4 g; f1 v8 t   case 1:
: v* d* V! k1 c- [- D    for(j=1;j&lt;=num_x;j++)
    {
1 Y. E* i, K3 N3 v4 }     printf("a[%d][%d]=%lg\t",i,j,a[j]);
( r2 ?! R+ m/ B2 y8 G1 N; J' \; d8 f    }
    printf("\n");
) t: `6 F4 C! ?$ W" E    break;
3 Y1 I9 W! r' L" \2 h+ H8 x. ]* _   case 2:
7 O  A- F5 M( S7 i    for(j=1;j&lt;=num_x;j++) {
/ b$ x* Y+ v! y, I- D     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]);
    }
    printf("\n");
    break;
. Z/ r' e) a+ Z' O9 l- V/ h7 ]( u/ W  }
}
//求检验数sigma
for(i=1;i&lt;=num_x;i++) {
  sigma=c;
) N% v% `1 l$ Z  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;
( g( h) w$ n0 w$ t+ a  switch(step) {
: C; q( c5 @  v' I* [   case 1:
    printf("sigma[%d]=%lg\t",i,sigma);
    break;
   case 2:
6 p; @/ m, O. L: ^: X& I    k_a=0;
    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
8 W# ~6 A. t8 u1 M' {! }    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
  c( ?0 P2 x* n5 {/ n    break;
7 h# W- V% N) q8 s) S0 A4 a  }
3 E; k4 q* X' N3 w/ |2 ^3 Q }
6 A# Z3 r4 X- x putchar('\n');
6 K2 T% \! {; w& G//检验检验数sigma是否全小于等于0
k=0;
1 d) K* y5 m0 C for(i=1;i&lt;=num_x;i++) {
  if(sigma&gt;0)
  k=1;
}
if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
& f5 K" V% L% j* g- @- m* S4 \  for(i=1;i&lt;=num_x;i++) {
  n7 `  ~, R7 v0 W5 U   k_f=k_a=0;
5 D5 W  r1 q$ n! O! t   for(j=1;j&lt;=num_ar;j++)
3 N, F5 f  ]. K! _6 b   if(i==arti[j]) k_a=1;
   if(sigma==0&amp;&amp;k_a!=1) {
" r2 q  m* |: B4 D; k    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
/ ?3 ^) ^- `- H  }. ~   if(k_f==0) {status=-1; return;}
/ W8 O6 q  {) r5 D; v- Z3 F, `6 c   }
4 f# T# a1 ]7 D( m  }
  status=1;
/ P- b$ U4 e9 F7 J0 I  return;
2 S7 ]; v: W8 D0 |* f2 R }
* @* {2 ^: X8 G+ `+ j/ o! E//检查是否为无界解
/ y+ _( i7 w$ U" X8 I7 _ for(i=1;i&lt;=num_x;i++) {
2 L4 u; }3 [7 M; q  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;}
0 B6 d9 J6 D0 p  }
}</P>
<>//确定换入变量
3 x, d# f% p/ o8 B  |) b for(i=1;i&lt;=num_x;i++) {
$ B# n9 M$ d/ i& \$ O  k=0;
  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赋初值
5 d* T$ I6 i4 K/ y- R* [8 t for(i=1;i&lt;=num_x;i++) {
2 d. Z5 ^  z% p. p  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
, |' n( g6 A( ], R1 a! N  if(k==0)
   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
3 x; A8 a0 J' w8 _6 E4 u: C; V    base_in=i;
- X* p( i) K2 z+ C2 _* |    temp=sigma;
   }
}</P>
% b& r9 K* Q2 T8 V& ~<>//确定换出变量
for(i=1;i&lt;=num_st;i++)
  if(a[base_in]&gt;0) {
( k4 j( }4 P; Z' i/ r2 i   temp=b/a[base_in]+1;
: a4 V* Q- h$ l. S* g   break;
2 H8 _/ R5 T9 _) ?4 L( Z" R( Q  }//temp赋初值
% E( o7 @, V% [& l4 Q9 i7 B# j for(i=1;i&lt;=num_st;i++) {
0 p" O8 q. F' {8 h+ ]  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
, \( {) Y! K% T9 j+ l/ s4 E1 W   for(j=1;j&lt;=num_ar;j++)
4 p9 O* c! \- P" r    if(base==arti[j]) {
     base_out=base;
     base_elem=i;
     temp=b/a[base_in];
/ ]9 f" G3 o: G1 t  x+ F+ g     break;
  s: j9 u$ z6 Z5 \7 F$ J    }
! J% F6 R5 ~7 @2 I" N' }( ^# D  }//人工变量优先换出
  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
+ e0 [% |. A) F2 W- K2 e) W   base_out=base;
   base_elem=i;
- A0 l# T; e+ k5 ?( }. f   temp=b/a[base_in];
  }
9 M; A) o, B( y. c/ b }</P>
<> printf(" 基变量：");
for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
) M1 X; _; ~  J9 n2 ?//基变量变换，进行新方程初始化后迭代
for(i=1;i&lt;=num_st;i++) {
, y8 m9 [1 \1 q; X& P' V  ]  if(base==base_out) base=base_in;
3 C+ ]6 b; z( q- O7 a. R% O9 |) H }
4 J( b& C: D- J' ^' I//初始化主元素行系数
( u8 |1 z2 {7 x  {) O5 t' V value_be=a[base_elem][base_in];
b[base_elem]/=value_be;
9 X  y3 P5 E5 {* F- E for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
<> for(i=1;i&lt;=num_st;i++) {
  if(i!=base_elem) {
   b-=b[base_elem]*a[base_in];
0 E& B4 u* ]  D" b  F1 C   value_be=a[base_in];
   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  }
) a8 Y' \& z. W }
stop++;
* J% l; y: O6 R4 x if(stop&gt;STP) {status=-2; return;}
5 U3 ]. h: ^$ o* i: o iterative();
}</P>
<>void output() {
int i,j;
double X[MAX];
printf("\n结果如下：\n");
5 A- z4 ]; g& d printf("\nX=(");
for(i=1;i&lt;=num_x;i++) {
" |% X' n, R  Q  i" `  for(j=1;j&lt;=num_st;j++)
  if(i==base[j]) {X=b[j];break;}
6 r. e) ]  U+ Y; A; |3 j5 g# m  else X=0;
  printf("%lg ",X);
}
printf(")");
! r$ E+ T. B/ Z* ^, `& F for(i=1;i&lt;=num_x;i++) max+=c*X;
if(ma_mi==1) printf("\nMax z= %lf\n",max);
+ X5 U; m% m% ?: z else printf("\nMin z= %lf\n",-max);
}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了