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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
: B& Q7 f" i4 J3 C#include&lt;math.h&gt;
+ a0 Z3 n! h3 b+ K  z! m#define X 5
0 n& Y) C; @/ B5 k#define Y 7
. f: N8 z3 K, y' [  R  M" _0 P) Fvoid xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
{
  int j;
( ^/ ^3 [3 y4 |9 ]$ U  *c=0;
$ f9 F* U# n. H$ b: Z3 _  for(j=1;j&lt;=*is;j++)
) {& B* m' y, i+ i    {
2 G/ V( \! J7 ?/ S# y# c' _3 h! f      if((*a)[*ir][j]-*c&gt;0)
    {
      *c=(*a)[*ir][j];
. T( E4 c1 t; i4 V      *j0=j;
9 P  X1 G- ?6 M# i+ z9 P    }
     }
# i% [8 G" m: j& Z" [/ q}</P><>      /***************** 参数说明 **********************/
      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
6 u8 t2 i. p' L& K+ ]6 r0 }* O      /********* <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],
! J  t" F; q  [! t* [( t       int(*k)[],float(*x)[])
% G7 @, ?* K' k5 m/ _. `{
  int i,m1,mn,j0,i0,j;
  float c,g;
  m1=m+1;
  mn=m+n;
  for(i=1;i&lt;=m;i++)
    (*k)=n+i;
1 j9 Y: B# a- @3 Y. V" W+ Lleap1:
  if(l1-1==0)
" l1 O/ q8 m7 I     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  else
     {
0 B5 s5 L' O& u9 n8 Nleap2:   if(l1-50==0)
      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
# k& K5 `: f  j0 r$ T1 H; \       else
      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
     }
9 @6 q4 A6 `# M0 {9 j! s   c-=1e-8;
   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
      {
' ?8 w0 k3 e, g0 L2 {    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
    return(-1);
      }
  ~5 ~  G( G; L( {" D   if((c&lt;0)&amp;&amp;(l1==1))
! ~8 T1 Q- y' \1 W. Y5 @       {
     l1=50;
     goto leap2;
' j0 q5 {$ i7 z8 V( v3 [( t       }
8 ^9 c4 @/ p: l0 O( M# y) Q; C    if(c&lt;=0)
1 s2 D  k0 j8 a9 }/ C       {
      for(i=1;i&lt;=m;i++)
; J6 @4 ?- n7 z0 f         (*x)=(*a)[mn1];
4 ?  h( M9 x- V5 g" |, K- y7 N" r% w$ X         printf("\n\t********Optimal solution********\n");
. c, K4 Z& x  O+ d      for(i=1;i&lt;=m;i++)
3 o: K; X+ ?  c1 p2 X4 k' g         printf("\ni=%d,x=%f\n",(*k),(*x));
      printf("\nF=%f\n",(*a)[m2][mn1]);
      return 0;
    }
     c=1e8;
     for(i=1;i&lt;=m;i++)
       {
     if((*a)[j0]&gt;1e-8)
        {
          g=(*a)[mn1]/(*a)[j0];
5 f4 M- N( ]5 Y$ [" S' u& R          if(g-c&lt;0);
         {
           c=g;
6 u# P8 a6 r1 u4 z1 Z% N           i0=i;
          }
         }
& {7 |" \$ ?$ F5 d  V- M5 i2 o& a    }
9 W1 ^& ~: ?% S) u  E4 v      if(c==1e+8)
; U+ @9 \+ W4 w7 N9 \    {
/ S/ K; m/ }+ N- m; E0 Y  z      printf("\n\t*******Lp no solution********\n");
      return -3;
5 k: Q: ~3 h! x! o5 c7 ?$ `6 M0 \     }
+ S" }, P# L: v" ?' p: e4 n+ w# M) H0 u" h       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
% h5 w+ P0 L6 p7 U      {
        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
7 p4 f; ~+ e0 b6 M. C5 [           continue;
  c" R6 F- u! [1 M         g=(*a)[i0][j]/(*a)[i0][j0];
" S, w0 u& j" L+ S6 N7 ^& n/ @           (*a)[i0][j]=g;
         for(i=1;i&lt;=m2;i++)
3 C# J1 ]$ |/ l/ a  n0 f           {
0 n2 p6 t, m' z/ l2 t" U' u2 ]         if((i==i0)||(!(l1==1)&amp;(i==m1)))
            continue;
         (*a)[j]=(*a)[j]-(*a)[j0]*g;
           }
       }
       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
       (*a)[i0][j0]=1;
$ u$ l+ A5 ~: e3 n# q* T5 y       goto leap1;
}  
main()
. Q* w% Y4 b8 o1 y( B+ G{
9 ~9 n1 K8 Q. y& k% e  float a[5][7]={{0,0,0, 0,  0,0,0},
                 {0,3,-4,3,  1,0,12},
7 H2 {6 R9 ?8 m                 {0,3,0, 6,  0,1,12},
                 {0,0,0, 0,  0,0,0},
) K5 T# g. W$ T5 C2 R  j& j                 {0,69,0,144,0,0,300}};
% R8 J. V7 b! Y/ P0 @6 v( Q  int k[3];
7 V7 l' g- b* ^2 p8 s. Y9 v5 G  float x[3];
: l" M* u6 L7 c/ d( R  clrscr();
! }+ d! |* P- F$ J9 x7 o8 s9 [  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
}
</P>

plgatc

<>/*************************************************************************
2 E; H9 J9 W+ C: U                    单纯型法解线性规划问题（两阶段法）        
; I! g' _( M" l9 X, X: S
编程环境：VC++6.0     
方程组输入说明：
变量非负，按提示输入相关参数。
' Z: \* g. L6 s: l" _ *************************************************************************/
#include &lt;stdio.h&gt;
#include &lt;stdlib.h&gt;
#define MAX 100
#define STP 100</P>
<>int stop=1; //迭代记数变量
5 S7 D! `% c$ o4 B0 I+ v, s$ {int status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
- t0 A, Q- Q$ ]8 mint step=1; //目前阶段</P>
<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
, h; K. o! R7 F- \/ I3 ^int num_x; //变量个数
int num_st; //约束方程数
int num_ar=0; //人工变量个数
9 X$ m& q! M. dint arti[MAX]; //人工变量下标
- Z) y5 e* l. c/ L3 C4 F8 pint base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
8 o& r1 r$ Y3 t2 i$ e4 w, l4 ^; O<>void create(); //建立方程组
void iterative(); //单纯型法迭代
void output(); //输出结果
' R# H: Z# n( @void banner(); //打印程序标题
void exchange(); //交换两阶段价值系数
# M+ m; M* M& W. s( f+ @  Kvoid show(); //输出方程组</P>
<>void main() {
int i,j,k;
8 c) l# D% F, u( L banner();
0 j% I9 V( _: Q) g. P7 ~8 n) D create();
//保存原价值系数，转换为第一阶段价值系数
for(i=1;i&lt;=num_x;i++) {
, s0 j. Z: D; ~  k=0;
( {& x  g5 m8 b  Y; [  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>
<> printf("\n\n第一阶段问题为：\n\n");
show();
step++;
/ S4 I. @% I/ S1 } printf("\n\n按回车开始第一阶段迭代");
' s( J% h( C/ r4 ?( F getchar();
& I* O% L7 L5 b7 W getchar();
: V5 w4 F' \$ _5 d- k iterative();
; r. T, P' }9 d/ D( R+ L- N if(status==-2) {
  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
( n: X! g' Q& ?& _# O* l9 B }
output();</P>
<> if(max!=0) {
4 E8 f% _5 o' l! x: \, I  puts("\n\n原问题无可行解。\n");
  puts("\n按回车结束");
' ~7 x. {( j7 L+ J$ A  getchar();
  exit(0);
}</P>
: N5 [4 E3 c# j6 u, {<> //转换为第二阶段价值系数
, a) L. I' F1 n exchange(c,temp_c);
0 G: n( o( i* L9 p6 b; B //把人工变量列全设为0
1 }  P+ t$ h" g5 L8 R1 g/ x/ T$ N for(i=1;i&lt;=num_ar;i++) {
  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
* |+ \2 P; `* N8 ]4 `! Q# w }</P>
4 \' f# S) n6 ?7 R: V<> puts("\n\n第二阶段问题为：\n\n");
show();
puts("\n\n按回车开始第二阶段迭代");
- }$ O# I; g4 o  ^/ F) U getchar();
iterative();
switch(status) {
case 1:
  output();
  puts("\n\n原问题有唯一最优解。\n");
" U  z' |( }) }  N9 M3 o2 g  puts("\n按回车结束");
! D( a! }1 `9 r/ M  \8 B  getchar();
  exit(0);
case 0:
  puts("\n\n原问题为无界解。\n");
7 W5 P9 \/ ~3 V" ]& N- B( j* O  {  puts("\n按回车结束");
  E, }: l  z+ U& U! W  getchar();
/ y1 ]0 a* g% @/ I# _6 _5 t& G  exit(0);
/ E. r7 i8 [0 G  a) [' J case -1:
  output();
  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
case -2:
$ U4 T$ c7 t6 s8 J  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
; J3 ?5 Y6 y4 u9 B1 M5 X }//switch
}</P>
<>void banner() {
& E4 R9 b, r: y; O printf("\t\t****************************************\n");
printf("\t\t         单纯型法解线性规划问题\n");
& P7 |- D, B0 L3 P6 [( @ printf("\t\t                         作者：Thunder\n");
/ P* k3 i$ m: a* P( ]7 w+ o printf("\t\t****************************************\n");
printf("\n");
}</P>
<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
//为1的不显示系数1，-1系数只显示负号
  r5 V& [# Z1 _ int i,j,k;
  W( l% y" a, P. @' f$ _1 w switch(step) {
% ^/ L: H; }. N2 ^5 | case 1:
  printf("min z= ");
# ]( Y/ Y$ l( N  printf("x[%d]",arti[1]);
# K! h7 S  b: l9 ^  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
1 g- o8 i5 h: i; t% d' f. R  break;
case 2:
' F) c4 B8 m. U7 ~* S printf("max z= ");
printf("%lg x[%d]",c[1],1);
! `$ y6 Y2 ?8 l' N9 a# z2 J* l for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
+ f5 `+ F1 S. X- L6 n0 c& z  else if(c==-1) printf(" - x[%d]",i);
  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
, v% \  b2 L* v  f3 }8 `+ i  else printf(" %lg x[%d]",c,i);
}
break;
- n+ _0 a) |! _ }</P>
) W" ]. ]5 K7 `1 E9 l$ e<> printf("\nst:\n");
: J2 U9 C6 E9 w. c for(i=1;i&lt;=num_st;i++) {
- y% p! p/ U0 p1 D  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);
1 d/ D" Y2 c7 \8 M, e2 S$ Q    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
" ~/ \: K: W% n: E; v    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);
" C7 R9 V, ?( ~" p0 w& |% v! a    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
" H, ?% @( u% f. P& \9 d    else printf(" %lg x[%d]",a[j],j);
4 v( d! h7 d8 h! V5 ?% Q    k=1;
   }
/ d- o, r4 m" u+ s  }
, L/ a3 l: f+ l5 ^) Y4 Y  printf(" == %lg\n",b);
}
printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
<>void exchange() {
int i;
double temp[MAX];
$ W) W5 R" E6 m9 n( O" r* U for(i=1;i&lt;=num_x;i++) {
  temp=temp_c;
/ R4 u" h" _) Y  temp_c=c;
# N+ }& z6 l; F- |1 @  c=temp;
}
/ F( Y8 m1 J/ e4 ]8 U1 n' {}</P>
3 h' u2 l9 u% k. q: x% n5 O/ ]<>void create() {
0 C5 p1 W4 B2 M//输入方程组系数，每个方程输完后回显确认
! Q* M% V9 l2 c) N, W int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
. }& P9 {, X& ?3 F# l. ^7 G; b char confirm;

while(1) {
  printf("请选择：1、求最大值，2、求最小值：(1/2)");
' s6 B3 T6 @$ S4 A  scanf("%d",&amp;ma_mi);
1 ]$ g7 q% f. V' h  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
% j( ]% C, k& q' |! N  else break;
1 ?7 \. g( x9 E9 l0 |4 N }
3 V; E. R- r- i
while(1) {
  printf("指定变量个数：");
  scanf("%d",&amp;num_x);
: S; D. Y( |- g! W6 y8 ?  printf("输入价值系数c1-c%d:\n",num_x);
  for(i=1;i&lt;=num_x;i++) {   
: W8 `$ Y/ v4 M6 ~& L! M4 O   printf("c%d=",i);
. o1 ]9 B# E4 K4 V! V   scanf("%lf",&amp;c);
* {0 C5 B+ r0 Y  }
, z+ n& u- T2 D- ~  if(ma_mi==1) printf("max z= ");
2 i8 x3 C3 s3 p  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);
  [/ h5 m1 L/ Z# A   else printf("  %lg x[%d]",c,i);
  }
printf("\n正确吗？：(y/n)");
getchar();
9 }5 A( _. `; S4 x6 \  K confirm=getchar();
/ V( N4 F) @2 u9 b% n9 Z* V% e2 _ if (confirm=='y') break;
else if(confirm=='n') continue;
}</P>
<> printf("输入约束方程组个数：");
scanf("%d",&amp;num_st);
for(i=1;i&lt;=num_st;i++) {
  printf("st.%d:\n",i);
- C& q$ u& b7 G( P/ [  while(1) {
   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("输入错误，请重新选择。");
- H* r, V) R6 c  N4 `7 y* P- M   else break;
6 K+ A/ N4 @  [$ z8 M, f- ^+ d4 v  }
  printf("输入技术系数：\n");
  for(j=1;j&lt;=num_x;j++) {  
% V: s- M+ `0 ^% v   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
  }
  printf("输入资源拥有量：\nb%d=",i);
4 n) Q9 Z& F  r: U7 U' e) u  scanf("%lf",&amp;b);
9 \0 B( l/ X2 X6 L* L: i+ S0 X! m  
  printf("st.%i:\n",i);
; _( Z7 i- L7 U% P- H* I7 y  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);
( a9 j& M6 s6 r6 p# C4 E: P   else printf(" %lg x[%d]",a[j],j);
" q; E. i0 N# j2 ^  }
1 }3 G" b+ G! }/ y+ Y$ I  switch(re_st) {
6 o+ C/ V# r' N$ U/ R7 ]: j   case 1: printf(" == %lg",b); break;
   case 2: printf(" &gt;= %lg",b); break;
; w* ?3 ?" y$ U1 z5 o' r1 B6 U7 A   case 3: printf(" &lt;= %lg",b); break;
$ [0 [% q. `% V! B" G  }</P>
<>  while(1) {
   printf("\n正确吗？(y/n)");
   getchar();
+ V5 t0 e4 i& z   confirm=getchar();
   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
  }
}</P>
<>//显示输入的方程组
8 C" @$ ?$ R$ }9 Z1 H printf("\n原问题为：\n\n");
! p  X  K% x+ l& _- f if(ma_mi==1) printf("max z= ");
6 x9 ~9 y7 R" \4 X else printf("min z= ");
* P! Q, |: x7 K printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
* ]& L! Q( f$ C' _& ?0 L  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);
! Z- y/ u5 A. C& ~. E# |  else printf(" %lg x[%d]",c,i);
}</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
- p- q  [) y: X8 B: ]9 {  k=0;
% J7 s' W$ f1 g. {1 G' B7 v  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
$ @# R+ D8 o  Z3 I, H3 [    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
' L, H/ S" m* h    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);
    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
+ \3 T+ y. ^0 v' f; q( ~' z    else printf(" %lg x[%d]",a[j],j);
- I5 Q" p6 v9 h7 ?0 f/ R    k=1;
8 F. U+ [" n9 ?, v" R, d9 r   }
  }
4 j8 d6 t7 B; s  switch(re_st) {
   case 1:
    printf(" == %lg\n",b);
3 Y7 H6 |8 e" j; L    break;
5 ]! G' W9 w+ L( t   case 2:
    printf(" &gt;= %lg\n",b);
    break;
   case 3:
    printf(" &lt;= %lg\n",b);
" R6 ?6 g* h: A3 _! L% A    break;
2 z5 K! M$ @- w( }  W: |  }
}
printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
: D+ s7 i: ?0 M" {5 {1 \<> tnum_x=num_x;
for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
  case 1:
  case 3:
, y* T3 r. y! w1 D/ D5 ^   num_x+=1;
8 F: k/ w9 x' F; {: X2 V! e! t4 g2 B   break;
! X# a% X4 j* `  d' t$ z" t  case 2:
   num_x+=2;
   break;
  }
}</P>
<>//化为标准形式
$ \3 H1 }* u4 b! w! F4 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:
- A/ z$ Q5 V( ~( O5 I' e# B    num_addv++;
    num_ba++;
+ H1 t( ?. R5 N5 n" s5 `    num_ar++;
    c[tnum_x+num_addv]=0;
! m; w  B; `& I& K; \    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;
5 o: K' R  k' f* i9 X, M# Y     else a[j]=0;
$ }2 n2 n* e& p" ^    break;
& l& m7 W+ v: D/ V* O# |) g' D   case 2:
    num_addv++;
    c[tnum_x+num_addv]=0;
    num_addv++;
    num_ba++;
    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++)
# z6 Z6 b1 F( h+ p     if(j==tnum_x+num_addv-1) a[tnum_x+num_addv-1]=-1;
; d9 M* Z0 C0 v* L- ]     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
    break;
   case 3:
' `9 b& s6 L9 k# g' A4 O& c4 Z2 W    num_addv++;
3 [" y% q# \. }7 K" W& k' C  h    num_ba++;
! g' m4 `  F. w" G' c    c[tnum_x+num_addv]=0;
    base[num_ba]=tnum_x+num_addv;
1 j7 s9 A- v1 c" g. M    for(j=tnum_x+1;j&lt;=num_x;j++)
     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
  M# x( Y/ r9 [3 {, }; U0 D     else a[j]=0;
6 l- Y- y7 d' }8 g    break;
  }//switch
. Z  P) C3 J! I. d2 Y( G  S }//增加松弛变量、剩余变量、人工变量、确定基变量</P>
- B$ n2 @6 C) g; d<>//显示标准化后的方程组
$ R4 e0 z( B) t$ l! F printf("\n化为标准形式后：\n\n");
if(ma_mi==1) printf("max z= ");
" F; Q5 k/ R& ~- {% T else printf("max z'= ");
: A  W: I) b: ]. j/ h 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++)
$ E' f: F! l; s" R% Y: C1 c! M0 u  if(i==arti[j]) k=1;
  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);
2 N' g  `$ ]& T  L4 d" R8 f7 {% W  else printf(" %lg x[%d]",c,i);
}</P>
  q* f( e" n1 W6 V/ v1 E" t1 B, ]" _<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
