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

lauqt1980612

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

winzipftp

是你吗？thunder？帅呆了

tw1982

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

tw1982

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

yqm10507

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

sg47

好详细！

plgatc

<>/*************************************************************************
                    单纯型法解线性规划问题（两阶段法）        
, n6 r$ h: X& A0 v; H) R
9 l9 H, _3 G; P 编程环境：VC++6.0     
方程组输入说明：
! o4 N) a8 \  C 变量非负，按提示输入相关参数。
*************************************************************************/
7 Q/ \+ Y* G+ I; y$ i6 |#include &lt;stdio.h&gt;
#include &lt;stdlib.h&gt;
0 ]6 e/ _' ~8 y8 _3 u2 V#define MAX 100
#define STP 100</P>
<>int stop=1; //迭代记数变量
int status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
int step=1; //目前阶段</P>
/ ~: |( J7 h9 Q- i; I# b3 x<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
int num_x; //变量个数
3 R1 f( r7 n$ fint num_st; //约束方程数
+ z, W1 F2 o& k" eint num_ar=0; //人工变量个数
int arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
8 H( C3 y' r+ B" i0 H: Aint ma_mi; //1为求最大值，2为求最小值</P>
$ [# \2 j3 S; B6 D<>void create(); //建立方程组
- Y% U; M: E# W  M, ~# Kvoid iterative(); //单纯型法迭代
void output(); //输出结果
void banner(); //打印程序标题
6 ^6 `( a6 @* |3 avoid exchange(); //交换两阶段价值系数
void show(); //输出方程组</P>
<>void main() {
( Z9 N, M2 m; y0 A0 w int i,j,k;
banner();
create();
7 T  P. R1 R- A" h' X) W( t! M //保存原价值系数，转换为第一阶段价值系数
for(i=1;i&lt;=num_x;i++) {
  k=0;
- j: e: H: L6 B& P& o1 M7 x; i  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
  if(k==1) temp_c=-1;
: }0 j. p1 T5 ?& c  else temp_c=0;
: Y, B6 U8 N. z# M }
+ Q; A2 o' v2 F0 |. D; q7 H: X exchange(c,temp_c);</P>
# ^& }" t" T, |<> printf("\n\n第一阶段问题为：\n\n");
7 o4 T8 j7 x# x8 L0 ~6 H0 X6 i. s show();
# p3 T5 i1 L$ K2 F7 P/ x step++;
printf("\n\n按回车开始第一阶段迭代");
getchar();
getchar();
$ f1 f1 h& u# h" ^' C- q4 [* ] iterative();
7 S$ d/ ], I. h; {8 m9 G) D% A if(status==-2) {
  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
( X* f( C/ B! k  getchar();
  exit(0);
) {0 T) L  \( }: W }
output();</P>
& }- K0 ^/ v) K7 s" b) }3 k: I<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
9 @  R0 `+ v6 I2 U7 }" i% O }</P>
* _& Q* {3 J4 Y, y+ i6 L0 \<> //转换为第二阶段价值系数
exchange(c,temp_c);
//把人工变量列全设为0
for(i=1;i&lt;=num_ar;i++) {
+ U3 X4 }' d( Q9 S) K% y4 ?  c[arti]=0;
! J% B; ]  U6 O' s  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
, M8 V1 y9 g, A  W. }- d" j }</P>
<> puts("\n\n第二阶段问题为：\n\n");
show();
) w5 D' L, w  O puts("\n\n按回车开始第二阶段迭代");
getchar();
iterative();
switch(status) {
  v$ A% n! f0 ] case 1:
/ k; K: F1 t/ k8 i$ K8 E9 @; N  output();
  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
% g3 z& k' ^4 h  getchar();
  exit(0);
8 i( j$ l5 t" c$ ~" Z2 V case 0:
1 d; F+ _) L" u6 I4 s$ ?6 U# b9 I. B, o  puts("\n\n原问题为无界解。\n");
  puts("\n按回车结束");
/ d- ?9 C7 y; H. R! T5 C  getchar();
5 U4 A0 h4 _* y% r! m  exit(0);
case -1:
  output();
  puts("\n\n原问题有无穷多最优解。\n");
6 @! v2 ^) A# u' M6 q5 F  puts("\n按回车结束");
  getchar();
