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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
#include&lt;math.h&gt;
' G8 Y9 E; p& {#define X 5
#define Y 7
8 @: q: _: Z6 u/ C# J" yvoid xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
( O4 Y' K" b) I# W& C  ?' s" L{
$ K& }$ q0 g6 [  int j;
  *c=0;
  for(j=1;j&lt;=*is;j++)
# j! k$ O6 @/ n- \/ D    {
/ W: {! V* D" O) o) t0 ]      if((*a)[*ir][j]-*c&gt;0)
    {
1 C' X( r: C- u. E, j      *c=(*a)[*ir][j];
+ m7 z% H6 k2 \. q8 b# N  n; \      *j0=j;
    }
     }
}</P><>      /***************** 参数说明 **********************/
      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
; E$ M8 [1 j2 a! ?      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
! U+ V# U5 A4 m* W: b) B+ G: A      /********* <a href="http://happyyangxu.home.sunbo.net*******/" target="_blank" >http://happyyangxu.home.sunbo.net*******/</A>
+ d3 r7 {1 E- n! y6 n0 I- Eint xi_sm(int m,int n,int m2,int mn1,int l1,float (*a)[X][Y],
       int(*k)[],float(*x)[])
1 b; y& S% K- e% ~1 m) W0 p{
- s* e8 k3 r5 c9 A! w$ H  int i,m1,mn,j0,i0,j;
  float c,g;
  m1=m+1;
  mn=m+n;
# D+ r9 P8 n" x6 Y  for(i=1;i&lt;=m;i++)
2 }/ J( I; x1 J) j    (*k)=n+i;
leap1:
  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
# G2 U! Z  ^: `3 [1 H  else
     {
leap2:   if(l1-50==0)
      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
+ z0 ~7 k+ |+ J, F2 P' v5 y       else
% L8 Z2 O# w7 ?      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
. U# ^  e3 S; x( p: R9 }% E     }
/ C& v. b$ U; N3 Q$ z7 B1 d. E   c-=1e-8;
) @2 g: K$ T4 v. E7 ]   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
* N9 k: k9 Q* W5 m% p# ?+ @  G3 g$ E      {
0 G: @8 m1 E+ V0 {3 |- n+ l/ O+ L' i    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
    return(-1);
2 |/ a6 e) l9 z% U% R      }
) l' {) a! _, ^. z+ ?( A2 L( @/ P   if((c&lt;0)&amp;&amp;(l1==1))
0 y- p. }9 p- R% {8 `       {
     l1=50;
9 b* V# \4 U; A2 e+ g+ K  `% x     goto leap2;
       }
