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

lauqt1980612

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

winzipftp

是你吗？thunder？帅呆了

tw1982

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

tw1982

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

yqm10507

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

sg47

好详细！

plgatc

<>/*************************************************************************
                    单纯型法解线性规划问题（两阶段法）        
4 v; ~  U' c& }- w2 k; n
6 z3 Q5 r; J1 ~% y/ V7 e; Z 编程环境：VC++6.0     
方程组输入说明：
/ W" o; V: ~. [7 F+ w3 U4 Y 变量非负，按提示输入相关参数。
+ ^% |9 P8 N( {: q4 m *************************************************************************/
#include &lt;stdio.h&gt;
#include &lt;stdlib.h&gt;
#define MAX 100
#define STP 100</P>
+ F1 E. W, P0 H* s3 ?! v+ @; U<>int stop=1; //迭代记数变量
int status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
int step=1; //目前阶段</P>
  o/ y& S. v3 B3 r! _* N5 T<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
9 d) R- \2 ?( X2 X+ Z. P6 Q' xint num_x; //变量个数
. }! V+ s. D8 p+ b# q+ F2 p: `! Gint num_st; //约束方程数
; P. i. c8 H5 U1 x8 a( e: e! _- Bint num_ar=0; //人工变量个数
int arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
<>void create(); //建立方程组
void iterative(); //单纯型法迭代
, s6 ^3 u3 i' e3 k9 Lvoid output(); //输出结果
  Q* K* K- T( o% @: h' w0 \void banner(); //打印程序标题
: I7 C5 X* m2 B/ ivoid exchange(); //交换两阶段价值系数
( A/ {' @  v# Tvoid show(); //输出方程组</P>
<>void main() {
int i,j,k;
banner();
create();
% g0 S# [, R1 h  V //保存原价值系数，转换为第一阶段价值系数
# }$ ~3 X7 G& {  G2 E6 V8 m& V for(i=1;i&lt;=num_x;i++) {
  k=0;
* \& H) H( K: U, Q4 }  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
$ A5 w1 D+ D% I7 W& f  if(k==1) temp_c=-1;
, Q7 z6 o! s/ X" o; l4 C8 @  else temp_c=0;
2 `* X8 V2 ~& L1 \+ |" ?1 t }
exchange(c,temp_c);</P>
<> printf("\n\n第一阶段问题为：\n\n");
show();
step++;
' r2 j0 X, a! J6 { printf("\n\n按回车开始第一阶段迭代");
getchar();
getchar();
3 q1 \1 t) i% l9 d) Z  t( u, f3 Q iterative();
if(status==-2) {
  puts("迭代超过限制次数强行终止！\n");
- `8 H; J3 M) `- D  puts("\n按回车结束");
  getchar();
  ?! G& P  q' s  exit(0);
}
. {9 P7 O3 v. O" A output();</P>
+ e, w1 ^& ]( i- N; V7 W" j3 A( P<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
4 ]# ?. n  q, `, G4 X: L+ k  puts("\n按回车结束");
  getchar();
  exit(0);
# @. T$ u- `2 ^2 l }</P>
5 D( o0 ~1 O" `. v<> //转换为第二阶段价值系数
exchange(c,temp_c);
//把人工变量列全设为0
/ x# H0 n4 _5 ~) J6 D; P for(i=1;i&lt;=num_ar;i++) {
2 W3 k; F4 n# A$ x/ z2 H! W  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
5 u# C! j  A1 H: T# H- @) H6 O }</P>
# n( j  t' J: x" V& h% G( }<> puts("\n\n第二阶段问题为：\n\n");
$ _. M, V/ k# ]7 i show();
$ m9 e" Y. y' f* p puts("\n\n按回车开始第二阶段迭代");
  @. b1 J8 e! v' B8 x  n/ f getchar();
iterative();
9 E9 X) E& g  x: Z" S& d switch(status) {
" v/ p& D; g% ?' t+ q# D case 1:
  output();
  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
case 0:
  puts("\n\n原问题为无界解。\n");
  puts("\n按回车结束");
: w: N+ P/ {1 `; Q8 ]. g  getchar();
  exit(0);
case -1:
* m# g; q$ {4 D0 H1 @: n9 I0 U  output();
# ~3 ~8 b2 [0 {0 Q' _3 L! n# W  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
7 r7 k9 Y! r9 k! z  getchar();
  exit(0);
" y. C' Q# a5 ?: w" T case -2:
  puts("迭代超过限制次数强行终止！\n");
# m* W( f& I  @' Y" S  puts("\n按回车结束");
% x7 e3 B. w* D/ c) G4 L  getchar();
  exit(0);
}//switch
+ |8 u9 i% J# S0 s$ ?- B3 {4 I  j}</P>
<>void banner() {
printf("\t\t****************************************\n");
% l' C6 B2 L7 D: g( b6 ` printf("\t\t         单纯型法解线性规划问题\n");
% Z& P- h# y6 x+ c0 F% _7 @. s" f printf("\t\t                         作者：Thunder\n");
% l& h( E& A9 a1 m+ s- L& @* l# q printf("\t\t****************************************\n");
. |" _, p3 I( j- _* [ printf("\n");
}</P>
, W" e2 v6 @8 b2 l7 \' ]! O8 ?<>void show() {
* t: l3 `. ~% w//对方程组以自然的格式输出，系数为零的x不显示
0 }* j; c8 X+ q1 D//为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);
- z4 v3 R: ?. \. a  break;
' [' Y7 w- X0 I case 2:
. E+ w# l0 C# o printf("max z= ");
' I1 ?7 [3 D" Z* _- A: @) r printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
  B* y* o2 I7 @+ C& f  if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
5 b- O3 Q$ M/ S# ^1 `! b8 r# l  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
. Y: G$ H$ S7 O  else printf(" %lg x[%d]",c,i);
. {" g8 G1 q9 ]) X2 \5 X2 T5 } }
" O  O" U) r9 r! |: j break;
+ G$ L1 a* ~! C3 t8 U" A }</P>
3 ^1 Q  k- q0 s. d<> printf("\nst:\n");
1 f+ H9 R" i3 Y/ p for(i=1;i&lt;=num_st;i++) {
/ o6 \% @5 m2 n  k=0;
  ^/ F, \& K0 ~" M4 N  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
: M. f( i9 v/ F' V# y! M    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);
5 u, [; ]0 y. @. L  ~    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);
    k=1;
* g* c2 F+ a. ?   }
  }
7 v# E( t; t3 ^. v$ H  printf(" == %lg\n",b);
}
/ ?( a% B+ U$ S printf(" x[1]~x[%d]&gt;=0",num_x);
4 j/ B" A& J0 L8 v}</P>
<>void exchange() {
% t) L1 L) @: c' ?% u& ]6 r int i;
. P4 G4 K2 l+ M double temp[MAX];
for(i=1;i&lt;=num_x;i++) {
  temp=temp_c;
  temp_c=c;
3 _. \5 K, R4 o) w  S  c=temp;
}
( n# q" t5 Q- x}</P>
<>void create() {
$ t: G/ Z+ M8 u& G$ e//输入方程组系数，每个方程输完后回显确认
int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
char confirm;
% G* |1 V% {$ [- z4 n8 }, z
$ K7 ?0 G7 w1 }8 C$ {, z7 X; @ while(1) {
* _# s4 f7 o/ c" d, u0 S  printf("请选择：1、求最大值，2、求最小值：(1/2)");
. E4 N9 T) ?1 t: |: _  A' A  scanf("%d",&amp;ma_mi);
0 Q+ ]" ]  E6 }' z  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
4 ]! g5 G  }; V9 a& H8 {  else break;
2 X1 ]1 E, O5 |8 k) F }
( I, }( p9 w/ U9 y  R
. f2 w  K9 c7 X$ L  }( A while(1) {
  printf("指定变量个数：");
  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);
  }
