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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
#include&lt;math.h&gt;
#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])
9 Y" U6 P1 y! R: @& L{
7 i6 X4 ^7 A9 n# ^# w8 p$ I/ d2 z  int j;
9 S: H& W8 [8 p& ~  *c=0;
9 ~) e" \' C) L) e) L; e$ [  for(j=1;j&lt;=*is;j++)
    {
) E3 d. T" P/ g      if((*a)[*ir][j]-*c&gt;0)
    {
      *c=(*a)[*ir][j];
/ p# A( h) J! _, u4 Y      *j0=j;
6 j, _8 Y, j; T; _  W5 S, C    }
     }
8 ^9 e, ?. N' g( B}</P><>      /***************** 参数说明 **********************/
) h9 x- }& Z. z9 J9 s: ^1 Q0 F* N" o      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
; x- ]* W$ ^& B8 G! ~7 K      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
2 }( N" C2 {, ]% Y      /********* <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],
, L' v. C! P0 c8 c       int(*k)[],float(*x)[])
{
  int i,m1,mn,j0,i0,j;
  float c,g;
  m1=m+1;
' ]* M: f9 b6 O  mn=m+n;
  for(i=1;i&lt;=m;i++)
    (*k)=n+i;
* q( ^9 d9 P+ I2 w: U& {5 A# |leap1:
  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
0 c0 ?! y, U& n# W9 V2 w  else
     {
1 I6 V7 L& n: p* N. D* H; tleap2:   if(l1-50==0)
2 x- w9 h- F6 Q1 Y- W: z9 q, x      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
       else
) o4 B/ }* o0 ?9 P  j      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
: i# D' Z: W$ W5 k! ?8 s0 C6 I     }
   c-=1e-8;
9 ^5 C! y6 |6 |+ u5 }: P; e   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");
    return(-1);
      }
6 q* I+ J5 C; ^& }4 _9 W0 j* ~   if((c&lt;0)&amp;&amp;(l1==1))
8 Y: D/ q! L  }. ?       {
     l1=50;
3 }7 A* H- N1 [$ b     goto leap2;
3 @9 E' }, b6 _2 k9 O0 L# s& L, S       }
/ [9 e! l- W5 i4 d- @    if(c&lt;=0)
, m, K( h7 D4 P$ E       {
      for(i=1;i&lt;=m;i++)
+ i$ z$ V4 W$ _9 q& }         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
& c7 o9 U8 c5 k- ]      printf("\nF=%f\n",(*a)[m2][mn1]);
      return 0;
0 f, X- v2 X! J" c& r; c    }
2 f( H7 i' _3 d     c=1e8;
4 A# m  y1 i' K* v+ c     for(i=1;i&lt;=m;i++)
       {
1 H9 L. ^0 X; K0 I5 d& {     if((*a)[j0]&gt;1e-8)
& i- L# l) Y. P' z9 S7 S5 ]        {
- x+ T- b+ y- X' a; d% c8 p          g=(*a)[mn1]/(*a)[j0];
          if(g-c&lt;0);
         {
           c=g;
           i0=i;
          }
4 G* b% V: `  V" [* n% z         }
5 |. _. E$ Z  `2 x  w* E! V; b& x    }
) f- w4 I- m& n6 e      if(c==1e+8)
. [  I- S, H4 z8 Q8 B    {
      printf("\n\t*******Lp no solution********\n");
0 F( P8 I! y) f2 g- l      return -3;
7 }0 U7 Q6 o; ]& {4 g* z     }
       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
  T# V  K- c  w9 f8 H' F      {
. h0 `% P; \2 u" Z) F6 Z        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
           continue;
* |. E3 @, y5 ~# C  }( c         g=(*a)[i0][j]/(*a)[i0][j0];
           (*a)[i0][j]=g;
* u' m7 @- s  {9 y- g  f         for(i=1;i&lt;=m2;i++)
           {
- d2 D5 q( [( G# y' h3 F- X         if((i==i0)||(!(l1==1)&amp;(i==m1)))
            continue;
         (*a)[j]=(*a)[j]-(*a)[j0]*g;
           }
       }
       for(i=1;i&lt;=m2;i++)
7 z' P2 [2 u/ H) T/ u* R( z      (*a)[j0]=0;
       (*a)[i0][j0]=1;
       goto leap1;
4 x1 r/ J' n6 S! t; X}  
main()
% E; ]0 {$ q( E{
  float a[5][7]={{0,0,0, 0,  0,0,0},
: H; q, ?: }9 R& a                 {0,3,-4,3,  1,0,12},
                 {0,3,0, 6,  0,1,12},
+ t' d+ g) A0 r9 `6 D( A                 {0,0,0, 0,  0,0,0},
5 @3 M- G8 c+ l                 {0,69,0,144,0,0,300}};
1 x+ Y. y% H$ }  int k[3];
! }: @" q1 j2 ?" d  float x[3];
  clrscr();
  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
. y5 B4 R* B6 q' T9 a$ d6 y/ K7 ?}
9 D$ T/ w& L. u/ a! k8 T2 Z</P>

plgatc

<>/*************************************************************************
3 E' x4 ~) r1 m6 ~1 j* H2 Q* H1 A6 _                    单纯型法解线性规划问题（两阶段法）        
: l4 e0 }5 J5 r  X3 I; ?
0 T, k# R  y! C 编程环境：VC++6.0     
方程组输入说明：
变量非负，按提示输入相关参数。
$ O) c) k7 O3 m2 c6 _ *************************************************************************/
7 F1 Y  q% A1 v% W! F#include &lt;stdio.h&gt;
#include &lt;stdlib.h&gt;
# X3 A( }' a9 u- j2 o#define MAX 100
% z% B8 K, V6 B% M0 j$ ]#define STP 100</P>
8 V" l5 {  S# h* P: K- c# a! `<>int stop=1; //迭代记数变量
int status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
int step=1; //目前阶段</P>
<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
1 k6 t' c' L1 Q7 jint num_x; //变量个数
$ c" t. h1 f3 O) a' xint num_st; //约束方程数
int num_ar=0; //人工变量个数
( v9 _/ A) v: X8 C8 ^int arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
1 i9 b2 C9 u3 Y. E* C; l$ Z<>void create(); //建立方程组
void iterative(); //单纯型法迭代
- d! Z) ?. Q3 L- }. `$ o' hvoid output(); //输出结果
void banner(); //打印程序标题
void exchange(); //交换两阶段价值系数
. ?8 e" C- t, F' ?2 E0 q5 yvoid show(); //输出方程组</P>
<>void main() {
int i,j,k;
' t/ g% w* D1 G) C3 M banner();
& `$ c* t0 z, l# s- @! s create();
//保存原价值系数，转换为第一阶段价值系数
for(i=1;i&lt;=num_x;i++) {
  h& f$ ~1 x- ~1 a  k=0;
1 l% O" N8 Z$ K; v. F  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
  if(k==1) temp_c=-1;
- D1 d9 ]: c3 u: x1 |  else temp_c=0;
0 y0 a1 U, ?4 F }
; B& ?4 q5 @" V; P3 F exchange(c,temp_c);</P>
* Y9 |4 V$ y& j+ {% O<> printf("\n\n第一阶段问题为：\n\n");
* ?1 H! p3 F3 W show();
step++;
' P6 T. H8 ~- p0 ?1 k printf("\n\n按回车开始第一阶段迭代");
/ C* b# y0 z. O5 c; I getchar();
! h$ \5 I: f5 [ getchar();
# y; R! C; x# f5 X9 X iterative();
) a5 w* x& L4 U. n* I0 ?9 R if(status==-2) {
' A( E; w2 Z/ |. L& g/ Z  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
- O% b0 d  O; {* s/ [9 b; D" J3 T4 J  getchar();
  exit(0);
) ]! t" G7 r$ \, E  U3 W, p }
7 O/ {1 g! z9 t' R2 ]6 T& n5 v output();</P>
<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
' I0 {! F7 ?+ a. N; D  puts("\n按回车结束");
  getchar();
  exit(0);
}</P>
* k2 v6 ^- O1 N0 y& k, e# g! U9 V<> //转换为第二阶段价值系数
! w6 A. u: ]  }- j3 X exchange(c,temp_c);
//把人工变量列全设为0
" z1 k1 v# m, {+ o for(i=1;i&lt;=num_ar;i++) {
  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
}</P>
2 U3 m/ {% p0 \' J<> puts("\n\n第二阶段问题为：\n\n");
show();
puts("\n\n按回车开始第二阶段迭代");
getchar();
& ~7 \1 \0 ~1 H; r' D iterative();
switch(status) {
case 1:
  output();
) ^* i! w- b3 g- Y( h  puts("\n\n原问题有唯一最优解。\n");
2 T- m9 @. r" ?' o# b- [) x. W0 }; u  puts("\n按回车结束");
  getchar();
" U5 r1 t' a- q+ \  exit(0);
3 M( B0 j6 ]' o case 0:
  puts("\n\n原问题为无界解。\n");
  puts("\n按回车结束");
) Y+ ?# V% o* K2 D3 O8 g6 m% e  getchar();
  exit(0);
6 Z) @6 x& ^3 _) X" ~ case -1:
  output();
  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
  getchar();