4 [9 H4 L, r- P+ \/ V% J  exit(0);
7 }) `) X$ M* l/ i case -2:
  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
  getchar();
6 }  A4 _9 q& U5 K3 ^6 }4 c  exit(0);
$ z! _9 D" Z5 ~4 Y% w( L) X }//switch
- l" M5 q* V6 s6 p! M. ]! c1 J  ~/ _  y}</P>
<>void banner() {
printf("\t\t****************************************\n");
; n# B) C3 p7 [5 c) I4 } printf("\t\t         单纯型法解线性规划问题\n");
0 E3 d% w- J2 n' W printf("\t\t                         作者：Thunder\n");
1 c4 }* S! x, _+ l printf("\t\t****************************************\n");
printf("\n");
. _0 e: `) U) D8 m; \2 Q3 j}</P>
<>void show() {
. K3 ?3 ^, W# A# t; M6 E! Z//对方程组以自然的格式输出，系数为零的x不显示
//为1的不显示系数1，-1系数只显示负号
$ z: e+ W' J7 W8 `( P" I int i,j,k;
switch(step) {
case 1:
9 N: c' t, M# `) @: e  printf("min z= ");
  printf("x[%d]",arti[1]);
& P+ r9 m( }0 e, N  b, g+ ^  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
/ \! V7 R: @) I  break;
case 2:
printf("max z= ");
printf("%lg x[%d]",c[1],1);
/ G1 g9 v: R" H1 s* ^' ? for(i=2;i&lt;=num_x;i++) {
8 l) y% R8 ^4 U( E* p( K5 x  if(c==1) printf(" + x[%d]",i);
4 {" m* G+ a! h& o& {  else if(c==-1) printf(" - x[%d]",i);
+ l( ]7 K6 }+ v  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
- n. U4 z' m. r. J" H0 b2 G# M4 x  else printf(" %lg x[%d]",c,i);
}
break;
# K7 \( x& v8 C" n. K; S }</P>
! w" k, p- S( _% \<> printf("\nst:\n");
1 H' c3 d, a! O4 U) h for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
5 u( K6 h. f! ]* s    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
    else if(a[j]==-1) printf(" - x[%d]",j);
, ?' x2 ?0 J" m3 @    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
$ l6 o" I( {% p    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
! N  I/ Q0 p) g9 I0 T4 ^    else printf(" %lg x[%d]",a[j],j);
$ ~/ [0 m7 S1 u1 U7 X    k=1;
5 u6 t/ y+ [7 X1 s- ?  S" G* ?   }
  }
7 y3 A$ v% B  h5 d- F  printf(" == %lg\n",b);
}
0 |5 H7 o0 n* e# |/ D printf(" x[1]~x[%d]&gt;=0",num_x);
6 j! v; b7 F( b* ?3 g/ f; @}</P>
<>void exchange() {
int i;
double temp[MAX];
for(i=1;i&lt;=num_x;i++) {
0 f& V  K3 N2 ~. X  temp=temp_c;
  temp_c=c;
  c=temp;
  D5 {& ~- D& F }
}</P>
<>void create() {
//输入方程组系数，每个方程输完后回显确认
int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
char confirm;
2 g. Y/ _; D: S1 j9 w
( x2 I& F) R4 j( ~ while(1) {
9 L1 x+ X& F5 G  printf("请选择：1、求最大值，2、求最小值：(1/2)");
: `, c" z% @2 W1 P! O  c  scanf("%d",&amp;ma_mi);
5 T! N( W$ q' a3 d: Z5 S  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
  else break;
}
8 a% U7 T$ v, C# T0 I
while(1) {
  printf("指定变量个数：");
  |; N$ Z  K7 a1 a  scanf("%d",&amp;num_x);
+ K2 `8 j: P$ X1 e  printf("输入价值系数c1-c%d:\n",num_x);
  for(i=1;i&lt;=num_x;i++) {   
   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);
& ?( y, P8 p- T* w   else printf("  %lg x[%d]",c,i);
. A9 x  w3 n& A. B  }
printf("\n正确吗？：(y/n)");
getchar();
0 `+ J% R6 N7 ?) m, e confirm=getchar();
$ W5 m8 R& R/ S4 V7 F- y9 W' L if (confirm=='y') break;
+ g# Z' [6 A: A+ k9 k else if(confirm=='n') continue;
. I( w3 w# T, q* k }</P>
' f+ V2 p* _4 j2 c) F- ^) C<> printf("输入约束方程组个数：");
scanf("%d",&amp;num_st);
) G$ I1 O! W' d; A7 O for(i=1;i&lt;=num_st;i++) {
  printf("st.%d:\n",i);
8 Q3 T' f4 @6 U  while(1) {
- v+ J% G$ t  t% G; c   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
" {, i5 \7 \7 u4 e7 V) _% ?2 T   scanf("%d",&amp;re_st);
   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
   else break;
2 ^! |" Z/ G2 v/ \9 Q  }
  printf("输入技术系数：\n");
  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
$ D7 I- e4 L9 x8 z* |; t# @0 w5 X   scanf("%lf",&amp;a[j]);
  }
" D. q$ B" ]- {- ?* h  printf("输入资源拥有量：\nb%d=",i);
  scanf("%lf",&amp;b);
  
  printf("st.%i:\n",i);
5 z5 Z; C+ X$ V5 M2 {2 P. c: ~- V  printf("%lg x[%d]",a[1],1);
2 B# H% O; K. u- o2 q  for(j=2;j&lt;=num_x;j++) {
( w/ h, N$ l2 k, T   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
   else printf(" %lg x[%d]",a[j],j);
  }
( o" k# L# u. {  switch(re_st) {
   case 1: printf(" == %lg",b); break;
   case 2: printf(" &gt;= %lg",b); break;
   case 3: printf(" &lt;= %lg",b); break;
  }</P>
<>  while(1) {
   printf("\n正确吗？(y/n)");
   getchar();
; K/ _  }' I! P& j5 l( X/ g5 |   confirm=getchar();
5 u6 g; s9 Z4 }1 v   if (confirm=='y') break;
6 D2 N8 J6 @% \" S/ Y0 r5 u. @   else if(confirm=='n') {i-=1; break;}
  }
' C# y1 L5 A  S9 w }</P>
<>//显示输入的方程组
2 x; K/ B% j( n$ i, U5 ?. a printf("\n原问题为：\n\n");
if(ma_mi==1) printf("max z= ");
else printf("min z= ");
! q% u% p+ [7 l. \ printf("%lg x[%d]",c[1],1);
5 A- K6 t+ v" J1 ]; ^ for(i=2;i&lt;=num_x;i++) {
' e9 a& r$ r$ X, b  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);
+ C$ c! ~+ l0 X& K5 }% c( y  else printf(" %lg x[%d]",c,i);
; [6 i- O2 D7 @( T }</P>
, W; B# ]% Z! x6 A/ q, Y3 g6 w<> printf("\nst:\n");
3 M- s0 [2 E+ v. ^5 { for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
5 a& b  Q* O+ `   if(a[j]!=0) {
; B6 W2 I/ R. s, C- q9 |    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);
$ R/ H: N- r. U' A3 Z' X/ l& A% Q    else if(a[j]==-1) printf(" - x[%d]",j);
- G* Y; b+ O2 l8 R$ f    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
" n2 _; G" A+ J: ?& T: }3 Y    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
3 m& {1 J* l  [% l    else printf(" %lg x[%d]",a[j],j);
9 f' i: d1 S8 U2 o- M1 q. y    k=1;
   }
  }
8 I+ A, x) z' o$ W0 s/ Y( }  switch(re_st) {
   case 1:
/ F# ~) _: Q" b9 o6 L    printf(" == %lg\n",b);
* Q# c2 }" Q* Z( Z    break;
   case 2:
% F# k; c; R. g6 y/ M    printf(" &gt;= %lg\n",b);
% g7 w+ D6 O/ j2 N* l# o" I5 C    break;
   case 3:
    printf(" &lt;= %lg\n",b);
    break;
  }
}
printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
<> tnum_x=num_x;
" d+ d$ }/ M) H& C: q% x for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
' }( {- a0 b6 O2 ^: Z& ]; O% N  case 1:
  case 3:
   num_x+=1;