. g! t  {: D! ]8 G5 n) S  if(ma_mi==1) printf("max z= ");
  else printf("min z= ");
  printf("%lg x[%d]",c[1],1);
+ Z  A( a7 J  P  for(i=2;i&lt;=num_x;i++) {
   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
. H, X% O: V7 g' {6 H   else printf("  %lg x[%d]",c,i);
' |/ E2 a- R; Y! ~2 {  }
printf("\n正确吗？：(y/n)");
' G2 Y$ V8 p& q getchar();
# [* e* }' {$ V, S confirm=getchar();
if (confirm=='y') break;
, E% x5 x% X9 n0 j- U1 v else if(confirm=='n') continue;
5 O8 u. {( u( L' q7 ? }</P>
<> printf("输入约束方程组个数：");
8 c# X" M+ |2 M2 |2 S' l scanf("%d",&amp;num_st);
for(i=1;i&lt;=num_st;i++) {
1 F; r1 B/ ^# {2 y/ A; M7 H+ T  printf("st.%d:\n",i);
( n4 o5 E& a& k% Y" \  while(1) {
   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
$ z+ Y! R8 A1 f( M: Y0 v6 S   scanf("%d",&amp;re_st);
   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
2 F7 J$ p& z$ D0 {# m# H   else break;
  }
  printf("输入技术系数：\n");
  for(j=1;j&lt;=num_x;j++) {  
+ O% U8 G+ Y+ v. w2 J" P  o   printf("a%d=",j);
) Z+ N) X3 W" `- ~0 ]7 _   scanf("%lf",&amp;a[j]);
  }
  printf("输入资源拥有量：\nb%d=",i);
  scanf("%lf",&amp;b);
# I7 K9 F: t1 `. V  
6 W* J' V# p$ }3 z  printf("st.%i:\n",i);
  printf("%lg x[%d]",a[1],1);
  for(j=2;j&lt;=num_x;j++) {
   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
- E1 q; t/ w6 Y4 |+ n   else printf(" %lg x[%d]",a[j],j);
: w, P) e1 o0 O  }
  switch(re_st) {
   case 1: printf(" == %lg",b); break;
   case 2: printf(" &gt;= %lg",b); break;
   case 3: printf(" &lt;= %lg",b); break;
  }</P>
1 o0 ~, K- r3 _" _<>  while(1) {
   printf("\n正确吗？(y/n)");
( c1 v* g7 g9 O  n' c' J+ x$ J& V   getchar();
   confirm=getchar();
   if (confirm=='y') break;
$ k5 r. O" ?3 l+ k. b/ D& E   else if(confirm=='n') {i-=1; break;}
+ g; K/ U' g0 d: q  }
3 \* r- y) ^) R# s, r, V$ F }</P>
3 v& L, k$ c- a) t<>//显示输入的方程组
" T! Z# x9 a& \* w printf("\n原问题为：\n\n");
; U# \  a* `" P  N& J: a if(ma_mi==1) printf("max z= ");
else printf("min z= ");
* G: B2 ^2 I0 C" v0 ~ 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);
, M0 r* N7 P# D+ C$ c  M+ k/ y  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
# v( A( k3 i1 Z/ D  else printf(" %lg x[%d]",c,i);
}</P>
<> printf("\nst:\n");
3 w$ ^! i; E3 N6 f- U5 ^5 W for(i=1;i&lt;=num_st;i++) {
8 n8 l- h8 C+ v3 c  k=0;
! P4 U, L# m0 k% Q0 Y1 C  for(j=1;j&lt;=num_x;j++) {
   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);
+ X6 b9 x  U4 M; O    else if(a[j]==-1) printf(" - x[%d]",j);
4 f! R3 Q/ n/ W1 K1 M: e9 F    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
" a# Y& U& |8 N    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
& f9 e: y, i, a2 K1 I- s    else printf(" %lg x[%d]",a[j],j);
    k=1;
/ q- D1 a! g: c; p   }
  }
