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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
#include&lt;math.h&gt;
#define X 5
  Q* `; s6 Q9 k5 _#define Y 7
void xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
{
' Y; g3 W2 `) K4 E5 [7 O  int j;
% `) J7 G* ]* [% L  y! a  *c=0;
5 B6 K, J7 s7 h4 a  for(j=1;j&lt;=*is;j++)
    {
      if((*a)[*ir][j]-*c&gt;0)
    {
& \' W1 n9 Y2 j& L      *c=(*a)[*ir][j];
      *j0=j;
2 F: W$ }0 N; j* S6 {  b8 U& c    }
     }
# a. |3 {  s1 b' r6 K8 @* w}</P><>      /***************** 参数说明 **********************/
      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
# f' o6 Z7 g4 b7 u+ c3 l9 V      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
2 }0 M4 ^, V, @& p, ^0 S; q& Y      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
- W6 d) C* ?* q* ^      /********* <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],
% i; ~% J) P6 a& |       int(*k)[],float(*x)[])
{
  int i,m1,mn,j0,i0,j;
  float c,g;
5 o* B6 A3 I/ E" z8 t  m1=m+1;
  mn=m+n;
  f# ?' l1 g; \5 Z5 v0 a  for(i=1;i&lt;=m;i++)
3 b1 @& x( l# g! A9 T0 k1 J9 E) P% k7 n    (*k)=n+i;
leap1:
$ J, E, G6 q( v' @& K+ J- k  if(l1-1==0)
! d- w7 Q& M& `6 t4 |     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  else
     {
leap2:   if(l1-50==0)
      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
       else
      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
     }
9 l* O) [* l# n! ~8 \/ U0 J   c-=1e-8;
, j0 l: p  \) F* G1 I   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
. I1 P  W/ m+ |! G9 Q      {
% _) D8 B* I% B4 P  ?/ P    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
    return(-1);
6 L$ O0 x: A, T- l. i) W/ r      }
7 q% b0 C" ?' v$ ?- V2 H4 W: P- U   if((c&lt;0)&amp;&amp;(l1==1))
& g7 ]1 z% s! |       {
% z) H- G: R" |& \& `8 V. L     l1=50;
     goto leap2;
0 X  i9 b( O) r       }
4 y& ~* k2 g& E1 I    if(c&lt;=0)
       {
* V0 @( Y$ Z) h4 o5 q1 s/ @2 J      for(i=1;i&lt;=m;i++)
4 n) ]& t: H3 J         (*x)=(*a)[mn1];
2 L5 i, A  v: y4 V% g         printf("\n\t********Optimal solution********\n");
      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
/ B& Y+ @: ]  w: U( O& f      printf("\nF=%f\n",(*a)[m2][mn1]);
      return 0;
4 I' {9 G4 @% |1 U8 g& M    }
     c=1e8;
% x) N! \; I5 S& B     for(i=1;i&lt;=m;i++)
& H, d9 l- w: x  S0 ~9 H/ b9 c3 G       {
$ {3 }: a. F6 h# [. \! T     if((*a)[j0]&gt;1e-8)
        {
% ~5 l! Z' X0 P+ }5 V2 b          g=(*a)[mn1]/(*a)[j0];
% a. a: B$ q; M* B$ C" D  k          if(g-c&lt;0);
         {
" [6 b$ H8 R7 B5 c2 I9 O           c=g;
           i0=i;
# L, u$ ~; @- l+ [) h          }
+ }3 p0 S, J/ v& p1 N& G* v         }
    }
      if(c==1e+8)
    {
4 c9 t3 I" s4 s% N- ?" ^      printf("\n\t*******Lp no solution********\n");
" Q, T$ y- e" X  E      return -3;
% t# l* |3 Y0 V. u0 x     }
9 K9 o6 @; h  o" @! J       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
      {
2 ^" V- N0 g7 w3 x        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
9 @1 ?5 F; z6 u9 y           continue;
         g=(*a)[i0][j]/(*a)[i0][j0];
           (*a)[i0][j]=g;
$ W  p0 j+ L9 w, E) C/ o$ G6 Z# d         for(i=1;i&lt;=m2;i++)
+ j% m6 l* L- u7 f3 B9 ]# U           {
         if((i==i0)||(!(l1==1)&amp;(i==m1)))
