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

lauqt1980612

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

winzipftp

是你吗？thunder？帅呆了

tw1982

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

tw1982

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

yqm10507

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

sg47

好详细！

plgatc

<>/*************************************************************************
                    单纯型法解线性规划问题（两阶段法）        

  Q2 [6 q- s; U, l& u, h  B: d 编程环境：VC++6.0     
* e, u5 `2 i" V5 ^" h/ M  L 方程组输入说明：
变量非负，按提示输入相关参数。
7 i) E4 a: a7 i# f *************************************************************************/
. _- m; h5 U4 F2 Q  q: f( C) f#include &lt;stdio.h&gt;
: C8 b2 `" ~4 ?8 E5 z#include &lt;stdlib.h&gt;
4 ~- p. b/ ~3 B" X4 P8 X6 j; V5 w) o#define MAX 100
#define STP 100</P>
( y; {0 }: \: z. o# j<>int stop=1; //迭代记数变量
! e- E% }6 R9 V3 n* E) pint status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
+ E$ w2 C5 a# B; i  @; }int step=1; //目前阶段</P>
# u) c# D  }/ [7 {- G" Y9 ~<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
2 q' Y& ^& O* F) e# J/ @int num_x; //变量个数
8 F2 n' `2 i  [& W: ^& j' Gint num_st; //约束方程数
/ V4 A1 k# k4 ~" [5 I, [- dint num_ar=0; //人工变量个数
% T3 X+ q. q# n$ A# bint arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
, X; K2 _' W2 T% q<>void create(); //建立方程组
void iterative(); //单纯型法迭代
+ y% M6 X5 }3 k$ Wvoid output(); //输出结果
void banner(); //打印程序标题
5 _6 T. g! B4 C0 [. {void exchange(); //交换两阶段价值系数
4 v; z  b0 o. e- Gvoid show(); //输出方程组</P>
. P# v5 G: W8 g<>void main() {
& t* q8 x( A/ r- q. q* u: |# y int i,j,k;
banner();
create();
4 f! d7 I1 ~: X; I. V //保存原价值系数，转换为第一阶段价值系数
# x8 V3 i, Q- z' \# K for(i=1;i&lt;=num_x;i++) {
7 c% H, D# ]3 B% b) g2 A  k=0;
+ \/ T, T' B* Y/ Y6 d  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>
5 x& i2 U2 J0 t3 b8 t! V* `2 [<> printf("\n\n第一阶段问题为：\n\n");
show();
( \( m6 n& G3 b- h7 O* O1 q; K step++;
. a& a9 I) w: D8 c1 t printf("\n\n按回车开始第一阶段迭代");
getchar();
3 I) S6 H- q- u. s" \) R getchar();
iterative();
if(status==-2) {
3 k/ I; @* Q, p! H  n  puts("迭代超过限制次数强行终止！\n");
" S+ _5 X' q3 B' H, o9 |  puts("\n按回车结束");
  getchar();
  exit(0);
! h. \" F: s. u4 M. B+ o- K" f }
output();</P>
<> if(max!=0) {
. y% x* H: H; p! g6 b" _# o  puts("\n\n原问题无可行解。\n");
0 b2 V" a/ H6 s% V( U& q/ C  puts("\n按回车结束");
  getchar();
* ~8 r3 f- z2 R" x  exit(0);
, c0 W. @% X- q1 j }</P>
<> //转换为第二阶段价值系数
exchange(c,temp_c);
  h5 Y8 S9 h  i //把人工变量列全设为0
for(i=1;i&lt;=num_ar;i++) {
) G/ H* y8 ^: T) f5 O$ v  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
% [) V. ?7 O! Z$ F. f1 e' b }</P>
  h2 Y" q  K2 s1 ^: ~<> puts("\n\n第二阶段问题为：\n\n");
show();
puts("\n\n按回车开始第二阶段迭代");
5 S/ ^& u, m% _$ ]3 D0 y2 N* M getchar();
* o9 b$ Z* T7 _) }9 n3 X3 ^4 [. Q iterative();
switch(status) {
case 1:
  output();
) z! {+ l& R  ?; M( s  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
; m  M3 |' v  _3 ?' m# ~3 r  getchar();
  exit(0);
case 0:
$ i7 T) e6 K* O  puts("\n\n原问题为无界解。\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
9 W0 R, w( H, H& h: U case -1:
  output();
  puts("\n\n原问题有无穷多最优解。\n");
8 t5 v& l: T/ t; P$ D. N! B  puts("\n按回车结束");
  getchar();
  exit(0);
case -2:
# ]# o- C9 Y/ Y1 F& `( ~" d# ]1 K  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
& F. \5 J- W# A, v8 ?2 u* e. ]  getchar();
  z. q% N# z3 E& v+ m1 Q: c  exit(0);
" Q- r4 S5 D' S8 v9 H7 w# n }//switch
}</P>
- L# z/ ^6 S) \3 e<>void banner() {
; I7 z+ y* q0 \9 _7 t printf("\t\t****************************************\n");
printf("\t\t         单纯型法解线性规划问题\n");
* _, `- G9 Z! s. H0 ^, i printf("\t\t                         作者：Thunder\n");
" M3 P9 X3 J! ?: t0 z) I printf("\t\t****************************************\n");
5 j: F+ C+ C  ^4 A; V( K printf("\n");
}</P>
1 R% s" Z4 ?8 v0 z* A<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
//为1的不显示系数1，-1系数只显示负号
! `- i3 D0 E+ ? int i,j,k;
switch(step) {
; A/ ^. X1 s& ^$ l" q8 s case 1:
  printf("min z= ");
  printf("x[%d]",arti[1]);
$ W" ]( u' ^( ]8 \  w8 b1 k  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
9 z1 V- p: I+ V: {, K  break;
) Q. T( m/ j& u1 G& U0 S! C case 2:
printf("max z= ");
+ j, W: A; x- B- [ printf("%lg x[%d]",c[1],1);
* z$ v2 ]' Z+ u" z8 u for(i=2;i&lt;=num_x;i++) {
& S6 O, m: H# ^8 n  if(c==1) printf(" + x[%d]",i);
1 G4 Q6 S  _& L- e5 M. j0 O" a5 S  else if(c==-1) printf(" - x[%d]",i);
" y+ B6 L3 {$ h4 S  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
  else printf(" %lg x[%d]",c,i);
}
  }* a8 m. V1 d5 k  w0 B break;
}</P>
  K8 `& r$ k' [<> printf("\nst:\n");
+ [0 F& H2 ?) J for(i=1;i&lt;=num_st;i++) {
* j! G" F, R9 @. [' [  k=0;
# c) V; c4 t. j- U- z0 ^6 {+ S. }  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
; X+ W3 x8 _" E* N' r    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
& k1 s9 M0 h2 g0 A0 L: c  s    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);
8 ]: T+ r) D7 W1 Y0 \    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;
   }
  }
) t' N$ D" {# C6 |; {% e6 a, G! {  printf(" == %lg\n",b);
}
: b  U4 O0 o) Z" q" l) i8 O8 e* ?6 r printf(" x[1]~x[%d]&gt;=0",num_x);
$ k4 [! D5 W! }( f}</P>
<>void exchange() {
4 w2 H% m5 h& y0 a int i;
double temp[MAX];
for(i=1;i&lt;=num_x;i++) {
  S& I( N  a/ M; o2 e+ n  temp=temp_c;
  temp_c=c;
3 v- n6 {4 C# ?/ {, w9 V  l! \  c=temp;
* G; y& o9 l) ?& p; c }
3 L; O$ f$ V) w}</P>
0 `7 \* }2 c, h$ k$ Q0 s- y<>void create() {
//输入方程组系数，每个方程输完后回显确认
9 H# |; ?' F6 V  q5 e, x int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
  D8 x6 Q7 J' ? char confirm;
% f5 x: F4 s  C+ n
while(1) {
  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  scanf("%d",&amp;ma_mi);
9 C3 c6 j2 E! Z! V$ x' F% W+ Y  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
  else break;
}

while(1) {
: `; L* Z; n* M" |7 k  printf("指定变量个数：");
5 L: ~- Z9 ~- K! D. Z7 x9 A" }  scanf("%d",&amp;num_x);
  printf("输入价值系数c1-c%d:\n",num_x);
# m  X) i% H! J" ~! D1 O  for(i=1;i&lt;=num_x;i++) {   
   printf("c%d=",i);
   scanf("%lf",&amp;c);
+ s/ Y# `$ l% x5 o# o0 p  }
  if(ma_mi==1) printf("max z= ");
% |$ E: m. L8 ?: z. W  else printf("min z= ");
  printf("%lg x[%d]",c[1],1);
# R5 C0 s& F( y4 l2 ~* e  for(i=2;i&lt;=num_x;i++) {
   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
8 a* g9 O' T8 h; R7 `7 }5 z1 {5 D+ P   else printf("  %lg x[%d]",c,i);
  }
printf("\n正确吗？：(y/n)");
/ j6 n3 }( C% G getchar();
  |+ L& g' e9 @% j confirm=getchar();
5 ?+ j$ A1 u; @ if (confirm=='y') break;
else if(confirm=='n') continue;
6 \$ |# d0 [% P$ o* a }</P>
<> printf("输入约束方程组个数：");
6 u* D8 x& K' R$ G4 g2 w' M scanf("%d",&amp;num_st);
7 D$ |/ K5 I' q1 i; z8 q for(i=1;i&lt;=num_st;i++) {
! O; G' o8 I7 B' ?) D  printf("st.%d:\n",i);
  while(1) {
5 \8 n+ f6 ]1 m7 C3 \8 O   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
- G# b. a: e; ]$ D; O* f* ?   scanf("%d",&amp;re_st);
. V7 I4 p+ ~; ^, u0 @0 t   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
4 ]- n, Z# \. }! q% M   else break;
  }