7 o- ?0 b3 y, f& k& V+ j/ Y/ t  switch(re_st) {
2 F* Q' {: z" b% n* f& y' l! o   case 1:
    printf(" == %lg\n",b);
2 N9 z0 w6 M2 `& y* t+ D    break;
   case 2:
( _) {! N' ?* d    printf(" &gt;= %lg\n",b);
) v1 ~- `; G: w: K/ k) K1 r    break;
% a% M3 Y4 Q6 |2 P3 P   case 3:
    printf(" &lt;= %lg\n",b);
, O7 B3 U4 x0 L6 i7 n2 ~  X    break;
; k( l$ _/ n6 P% W/ V  }
}
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;
( x/ ^* z& Y' G" _: F   break;
  case 2:
   num_x+=2;
   break;
3 h" O5 @! E7 F  }
}</P>
<>//化为标准形式
if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
for(i=1;i&lt;=num_st;i++) {
4 v9 ^. {! B& l. ^; \% P9 d  switch(re_st) {
9 j4 r! q: k5 j9 g8 a1 }0 T   case 1:
    num_addv++;
) G( `, O+ O1 L7 n* J& D; y- Y. O+ @    num_ba++;
    num_ar++;
    c[tnum_x+num_addv]=0;
) M9 S8 B4 O$ V* M9 V& q' f$ V+ S    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;
/ s7 Q  G8 H% e/ P* E$ F% `     else a[j]=0;
    break;
) R; }) v: Q7 Q: t1 v  n: b   case 2:
1 o$ `/ ]+ S' x+ b( X. e! u, e) C% a    num_addv++;
8 l' n9 }$ a% `1 ^/ j5 E% h3 E    c[tnum_x+num_addv]=0;
    num_addv++;