& e! L9 P* ^- D7 [            continue;
; B7 J* E0 F( J4 W$ u$ N         (*a)[j]=(*a)[j]-(*a)[j0]*g;
           }
       }
( R8 c4 w1 q$ v1 Y# ?       for(i=1;i&lt;=m2;i++)
2 @+ U! H/ ~! V+ w$ R5 W* l      (*a)[j0]=0;
9 l/ E) a: }) q4 ~, b       (*a)[i0][j0]=1;
! h7 U2 P1 v' w& Z% }8 J       goto leap1;
* I( E# I5 Q& z% }9 Y3 [}  
main()
{
% F4 s- z7 I& K- k. O# {  float a[5][7]={{0,0,0, 0,  0,0,0},
0 ]' q- e7 \, ~/ f( z                 {0,3,-4,3,  1,0,12},
                 {0,3,0, 6,  0,1,12},
! y& f& a7 `. n, _. X9 I7 f                 {0,0,0, 0,  0,0,0},
                 {0,69,0,144,0,0,300}};
  int k[3];
  float x[3];
  clrscr();
  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
}
</P>

plgatc

<>/*************************************************************************
6 R- G( U1 y  r% p* I                    单纯型法解线性规划问题（两阶段法）        

编程环境：VC++6.0     
6 B4 X9 B" \# F: M. b 方程组输入说明：
0 y# _) d" T& }0 F6 s: V 变量非负，按提示输入相关参数。
* g& S5 f6 j: O) l9 U2 Z" {3 X; J; u* J& v *************************************************************************/
7 o8 v5 T, I9 E* p5 E  A#include &lt;stdio.h&gt;
& e0 A4 k. y7 b( F#include &lt;stdlib.h&gt;
#define MAX 100
#define STP 100</P>
<>int stop=1; //迭代记数变量
int status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
) J- F! a$ }  G+ Sint step=1; //目前阶段</P>
<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
5 A3 U% x2 n* y0 j) ^! Aint num_x; //变量个数
1 ]+ o% \& G9 J2 Hint num_st; //约束方程数
6 \, o  I3 f( f8 U$ Gint num_ar=0; //人工变量个数
, _, D- }+ X* K4 @) b' ^int arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
9 D% v  ^. Q' w; `6 n% O8 K$ \- ]<>void create(); //建立方程组
void iterative(); //单纯型法迭代
void output(); //输出结果
+ N$ [- j2 O$ `2 q0 k+ Vvoid banner(); //打印程序标题
/ d! ?3 C$ ]- Rvoid exchange(); //交换两阶段价值系数
void show(); //输出方程组</P>
<>void main() {
+ I' }2 T1 G: @) }2 L int i,j,k;
) _/ ~; g5 z+ x5 c; O banner();
+ _, V7 y8 D* Y. l- p: ?# O: ?& j& D create();
% b0 ~/ W7 |5 Z2 v5 n //保存原价值系数，转换为第一阶段价值系数
for(i=1;i&lt;=num_x;i++) {
% @! K' M4 u* g3 D  k=0;
! G* j5 k8 V2 W  p) }# H, Q5 T  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
  if(k==1) temp_c=-1;
0 }/ r& D% u, p  else temp_c=0;
}
exchange(c,temp_c);</P>
<> printf("\n\n第一阶段问题为：\n\n");
show();
step++;
- T1 }' T3 {1 W& Y9 h printf("\n\n按回车开始第一阶段迭代");
getchar();
6 C, u, u. S7 z: y2 J9 n getchar();
iterative();
if(status==-2) {
  puts("迭代超过限制次数强行终止！\n");
) E2 f. S0 D, b3 X! Q  puts("\n按回车结束");
+ ~0 W3 d1 ?# ?  e  getchar();
  exit(0);
}
output();</P>
. _9 t- _& a* }- V<> if(max!=0) {
- M; i$ ^3 L& Q/ w, h+ ^  puts("\n\n原问题无可行解。\n");
8 v1 t9 v4 C/ x3 h; {& m$ D  \  puts("\n按回车结束");
$ w6 w4 I2 d% p2 m  getchar();
  exit(0);
}</P>
<> //转换为第二阶段价值系数
exchange(c,temp_c);
//把人工变量列全设为0
9 S4 w+ q- E% C; L for(i=1;i&lt;=num_ar;i++) {
  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
) ^% k# w" f1 |5 [% z }</P>
<> puts("\n\n第二阶段问题为：\n\n");
" A' q0 p1 Z" {6 P8 I show();
puts("\n\n按回车开始第二阶段迭代");
$ q4 \, s& A5 c0 ?: v3 _2 ~ getchar();
; Z; [& q( h- V. j, H iterative();
switch(status) {
case 1:
  output();
  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
$ x( x! @& ~8 ~1 t  getchar();
' M9 [: Q  U5 l4 S* O  exit(0);
case 0:
9 a9 q! N+ u( P7 t) W7 P$ F3 c  puts("\n\n原问题为无界解。\n");
, f/ r  C# |3 A2 F4 B  puts("\n按回车结束");
  getchar();
  exit(0);
case -1:
  output();
0 ~' S( t! G. k# }& U$ w2 u" ?  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
% @- J+ k2 O5 I& D& N* ?7 }: C0 l case -2:
2 F1 K% }0 p; B; z  puts("迭代超过限制次数强行终止！\n");
/ D1 i* ^4 b! l6 Y( |  puts("\n按回车结束");
( p, U3 z( z6 t+ f) @. x! Z& |% c& R  getchar();
& j! J# U3 K6 |# c( z  exit(0);
}//switch
}</P>
0 X6 E- g# r) @$ j0 h  Q<>void banner() {
/ a6 t$ v8 s! @$ ]) V5 V5 {6 o printf("\t\t****************************************\n");
printf("\t\t         单纯型法解线性规划问题\n");
printf("\t\t                         作者：Thunder\n");
printf("\t\t****************************************\n");
5 J  L  L5 T+ I- ?0 W' M2 ~$ p printf("\n");
; D( Y; D  P+ h! G' r) c}</P>
" S: T" B0 Q) i) [0 u<>void show() {
. s! ~( y! E" h( M% F# J//对方程组以自然的格式输出，系数为零的x不显示
3 x9 R9 l" f7 I! ]/ x//为1的不显示系数1，-1系数只显示负号
int i,j,k;
switch(step) {
- W7 A+ e$ \0 b7 o+ y* R3 g case 1:
  printf("min z= ");
3 ~! w& R; @$ C1 [& j  printf("x[%d]",arti[1]);
  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
4 i# G3 B, M" V case 2:
/ y) P& Z3 I) @0 n; ~$ y printf("max z= ");
" \9 e1 l' l- A; K" P) I0 Y2 Y printf("%lg x[%d]",c[1],1);
- c  V- i1 w) F/ ]4 a for(i=2;i&lt;=num_x;i++) {
  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);
) D  F" G& Q6 j6 e2 _ }
9 A) z9 L6 t9 e3 O( o1 v3 ^ break;
}</P>
4 `( Y/ E* I  P( _3 `% D# I<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
8 C; i( M5 F3 Q% l  k=0;
  for(j=1;j&lt;=num_x;j++) {
& d, v) v) u6 ]. U- G  u5 s   if(a[j]!=0) {
$ q* Z4 M1 X4 p# q) Q. x    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);
7 {8 ~0 N$ b# ?" S; B$ K: w6 j    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
+ R; H  C& b. _1 i    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);
! C6 `9 F: G: R% c% q    k=1;
( v, ~  X1 O, @; Y   }
  }
  printf(" == %lg\n",b);
}
printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
<>void exchange() {
/ E- \' {- A9 C" K$ h int i;
+ W: b" k/ w& s  V double temp[MAX];
& `: l2 D$ @( O( p6 {' S) | for(i=1;i&lt;=num_x;i++) {
  temp=temp_c;
  temp_c=c;
  c=temp;
}
}</P>
<>void create() {
0 m+ y7 l9 W& n. Q2 H* f  W' A//输入方程组系数，每个方程输完后回显确认
: v8 ]. W/ A+ B: k int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
char confirm;

7 }" U  V& y# |8 _2 A; g' ` while(1) {
  printf("请选择：1、求最大值，2、求最小值：(1/2)");
0 w/ N+ S; v1 x+ G$ G1 p2 i  scanf("%d",&amp;ma_mi);
* \6 F3 d* T7 h3 P# a& ^+ ~  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
  else break;
}

while(1) {
  printf("指定变量个数：");
  scanf("%d",&amp;num_x);
  printf("输入价值系数c1-c%d:\n",num_x);
3 V' g5 H0 V/ ~$ p  for(i=1;i&lt;=num_x;i++) {   
   printf("c%d=",i);
   scanf("%lf",&amp;c);
  }
- a/ C2 V* ~2 T  if(ma_mi==1) printf("max z= ");
8 o' b, Y) O! T( z+ \  else printf("min z= ");
  printf("%lg x[%d]",c[1],1);
  for(i=2;i&lt;=num_x;i++) {
) w/ \# d+ A0 m: x4 V   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
' H, Q) L' @* p   else printf("  %lg x[%d]",c,i);
4 R) S+ E* ?3 A6 i0 r  }
. q8 G' b$ V/ x2 u2 r) z$ [6 a printf("\n正确吗？：(y/n)");
9 X2 _# N3 C; u$ G4 Y getchar();
confirm=getchar();
if (confirm=='y') break;
6 o2 Q% K8 a! Z& h else if(confirm=='n') continue;
}</P>
$ ^4 k1 Q/ E2 b" C3 e5 F<> printf("输入约束方程组个数：");
1 {1 W. O* V+ g9 t  S' B scanf("%d",&amp;num_st);
5 n! n( b. i. D. E7 v* M" ~ for(i=1;i&lt;=num_st;i++) {
  printf("st.%d:\n",i);
. @7 @7 @1 I- j) Z, Y# z  while(1) {
! Z: e1 Z4 [; e  h$ C* P   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("输入错误，请重新选择。");
7 _+ F( X5 e$ {, V0 U$ \   else break;
: S' _9 f* S- P+ C4 I3 c  }
7 `5 W8 Z9 U, B: m( K  printf("输入技术系数：\n");
5 j; c  O* R1 \  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
! \' d; M4 V7 A- s6 h9 @  }
  printf("输入资源拥有量：\nb%d=",i);
: n# f7 r8 j1 G  scanf("%lf",&amp;b);
' l1 C/ @$ J  ^1 V, h+ \  
  printf("st.%i:\n",i);
) Z  ~9 O+ R$ ~, |; K3 q5 @3 a  printf("%lg x[%d]",a[1],1);
+ W! V& w$ D6 Y8 B! `  for(j=2;j&lt;=num_x;j++) {
4 r5 y' w- a% s   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
7 ~0 K! {2 t& h/ F   else printf(" %lg x[%d]",a[j],j);
* X3 ?; u1 z0 d/ a/ r  }
  switch(re_st) {
' ?" M1 _9 V& u   case 1: printf(" == %lg",b); break;
   case 2: printf(" &gt;= %lg",b); break;
   case 3: printf(" &lt;= %lg",b); break;
  }</P>
1 F$ N6 Z. y' g<>  while(1) {
   printf("\n正确吗？(y/n)");
% }# `1 {  z  ^8 j  g   getchar();
# R4 m" L5 C5 V+ n# D6 W% A2 y5 _: s   confirm=getchar();
; O. w) ?9 Q. L- W, {$ x   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
  }
}</P>
: ~0 D5 _2 y% ~: |. q, L( P- c6 e<>//显示输入的方程组
' U. C8 |( f; h. { printf("\n原问题为：\n\n");
if(ma_mi==1) printf("max z= ");
else printf("min z= ");
$ }0 L$ n2 v: F% Y, T/ ?/ h 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);
  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
5 z7 I% x; X8 p; K3 F/ V  else printf(" %lg x[%d]",c,i);
: m* m$ M5 T4 q" N: ?3 M }</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
/ h9 V7 O0 P" g& E! k  E  k=0;
  for(j=1;j&lt;=num_x;j++) {
# g: B5 i/ H  K0 Y* V/ P5 ?   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
+ S4 `9 Q$ |, ~/ ]" Z    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);
; l; X/ g9 Z4 J1 Y' V* \    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
& Y3 R  `8 q$ h    else printf(" %lg x[%d]",a[j],j);
; _9 A6 C9 A% S# K    k=1;
) P  A3 S7 \# Z   }
  }