4 W, h7 W, o+ z. x7 p5 X    if(c&lt;=0)
1 b5 A6 n6 F, C8 t! @* U$ t       {
      for(i=1;i&lt;=m;i++)
$ |5 a) `  F. b9 r5 K( w' R$ E# b         (*x)=(*a)[mn1];
% p4 t; B$ d) {; i         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]);
      return 0;
: [7 |" U( u% F5 ]    }
- \2 ^: `4 o! X# h! a* N     c=1e8;
     for(i=1;i&lt;=m;i++)
       {
! u' U: W3 m0 X4 \     if((*a)[j0]&gt;1e-8)
: J0 {+ c7 X, u3 V        {
$ v4 o, V6 R, L          g=(*a)[mn1]/(*a)[j0];
4 ], P  P: Y, G  K: B          if(g-c&lt;0);
" Z. ?3 Q# r7 c* E) M         {
           c=g;
4 T$ b; \2 ^, T           i0=i;
, S! d5 J2 N+ ?( ^          }
( A1 C/ B; p; i+ I1 ^         }
/ L% P( q/ M4 U2 p) [7 {    }
2 r* b4 n% _) J8 z+ I" B# q; v2 ]4 Q      if(c==1e+8)
    {
5 x3 q& l' k/ B1 L' a      printf("\n\t*******Lp no solution********\n");
2 H/ d2 a  H: A  V# S& E3 Z9 z      return -3;
& z# G5 i& u/ Q+ N% N( o     }
4 v( O8 a+ m  }; u, l: u       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
      {
        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
           continue;
         g=(*a)[i0][j]/(*a)[i0][j0];
" e) O. T1 v6 n; f           (*a)[i0][j]=g;
  R6 D% h" K) ^  _) q6 ~7 J         for(i=1;i&lt;=m2;i++)
           {
2 ?! Q, x4 M( o6 g2 c         if((i==i0)||(!(l1==1)&amp;(i==m1)))
1 M5 r$ F) _( d6 p0 y! v            continue;
         (*a)[j]=(*a)[j]-(*a)[j0]*g;
           }
/ Z+ f% }1 s  S* E6 K) Y2 U       }
# U5 ~% ?8 P  b6 U: C( n       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
       (*a)[i0][j0]=1;
3 F! p  X& e8 {5 {. z       goto leap1;
5 R* @/ C  F  Y; X  g% J/ k}  
main()
{
  float a[5][7]={{0,0,0, 0,  0,0,0},
! W* q2 z1 G0 `8 c; E5 C                 {0,3,-4,3,  1,0,12},
9 \& L, ~, Z' M. r- |                 {0,3,0, 6,  0,1,12},
                 {0,0,0, 0,  0,0,0},
                 {0,69,0,144,0,0,300}};
  int k[3];
" ^) g3 k3 e7 ^* m$ e& L" |; G! j  float x[3];
  D4 x* f! Q4 g* ?) ^3 T  clrscr();
2 C' D' b. k. h  ^  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
! @+ I$ t0 O) m$ t9 l# W. C* `}
1 Y1 E  n/ b6 Q# {6 k</P>

plgatc

<>/*************************************************************************
& o, e; a& b# [& M5 G                    单纯型法解线性规划问题（两阶段法）        

* q# l# H8 A: U' f# l8 n% e( i 编程环境：VC++6.0     
; B9 P2 e* R- F: G* D9 x 方程组输入说明：
变量非负，按提示输入相关参数。
*************************************************************************/
/ a  ?8 U# e1 U* t/ m#include &lt;stdio.h&gt;
* E0 }) t3 i9 x' q+ C: ^0 [#include &lt;stdlib.h&gt;
#define MAX 100
, w9 ?, x7 ?  H" G. _#define STP 100</P>
8 \2 s3 ^/ l2 d& }2 L9 {! O" F$ W<>int stop=1; //迭代记数变量
int status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
: a3 R% b- p1 _5 J( f. |int step=1; //目前阶段</P>
+ k# N/ p* z  N<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
5 \  O. ]  a) o6 H. [int num_x; //变量个数
+ D( L# d, k: ^" M5 `* n5 N8 Oint num_st; //约束方程数
4 \: k& E, V; rint num_ar=0; //人工变量个数
! X8 e9 |% W6 S$ n0 ^! \( ]( tint arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
2 Q) P9 w4 d. b+ y. Y<>void create(); //建立方程组
void iterative(); //单纯型法迭代
void output(); //输出结果
( `; B& P$ [, Z! Z/ uvoid banner(); //打印程序标题
+ P- F7 c# m+ n; A5 k# Ivoid exchange(); //交换两阶段价值系数
. `' [% c- z% U4 n; C4 t) cvoid show(); //输出方程组</P>
<>void main() {
int i,j,k;
6 g& ?' `( a3 @9 t* S1 w# A2 e banner();
- O" g/ P* P4 Z create();
//保存原价值系数，转换为第一阶段价值系数
# p( c  O3 ]# n5 i for(i=1;i&lt;=num_x;i++) {
# j9 {! F! X( |+ _  k=0;
  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
  if(k==1) temp_c=-1;
" _; f% s9 P  G7 V5 O& `' U  else temp_c=0;
}
: h* ^5 q$ Y' A- x6 ?0 [ exchange(c,temp_c);</P>
<> printf("\n\n第一阶段问题为：\n\n");
show();
  x. k0 u4 n: M( K- ]) S step++;
printf("\n\n按回车开始第一阶段迭代");
getchar();
: y) z$ J' C  |' I getchar();
iterative();
/ e4 \0 B9 l, y# T4 y if(status==-2) {
" i. }: }; k, \, k! ]" e0 j6 Q* F  puts("迭代超过限制次数强行终止！\n");
" o" C7 q/ U! A# w  puts("\n按回车结束");
: d1 B' w6 v8 b" v) N9 A2 \: ^  getchar();
  exit(0);
/ M: J; ~# m' N5 \ }
# m# r" L: k6 P% b  M8 i output();</P>
<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
  puts("\n按回车结束");
5 x6 D0 u! i9 l. f  getchar();
, ]5 I7 E0 d5 W) J7 N2 L  exit(0);
}</P>
) }$ j( o2 C+ I, d<> //转换为第二阶段价值系数
exchange(c,temp_c);
/ W" s) @/ G* ?3 l- b: Y //把人工变量列全设为0
+ u7 w( T$ C& u7 n4 }# a' @% Y for(i=1;i&lt;=num_ar;i++) {
. O- u$ h! {, B1 X  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
& }- I, K8 G1 ~) }6 V }</P>
* q4 X1 n' y* `$ l0 P) L<> puts("\n\n第二阶段问题为：\n\n");
show();
puts("\n\n按回车开始第二阶段迭代");
getchar();
) s, I/ C2 K/ @9 T iterative();
switch(status) {
$ b/ C4 o3 ]9 Z" N8 K case 1:
  output();
  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