6 d. S: S2 Q' s5 r% @4 z# O% }  V    num_ba++;
! Z* ?. H* k/ t2 e! o1 y0 ?    num_ar++;
    c[tnum_x+num_addv]=0;
5 s! H% ?% D! n, e. s3 s* s6 e    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
    for(j=tnum_x+1;j&lt;=num_x;j++)
! b' S! d6 ]) ^3 U5 k     if(j==tnum_x+num_addv-1) a[tnum_x+num_addv-1]=-1;
5 j  a8 Z; y  _     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
% ~! x5 H9 N! g, J    break;
6 K# D# ]7 _! J6 x. O: g" \/ d   case 3:
    num_addv++;
& O" |& C* t3 X) z2 A( s. w/ [$ l    num_ba++;
% n, c" W' K* E) g    c[tnum_x+num_addv]=0;
: D4 v; {9 z1 Q. F  J% v& _    base[num_ba]=tnum_x+num_addv;
    for(j=tnum_x+1;j&lt;=num_x;j++)
, g% D9 u! y( D5 v  r% t     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
9 G6 J8 ]) g- B" }3 @     else a[j]=0;
4 k8 z7 e1 I6 t* V    break;
' t: D  k# d. [& P- n  }//switch
}//增加松弛变量、剩余变量、人工变量、确定基变量</P>
% F5 M, k, o5 `5 d% A<>//显示标准化后的方程组
" v$ A1 P% M7 j% E printf("\n化为标准形式后：\n\n");
if(ma_mi==1) printf("max z= ");
# V/ O) ]* Z/ {" R0 ?- f else printf("max z'= ");
5 P  w2 B1 T  Z4 q/ d printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_ar;j++)
  if(i==arti[j]) k=1;
& C8 y$ w- N* p, A  if(k==1) printf(" -M x[%d]",i);
. J: m; @; Q) B  R7 `; M8 K- ]  I$ j  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);
* J, E) K/ a( y" w  else printf(" %lg x[%d]",c,i);
}</P>
<> printf("\nst:\n");
/ N  l) C1 ?) b" I2 w& r* y for(i=1;i&lt;=num_st;i++) {
  k=0;
7 _0 Y/ p# [9 W  for(j=1;j&lt;=num_x;j++) {
   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);
    else if(a[j]==-1) printf(" - x[%d]",j);
& ^3 R! R6 n$ a8 J2 [& X  q, ~    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;
   }