8 f$ o# B& Z! T+ i# [   break;
1 j0 \: Y0 R. Q  case 2:
/ {6 H' J" x% o% ~2 I   num_x+=2;
   break;
7 ?' M9 Q; S& {0 o- C: E  }
}</P>
<>//化为标准形式
# E. z1 K: n( r# l3 A if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
' r8 |% t$ z& `( e0 o7 } for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
5 b! C. R2 ]% u- l   case 1:
. K/ J% i# m- h  Q6 s8 o: G    num_addv++;
& J( d! L4 \# _4 S    num_ba++;
+ x  |) r5 R, K& z    num_ar++;
; T8 h/ }* C$ v8 Q- a    c[tnum_x+num_addv]=0;
0 g  S3 K4 f+ L9 a7 p( K& [, F    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
, n- e+ D% P$ D2 @    for(j=tnum_x+1;j&lt;=num_x;j++)
9 g- o! f. G' H: n3 h3 q     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
    break;
' W1 V1 J0 f4 ~- g2 e   case 2:
    num_addv++;
    c[tnum_x+num_addv]=0;
" S" l: U5 ]8 C) ?: O    num_addv++;
    num_ba++;
* P3 f1 Q: i8 N6 `8 a9 L/ C    num_ar++;
    c[tnum_x+num_addv]=0;
    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
9 ?; U7 \7 e6 e# K5 a9 X. \; d: W    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;
+ r3 d( c5 s2 a4 v7 ?& ?     else a[j]=0;
    break;
   case 3:
- q% b" e3 k  U: ?/ r; [3 ?    num_addv++;
    num_ba++;
    c[tnum_x+num_addv]=0;
    base[num_ba]=tnum_x+num_addv;
0 Q. _! S/ o' q+ v, H    for(j=tnum_x+1;j&lt;=num_x;j++)
( q! r0 d& E/ p     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
  t4 Y: @5 ^( R     else a[j]=0;
) h% ~$ I" U; o! \! y5 A8 V$ W    break;
& E; Y2 A) P6 B  G$ g+ T7 ^  }//switch
}//增加松弛变量、剩余变量、人工变量、确定基变量</P>
1 U3 w" L* t# J<>//显示标准化后的方程组
% [# S( k+ T; l! w2 ^ printf("\n化为标准形式后：\n\n");
" g. R3 u- P5 k4 M5 I7 h if(ma_mi==1) printf("max z= ");
else printf("max z'= ");
  l& v2 {  T, n! t( B" B9 w; W. G 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++)
8 r7 z# |, u( F/ ?, W+ K  if(i==arti[j]) k=1;
  if(k==1) printf(" -M x[%d]",i);
6 R- x9 Q& \6 S0 G5 W  else if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
; u4 k1 g+ j$ p& a% j6 B  Q6 G  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
1 c/ I8 W. T; ~5 I3 T  else printf(" %lg x[%d]",c,i);
1 S4 Y+ u6 A. T& x3 h# e: D }</P>
9 I: Y. I( V! Z  n% F/ l$ G( o& p' {<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
% `! ^% p6 K9 z% N) A1 G   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);
+ T) x6 B: g, l$ K    else if(a[j]==-1) printf(" - x[%d]",j);
+ ]1 z$ k/ ~8 ?0 e/ B    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;
8 E1 Q8 R1 g! `7 d* J# x/ U8 F   }
  }
) C8 e) o: m1 v9 i  printf(" == %lg\n",b);
; y( U8 o/ B3 f4 H: u }
3 k: }& ?. `( S9 r8 p printf(" x[1]~x[%d]&gt;=0",num_x);
- o1 e% ~- s1 I$ o' ]}</P>
3 `! ^( m& _- _! B. P1 A<>void iterative() {
7 Q0 H7 w; g# M1 x$ ]" \8 Q" `# a$ j int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
& i4 A. V2 ^0 Q, c0 V' p int base_elem;
5 Q0 W3 y. ^3 c4 V/ o$ I- i  u int base_out,base_in;
double sigma[MAX],temp;
0 q/ Q8 i. {2 }; T* |: P% K( m double value_be; //高斯消元里保存主元素值</P>
<> printf("\n\n第%d次迭代：\n\n",stop);
+ L$ [8 s6 Q2 c# x3 D% ? for(i=1;i&lt;=num_st;i++) {
  printf("c%d=%lg\t",base,c[base]);
8 U1 c8 P* p( |  printf("b%d=%lg\t",i,b);</P>
<>  switch(step) {
0 b3 I8 g! n+ r# @) T/ \+ b   case 1:
    for(j=1;j&lt;=num_x;j++)
    {
     printf("a[%d][%d]=%lg\t",i,j,a[j]);
! G, L* p: o4 R! w* ?% J    }
    printf("\n");
9 R# P2 M" j: W$ J& H    break;
   case 2:
    for(j=1;j&lt;=num_x;j++) {
     k_a=0;
* d- ^: p5 t1 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]);
/ N6 R$ r- \* B' j% A9 f4 J    }
1 P! ~! k7 ]6 v; \% t# r( `    printf("\n");
    break;
  }
% q( D; w2 n1 A! p6 i% Q& d7 L }
0 w4 f4 V; X" n  U" m//求检验数sigma
for(i=1;i&lt;=num_x;i++) {
  sigma=c;
  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;
  switch(step) {
3 f+ O/ |9 Y: R* h6 }   case 1:
3 C$ Q8 Y, T( \0 @. N    printf("sigma[%d]=%lg\t",i,sigma);
    break;
; r& e9 ]* U. G6 W7 t, G7 n! n$ ~. e: d$ S3 B   case 2:
    k_a=0;
    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
" i" ]' o, T" _/ ~8 n    break;
% }  l) c( V2 @1 I. j7 P  }
4 s7 k7 b# m8 g! C* b. Q }
3 q' q9 D) R8 y& U! q) g( w3 p putchar('\n');
, l; }8 S: n: `+ z//检验检验数sigma是否全小于等于0
% v  A! `. N  l0 i. E5 K k=0;
8 a/ z1 D& T9 q+ z for(i=1;i&lt;=num_x;i++) {
6 E) |' T+ r4 u- L5 ?9 t  if(sigma&gt;0)
  k=1;
}
if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
  for(i=1;i&lt;=num_x;i++) {
   k_f=k_a=0;
( \# k, ]& z6 G3 n. S+ {   for(j=1;j&lt;=num_ar;j++)
) W& }, r  w( N2 m* V7 Y% T   if(i==arti[j]) k_a=1;
   if(sigma==0&amp;&amp;k_a!=1) {
2 ?/ W: E. o) U; b; B    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
* I$ w! y5 g3 \8 w: o   if(k_f==0) {status=-1; return;}
   }