: Q+ m7 k4 |3 C1 W5 _  getchar();
( C; U' K( q2 [' |- E: I* Q  exit(0);
case 0:
  puts("\n\n原问题为无界解。\n");
  puts("\n按回车结束");
) w1 ]% P; C- m; Z  getchar();
/ X" O0 v/ j( T, y% t$ k1 o: ~* G& F  exit(0);
2 u  d  {  c$ e& M7 G case -1:
! @; s+ j  q4 Q% ?/ z+ A7 E0 J  output();
# N+ Y+ t) _# t- n( a  puts("\n\n原问题有无穷多最优解。\n");
1 Q" Y: Y: `1 a) u$ F  puts("\n按回车结束");
9 T1 a% `& Z& k  q  getchar();
  exit(0);
case -2:
+ t( Q% @: T) r% N5 F& l6 v  puts("迭代超过限制次数强行终止！\n");
5 a' k4 p9 @6 W8 ^' {2 O  puts("\n按回车结束");
  getchar();
2 q6 x' K6 h! N+ Z- `8 U8 x  exit(0);
# l1 J2 j' r- _, t }//switch
& K9 m4 @7 ?! D( c% m( H) C' B}</P>
' a/ P5 {4 i7 f! A  U* T<>void banner() {
  I4 t, a( b- z- \- V! {# i printf("\t\t****************************************\n");
+ e8 r( j; d3 C( B, } printf("\t\t         单纯型法解线性规划问题\n");
- V' ]! ?9 {4 l7 ` printf("\t\t                         作者：Thunder\n");
printf("\t\t****************************************\n");
* b* [7 ?  o0 t printf("\n");
7 z- V- V" ?8 m( D( i}</P>
<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
//为1的不显示系数1，-1系数只显示负号
% X  O% U$ f- D- Q# u; i int i,j,k;
7 K# J( ^/ V  x$ a4 ~ switch(step) {
- R( ?! M* s' L9 t' L9 T case 1:
  printf("min z= ");
  printf("x[%d]",arti[1]);
% p# b% `$ \1 e  e% b: g  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
case 2:
) r6 ^2 z( e" c2 }8 C& e$ A# ^ printf("max z= ");
printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
4 E; ^) }( l8 V3 v8 F  V5 R( G1 K2 f  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
; \/ F: R( r- x) w) w* Y/ W* I  else printf(" %lg x[%d]",c,i);
- S2 H% F4 I# L7 U5 Q }
break;
}</P>
$ P! \: Y4 K0 W1 }& `7 v<> printf("\nst:\n");
5 d, T, A$ J' a. [. c( Z7 o' B! C for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
( w  M; M0 b4 q1 ^    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
7 q4 C+ ]5 @8 u0 C; ]# l- t. f    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
  I0 k6 t7 a/ b9 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);
  E7 Q0 H$ c/ L4 y- q, W8 A3 X: q* }    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
5 i/ X! u" [1 \9 ?; i# Y1 y    else printf(" %lg x[%d]",a[j],j);
    k=1;
   }
  }
  printf(" == %lg\n",b);
8 i( C6 U! ?3 y/ @* C4 m }
" Z9 r5 y, u! | printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
<>void exchange() {
; s& \7 N: {6 i  o int i;
5 e8 @* h7 q' O7 }- h double temp[MAX];
3 e) a6 v+ C7 @2 w% } for(i=1;i&lt;=num_x;i++) {
  temp=temp_c;
" p" k! `- l1 O  temp_c=c;
6 u7 l* y' ^" k* }. B$ N: w6 c  c=temp;
! |5 x1 ?  J1 E9 z$ a$ o }
7 P& Y: g* I; l: N) u& D% n}</P>
2 p. Q0 a9 A) R7 c8 C. V+ ?* A<>void create() {
//输入方程组系数，每个方程输完后回显确认
int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
- ?, F% w( V1 W0 S( `1 s1 v/ S char confirm;

; @' X* L2 B/ [ while(1) {
2 S; W3 ?+ [3 j  printf("请选择：1、求最大值，2、求最小值：(1/2)");
; }/ j7 @8 C6 m  scanf("%d",&amp;ma_mi);
, n, |" i% Z6 U- p. O$ G  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
  else break;
$ X9 w7 ]+ F4 G6 V& x0 K }
& Y* }: H+ v! i% D
) @5 N$ ]6 V# b6 l6 h while(1) {
  printf("指定变量个数：");
  scanf("%d",&amp;num_x);
  printf("输入价值系数c1-c%d:\n",num_x);
4 Y* G+ A: \) ]' E5 ~" w  for(i=1;i&lt;=num_x;i++) {   
7 B' x( A* D( L1 t   printf("c%d=",i);
) m. x1 S) z) s& j2 C   scanf("%lf",&amp;c);
  }
  if(ma_mi==1) printf("max z= ");
  else printf("min z= ");
0 b& F. M+ N& a3 N, {  printf("%lg x[%d]",c[1],1);
  for(i=2;i&lt;=num_x;i++) {
$ C) C( E6 x7 C! ]$ x   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
   else printf("  %lg x[%d]",c,i);
  }
printf("\n正确吗？：(y/n)");
: L) E. t2 K) |: T getchar();
confirm=getchar();
if (confirm=='y') break;
else if(confirm=='n') continue;
}</P>
! o3 s6 @8 o3 l  E* Q$ }9 y- x<> printf("输入约束方程组个数：");
scanf("%d",&amp;num_st);
for(i=1;i&lt;=num_st;i++) {
  printf("st.%d:\n",i);
7 {' Q# D( o5 Z# `  v  G  while(1) {
* l5 @1 D) i  M2 {5 Z   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("输入错误，请重新选择。");
   else break;
  }
+ C7 ]6 p- O, l$ @  printf("输入技术系数：\n");
  for(j=1;j&lt;=num_x;j++) {  
1 S/ D' g! U6 e' m- d! w3 z   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
  }
, {8 H1 i/ I4 y% ]0 {+ s  printf("输入资源拥有量：\nb%d=",i);
5 ~( Q( y; f  X7 O& ~! ^6 p  scanf("%lf",&amp;b);
  
, e( \' x6 x$ K* x& ]8 `( k  X) M6 E  printf("st.%i:\n",i);
  printf("%lg x[%d]",a[1],1);
  for(j=2;j&lt;=num_x;j++) {
4 k) C! E* i+ ]/ |   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
   else printf(" %lg x[%d]",a[j],j);
) k5 i) c  ^& Q  T$ M/ F* y  }
  switch(re_st) {
   case 1: printf(" == %lg",b); break;
- i5 R8 {  k! E9 }! ]" ~; B0 @' P   case 2: printf(" &gt;= %lg",b); break;
4 r7 U) ]6 G+ ?! i# Z. _   case 3: printf(" &lt;= %lg",b); break;
  }</P>