: w" I# d7 x  K! j* x  ?& f1 i  exit(0);
case -2:
% ^. a* U. {+ e+ b) q  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
; f4 ?* v. r! ]4 M5 ^  M$ b }//switch
}</P>
  S8 j  Z1 Q6 [  K) q' R. |<>void banner() {
5 X9 J) S1 K% L+ g& U printf("\t\t****************************************\n");
printf("\t\t         单纯型法解线性规划问题\n");
+ m) d* C: ]* Q2 b7 ~ printf("\t\t                         作者：Thunder\n");
2 O' U! a8 W  j) X7 Q$ P5 t printf("\t\t****************************************\n");
printf("\n");
}</P>
<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
, v) m8 ^$ [" D//为1的不显示系数1，-1系数只显示负号
int i,j,k;
switch(step) {
case 1:
% o& F6 }2 _  [( ~/ x  printf("min z= ");
  printf("x[%d]",arti[1]);
  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
) }  a  [3 _3 J8 r case 2:
, c2 [& M7 ^6 G$ [ printf("max z= ");
4 w3 X7 t$ D& n% D printf("%lg x[%d]",c[1],1);
. N# Q5 t; j! m" {/ r for(i=2;i&lt;=num_x;i++) {
8 l3 _* `# j8 M8 j! O4 k3 w  if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
# `6 B- K8 b. I8 V8 ]  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
  else printf(" %lg x[%d]",c,i);
4 S) ]1 }, f; w0 b# ]$ T0 b! c }
break;
. J$ P. u6 N8 j. Z" x }</P>
  R. ]- I( _9 h8 V+ Y0 E% `( r. i<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
( g+ r' E5 U0 W. F  k=0;
  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
/ d) l. X2 z9 Y2 `( A    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);
/ W% ?* m! K+ O) X+ p    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);
6 r% z% ]( `  i    else printf(" %lg x[%d]",a[j],j);
    k=1;
: V( w& J( e4 F' H2 m$ f   }
. \" |& {0 s2 \9 Q- p+ X5 x0 |6 x  }
  printf(" == %lg\n",b);
}
; `0 q  @( A; C6 K+ N, U; A printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
5 e9 t" ~$ c& u6 f<>void exchange() {
int i;
double temp[MAX];
' ]- `4 V! p5 I7 ]& Q% x for(i=1;i&lt;=num_x;i++) {
( _* V4 U2 {) A" m/ `  temp=temp_c;
  temp_c=c;
  c=temp;
8 u0 Z3 h' O" j }
}</P>
<>void create() {
//输入方程组系数，每个方程输完后回显确认
; M' U# {- ^' P+ S# m' B. {4 ` int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
char confirm;

while(1) {
  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  scanf("%d",&amp;ma_mi);
  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
' x8 R2 m5 G, I3 h2 V6 e# Z) o. H( e  else break;
- _4 g8 C5 r6 I# n3 \& e8 | }
  q. X+ V: B( t9 P5 ]) Z
while(1) {
) w/ R) D: o1 {  printf("指定变量个数：");
  scanf("%d",&amp;num_x);
: b: _& A" P- ^  q9 X/ i  printf("输入价值系数c1-c%d:\n",num_x);
  for(i=1;i&lt;=num_x;i++) {   
   printf("c%d=",i);
; n) h0 g( {8 O  v2 C1 H   scanf("%lf",&amp;c);
% A+ q$ a0 H- B, W9 ?  }
8 g6 |6 q9 u. Y9 E$ e4 S0 ?  if(ma_mi==1) printf("max z= ");
1 e9 [" V5 l5 s& _  else printf("min z= ");
  printf("%lg x[%d]",c[1],1);
  for(i=2;i&lt;=num_x;i++) {
; D' F9 N" Q$ }% w0 Y# q. L, r: |   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
   else printf("  %lg x[%d]",c,i);
  }
3 [5 o- s/ H/ I1 R! a" i: ? printf("\n正确吗？：(y/n)");
. E5 b  m# A* B getchar();
confirm=getchar();
' r* y1 a' W  `$ T% E if (confirm=='y') break;
else if(confirm=='n') continue;
; r! G: N8 j6 ?+ h/ V( o3 b. j }</P>
<> printf("输入约束方程组个数：");
! b& b) R5 y: n) M8 L4 [ scanf("%d",&amp;num_st);
for(i=1;i&lt;=num_st;i++) {
+ V; i4 z7 Z9 ?  printf("st.%d:\n",i);
0 D5 F5 F: J* ?9 i" N  while(1) {
   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
1 C" R, _% s1 }& E( j   scanf("%d",&amp;re_st);
   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
, B4 c8 P1 c) R, \   else break;
  }
  printf("输入技术系数：\n");
& L7 R& V( S" m! }  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
  }
  printf("输入资源拥有量：\nb%d=",i);
  scanf("%lf",&amp;b);
6 T) w3 m/ h- k+ C4 O  
  printf("st.%i:\n",i);
  printf("%lg x[%d]",a[1],1);
: Y9 d) |- y6 g% `  for(j=2;j&lt;=num_x;j++) {
   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
. ^8 \+ p& ^$ {7 ~, |   else printf(" %lg x[%d]",a[j],j);
  }
+ O9 K9 G% {# X  switch(re_st) {
; D3 W. Q" I: H- ^   case 1: printf(" == %lg",b); break;
/ _1 U8 [( }0 r- g3 F   case 2: printf(" &gt;= %lg",b); break;
   case 3: printf(" &lt;= %lg",b); break;
  }</P>
<>  while(1) {
   printf("\n正确吗？(y/n)");
   getchar();
   confirm=getchar();
   if (confirm=='y') break;
6 r5 f& N5 p% |9 ]; f& f1 J  x! l( U* q6 `   else if(confirm=='n') {i-=1; break;}
) G7 ?4 f& ^8 P$ g  T9 o$ V  }
) i4 |$ I( @* g5 y, j }</P>
<>//显示输入的方程组
; v% X1 j- W" v# ] printf("\n原问题为：\n\n");
% }: G: c2 _' }( s6 q' }3 | if(ma_mi==1) printf("max z= ");
else printf("min z= ");
1 A3 {+ V1 W. R# Y printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
% Y4 h( s& V" `. F! z  if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
$ B( k2 f- ~* V" f9 F0 u+ B  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
  else printf(" %lg x[%d]",c,i);
}</P>
* G& j- M1 r4 O<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
8 e" m2 l' _/ M9 }0 e  k=0;
  for(j=1;j&lt;=num_x;j++) {
7 w" ^, A* k2 X4 N7 W! Q) s. F   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
( H( c- X  _9 ?    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
' @" S( p# Q( ^$ \4 t) n4 h    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);
( o1 c" w* B- r8 Q" Q. Y    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;
   }
1 d. ~2 A6 G7 u5 \  N0 |  }
5 {! |9 F) C! B  switch(re_st) {
   case 1:
    printf(" == %lg\n",b);
; W) I! g! H' V* {. W8 ~) e    break;
6 c1 p' _( B- O) T+ o7 B/ T   case 2:
& j2 s0 \, {) W2 p. d3 u+ a    printf(" &gt;= %lg\n",b);
& ?2 K8 y+ d! P    break;
& w" S5 O9 u0 \4 h8 f( h   case 3:
0 \5 g4 Y9 P. \4 S# d. c    printf(" &lt;= %lg\n",b);
- k9 M# I; c: f    break;
- ^- ]2 l" Z' v( x  }
}
$ d3 c2 k' r+ M) M+ p" W 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) {
/ J3 i! _- `  A6 _$ U/ z  case 1:
  case 3:
1 b4 W4 i1 X1 g" |0 T   num_x+=1;
& Y4 U' k# E( t   break;
1 W# d& E  R! H) m! z2 s% U  g  case 2:
- H; B  z1 ^0 d8 s   num_x+=2;
   break;
2 E/ r+ v, q; A# ~* [6 i$ b  }
}</P>
# k: Y' q: L( A4 R( [8 N8 Z* P<>//化为标准形式
2 p- m; }7 V$ {8 M if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
' d( D% D: D7 o' E, ?+ P for(i=1;i&lt;=num_st;i++) {
, g1 @& Z' l5 x7 B1 A1 ~  i1 O  switch(re_st) {
; h  R- n' \! g6 X3 D9 X   case 1:
    num_addv++;
& {( ?. P" I! n# l1 T    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++)
     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
; S' A* |% B( i5 @" g. I8 z    break;
   case 2:
: u; p* {6 i' v0 K* s    num_addv++;
    c[tnum_x+num_addv]=0;
+ y' f9 B1 S# X1 ]( X    num_addv++;
- ^; n7 u, \2 S% R$ }: j- H    num_ba++;
( {! ?% ^. B4 E+ D5 X3 ?8 w    num_ar++;
    c[tnum_x+num_addv]=0;
    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
7 t/ O1 G7 d! b* d" H$ L) J    for(j=tnum_x+1;j&lt;=num_x;j++)
     if(j==tnum_x+num_addv-1) a[tnum_x+num_addv-1]=-1;
# W0 g$ N' u" g: N     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
' B1 Q" c2 ]! R1 t: Y9 Y6 Z! E- D3 r    break;
   case 3:
6 `2 [# \3 [: b% k    num_addv++;
    num_ba++;
+ H# Z# C7 c5 m) l9 U4 }: m; C    c[tnum_x+num_addv]=0;
    base[num_ba]=tnum_x+num_addv;
5 U8 ]* I) n4 f    for(j=tnum_x+1;j&lt;=num_x;j++)
8 l# a) B* Y# |     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
& @2 ^$ p9 \* t8 R6 i3 v     else a[j]=0;
    break;
* ~7 i, N. F. B/ f. c  }//switch
}//增加松弛变量、剩余变量、人工变量、确定基变量</P>
<>//显示标准化后的方程组
printf("\n化为标准形式后：\n\n");
+ m5 K! Y, i+ f; i; g. n3 _ if(ma_mi==1) printf("max z= ");
else printf("max z'= ");
printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
; y$ m( m0 Z- E9 P4 H  k=0;
  for(j=1;j&lt;=num_ar;j++)
4 u" U, ~5 a7 t) I2 M  if(i==arti[j]) k=1;
  if(k==1) printf(" -M x[%d]",i);
+ x5 [% ^/ f4 Q* ?4 L  r; P* e  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);
  else printf(" %lg x[%d]",c,i);
}</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
( k# B! {2 @4 [2 p: z4 D4 {  k=0;
  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
% M& g$ ~! l' P+ U8 t" B( G    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);
    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);
. m: X) z% K3 e    else printf(" %lg x[%d]",a[j],j);
    k=1;
( S* e) b2 V& V& [. ]( @, K5 z# f0 I) _   }
3 r; A, q: \* J' {4 h! p6 z3 O6 p  }
/ v) @3 A+ q' Z: S% m' `  printf(" == %lg\n",b);
+ k( h0 d0 Y7 L/ D, J3 H2 K" j7 N/ ~  ~ }
2 |( W" Z+ L, U9 z- B printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
8 o: I! ?' O% ]; x<>void iterative() {
int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
% k8 d- h: t, j; M1 O int base_elem;
: Y/ i1 F& P( q* _6 e% M6 | int base_out,base_in;
double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
8 I9 r3 t' j* ~" ?- A  u' k  c2 P6 F* [<> printf("\n\n第%d次迭代：\n\n",stop);
for(i=1;i&lt;=num_st;i++) {
  printf("c%d=%lg\t",base,c[base]);
7 O$ @! D- J: U  printf("b%d=%lg\t",i,b);</P>
; q7 G; W9 }- q$ t) H, D6 S, A<>  switch(step) {
   case 1:
    for(j=1;j&lt;=num_x;j++)
; D" e- t! _& O! c* Q0 g/ t    {
; R. f& ?+ u0 R$ ~     printf("a[%d][%d]=%lg\t",i,j,a[j]);
) h5 P" v  o5 t. [4 @# F    }
, w; W* C, ~) v& f, T    printf("\n");
    break;
5 k( I) u0 M8 s( O- M   case 2:
0 V  S2 N7 M0 I/ }* S7 V    for(j=1;j&lt;=num_x;j++) {
. a3 w0 {3 T3 |. G, @5 S     k_a=0;
     for(l=1;l&lt;=num_ar;l++) if(j==arti[l])k_a=1;
- S* R5 \' s4 j     if(k_a!=1) printf("a[%d][%d]=%lg\t",i,j,a[j]);
0 {( w: w0 A# z0 S8 `+ ~. _    }
0 K. L8 B5 a3 ?5 e6 z    printf("\n");
) I9 i4 F  J  D( Z    break;
  K* b' q% Q- N1 k  }
1 [  T' J8 T, n5 d' A9 h }
: T% A5 c& _# a. B( [. j5 ]0 ]+ D//求检验数sigma
for(i=1;i&lt;=num_x;i++) {
  sigma=c;
( ~7 f. j5 q; u9 z% N8 o: ~  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) {
" ?$ T+ o$ n/ r6 B   case 1:
    printf("sigma[%d]=%lg\t",i,sigma);
    break;
   case 2:
2 c# p0 k9 P' T: K0 y    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);
8 ]( x+ Y& ?4 Y0 V8 p/ {    break;
; |- M9 V  N$ @8 B8 I* p9 \  }
- z2 O5 y0 _: Y+ @ }
putchar('\n');
//检验检验数sigma是否全小于等于0
% F; D3 j. c6 q: L% |. J k=0;
for(i=1;i&lt;=num_x;i++) {
4 W& ~) f1 i- |- w& N' m# k4 ?  if(sigma&gt;0)
, a& S4 y  V2 ^* T) a  k=1;
8 @' t9 i; R5 C4 p, M }
if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
  for(i=1;i&lt;=num_x;i++) {
# e$ ^4 B( z7 |) a! v   k_f=k_a=0;
   for(j=1;j&lt;=num_ar;j++)
   if(i==arti[j]) k_a=1;
7 ]& _& H! j8 E: x   if(sigma==0&amp;&amp;k_a!=1) {
1 j9 O6 u- m7 p* o    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
! m8 A- d& o, U- n( m  M: h   if(k_f==0) {status=-1; return;}
   }
  }
9 p$ E- {5 c  ?* r  status=1;
  return;
}
//检查是否为无界解
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;
  r5 r: x% I2 w6 k/ r   if(k_f!=1) {status=0; return;}
  }
}</P>
. G0 o" z, Q- g; s, w) O0 e<>//确定换入变量
for(i=1;i&lt;=num_x;i++) {
& m- _6 _+ D9 ^- N' @4 S9 B; X  k=0;
! I9 ]9 G, `% @7 V  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
( I2 i) i0 h; f: p2 C8 d  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
1 _/ j1 h* R" x. {, k4 b1 q }//temp赋初值
5 G4 Z, a6 M/ p% X+ v1 g for(i=1;i&lt;=num_x;i++) {
  k=0;
( L$ d3 e- g9 p  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
+ L# `# m9 y# A; T1 [1 v  if(k==0)
   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
    base_in=i;
& @% J( a6 X# o" x: K0 a5 Q% |$ H- T    temp=sigma;
9 Q$ M$ ?  r4 L- i/ O$ K( K  {4 Z   }
}</P>
$ [, L- d0 G! T/ e& P4 n1 J<>//确定换出变量
for(i=1;i&lt;=num_st;i++)
+ k5 z: P& B2 G7 B2 g0 A  if(a[base_in]&gt;0) {
   temp=b/a[base_in]+1;
; |  R2 i+ }9 n9 f2 i% y2 r6 }   break;
5 c9 P8 B$ O. V& Z! P5 n% D  }//temp赋初值
) ~; y# c- T* O$ ?3 x2 ?5 ?" D for(i=1;i&lt;=num_st;i++) {
: h, ]+ O! i; \+ R- x. L  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
   for(j=1;j&lt;=num_ar;j++)