2 W4 ~0 C/ G7 X. ?  }
  status=1;
' @( x; h0 G& o% d  return;
}
4 w, M0 `7 `  v1 _+ R! ?% X: Q//检查是否为无界解
  K. E9 i1 z1 Z6 z+ Y1 U+ R for(i=1;i&lt;=num_x;i++) {
  }6 b; T) W. E2 h  k_f=0;
  if(sigma&gt;0) {
3 N# F- m7 \. }# J; [' x8 O   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
   if(k_f!=1) {status=0; return;}
; v; @+ w; j' S# i  P) |! H6 j  }
}</P>
. O, e( e) a  y* l# `  }<>//确定换入变量
for(i=1;i&lt;=num_x;i++) {
  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赋初值
for(i=1;i&lt;=num_x;i++) {
+ K$ e" x8 K' x1 E6 `; ]  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  if(k==0)
   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
- Y: F1 C& j1 v# {    base_in=i;
    temp=sigma;
   }
}</P>
8 n8 k5 A7 o0 d$ s/ l9 G/ ~<>//确定换出变量
7 p! g  w- T4 l6 Z for(i=1;i&lt;=num_st;i++)
  if(a[base_in]&gt;0) {
. |0 \: N+ q% B5 X. w. t   temp=b/a[base_in]+1;
5 d4 Q! F' P' N2 b4 q  I# V   break;
4 B- ]  h, {0 G  }//temp赋初值
% r% `3 p/ W8 @  @' Q9 H for(i=1;i&lt;=num_st;i++) {
  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
( `. ~3 b# @3 _6 y4 b   for(j=1;j&lt;=num_ar;j++)
! a" Z% d  c' c! A" K    if(base==arti[j]) {
     base_out=base;
     base_elem=i;
6 v, R# f" |* Q1 p) h& k     temp=b/a[base_in];
$ H0 ^3 \) m( A3 A9 l     break;
    }
  }//人工变量优先换出
  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
   base_out=base;