8 [1 `( _$ J  P) c1 {, h<>  while(1) {
& u1 u, R0 C& G% O3 x* y   printf("\n正确吗？(y/n)");
& {1 B- T  N/ P% J/ A8 w( _   getchar();
   confirm=getchar();
   if (confirm=='y') break;
- y2 y  \7 G; p* z3 T   else if(confirm=='n') {i-=1; break;}
  }
9 Y7 X. ]% E+ S- p" T+ Z4 f }</P>
7 \4 Y3 v$ f1 x1 c) B  ^: v1 ~<>//显示输入的方程组
  g9 a- o+ B  G4 u6 L printf("\n原问题为：\n\n");
# ?( z' z2 S; N. F% C3 w if(ma_mi==1) printf("max z= ");
* y4 g3 l  [' |: } else printf("min z= ");
  t* _5 a0 ], ~. Z# a5 k" y printf("%lg x[%d]",c[1],1);
! `0 O! v. x4 j& J  u/ Q1 e for(i=2;i&lt;=num_x;i++) {
$ \; v0 p* g% y! }) \* V  if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
: @7 g' ~& n" B& P; [7 G  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
3 Q* z# H& D" \. x: A! T& M  else printf(" %lg x[%d]",c,i);
}</P>
0 I8 t! \) ]2 v  i<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
7 S. O" o4 \2 r: `. u$ G   if(a[j]!=0) {
7 C; x5 j3 J& V& Z    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
0 n3 r5 I: D, t& ~" j* X1 y    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
    else if(a[j]==-1) printf(" - x[%d]",j);