+ ?, i2 T4 F" I; i# _3 [$ C+ O  printf("输入技术系数：\n");
  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
" R" m, @1 [# q1 o% X9 |   scanf("%lf",&amp;a[j]);
/ ]3 Z( m+ @+ s$ X  }
6 n4 o, l9 k5 {7 f/ h  printf("输入资源拥有量：\nb%d=",i);
: t$ F- [% Z& p; Y  scanf("%lf",&amp;b);
  
, h# W% B3 e2 Y( Q$ r  printf("st.%i:\n",i);
! y# D3 @6 D; ]$ P  printf("%lg x[%d]",a[1],1);
" G4 O' O. L& F- Z& }  for(j=2;j&lt;=num_x;j++) {
& Y+ K6 @" M0 y2 v9 X/ \   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
/ }; p7 a. G2 @& _4 x5 Y! z   else printf(" %lg x[%d]",a[j],j);
  }
  switch(re_st) {
* ?/ z; l  {1 ?' z5 d+ l# j8 h   case 1: printf(" == %lg",b); break;
   case 2: printf(" &gt;= %lg",b); break;
; e$ ~% k" ^# M   case 3: printf(" &lt;= %lg",b); break;
$ }: |; @6 R$ ~  T5 F( l2 ]  }</P>
' G$ j; j2 j3 ^2 b<>  while(1) {
0 o6 e; }+ s/ ]3 H6 K6 d* v  J( z$ \   printf("\n正确吗？(y/n)");
   getchar();
% W7 d6 @( m) A4 y6 O: f9 [   confirm=getchar();
0 p3 m/ S8 ]- ?/ z+ m) K: U" n   if (confirm=='y') break;
) V# {1 O+ k' R8 }$ R   else if(confirm=='n') {i-=1; break;}
  }
}</P>
<>//显示输入的方程组
printf("\n原问题为：\n\n");
' l% Q9 G1 n$ K) Q$ { if(ma_mi==1) printf("max z= ");
, B3 B% L' m/ R% }: r2 U- c else printf("min z= ");
. v7 b; f- E! ~9 m1 W  D1 f1 }, I printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
3 D: _& T" {' y  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);
' Q$ n7 H! w: v  else printf(" %lg x[%d]",c,i);
}</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
& d/ E/ c/ O* c" i9 }" r  k=0;
; b+ N( C+ b# Z- |  for(j=1;j&lt;=num_x;j++) {
7 a" b: k' @* A; G/ I8 c" t   if(a[j]!=0) {
5 B, T0 r3 H0 C    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
' S) I0 U7 S6 g; w* E4 V% d! C    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
- x3 O# r) m% [( m9 U' `3 Z    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);
1 _( e* _9 }& @: i! v    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;
2 G7 I2 X/ W8 k* d6 T- B   }
  }
  switch(re_st) {
   case 1:
    printf(" == %lg\n",b);
    break;
! `; D! |5 U" ?; l   case 2:
    printf(" &gt;= %lg\n",b);
9 A) D9 @2 t; N0 Z& @: }    break;
   case 3:
. P+ \3 `/ Q& A5 D8 K    printf(" &lt;= %lg\n",b);
6 R/ \  s  q; Z8 s; U    break;
  }
}
0 y; w; a/ ^! S8 k# Q, | printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
! X& s* d/ ]9 m5 r0 r<> tnum_x=num_x;
for(i=1;i&lt;=num_st;i++) {
" }3 k/ a8 D# k7 I& e% }  switch(re_st) {
  case 1:
( S: R/ h' F' z7 p- @+ X  case 3:
8 {1 H9 R5 C6 P5 D; @   num_x+=1;
   break;
5 u* I9 n# X! D- @. H  case 2:
) y# p0 f8 O& |1 K   num_x+=2;
   break;
- _" M1 \0 R; k% ^: r' V  }
( ^, L2 l0 D& ^ }</P>
" x5 [. u1 D# V: h" Q& X<>//化为标准形式
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) {
   case 1:
    num_addv++;
    num_ba++;
) h% G* |9 l" r5 ~    num_ar++;
    c[tnum_x+num_addv]=0;
    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
    for(j=tnum_x+1;j&lt;=num_x;j++)
* K) Q2 e/ y4 e     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
- f( D$ w7 `0 y5 D' Y6 \    break;
   case 2:
    num_addv++;
: N( y5 u9 D# O  f$ Q. z    c[tnum_x+num_addv]=0;
    num_addv++;
    num_ba++;
" f  v+ t0 N2 D0 P& y, y    num_ar++;
    c[tnum_x+num_addv]=0;
    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;
3 X5 o! i, h; l1 l9 u# z# x     else a[j]=0;
    break;
" r/ H8 N+ F2 q- w   case 3:
    num_addv++;
    num_ba++;
1 ~" `9 ~7 _5 J' F, c% u( Q    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;
5 u6 n3 R8 i* d0 |  `& L9 h     else a[j]=0;
    break;
- x: v. d( @1 l* E  F: G# w. a  }//switch
, l( ~" x, G- P0 \) ~6 i }//增加松弛变量、剩余变量、人工变量、确定基变量</P>
<>//显示标准化后的方程组
printf("\n化为标准形式后：\n\n");
, j  N# M9 S7 S7 e( y if(ma_mi==1) printf("max z= ");
7 A! D. [$ ~- a! _9 W' f else printf("max z'= ");
printf("%lg x[%d]",c[1],1);
5 M! k5 }. m6 [7 r$ L for(i=2;i&lt;=num_x;i++) {
  p$ H7 U) q# I. D# q/ |  k=0;
# k+ E, f% Z& h, y- N  for(j=1;j&lt;=num_ar;j++)
& W% m2 `; m4 [- }( N7 H5 x( Z0 v  if(i==arti[j]) k=1;
  if(k==1) printf(" -M x[%d]",i);
  else if(c==1) printf(" + x[%d]",i);
( H1 G; F" x' a. L3 h5 m3 J8 @  else if(c==-1) printf(" - x[%d]",i);
  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
5 X) T) T8 W* P) }5 ^+ o  D  else printf(" %lg x[%d]",c,i);
}</P>
% f/ e; u7 _0 h/ l/ `<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  k=0;
  v) r6 x1 j0 N& J8 `) ^: ^" A3 k) E  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