! R! a+ j& |+ k8 _   base_elem=i;
9 \! A0 i  e# C8 b" \% i& }- V   temp=b/a[base_in];
! o& V$ X/ N& t0 M$ F, ]  }
}</P>
<> printf(" 基变量：");
for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
! Z7 J  r: D# ~$ P printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
//基变量变换，进行新方程初始化后迭代
* V6 e: t: ~4 V' }/ X for(i=1;i&lt;=num_st;i++) {
  if(base==base_out) base=base_in;
}
//初始化主元素行系数
value_be=a[base_elem][base_in];
: K9 M: s/ _7 ` b[base_elem]/=value_be;
. ^' [0 Q* ?3 i8 `! s 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];
& |2 U  d' `  l9 x6 H; ^   value_be=a[base_in];
5 [+ ~. q; L% H6 }  t   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  q1 ~9 O3 m6 P# H  }
}
stop++;
+ [  V; b# u  c! t) [* @ if(stop&gt;STP) {status=-2; return;}
! S0 r' k- E- @9 z, M6 U iterative();
; M) a# J( j( `" a2 N$ X}</P>
5 N8 `( E, ~% E7 k+ e<>void output() {
int i,j;
6 c7 `* \0 `( q- q double X[MAX];
9 A5 M+ S9 ^8 N& ]  W3 i printf("\n结果如下：\n");
printf("\nX=(");
for(i=1;i&lt;=num_x;i++) {
% }, s+ i, G- t& S- B  for(j=1;j&lt;=num_st;j++)
8 D8 O3 {' y. p4 z* ~  if(i==base[j]) {X=b[j];break;}
9 v1 |* M) H; i% u! S  else X=0;
3 ^3 J# d, T4 K4 v  printf("%lg ",X);
7 |+ O! J" W' R  T2 t }
printf(")");
for(i=1;i&lt;=num_x;i++) max+=c*X;
$ `! [. q5 [( h$ y; d if(ma_mi==1) printf("\nMax z= %lf\n",max);
else printf("\nMin z= %lf\n",-max);
0 B8 `# }0 N5 v+ Y* O* e) l6 x}</P>

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
: E. x. m, j1 c2 Z- l& @/ L#include&lt;math.h&gt;
" P7 q: K$ y" M9 o6 x3 r3 q' M#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])
{
  int j;
( J, g0 Y% U; z5 N% A  *c=0;
  for(j=1;j&lt;=*is;j++)
; P- c' _5 ~5 J2 Z; N    {
      if((*a)[*ir][j]-*c&gt;0)
* `# B0 {* P) a/ W/ v    {
2 o" I- R( W& e  L7 B/ t+ @8 A      *c=(*a)[*ir][j];
# f& v5 J6 x( c, g( R1 E. ^      *j0=j;
4 ?* V) U4 M9 Z% k# V* t. d    }
     }
}</P><>      /***************** 参数说明 **********************/
      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