6 S9 w+ W6 \# Y    if(base==arti[j]) {
     base_out=base;
     base_elem=i;
( _$ ^2 Z! [9 |     temp=b/a[base_in];
3 M! C7 y' [: p$ v* L     break;
/ i% R) F4 e+ {' e    }
- G- o* Z9 }5 z. b( A' e; z  }//人工变量优先换出
  V% w+ R# S8 ]: U  l  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];
( Z! |" T( e$ z# y5 {5 H# `  }
}</P>
<> printf(" 基变量：");
for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
4 ]& d# H7 h8 X+ }. i4 m6 \//基变量变换，进行新方程初始化后迭代
3 Y0 Q2 Z+ o( H' {# I" }- r for(i=1;i&lt;=num_st;i++) {
. y( K0 ]% r6 `. I: z  if(base==base_out) base=base_in;
}
; p! j' F: v% m: b# m* M: B//初始化主元素行系数
value_be=a[base_elem][base_in];
b[base_elem]/=value_be;
$ e  o0 `. {. W5 C% Y, H4 N' { for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
<> for(i=1;i&lt;=num_st;i++) {
: S& J. E7 i& A& z  if(i!=base_elem) {
2 k. m+ v# [* u% w2 I# `- z* |   b-=b[base_elem]*a[base_in];
" {7 w* I# |: r" }  E   value_be=a[base_in];
- x$ R0 F3 b3 {" i+ i; r   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
, R3 Y0 d4 i( m7 g% m' d  }
& G3 |5 H) T: a) N" s7 U7 t* _2 M }
$ i4 Z+ U, z5 f, y! ]0 t4 D' p stop++;
if(stop&gt;STP) {status=-2; return;}
iterative();
}</P>
<>void output() {
int i,j;
double X[MAX];
- }: A0 u1 G3 N6 I4 I- G  l# [ printf("\n结果如下：\n");
printf("\nX=(");
) w( F% H  X2 r for(i=1;i&lt;=num_x;i++) {
' a# S9 n1 R# C. ~/ Y( Y  for(j=1;j&lt;=num_st;j++)
  t5 d! \* i  ?% z* \5 }  if(i==base[j]) {X=b[j];break;}
# T& S5 I# t) n3 E, [  else X=0;
  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);
7 \5 I' f1 {) r: n  ^! y& [* i9 [ else printf("\nMin z= %lf\n",-max);
}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了