* @$ F. k. f! k+ S, l/ m1 W  k=0;
  for(j=1;j&lt;=num_x;j++) {
9 B$ Q4 w0 s5 _. |- n   if(a[j]!=0) {
& S5 M; T% Q: J/ P2 t1 u+ ?    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);
8 N" m) R& T0 M9 @4 Z' V/ n' C    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;
   }
  }
* V% Y6 M- y2 x9 R4 l  printf(" == %lg\n",b);
( M0 K2 l3 g1 C' q8 Z }
, U4 a7 i2 g6 ~1 `$ z+ x printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
<>void iterative() {
7 R# z6 S9 ?5 @6 M0 e 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;
0 K& Z) `& X, Y  u2 L9 I% N6 t7 ] double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
% ]! b+ g$ v: t  w+ e  B  U3 N<> printf("\n\n第%d次迭代：\n\n",stop);
for(i=1;i&lt;=num_st;i++) {
  printf("c%d=%lg\t",base,c[base]);
  printf("b%d=%lg\t",i,b);</P>
<>  switch(step) {
: S6 @4 b/ v% }3 d   case 1:
8 S5 L% A5 M: Z) u    for(j=1;j&lt;=num_x;j++)
6 ^  Z1 |# R. y3 Q  ]    {
* b# j9 h1 j- [: `4 _     printf("a[%d][%d]=%lg\t",i,j,a[j]);
/ r1 ~  p$ d4 @5 e" I; F" c    }
8 V3 {' h% q8 E+ E* l    printf("\n");
    break;
   case 2:
    for(j=1;j&lt;=num_x;j++) {
5 y1 \) Z! @/ d& m     k_a=0;
1 n: _: F$ c0 R3 a9 n     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]);
" ~8 K3 e6 H: l" G0 k: U    }
    printf("\n");
    break;
  }