% ~! [1 C; q6 i      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
      /********* <a href="http://happyyangxu.home.sunbo.net*******/" target="_blank" >http://happyyangxu.home.sunbo.net*******/</A>
& y% M% y1 x, [1 `4 Oint xi_sm(int m,int n,int m2,int mn1,int l1,float (*a)[X][Y],
8 L: i; b9 G; W: p2 j: m       int(*k)[],float(*x)[])
: a- E# i& ?6 {5 _- r" z{
2 r% T) Y0 S2 W1 S5 a) m  int i,m1,mn,j0,i0,j;
8 M) T* ~0 q/ p- R( P8 J  ^( O' Y7 j  float c,g;
, e- O, w+ A, ]; t  m1=m+1;
) X1 {) A% G1 [# r& m# f+ i* H1 @: r  mn=m+n;
; z( H) E" K4 g- }  for(i=1;i&lt;=m;i++)
7 l" E6 i/ O( [    (*k)=n+i;
leap1:
  if(l1-1==0)
3 I; S; {! U& j! f     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  else
     {
leap2:   if(l1-50==0)
4 s5 [1 x2 P! r3 L      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
" ]) d0 H( \& E8 _& U; I0 |       else
* \# l  u( l- X: K; h      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
" s# t- s' ^1 k' U2 @7 _     }
, g8 x4 e% @, t' ~% _( R. _   c-=1e-8;
   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
      {
" `6 K! X# {1 G# t  t! O2 L    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
    return(-1);
; B  Y# v, [: @      }
   if((c&lt;0)&amp;&amp;(l1==1))
8 k0 T& U- I- H2 r& g; Q) C- j6 R       {
     l1=50;
5 A6 d7 i& K) S4 k- H* [6 t/ h! D     goto leap2;
, o+ H0 i2 }# z4 p5 T( X2 E6 v       }
- M+ _9 }# x& G    if(c&lt;=0)
5 F7 m$ i& q( B% _( l& b       {
      for(i=1;i&lt;=m;i++)
1 t& p) u* k; h$ N9 k" o         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
5 u1 B5 C$ Q$ N8 ^' ]- y: u: r" m      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
      printf("\nF=%f\n",(*a)[m2][mn1]);
2 u- N4 _, p1 F! j8 a- X# C, Q      return 0;
2 B1 d/ p: k* i3 X5 {; K( L2 d    }
     c=1e8;
' b+ P% [% P# J2 y. ?5 J+ r     for(i=1;i&lt;=m;i++)
- q2 z; b7 Z" c! e+ A. g       {
% `9 G+ N; Z# M, J+ A) w     if((*a)[j0]&gt;1e-8)
; t' \4 e9 D6 U$ \6 f        {
: r9 W, R( t' i9 j% a5 g: s          g=(*a)[mn1]/(*a)[j0];
& Z* J$ \( {  E3 U          if(g-c&lt;0);
         {
           c=g;
3 P+ y1 _7 r: O( g0 c+ @6 ^- z           i0=i;
$ {; q" J! M' r2 O  z: I2 j: w          }
         }
5 U* W3 J3 n6 ^; [  u    }
      if(c==1e+8)
) s" o- W8 ?6 Q$ h8 i* q& n4 J    {
      printf("\n\t*******Lp no solution********\n");
      return -3;
5 Q& ?* {) g$ X     }
+ c% A( i* d" M8 ^; ]/ F7 i       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
" \: H# l7 R: U: g  x& P      {
        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
           continue;
         g=(*a)[i0][j]/(*a)[i0][j0];
) O! b5 N) Q9 f; k7 T) h           (*a)[i0][j]=g;
         for(i=1;i&lt;=m2;i++)
- k8 P* O8 ~' K3 p. o           {
         if((i==i0)||(!(l1==1)&amp;(i==m1)))
% f4 D2 u. r7 r+ t2 W# N6 K# A            continue;
" U  n* A% b9 }7 G         (*a)[j]=(*a)[j]-(*a)[j0]*g;
           }
8 J0 n) a8 z4 j3 ?8 G* s# g       }
       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
( t: K: M+ g+ Y% v# T       (*a)[i0][j0]=1;
* \8 Z; B. g; R* `4 }# d2 g+ M       goto leap1;
}  
main()
7 d$ l9 y5 D# }. q- I{
  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},
1 V0 o9 S2 E9 v* V8 U                 {0,0,0, 0,  0,0,0},
                 {0,69,0,144,0,0,300}};
  int k[3];
  float x[3];
% w/ E' g' V' r( r/ W& o  clrscr();
; s# Y. D+ u1 v0 h  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
}
* a. M7 ^) K  W3 [( \# {* z+ z6 [/ j( D</P>

julysea

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

lauqt1980612

最好是C语言的!