0 H. ?+ U, h0 Q$ h    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);
    k=1;
   }
  }
& f; d* u, \/ l9 |# t" Q  switch(re_st) {
$ l" I) u. ^4 Y1 u   case 1:
    printf(" == %lg\n",b);
    break;
% h7 _' \- f' j- G/ b/ t/ N$ e  E   case 2:
    printf(" &gt;= %lg\n",b);
, r: U  C% e( F& b    break;
   case 3:
# I. W+ F) [, b; O/ `+ O$ m7 `    printf(" &lt;= %lg\n",b);
    break;
. N* f+ U& p- m5 d+ w: A  }
; a$ q8 W  D7 y  m6 ? }
# o" j) R/ z( N8 i7 p 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) {
8 {+ o/ y; x% [# C( C  case 1:
  case 3:
6 A) d% s: Z9 u2 E/ R! }6 X   num_x+=1;
7 i4 ~- a$ f' G3 ^6 g8 k   break;
  case 2:
' f' J8 U/ A( U+ f0 r% W/ V: B1 u- i   num_x+=2;
$ E7 B/ L% J. _; U. L* h9 g   break;
2 |6 J" ^& ?2 h  y, Q2 j* o; T' u  }
}</P>
<>//化为标准形式
6 ?' Y' P: q( u5 x( m if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
for(i=1;i&lt;=num_st;i++) {
2 i% U$ ^& `- C  switch(re_st) {
   case 1:
    num_addv++;
    num_ba++;
7 l0 m3 z) T6 _+ g9 U! ^/ g* E    num_ar++;
    c[tnum_x+num_addv]=0;
3 N$ x& a& t- A& L5 {4 r    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;
, M- N9 A& a8 T3 V& W     else a[j]=0;
    break;
5 ]$ f8 c0 z8 [" o   case 2:
) y, A3 t; Z& y) Q    num_addv++;
    c[tnum_x+num_addv]=0;
    num_addv++;
    num_ba++;
2 Y1 A+ J" L9 Q6 U. V1 `    num_ar++;
    c[tnum_x+num_addv]=0;
! m, s( [$ A( Q8 F5 p$ A    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
7 I0 E: o; m; [5 [2 G# E8 ^    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;
3 D, e' C, Z5 a& u     else a[j]=0;
    break;
   case 3:
1 x% M7 T. o, Z    num_addv++;
. {& w9 j" D7 v0 R5 X    num_ba++;
    c[tnum_x+num_addv]=0;
3 R( v. {# q5 o7 v) l. }    base[num_ba]=tnum_x+num_addv;
% B; [9 n, L! U, H+ z3 p    for(j=tnum_x+1;j&lt;=num_x;j++)
; {. W6 u& t; y% S     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
) _! s3 }" x! l$ f' i* P    break;
8 w( A+ }2 @) {; b  }//switch
}//增加松弛变量、剩余变量、人工变量、确定基变量</P>
<>//显示标准化后的方程组
printf("\n化为标准形式后：\n\n");
if(ma_mi==1) printf("max z= ");
' s( l, q& S! O; s3 r else printf("max z'= ");
% u; G; t% u2 i. o) ? printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
* m3 @- o) b) i. I5 ~  k=0;
  for(j=1;j&lt;=num_ar;j++)
  if(i==arti[j]) k=1;
- r+ U, f5 r- H7 e" S  if(k==1) printf(" -M x[%d]",i);
# R( c2 P* D! s8 {6 k- ?  else if(c==1) printf(" + x[%d]",i);
& g0 a) R9 A( X  c  else if(c==-1) printf(" - x[%d]",i);
6 \: Y0 e; _! t3 w$ P  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
( o2 M7 k: c+ `3 j" [' P  else printf(" %lg x[%d]",c,i);
7 }. k7 D* L2 K }</P>
3 r- Q# }5 l% y5 d<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
& h, Z9 l3 L& _( k. C1 W  k=0;
8 u! Q/ l7 G$ ~9 [( v& J% K  for(j=1;j&lt;=num_x;j++) {
  a2 m0 r/ o, l   if(a[j]!=0) {
0 Y' i* c5 s/ W, X9 J    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);
/ w0 m/ |# b, L# C- C( L/ z% I; L    else if(a[j]==-1) printf(" - x[%d]",j);
- v: l: J- ~  t' G" ]  c    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);
3 s( A  `! b7 l, X    else printf(" %lg x[%d]",a[j],j);
    k=1;
2 B/ i$ `$ e6 J   }
- \3 o3 e2 k9 U3 h# a. d  }
  printf(" == %lg\n",b);
}
( c1 z6 X9 q, ]/ J3 S; w printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
! \) a# U) O" ?* K4 C<>void iterative() {
& U1 `3 d! |/ H9 V) b; H& N7 Q7 ] 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;
! ~$ M* U& @! D- Z! X# u! j1 ] double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
<> printf("\n\n第%d次迭代：\n\n",stop);
. y/ r4 H: ?. u$ L; y for(i=1;i&lt;=num_st;i++) {
9 q5 `, q2 Y! H: [* J2 h2 M, N  printf("c%d=%lg\t",base,c[base]);
( i- l& K1 }+ |9 ~2 q$ C  printf("b%d=%lg\t",i,b);</P>
<>  switch(step) {
   case 1:
    for(j=1;j&lt;=num_x;j++)
9 c6 r8 n4 ?) `& }% s1 G4 }    {
# w& T: C* \$ b8 J5 e4 f8 _     printf("a[%d][%d]=%lg\t",i,j,a[j]);
- I5 x+ \& U, }    }
    printf("\n");
- ]0 r( ~9 [  N( [* J, Y8 M8 c' }    break;
! F" G0 \6 O* A, G. V, ]2 V   case 2:
    for(j=1;j&lt;=num_x;j++) {
     k_a=0;
0 n! r: a" ]- v; |2 n# d     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]);
    }
% [6 z* v) X; |6 T+ P    printf("\n");
1 p3 Y& \: p& ]# m$ Z3 K" @    break;
: {" Y4 {: I' T8 }# z% k' @2 L  }
}
//求检验数sigma
1 F7 W' v# U6 M$ c9 ^2 I  @ for(i=1;i&lt;=num_x;i++) {
- Q# Z# O; h  b  Y: y8 o* M$ W  sigma=c;
  for(j=1;j&lt;=num_st;j++) sigma-=c[base[j]]*a[j];
; i- M" V" c5 t) x  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
! P. _2 @/ I, Q+ u* X7 o9 ~1 O. K  switch(step) {
   case 1:
9 z& p; L' z5 E# I9 y+ b0 \    printf("sigma[%d]=%lg\t",i,sigma);
    break;
( }; [" t' w- |* F! M   case 2:
    k_a=0;
    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
- a( _; q8 n7 \- A; W$ Z8 w+ W    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
, D) m6 c: X% |) H    break;
7 ]) G  K' n# f  E9 q2 f5 K. t" ^  }
2 M8 }7 F0 m1 M8 D }
, J3 ^: o& x9 W9 ^ putchar('\n');
- e9 I  \: p# l! A5 K0 n5 J3 i//检验检验数sigma是否全小于等于0
k=0;
for(i=1;i&lt;=num_x;i++) {
2 s0 U( S4 Y+ J' B3 Y7 c% o( x  if(sigma&gt;0)
  k=1;
}
if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
) f# [8 |1 {7 L2 i  b; C  for(i=1;i&lt;=num_x;i++) {
4 t4 Z9 s% r! }; O# m" @   k_f=k_a=0;
  g7 B8 l5 |- g; j$ [  _   for(j=1;j&lt;=num_ar;j++)
   if(i==arti[j]) k_a=1;
   if(sigma==0&amp;&amp;k_a!=1) {
    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
   if(k_f==0) {status=-1; return;}
   }
( _2 E4 w6 L2 Q: l9 f2 B/ i  }
  status=1;
- w: Y. c. D4 c5 {  return;
; I% u7 w+ U* e7 f4 e/ [ }
//检查是否为无界解
$ h* F9 k  z! a for(i=1;i&lt;=num_x;i++) {
, h- w7 ]2 h  A( `  k_f=0;
  if(sigma&gt;0) {
   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
& B8 m1 f4 m2 `   if(k_f!=1) {status=0; return;}
  }
5 e% T1 r2 b! w! c8 q1 V" T0 _- Y }</P>
$ h8 a& [5 L; H  d; z<>//确定换入变量
" k2 h: ?+ X4 v: i; K% P6 N5 k for(i=1;i&lt;=num_x;i++) {
( ~# ?2 x, N: G# o, s1 v  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
' l) Z2 d6 m' V$ S+ Q' ]9 C  E$ s  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
}//temp赋初值
for(i=1;i&lt;=num_x;i++) {
  k=0;
4 C( g: W: i, V: l" H% q/ Q  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
; l/ p9 D7 {7 H7 g/ {  if(k==0)
   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
8 k1 u4 p8 _% u' x1 x    base_in=i;
- Q# q3 A9 f- S, O" r( Z    temp=sigma;
   }
}</P>
9 D' b; Y1 ]( V; a/ T8 k<>//确定换出变量
' B% z4 j% ^6 `9 o0 R) [+ I# R5 ~ for(i=1;i&lt;=num_st;i++)
! f- `1 Y/ \4 q4 O- T  if(a[base_in]&gt;0) {
   temp=b/a[base_in]+1;
   break;
" R) Y$ A$ n+ h6 t  }//temp赋初值
. E& s" g6 A7 ]! F* k, X/ S 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++)
5 I; i5 j# w4 g) O) z: i% S    if(base==arti[j]) {
  T7 u0 l6 k& S     base_out=base;
2 f& c! f$ |# r! C" c     base_elem=i;
  T5 ?: o% Z. c/ ^" ?( l& o. o     temp=b/a[base_in];
/ ^+ Q5 A- s/ ^% J# M     break;
    }
- Z4 _- L) P3 \! o8 Q2 H: }  }//人工变量优先换出
. \  J; P* l: x" _  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
- S6 b2 b8 i* I' l, e+ ~# V7 q/ H$ H   base_out=base;
$ f. Q- \4 t% i& J3 A+ n   base_elem=i;
   temp=b/a[base_in];
7 E( U* a: d) u  }
- w+ B2 A$ {+ Z+ k; x/ Y4 w }</P>
<> printf(" 基变量：");
; x3 Y2 \- b$ H- m7 q for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
. p2 x# u! d; X2 w- [9 j printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
! ~$ _1 `4 y! j6 p& |$ q//基变量变换，进行新方程初始化后迭代
4 G) ?, l6 X( W for(i=1;i&lt;=num_st;i++) {
9 N# D5 D, w, P5 h  if(base==base_out) base=base_in;
}
+ l: T" D9 x* E: X; W//初始化主元素行系数
* i) c; U( R: ^ value_be=a[base_elem][base_in];
b[base_elem]/=value_be;
; u6 X5 g0 C$ \/ F) E1 H2 S for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
  A# U; p6 j, I<> for(i=1;i&lt;=num_st;i++) {
' n/ Q( q* ]& v  if(i!=base_elem) {
   b-=b[base_elem]*a[base_in];
   value_be=a[base_in];
   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  }
+ {( c0 N2 D/ F, {3 q" f( @ }
( p# H. K+ q7 @# @ stop++;
% h, \5 o5 A3 z  x/ r/ B. z if(stop&gt;STP) {status=-2; return;}
9 l' T7 O) Z0 h& b! R) C iterative();
}</P>
<>void output() {
. C+ E" U; d4 @" E  g2 J int i,j;
- G0 W2 N$ V, n& x* U- n double X[MAX];
  l' m9 }. Q3 z3 U4 i8 A: V printf("\n结果如下：\n");
4 P# j3 K4 W4 i; o; T" c printf("\nX=(");
for(i=1;i&lt;=num_x;i++) {
  for(j=1;j&lt;=num_st;j++)
, k% E, P6 W1 g) _  if(i==base[j]) {X=b[j];break;}
  else X=0;
4 v2 ^7 q! k4 x2 L2 U% ?0 s  printf("%lg ",X);
}
printf(")");
7 g1 M: |. ^5 n. o) p for(i=1;i&lt;=num_x;i++) max+=c*X;
' T0 P( M4 Z' T2 b9 F, d if(ma_mi==1) printf("\nMax z= %lf\n",max);
else printf("\nMin z= %lf\n",-max);
. k; p& x3 c7 p/ E. e  O+ W1 G}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了