+ l7 \+ C% O: U: L, F7 t  }
5 _1 \& C  Z4 Y  printf(" == %lg\n",b);
}
2 W( m( |0 ]/ a6 n printf(" x[1]~x[%d]&gt;=0",num_x);
# ~( i- X9 `; d8 m3 B8 v}</P>
<>void iterative() {
/ k- z: o# O1 y4 r' J3 l int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
0 B9 u2 H9 {4 s int base_elem;
int base_out,base_in;
double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
5 J+ Q! d! J5 Z6 c0 N% J<> printf("\n\n第%d次迭代：\n\n",stop);
& g4 p8 K) Y5 I! j' k for(i=1;i&lt;=num_st;i++) {
  printf("c%d=%lg\t",base,c[base]);
  printf("b%d=%lg\t",i,b);</P>
6 D$ `* t. [! P* K, _5 \- x<>  switch(step) {
; ~" m! J: L9 a) K0 F   case 1:
    for(j=1;j&lt;=num_x;j++)
    {
     printf("a[%d][%d]=%lg\t",i,j,a[j]);
, ]# s  C/ h# Q2 o; S    }
' e8 p$ I0 n) c; {    printf("\n");
. {" m2 W0 ?  O4 j5 `: D7 j& H- w    break;
   case 2:
    for(j=1;j&lt;=num_x;j++) {
6 D% d" ]5 x4 I) o     k_a=0;
& G+ M7 U; d2 V2 V     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]);
    }
& I+ i2 D7 j5 I0 d! q! Q% H    printf("\n");
    break;