* R0 }/ L, L5 R* y# I6 Q }
//求检验数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];
( G8 A& M' T" o7 q: E1 i5 \* n  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
1 b6 U" `8 g" i7 c5 \; z  switch(step) {
   case 1:
  u$ ^: w( g. y    printf("sigma[%d]=%lg\t",i,sigma);
    break;
   case 2:
    k_a=0;
    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
) g& H* n: K% A    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
    break;
  }
}
+ _6 h& {9 G( ?1 X putchar('\n');
% e7 a: v& m2 B8 [//检验检验数sigma是否全小于等于0
k=0;
for(i=1;i&lt;=num_x;i++) {
  if(sigma&gt;0)
  k=1;
}
! A+ G. b) {# x2 w0 Y4 k if(k==0) {
. ^4 g' E4 g! w! I+ i3 w& `  //sigma是全小于等于0时，检查是否为无穷多最优解
+ Y- {1 B$ z4 k' D  for(i=1;i&lt;=num_x;i++) {
   k_f=k_a=0;
6 _7 `& }$ y4 o4 p5 O   for(j=1;j&lt;=num_ar;j++)
   if(i==arti[j]) k_a=1;
   if(sigma==0&amp;&amp;k_a!=1) {
2 t8 q7 T2 S. r) Q! o    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
; y+ y# h7 e* P$ P% V, ?/ X   if(k_f==0) {status=-1; return;}
   }
  }
/ |- B/ N! N: ^2 O( _+ L% K) p  status=1;
  return;
, E7 O3 d  j9 m: W$ G }
0 J7 F, u- _# d# d+ B5 L4 n2 M//检查是否为无界解
4 ?+ D- X+ m6 m for(i=1;i&lt;=num_x;i++) {
& \2 k8 Z/ @8 y% v  k_f=0;
3 v! c3 T* {# c9 C  if(sigma&gt;0) {
0 q  D# ?/ q" \   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
   if(k_f!=1) {status=0; return;}
, b8 L6 D1 `: A9 F  }
% p7 _5 c* D# q0 m0 c+ ] }</P>
! f/ H0 t+ t' y9 i<>//确定换入变量
for(i=1;i&lt;=num_x;i++) {
, `/ d: F9 A: F4 D1 y  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;
3 Q, ]* z' J6 `. i' u }//temp赋初值
  A( m8 r1 ^2 }  v7 i; C4 j 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)
5 N# K1 f9 D4 \) l   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
# v' N/ F5 t( r8 w* m) I    base_in=i;
9 T- U1 ~/ N7 J: j! E; |' ^    temp=sigma;
9 Z8 Z6 o0 {! F+ @+ Q  m* y   }
8 f5 A/ D. G, s& ` }</P>
7 T4 W4 L0 ~$ q0 Z& Z<>//确定换出变量
# Y$ c, v( s1 x/ e; y# T5 n for(i=1;i&lt;=num_st;i++)
3 [" X6 J& H3 h7 H; P  if(a[base_in]&gt;0) {
   temp=b/a[base_in]+1;
; I2 {9 ?" X. o% R3 F9 Y4 a) o   break;
- s" U% t6 b% Z7 g& ?% J; g9 h  C  }//temp赋初值
for(i=1;i&lt;=num_st;i++) {
# a. S6 K9 V& p' ~1 p" R  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
& \0 R6 f) r* R% v" O   for(j=1;j&lt;=num_ar;j++)
( p. }8 B0 n, I6 b; \: w( G# B    if(base==arti[j]) {
5 b1 c3 t" G) [1 _6 J+ r+ G     base_out=base;
     base_elem=i;
     temp=b/a[base_in];
     break;
    }
  }//人工变量优先换出
5 q) @( W) {; }6 q+ r" _* d; K  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
   base_out=base;
   base_elem=i;
' r8 W+ j, l- L% w7 \   temp=b/a[base_in];
( n& i9 Q. D' X/ Q0 x: y  }
}</P>
<> printf(" 基变量：");
# k  n3 g; A8 \2 H' _ for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
& M' w0 U% t8 _( j1 `9 m//基变量变换，进行新方程初始化后迭代
for(i=1;i&lt;=num_st;i++) {
" E8 k% T+ {, v6 \' l$ R( c$ N  if(base==base_out) base=base_in;
5 t2 U# d9 n& C5 G }
//初始化主元素行系数
3 m5 i- E1 R9 d/ w/ J/ _, B, s5 n value_be=a[base_elem][base_in];
b[base_elem]/=value_be;
& Z# ^* h' F# e0 e  [: H% J0 | for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
<> for(i=1;i&lt;=num_st;i++) {
) _6 Y* X% s7 q; g  if(i!=base_elem) {
   b-=b[base_elem]*a[base_in];
: O2 Y' {; U7 O: u! f$ [   value_be=a[base_in];
   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  }
}
  `( D0 ]' ^( P* M% O1 f4 O5 z stop++;
if(stop&gt;STP) {status=-2; return;}
3 _) M& C. \4 @, ^. u iterative();
}</P>
<>void output() {
int i,j;
4 u+ \) ]$ J; [+ U5 n3 h double X[MAX];
printf("\n结果如下：\n");
0 ?, \/ p/ J* s! c+ }/ d; r  S' U1 v' i printf("\nX=(");
for(i=1;i&lt;=num_x;i++) {
  for(j=1;j&lt;=num_st;j++)
  if(i==base[j]) {X=b[j];break;}
. @; Q- E- H4 _# m. h  else X=0;
  printf("%lg ",X);
}
printf(")");
7 F7 s) _# d9 T2 \/ S) M for(i=1;i&lt;=num_x;i++) max+=c*X;
if(ma_mi==1) printf("\nMax z= %lf\n",max);
$ M& _& |- l0 k7 L6 N! b9 I else printf("\nMin z= %lf\n",-max);
}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了