- ?: H1 S( J  o8 ~$ H    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
: i6 I& r: U8 O7 v# n2 @9 A    else if(a[j]==-1) printf(" - x[%d]",j);
+ g5 j  G( [6 w: h9 s& ~( v    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
* o4 m( I  `1 i" q. z7 n    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
3 |+ `+ s& u- M3 z3 B    else printf(" %lg x[%d]",a[j],j);
  M3 l2 J# h2 p& f! ^    k=1;
   }
  }
: ^% A7 ]. H3 i  printf(" == %lg\n",b);
: `5 J# V: {& Y- ~: x }
printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
. o% Z8 T9 I- D& V" @3 b<>void iterative() {
int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
; Z3 L9 I4 ~, h! B$ r$ C  _ int base_elem;
int base_out,base_in;
double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
) u; L5 E( A+ y0 I+ `<> printf("\n\n第%d次迭代：\n\n",stop);
for(i=1;i&lt;=num_st;i++) {
- F9 ?9 P3 E! [' z# b1 k  printf("c%d=%lg\t",base,c[base]);
  printf("b%d=%lg\t",i,b);</P>
/ C8 P* m& W5 j5 C: P<>  switch(step) {
" W% u$ P3 ]6 @   case 1:
    for(j=1;j&lt;=num_x;j++)
    {
     printf("a[%d][%d]=%lg\t",i,j,a[j]);
    }
    printf("\n");
    break;
$ o2 R: F5 A8 I7 z+ F: }& F2 d   case 2:
    for(j=1;j&lt;=num_x;j++) {
     k_a=0;
     for(l=1;l&lt;=num_ar;l++) if(j==arti[l])k_a=1;
' ]& y$ d. L) _. K! N: Q2 T) k+ n     if(k_a!=1) printf("a[%d][%d]=%lg\t",i,j,a[j]);
* b' O3 B2 M0 q3 ~; V, b* P    }
    printf("\n");
    break;
% m8 A- Z' g) p( a  }
0 ~' Y! ^/ z8 F, a  j5 T }
; j& n: ]2 w) ?0 f//求检验数sigma
* y- v* Z' s7 X2 V) P2 y& F for(i=1;i&lt;=num_x;i++) {
  sigma=c;
" Y. S: i1 i+ C; q7 x4 r3 n  for(j=1;j&lt;=num_st;j++) sigma-=c[base[j]]*a[j];
( u8 J+ l  q/ ]$ i. n  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
" |; K% _. P* Q3 Y: y  switch(step) {
: w, v2 Z0 e% R1 v7 z+ z2 ?   case 1:
    printf("sigma[%d]=%lg\t",i,sigma);
1 S; N9 g0 {" D  @8 o6 [3 |2 S    break;
   case 2:
6 B9 i/ b5 E' c. K( N% {    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);
& S$ N/ R+ F4 [    break;
4 s1 D! @. r( e* B. V. c; |7 w' e+ H  }
}
putchar('\n');
* s3 P/ V( w; @. |$ f6 Q//检验检验数sigma是否全小于等于0
k=0;
7 i5 M$ P2 B6 T) g for(i=1;i&lt;=num_x;i++) {
  if(sigma&gt;0)
  k=1;
}
0 j4 @+ W4 |& j2 o( x3 K if(k==0) {
& ]# `  @' A5 y, u  //sigma是全小于等于0时，检查是否为无穷多最优解
  for(i=1;i&lt;=num_x;i++) {
   k_f=k_a=0;
- \' z# l0 }5 z   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;}
   }
; q3 L6 j$ i* k* [  }
5 S% K( p* H3 J# o' f. l4 c' d  status=1;
  return;
}
6 ^$ b5 |: Q; g/ `//检查是否为无界解
4 B; O& J" L+ V/ A* R1 e& X for(i=1;i&lt;=num_x;i++) {
  k_f=0;
  if(sigma&gt;0) {
   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
4 n5 y; Q, u8 |( |; m8 D: y   if(k_f!=1) {status=0; return;}
  }
}</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;
  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
  o& N  Q3 b! W7 }' b: k: t }//temp赋初值
4 _  w4 W- q4 U6 B+ q  P  N+ l4 \4 i for(i=1;i&lt;=num_x;i++) {
  k=0;
. _& p, u5 ?. F6 F9 r  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
; J8 }5 b9 @9 D1 W, ^! t3 Q0 c7 d  if(k==0)
   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
    base_in=i;
2 x- j, M$ S" J    temp=sigma;
   }
}</P>
) y' f; Z" T; E+ F<>//确定换出变量
for(i=1;i&lt;=num_st;i++)
/ d" e+ R# @2 q. O* `8 |& e$ U  if(a[base_in]&gt;0) {
8 [+ o. Z  S9 x   temp=b/a[base_in]+1;
. Z+ X2 p& r9 {   break;
  }//temp赋初值
- h7 C; h) z! s" h for(i=1;i&lt;=num_st;i++) {
8 J' s' w4 F. }" V, H* R* _  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
   for(j=1;j&lt;=num_ar;j++)
    if(base==arti[j]) {
     base_out=base;
( }* Z* \9 Z8 I$ E* ^     base_elem=i;
     temp=b/a[base_in];
     break;
    }
9 ]5 W: z1 R' N8 Y" P# V( `  }//人工变量优先换出
1 w3 L' v+ Z9 H/ g  ]  c$ U5 j  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
   base_out=base;
   base_elem=i;
# Q2 K$ V& N3 @+ R   temp=b/a[base_in];
  }
}</P>
( z+ C1 b& `0 L' S( N<> printf(" 基变量：");
for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
6 _: K$ p/ V' X printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
//基变量变换，进行新方程初始化后迭代
: |: D$ E4 b% R; y for(i=1;i&lt;=num_st;i++) {
2 L4 ?) d8 W/ [. L3 \  if(base==base_out) base=base_in;
}
+ x* H& d% O2 X9 ]//初始化主元素行系数
value_be=a[base_elem][base_in];
b[base_elem]/=value_be;
for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
8 `  o) J& l- [! [% ]6 o) a6 i<> for(i=1;i&lt;=num_st;i++) {
  if(i!=base_elem) {
   b-=b[base_elem]*a[base_in];
   value_be=a[base_in];
2 Z% ^6 R4 Y# T2 X: @1 Q" a   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  }
3 L# p/ c4 p+ S4 ~2 s+ o7 | }
stop++;
4 ^( A/ i6 W9 k( g9 [2 N if(stop&gt;STP) {status=-2; return;}
iterative();
* O* g8 V  M7 s, V}</P>
2 v) l! A" p0 o; o* a) r<>void output() {
4 M6 w9 U) f' C, x2 W9 x! d6 W int i,j;
double X[MAX];
$ h2 ~6 W$ X  p2 @5 @) L printf("\n结果如下：\n");
printf("\nX=(");
for(i=1;i&lt;=num_x;i++) {
: b. M6 x1 b' L% _7 Z" P  for(j=1;j&lt;=num_st;j++)
  if(i==base[j]) {X=b[j];break;}
  else X=0;
8 Z$ F9 Z4 ^, r  l/ W  printf("%lg ",X);
& \) ^7 U* P/ f4 A; X$ C! 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);
' R( ^; r/ Y; y8 |. s0 u! s8 D}</P>

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
1 l+ L2 n( D6 J( c; ?#include&lt;math.h&gt;
#define X 5
4 \) h. d; V" \+ \! M#define Y 7
! D' c% Q* `* ^, Yvoid xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
{
/ H! e) d6 j5 d$ T5 S) v  int j;
  *c=0;
% B- U: _# v8 ^* X  for(j=1;j&lt;=*is;j++)
    {
! ^# V5 `4 r: @5 L      if((*a)[*ir][j]-*c&gt;0)
    {
2 V+ x  ^, H7 q/ S: Q+ t      *c=(*a)[*ir][j];
9 J0 u* U/ [6 @; T$ j. ~* P& ^      *j0=j;
6 X$ D8 M1 S2 w; }- y1 w    }
& g/ ~) v, q+ s( R# P+ g" d     }
}</P><>      /***************** 参数说明 **********************/
$ \7 _9 |" I8 c1 W" k, b      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
$ v- ~- t7 W: s0 P4 Q5 O      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
      /********* <a href="http://happyyangxu.home.sunbo.net*******/" target="_blank" >http://happyyangxu.home.sunbo.net*******/</A>
- Z6 C4 Q- [6 Y2 O6 [int xi_sm(int m,int n,int m2,int mn1,int l1,float (*a)[X][Y],
( ?1 O1 Q, R' v; B: Z+ b" E       int(*k)[],float(*x)[])
{
  int i,m1,mn,j0,i0,j;
  float c,g;
, Q% g; F1 q% X$ X1 m' M( |  m1=m+1;
% ]' z; p2 x. c4 T* o  mn=m+n;
) ~. x2 Q" O7 w( D3 r' w2 M  for(i=1;i&lt;=m;i++)
    (*k)=n+i;
7 b* j7 E8 a, Y' M0 J0 i* kleap1:
  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  else
3 P! }% X% F- |" V     {
' K6 L1 V, _4 W0 C( @leap2:   if(l1-50==0)
      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
# U. ?! k: X" _* p  M       else
      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
     }
# t8 [+ B" t4 @0 K' G- u   c-=1e-8;
   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");
1 m" g7 ], t# z& Y- b; ^% {    return(-1);
% N) E  ]. t1 H4 ^/ n      }
9 ^  O& Z" `# t6 `% z4 a4 u   if((c&lt;0)&amp;&amp;(l1==1))
) O  n0 ^; y! U# h- @! q       {
) {5 r& U, J7 ^, [& M3 W1 B     l1=50;
     goto leap2;
, m  v( _" Q- D, X( V  @/ S3 u  {       }
    if(c&lt;=0)
# N5 d) H$ B# A; h2 _9 F" v       {
      for(i=1;i&lt;=m;i++)
         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
2 @/ x2 D3 o+ a! _" j* K      for(i=1;i&lt;=m;i++)
' Y; R/ j4 g  G- W3 S, C* F         printf("\ni=%d,x=%f\n",(*k),(*x));
7 a5 ^, S( Z! p5 E( E6 l& l      printf("\nF=%f\n",(*a)[m2][mn1]);
      return 0;
+ L. g+ l! @' D  z6 P    }
) L6 P2 T0 O! y     c=1e8;
     for(i=1;i&lt;=m;i++)
6 u* J7 E/ W* p# m" _       {
: u) A* O& E6 ?  K     if((*a)[j0]&gt;1e-8)
# n" Z- u) Z# K        {
          g=(*a)[mn1]/(*a)[j0];
          if(g-c&lt;0);
         {
  P( N: l* A' j/ |" A, {           c=g;
           i0=i;
9 J0 R' t  m! z  Q4 N          }
. k: f8 D2 ^' n8 K  ~  h0 G/ Y; r         }
    }
( S* e. j" d. b3 D3 a      if(c==1e+8)
, Z6 ^6 r, X3 P* o) u    {
9 F8 y; Y2 K& Z8 {      printf("\n\t*******Lp no solution********\n");
7 d2 C9 |, g: c      return -3;
     }
* V- E3 S3 U6 l1 ~# G       (*k)[i0]=j0;
; L: B3 C; V) j4 D8 _) M       for(j=1;j&lt;=mn1;j++)
; Q$ ~2 c9 M1 K$ n7 v      {
        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
           continue;
" E- s- I* y8 w1 u7 j9 F& z! \         g=(*a)[i0][j]/(*a)[i0][j0];
           (*a)[i0][j]=g;
         for(i=1;i&lt;=m2;i++)
           {
         if((i==i0)||(!(l1==1)&amp;(i==m1)))
/ ^6 W! X6 S) u% }4 ~            continue;
2 t3 [1 r; E+ ^: T         (*a)[j]=(*a)[j]-(*a)[j0]*g;
# E( F5 l# N4 H           }
       }
) y" O$ F2 O  w9 `0 I! W2 U       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
0 l$ Y  C2 t- |       (*a)[i0][j0]=1;
       goto leap1;
}  
main()
{
  float a[5][7]={{0,0,0, 0,  0,0,0},
, Y9 ]) v2 F  M* F$ a  j* ]6 i, p1 D                 {0,3,-4,3,  1,0,12},
                 {0,3,0, 6,  0,1,12},
$ p- D4 m9 n$ e9 c                 {0,0,0, 0,  0,0,0},
                 {0,69,0,144,0,0,300}};
  int k[3];
  float x[3];
9 a. g3 z$ X6 B6 m. c0 C  clrscr();
  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
}
</P>

julysea

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

lauqt1980612

最好是C语言的!