* o( a1 ?9 M6 n6 |1 o2 i1 ^  }
9 e. W" r: v  \ }
//求检验数sigma
for(i=1;i&lt;=num_x;i++) {
2 Q5 O  O* ?; m6 p  sigma=c;
2 O" Y: i; i+ r/ _& B2 L  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;
1 J5 p, y( ~) i# H% f* _  switch(step) {
   case 1:
    printf("sigma[%d]=%lg\t",i,sigma);
1 s6 A9 v" D0 u) u1 V1 c    break;
) j$ y( E. \! k0 X- ~   case 2:
    k_a=0;
5 L: z+ X% H3 `: A    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
/ B7 l* e4 q9 q    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
    break;
2 z: ?2 Y2 O$ U- r. S4 k4 j1 U  }
}
putchar('\n');
//检验检验数sigma是否全小于等于0
' K, }; M3 S  x$ C k=0;
for(i=1;i&lt;=num_x;i++) {
  if(sigma&gt;0)
. C0 H9 S# v9 O. ~  Z3 i6 v  k=1;
7 b5 D0 j, l. v8 V! k }
2 Q* I* G$ B% L if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
  for(i=1;i&lt;=num_x;i++) {
$ Q, ~4 ]" \9 g3 i8 v   k_f=k_a=0;
   for(j=1;j&lt;=num_ar;j++)
   if(i==arti[j]) k_a=1;
8 K* K" l. g6 j  X: V5 X   if(sigma==0&amp;&amp;k_a!=1) {
0 x+ |1 J( q5 Z+ g4 V    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
. R6 b5 k& r! m- v. @1 r   if(k_f==0) {status=-1; return;}
   }
  }
  status=1;
+ ]5 J7 n/ B0 J* O  return;
}
8 J  X1 O9 J" n2 ?7 p. J( B& v, E//检查是否为无界解
for(i=1;i&lt;=num_x;i++) {
  k_f=0;
$ h# Y+ V' b5 S# y# E3 p7 D  if(sigma&gt;0) {
5 b  S  P; Z2 \   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
   if(k_f!=1) {status=0; return;}
  }
9 H, V) }# H( s7 L }</P>
<>//确定换入变量
7 e& r4 m: }2 h5 u1 |3 j2 x6 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;
2 c; N: f4 j/ |4 [3 m' p' ]: T  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
+ y6 v+ g6 F6 ]! J1 M }//temp赋初值
for(i=1;i&lt;=num_x;i++) {
  k=0;
. a# A3 I* u( J1 o% q7 [' R0 a  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
3 E5 M# M; o/ F2 J  if(k==0)
+ X2 I) A& i0 D0 }   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
    base_in=i;
    temp=sigma;
   }
; P1 `3 t7 U5 i& I3 w$ c }</P>
<>//确定换出变量
! E. J1 u  N5 e; K$ Y" i6 r for(i=1;i&lt;=num_st;i++)
! Z) }1 J! Z. T) Q8 O9 O* z; @$ `  if(a[base_in]&gt;0) {
   temp=b/a[base_in]+1;
   break;
  }//temp赋初值
for(i=1;i&lt;=num_st;i++) {
/ D7 {" k0 V8 {/ B" l) t  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
4 g. }) B; U8 n   for(j=1;j&lt;=num_ar;j++)
    if(base==arti[j]) {
- q2 H, W/ z  ~     base_out=base;
     base_elem=i;
     temp=b/a[base_in];
     break;
    }
$ v3 |/ G; ~* |  }//人工变量优先换出
  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];
' r) d: O; j8 U! \  y; \  }
}</P>
<> printf(" 基变量：");
for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
- g$ g  I, Z  R6 e printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
//基变量变换，进行新方程初始化后迭代
8 r& B  b1 F1 [5 I, y. X) a- j for(i=1;i&lt;=num_st;i++) {
4 x3 p5 [# Y( {4 G1 I# j$ x  if(base==base_out) base=base_in;
8 y# [/ \5 n- f- _& x; O* F }
2 `+ E4 @: a- q" ?2 n  N) \//初始化主元素行系数
value_be=a[base_elem][base_in];
b[base_elem]/=value_be;
9 z1 W; C7 S. q  `7 U7 p 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];
   value_be=a[base_in];
4 y% `7 W; O! C   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
: ^. S. P: {! X( I3 B5 H0 \/ y  }
( G! o8 b. }# r. W) a1 I0 u2 L5 g3 n }
; L. O# R% Z/ s+ i- c$ e stop++;
  J% O7 E: T/ h- ? if(stop&gt;STP) {status=-2; return;}
# E% d% n/ P" `0 p& |! q iterative();
2 i5 `2 F4 ?+ F  p$ |8 r}</P>
<>void output() {
int i,j;
; Y2 _# w' H/ c double X[MAX];
printf("\n结果如下：\n");
0 F/ ^- g4 v: W" p3 @9 }( p printf("\nX=(");
$ R% ]  d( Q4 ^. ^( L7 Z for(i=1;i&lt;=num_x;i++) {
  for(j=1;j&lt;=num_st;j++)
$ D6 w  a! T& J5 Q! e; c  if(i==base[j]) {X=b[j];break;}
- [5 Q* X- S/ a+ U% f( J  else X=0;
; |1 N4 R, w6 K) O; S  printf("%lg ",X);
}
printf(")");
for(i=1;i&lt;=num_x;i++) max+=c*X;
if(ma_mi==1) printf("\nMax z= %lf\n",max);
else printf("\nMin z= %lf\n",-max);
}</P>

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
, L) Z( L  B8 g4 [4 v8 F#include&lt;math.h&gt;
#define X 5
#define Y 7
1 m( M& r  v) u  s- n9 bvoid xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
{
  int j;
  *c=0;
/ L* G! z% v) U- _+ H/ r& e  for(j=1;j&lt;=*is;j++)
    {
      if((*a)[*ir][j]-*c&gt;0)
    {
      *c=(*a)[*ir][j];
      *j0=j;
    }
& m( S; [7 b# j* P5 V9 _+ {+ b( P2 Y     }
5 ~+ d: N* Q, E# L* q6 @}</P><>      /***************** 参数说明 **********************/
) W+ x: Y7 G7 \      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
. F* N, m; M- F      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
1 Z+ a+ ~" p9 t1 O      /********* <a href="http://happyyangxu.home.sunbo.net*******/" target="_blank" >http://happyyangxu.home.sunbo.net*******/</A>
9 H0 T1 b9 }6 F) N6 \, j6 Fint 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;
  float c,g;
. h! ~$ {8 U+ v% g! c7 Q0 F; \# v  m1=m+1;
  mn=m+n;
  for(i=1;i&lt;=m;i++)
    (*k)=n+i;
8 o0 @% D5 {2 L5 O& w3 Y1 oleap1:
- u8 j, H5 t  J, W  if(l1-1==0)
9 Z$ I9 I; h) G" v( ?4 r8 D9 g( s+ G     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
0 ~, R# `/ B! Z: ~' ^- R  else
4 ~& ?2 m, g( F( ~! O     {
$ y% {7 `- a9 k: W. v$ wleap2:   if(l1-50==0)
" j) \6 d# s' M$ n4 c" F! I      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
7 g; c7 c( m( D, [, A       else
      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
     }
   c-=1e-8;
   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
  I- {$ K, c+ Z& j      {
    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
- j0 j/ V: e$ o4 G% L5 M" F    return(-1);
      }
: j0 C1 e7 O3 g" C/ w7 }) \1 Q9 [   if((c&lt;0)&amp;&amp;(l1==1))
       {
4 J2 g9 g7 ^+ c+ z2 \$ R     l1=50;
, C- S8 z' {) D* p# m     goto leap2;
& u7 J* q  S2 J7 _) `2 W) V       }
    if(c&lt;=0)
/ Z- v6 z+ I6 [0 e' Y" R3 }  E' u2 m       {
      for(i=1;i&lt;=m;i++)
7 t% Q$ H7 |  h+ o/ F         (*x)=(*a)[mn1];
4 K) a# x' k( p' Y5 V+ G* ]# y         printf("\n\t********Optimal solution********\n");
      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
: F3 J# H, B8 }! n+ b+ w4 b1 c/ m      printf("\nF=%f\n",(*a)[m2][mn1]);
      return 0;
    }
, ?( D+ B, ~5 p& z( Y$ R3 O     c=1e8;
% _" R. h/ x& [; |     for(i=1;i&lt;=m;i++)
       {
/ Z1 V9 J) Q9 W) P& @1 j     if((*a)[j0]&gt;1e-8)
/ a$ ]( m; }4 ^& f: G        {
# F7 H; q  Z$ J* N7 G          g=(*a)[mn1]/(*a)[j0];
          if(g-c&lt;0);
         {
: b0 }0 @3 e7 l8 P, b0 N. _/ n- f           c=g;
& ^6 {6 U1 b2 f- A           i0=i;
          }
( I* [. x9 W# ~% O+ ~         }
6 o% ^) p3 v) M  [! o$ ^    }
      if(c==1e+8)
    {
      printf("\n\t*******Lp no solution********\n");
      return -3;
' B6 M: l& d5 i" G     }
       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
! d3 A; E1 m+ u9 G3 r      {
        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
           continue;
% [  d! P, S! L8 f         g=(*a)[i0][j]/(*a)[i0][j0];
           (*a)[i0][j]=g;
         for(i=1;i&lt;=m2;i++)
3 S6 n7 M; u; A' v6 V( `; B2 M           {
* l( R0 C1 C! J# {6 g         if((i==i0)||(!(l1==1)&amp;(i==m1)))
1 ~+ k: E% t, R4 W            continue;
         (*a)[j]=(*a)[j]-(*a)[j0]*g;
           }
2 z: t( y4 E& a       }
       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
       (*a)[i0][j0]=1;
       goto leap1;
3 ]1 y+ u) j+ Y" _2 |2 }& Z/ W}  
& D5 C( T4 H# _; F; H0 g8 wmain()
{
  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},
                 {0,0,0, 0,  0,0,0},
                 {0,69,0,144,0,0,300}};
. W1 s$ I9 e' F5 y8 \  int k[3];
  c- @! p/ j- I0 \& _" Q- @  float x[3];
  clrscr();
  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
}
- S+ G0 A. C! I</P>

julysea

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

lauqt1980612

最好是C语言的!