* c5 {; [5 P1 A; B3 D  switch(re_st) {
   case 1:
# _4 M. x( a' A2 k) p9 @. k2 h    printf(" == %lg\n",b);
    break;
! @+ P1 l' H/ U# P7 c0 M   case 2:
6 G' h3 E$ x/ z( J0 c; t9 \" D    printf(" &gt;= %lg\n",b);
    break;
   case 3:
    printf(" &lt;= %lg\n",b);
5 ?) l" Q% g7 o" G    break;
  }
" j9 P' B. a7 w. ?$ z- r }
printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
<> tnum_x=num_x;
; o( U) ~0 T* J2 @+ ~  Z for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
* ^  e' Q( g1 r/ K' r4 k" J' X- O  case 1:
  case 3:
   num_x+=1;
   break;
% ^: `6 y- k+ l$ r5 A5 s  Q  case 2:
   num_x+=2;
   break;
  }
$ u. U5 j% x3 Y2 ?: Z/ n }</P>
  p5 R) R+ R) C9 t<>//化为标准形式
if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
; i  J7 P) U; }( y8 @8 C8 @! K for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
* I) R  V1 L! v  g$ a   case 1:
% Z9 @$ E+ Z4 ?5 T  Q    num_addv++;
    num_ba++;
$ J. f/ G* g2 x* k, N& U7 N0 |    num_ar++;
( h6 D7 n9 _: R8 K2 l0 v. l  A! H0 i    c[tnum_x+num_addv]=0;
! C  C3 r0 C9 K$ l& B! C- l$ B    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
4 y: U- {& h5 i( @- Y/ J5 ]. ~    for(j=tnum_x+1;j&lt;=num_x;j++)
     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
/ S. N0 b8 E) Z; @9 m) j# q     else a[j]=0;
2 \  ?/ o" L% \# s- r    break;
   case 2:
. k7 J& W4 f8 w' W3 a+ Q2 N    num_addv++;
3 O  W  m5 [% w- W& C    c[tnum_x+num_addv]=0;
    num_addv++;
    num_ba++;
  b3 v! g4 G. H    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;
     else a[j]=0;
    break;
   case 3:
    num_addv++;
    num_ba++;
" F) e& U# B1 J% |    c[tnum_x+num_addv]=0;
    base[num_ba]=tnum_x+num_addv;
8 y) v' W, I$ V  H1 a# Z+ R, ^    for(j=tnum_x+1;j&lt;=num_x;j++)
     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
& ^1 L1 k% j/ q* b3 m0 Z     else a[j]=0;
    break;
  }//switch
}//增加松弛变量、剩余变量、人工变量、确定基变量</P>
7 `1 X3 G8 [9 t& K1 S& t<>//显示标准化后的方程组
3 n! v( |9 h$ w, { printf("\n化为标准形式后：\n\n");
5 f& U7 s, J/ b& B. N/ z  I if(ma_mi==1) printf("max z= ");
else printf("max z'= ");
: _, }# G7 @2 I: a/ f printf("%lg x[%d]",c[1],1);
* o/ F: H0 l3 P for(i=2;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_ar;j++)
  if(i==arti[j]) k=1;
  j+ {/ [) @: N  if(k==1) printf(" -M x[%d]",i);
  else if(c==1) printf(" + x[%d]",i);
( p/ U, J/ h, J7 q8 Y9 c% Z  else if(c==-1) printf(" - x[%d]",i);
$ p: m! b# a+ G7 K  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
  else printf(" %lg x[%d]",c,i);
* `: u; J2 X5 \7 s/ D( d! i/ H }</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  k=0;
1 H5 W* w  n4 E& ]" ]! s$ v  for(j=1;j&lt;=num_x;j++) {
8 N9 V3 f7 N- V   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);
: ?% ]$ V, w* L7 E. h6 M    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);
/ p$ R2 Z) x2 D( G, @. k, a- O    else printf(" %lg x[%d]",a[j],j);
    k=1;
3 J; b6 q8 c6 A/ f- V   }
3 e+ C/ k4 q5 `! f  }
  printf(" == %lg\n",b);
}
$ x* T: p9 c/ } printf(" x[1]~x[%d]&gt;=0",num_x);
6 g7 F" E3 ]$ r3 Z}</P>
- T$ u7 D9 u1 r- p/ i<>void iterative() {
& X8 ~6 ~, g  W. I, Z int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
int base_elem;
) Q* E* d: S, n5 P( K4 v int base_out,base_in;
double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
! A& A+ B* _& f- @* S<> printf("\n\n第%d次迭代：\n\n",stop);
for(i=1;i&lt;=num_st;i++) {
. H+ }- J! g1 f  printf("c%d=%lg\t",base,c[base]);
  printf("b%d=%lg\t",i,b);</P>
<>  switch(step) {
   case 1:
8 S$ [6 t4 n2 [9 o4 q6 t8 d    for(j=1;j&lt;=num_x;j++)
    {
& U) u' D, q( |     printf("a[%d][%d]=%lg\t",i,j,a[j]);
: ^% [) H7 F; p7 Y! H- g    }
    printf("\n");
3 W# i# g# G/ m5 p- S. V  y    break;
   case 2:
- b6 E& C% V: `& {& u( o- s9 b- b    for(j=1;j&lt;=num_x;j++) {
: e4 y  u" b% w3 p     k_a=0;
' a% \6 ], s+ H7 o/ X6 V& t' W2 J     for(l=1;l&lt;=num_ar;l++) if(j==arti[l])k_a=1;
6 z; \+ D* Z! ~     if(k_a!=1) printf("a[%d][%d]=%lg\t",i,j,a[j]);
( J7 j  Z0 j/ G" I, D2 x( B2 B    }
! l. p% j4 P9 u  H    printf("\n");
    break;
( G- _7 V% J% S. W9 g3 q0 `+ d3 v  }
}
3 _; j( W$ F5 N. l//求检验数sigma
# E7 h/ L8 z9 a, T1 Y6 \ for(i=1;i&lt;=num_x;i++) {
  sigma=c;
1 d4 f, e/ w% k  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) {
   case 1:
; K1 e! y% N$ I7 H5 ?% P  \# {    printf("sigma[%d]=%lg\t",i,sigma);
. O( q; L3 A/ n    break;
   case 2:
    k_a=0;
    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
9 ^. n6 P3 n9 O8 I7 y    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
    break;
7 j$ j1 P  f, o. x  }
}
putchar('\n');
//检验检验数sigma是否全小于等于0
k=0;
, B% v$ {* c$ p! n) F for(i=1;i&lt;=num_x;i++) {
8 c+ I& s1 L- W; a  if(sigma&gt;0)
  k=1;
}
8 o  C" H* ]! @# g5 P if(k==0) {
& U% K* x" j; ~& d5 d1 I6 z( q  //sigma是全小于等于0时，检查是否为无穷多最优解
- |, J7 O. z; |  for(i=1;i&lt;=num_x;i++) {
   k_f=k_a=0;
6 O( x4 B$ F9 U# V) L7 d  B   for(j=1;j&lt;=num_ar;j++)
* ~9 ?; \: p  Y) B1 O- [& O   if(i==arti[j]) k_a=1;
- d& D: g6 }. o% x0 a   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;}
9 Z0 m+ I1 T: p, k7 y" n# j   }
. y+ S' W. ?6 S9 {) Y+ d3 ]( u" C  }
  status=1;
  return;
2 ]' D8 i9 ]) A }
( i& N8 X  k/ G/ f* r* q& p//检查是否为无界解
for(i=1;i&lt;=num_x;i++) {
" s* e) g/ B, u  k_f=0;
! A$ {* m8 q/ t2 [( y+ S  if(sigma&gt;0) {
   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
0 y6 T+ Y/ a1 s   if(k_f!=1) {status=0; return;}
' M; o! V& r/ W4 F+ O9 w) m. l  }
}</P>
<>//确定换入变量
for(i=1;i&lt;=num_x;i++) {
6 V; C1 ]5 b7 ]/ U% o* v  k=0;
# B" [5 r5 S0 Z6 R2 W1 }% O  z  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
5 k# i* Q0 _' `" k8 u+ [  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
}//temp赋初值
8 r! i2 m" b' a' ~* [7 s for(i=1;i&lt;=num_x;i++) {
2 G3 v; s2 |1 j  K0 S  k=0;
7 _3 ^5 X" b7 V2 ]  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  if(k==0)
1 Q% y+ a6 @9 o. v: V   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
    base_in=i;
    temp=sigma;
   }
}</P>
<>//确定换出变量
for(i=1;i&lt;=num_st;i++)
  Q; }- E; q) `- U# Q4 b4 w  if(a[base_in]&gt;0) {
4 U# r  K- t  A/ e& l3 o5 `  c   temp=b/a[base_in]+1;
   break;
' A- d" e) d+ x* r* a' [6 R+ O  }//temp赋初值
: G: {+ L+ A+ T- }! |5 ^( H, h9 r3 E& ? for(i=1;i&lt;=num_st;i++) {
  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
   for(j=1;j&lt;=num_ar;j++)
( o. `# N2 j4 f, ]& c" F& m& [    if(base==arti[j]) {
4 e5 ^' }* @9 z/ u, V; F     base_out=base;
     base_elem=i;
8 b2 u" l! p) ?7 f& Q% _. \     temp=b/a[base_in];
- x! v2 b1 |- X: c% \" H     break;
' l6 E  N. u4 s    }
  }//人工变量优先换出
  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
   base_out=base;
   base_elem=i;
4 P' ^6 T) L8 E5 f% |   temp=b/a[base_in];
8 j: a( |$ W9 J: D# H0 G1 B' {  }
4 \) w- ?0 X' j, Z& b# p# ?& F% |5 ? }</P>
/ Z$ ?" @- M" |<> printf(" 基变量：");
for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
//基变量变换，进行新方程初始化后迭代
for(i=1;i&lt;=num_st;i++) {
5 H& \; a( B$ t" G! Y6 w  if(base==base_out) base=base_in;
9 T$ y/ G$ X0 L9 ~5 z }
7 m6 d. ^; A4 S* f- f# o//初始化主元素行系数
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>
<> for(i=1;i&lt;=num_st;i++) {
0 v! e! d% g( T9 M. L2 [3 h  if(i!=base_elem) {
% U) z$ t) \4 D' L3 o. W   b-=b[base_elem]*a[base_in];
2 Q, I4 s6 }7 y2 f   value_be=a[base_in];
! r% o3 G5 A! Y9 A) b* j   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
4 g  @6 [4 H& X9 u! \5 o. [6 {  }
}
stop++;
/ W( X+ h8 m' c% a" {: v4 O. } if(stop&gt;STP) {status=-2; return;}
! {; q! e' Y6 O* [! d iterative();
3 E8 E8 Y6 c: g+ `- b3 E* F}</P>
<>void output() {
int i,j;
double X[MAX];
printf("\n结果如下：\n");
' D1 m  a+ x7 k* D3 ^) Z9 T& u printf("\nX=(");
- J6 b& y" I( c8 A: }  F for(i=1;i&lt;=num_x;i++) {
/ m8 Y9 u, s/ h# I) ]1 G  for(j=1;j&lt;=num_st;j++)
8 _2 M6 u/ m  G# s" p3 U  if(i==base[j]) {X=b[j];break;}
7 @/ n9 s- n; {% W) {( l! K1 b- `  else X=0;
  printf("%lg ",X);
}
/ y' e% \" Y  i  U printf(")");
* F; }5 Y+ r$ j* l" ` for(i=1;i&lt;=num_x;i++) max+=c*X;
2 k; [5 {# E% o) a if(ma_mi==1) printf("\nMax z= %lf\n",max);
else printf("\nMin z= %lf\n",-max);
}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了