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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
: Y: R% w" ?( x, B/ Z  Q#include&lt;math.h&gt;
% n' H- e  a9 M8 Q#define X 5
- G/ n" d! @7 b4 p% r* K#define Y 7
! V: }2 F6 C; [3 S% m9 m0 f0 ~void xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
{
  int j;
* y) m2 J9 \- D6 S5 [  V$ y& \2 I  *c=0;
4 u5 \* f! K7 O; J5 k1 a) E. _  for(j=1;j&lt;=*is;j++)
4 n8 q: f# J  S- r    {
+ [& h' Y. E/ f0 [5 ]      if((*a)[*ir][j]-*c&gt;0)
    {
4 T# F. V5 _/ H4 R      *c=(*a)[*ir][j];
      *j0=j;
    }
3 q- x2 F% z* ?/ j1 h     }
% O3 _8 {7 q/ O% C# v}</P><>      /***************** 参数说明 **********************/
) _  y! k! [. |/ l$ S9 W      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
) y$ L0 F0 n9 ^! k      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
8 n: i, h6 Y# n% E. B5 I4 H      /********* <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],
( U) f! l3 J* V8 w       int(*k)[],float(*x)[])
6 B  _( U' v- F0 r$ |5 j{
/ J, o3 @6 T: p  int i,m1,mn,j0,i0,j;
( ?9 B  M+ [/ i; J# \  float c,g;
! I& d1 f& h& ~- h7 L; K0 T  m1=m+1;
9 ?0 E# f$ J# @2 Q6 L2 T2 S  X; A  mn=m+n;
. D- v8 O) I8 K, m  for(i=1;i&lt;=m;i++)
    (*k)=n+i;
leap1:
' i* q; W3 X3 [6 u( @# ~) O6 A  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
* ^0 y, p8 n0 D  else
( l- A- S) K$ ?9 S- i- v     {
3 R2 z, A! M# z( j/ L+ v2 }# bleap2:   if(l1-50==0)
2 @% s, F0 @& Z3 a, }1 B      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
2 X5 O' H7 U+ ]+ g, M# k7 B4 Z" J       else
$ R6 C4 I( Z( ?/ j  W! \4 d( b      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
8 w" a8 j. E3 f. H6 a7 H% S% z     }
2 f3 j( j* I! T1 x% A   c-=1e-8;
( h* a# p2 x+ L& L$ Q1 r   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
+ {  L$ ]1 x" J8 f* L      {
    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
    return(-1);
      }
/ P5 J# X- U& }4 \; Y1 I# h7 \   if((c&lt;0)&amp;&amp;(l1==1))
- p) \4 i  Z% O1 T& _8 w1 _9 n       {
     l1=50;
& @5 f( S5 L# F0 D4 y! o     goto leap2;
       }
    if(c&lt;=0)
       {
$ k1 Q1 L$ q7 q, d& T, S      for(i=1;i&lt;=m;i++)
  m; d2 q$ J0 n- v& b  d# t$ E/ Z. Y         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
4 o" A) W; K- G: `) F% a      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
7 ]$ k2 M: K2 I( s: ~5 f) E; O8 m      printf("\nF=%f\n",(*a)[m2][mn1]);
7 m+ _/ N7 Q' X# l4 X7 @: B      return 0;
+ T; c' v7 P7 ~0 |- X' W, t    }
) I* i$ {+ l+ @' _     c=1e8;
     for(i=1;i&lt;=m;i++)
5 y5 t2 C0 i7 a2 J4 r! S       {
     if((*a)[j0]&gt;1e-8)
        {
0 E. m5 b8 b6 ^          g=(*a)[mn1]/(*a)[j0];
          if(g-c&lt;0);
9 w  X5 A4 i8 D: `" |         {
           c=g;
           i0=i;
          }
         }
    }
  @6 h' p& K( z" K      if(c==1e+8)
    {
* @) v) T; x8 q9 M% E1 [      printf("\n\t*******Lp no solution********\n");
      return -3;
     }
) ], ?' K6 O1 [  C0 i0 c; |       (*k)[i0]=j0;
1 n! F/ H; H0 a4 n       for(j=1;j&lt;=mn1;j++)
" Q8 s5 k( Z% ]4 Z- o: a% M% d% X      {
5 V6 [" H$ Z6 u) l8 ]7 X) |2 G        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
           continue;
         g=(*a)[i0][j]/(*a)[i0][j0];
           (*a)[i0][j]=g;
         for(i=1;i&lt;=m2;i++)
0 F. _4 s6 M) |; `' R           {
         if((i==i0)||(!(l1==1)&amp;(i==m1)))
            continue;
         (*a)[j]=(*a)[j]-(*a)[j0]*g;
  h4 G; Y* b; ^           }
; e. X, ~# u" @" a; f+ _" H! L       }
       for(i=1;i&lt;=m2;i++)
  ?. c6 I7 [6 Q$ n) Q4 V      (*a)[j0]=0;
+ _9 B; i& c, r) |3 L       (*a)[i0][j0]=1;
4 M5 S* G4 B5 _! E0 K" G+ X       goto leap1;
}  
7 |. F& i8 p# W: z  D  D3 b; cmain()
{
. r$ C- K$ y* j& L, i7 v  float a[5][7]={{0,0,0, 0,  0,0,0},
. ?  B( m, F" ]9 W- [. {& g0 C                 {0,3,-4,3,  1,0,12},
) f; ~, P9 A4 S  K" ^                 {0,3,0, 6,  0,1,12},
$ O4 Y- v: y5 O/ P9 L                 {0,0,0, 0,  0,0,0},
! O$ t4 R& [" J+ H8 l3 c; u4 p                 {0,69,0,144,0,0,300}};
, `  A7 |. G' \. E# n  [  int k[3];
9 `( o3 w2 T7 T% G  float x[3];
" s1 w; Y0 g% t/ z: j8 \1 O: c  clrscr();
8 g4 w" b+ M2 c) K  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
}
2 l4 E# e% ]5 U$ R</P>

plgatc

<>/*************************************************************************
# G- u6 P5 ^* ]1 x6 E% g& G                    单纯型法解线性规划问题（两阶段法）        

编程环境：VC++6.0     
方程组输入说明：
( ~; o) v2 _# f7 T; s) s, T 变量非负，按提示输入相关参数。
( @4 _; Q) m5 ^ *************************************************************************/
#include &lt;stdio.h&gt;
) l& C+ x4 k+ T#include &lt;stdlib.h&gt;
) e! b0 E" a8 j2 e5 ]  g#define MAX 100
, l% t1 p8 a7 q4 ?/ x#define STP 100</P>
) Y: g' q1 a; ?0 K<>int stop=1; //迭代记数变量
! a$ R7 F4 H3 ~" e3 bint status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
int step=1; //目前阶段</P>
/ N" X  Q: V2 C/ V! ~<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
4 D5 a8 m& U4 i/ sint num_x; //变量个数
int num_st; //约束方程数
int num_ar=0; //人工变量个数
, K) ~% S$ w$ D5 R; W! M# aint arti[MAX]; //人工变量下标
) k7 n& W; h  }- Pint base[MAX]; //基变量下标
8 @0 Z  O8 E" z6 D; L' I2 Bint ma_mi; //1为求最大值，2为求最小值</P>
  @3 P; }7 a& f$ o8 T  l<>void create(); //建立方程组
void iterative(); //单纯型法迭代
void output(); //输出结果
7 U* s; O% `$ w$ G) C$ fvoid banner(); //打印程序标题
void exchange(); //交换两阶段价值系数
) N2 ]+ O, k4 Z* }. xvoid show(); //输出方程组</P>
<>void main() {
int i,j,k;
banner();
( A2 y0 I$ n5 r! L) z create();
1 _3 l8 p8 ^! s! e+ x; L% C0 I( x; W //保存原价值系数，转换为第一阶段价值系数
for(i=1;i&lt;=num_x;i++) {
( T& q+ s; e5 k5 V! _  k=0;
  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
; e' P$ W6 P2 F* P6 T' j  if(k==1) temp_c=-1;
  else temp_c=0;
3 }$ k4 p. F8 {7 J; w( E }
exchange(c,temp_c);</P>
<> printf("\n\n第一阶段问题为：\n\n");
8 h5 ]5 x# N% b4 ~ show();
step++;
! d2 _" l: J/ U9 m* M) \ printf("\n\n按回车开始第一阶段迭代");
+ H) G0 |6 Z! i  g6 o6 [ getchar();
getchar();
iterative();
if(status==-2) {
  puts("迭代超过限制次数强行终止！\n");
9 A+ ~/ B# h1 i  puts("\n按回车结束");
5 D: A- I) p- }0 Y& d. {+ b& ~  getchar();
  exit(0);
7 z, g- J2 r' y/ ^- v! A }
  M9 ^3 i5 A( A+ U! n4 q, E* Y output();</P>
<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
  puts("\n按回车结束");
  getchar();
" _6 D& m( a! F, S; F9 Q& c7 p+ Z  exit(0);
0 g8 M# W; ?3 y0 p0 k& y }</P>
3 _3 V0 Y# e3 l' \, N<> //转换为第二阶段价值系数
. F0 m" h# B% H* z exchange(c,temp_c);
% V+ Z: ^7 b1 I& v, `( n //把人工变量列全设为0
. [. J* `* m8 _/ P: } for(i=1;i&lt;=num_ar;i++) {
6 m- K; s' Q4 M  l9 H  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
' g# j& G, o3 }0 o }</P>
<> puts("\n\n第二阶段问题为：\n\n");
  J7 }4 |$ U! m1 ?8 s show();
puts("\n\n按回车开始第二阶段迭代");
getchar();
9 D3 g0 `8 F: B- n, F/ ^0 |2 f iterative();
9 r! B# }  _% z6 [/ m: J switch(status) {
case 1:
7 D. ?+ o' Y4 t: ]- Z( R  output();
. i9 V0 J9 Y: N- X! ]7 L  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
1 M4 m7 k1 g7 I( H5 ~ case 0:
  puts("\n\n原问题为无界解。\n");
# i7 g% h+ Y* e' g2 i& l  puts("\n按回车结束");
4 N. U8 s1 E. L& ~4 z1 m1 w6 O  getchar();
  exit(0);
. ~) g$ Q2 |% h case -1:
  output();
  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
- y0 [# p% A& l- q/ @ case -2:
  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
}//switch
, i# u7 @- `5 X( L) B) {9 I4 }}</P>
<>void banner() {
printf("\t\t****************************************\n");
printf("\t\t         单纯型法解线性规划问题\n");
: e2 F; w: _+ m( ?. @. i printf("\t\t                         作者：Thunder\n");
printf("\t\t****************************************\n");
printf("\n");
}</P>
- d! D! N- K6 J. v& h7 b+ `; r" _<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
//为1的不显示系数1，-1系数只显示负号
& c- s5 q2 ~( M1 T. q: D int i,j,k;
switch(step) {
case 1:
  printf("min z= ");
4 U7 J1 d4 P7 Q+ G  printf("x[%d]",arti[1]);
  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
case 2:
& n2 E6 V1 j9 \$ r2 r- y$ P printf("max z= ");
printf("%lg x[%d]",c[1],1);
! J7 O1 S: [# S" \& \0 p for(i=2;i&lt;=num_x;i++) {
) E( f  h- ]3 T4 F- s5 m3 ?. h$ u  if(c==1) printf(" + x[%d]",i);
6 ?- b& h7 m7 k9 t9 z, z  else if(c==-1) printf(" - x[%d]",i);
  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
6 @' u2 N7 g8 `/ x" F, _  else printf(" %lg x[%d]",c,i);
}
' Q$ V3 x! v/ X& _. h break;
}</P>
' g2 ~3 w* y, R  b) l( n5 y<> printf("\nst:\n");
& [8 N3 v0 w9 Y* P& g for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
$ [0 F7 C2 `# J! |# L9 M( O. ]3 m# R   if(a[j]!=0) {
8 V" C' w8 d; T  E, M2 {    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);
1 B& S* [, g/ o    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);
" u* f6 E# e3 j  Z6 g. e" [+ ]  T! W    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);
2 B6 T( }; J/ ]' {    k=1;
   }
0 B4 d& v6 s! p0 v- n0 D  }
, H0 {6 Q1 n, l  printf(" == %lg\n",b);
+ \# a1 L( \3 Y, Z* a1 F3 \8 I& o }
6 e& u$ X  P3 ]$ n: ?# @$ e printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
3 I" M* m; h2 }<>void exchange() {
% ~! v+ p- U- @% E& Z+ w5 z( ? int i;
; r* O8 W+ @5 C" G* T double temp[MAX];
0 i+ A/ L9 u$ y. t2 x$ Q$ S for(i=1;i&lt;=num_x;i++) {
  temp=temp_c;
" d' G1 ~! R9 F8 K' s. G, J3 s  temp_c=c;
: g& t# Q9 J: E! ?8 |  c=temp;
}
# ?, s1 `4 v& C/ E  r}</P>
<>void create() {
) G9 W7 t, C* b0 _( i( }8 M; T//输入方程组系数，每个方程输完后回显确认
7 T5 k% g; U- z) W8 C3 [) |9 M8 K int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
6 ?4 e% o" m6 }6 v* Z! f. T& x7 w char confirm;

, l  k. A4 T' p9 e while(1) {
  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  scanf("%d",&amp;ma_mi);
  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
: S% M2 g+ Y% `3 i2 o2 Q- F  else break;
4 J2 R2 u! _. v5 E }

+ x5 T: ^7 `- C% l# w while(1) {
( a8 L; \1 Q& N, L: J( K  printf("指定变量个数：");
  scanf("%d",&amp;num_x);
  printf("输入价值系数c1-c%d:\n",num_x);
8 @5 h1 s# I4 H  for(i=1;i&lt;=num_x;i++) {   
   printf("c%d=",i);
9 {( k' A$ c$ p' m- O+ |/ G/ ^   scanf("%lf",&amp;c);
  }
  if(ma_mi==1) printf("max z= ");
  else printf("min z= ");
5 K2 w8 F% N3 W  c2 A  printf("%lg x[%d]",c[1],1);
  for(i=2;i&lt;=num_x;i++) {
! f( S5 r* y! ]9 i$ g: v8 }' _2 Q) V   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
+ P; C8 \9 M. o* }   else printf("  %lg x[%d]",c,i);
1 t9 a5 T( d7 _( l& x  }
8 j  M  z$ J. i printf("\n正确吗？：(y/n)");
getchar();
confirm=getchar();
( X/ L0 `1 N* k% x( H' B if (confirm=='y') break;
else if(confirm=='n') continue;
3 w. I5 J3 e# T) ^) P }</P>
<> printf("输入约束方程组个数：");
$ R1 U6 O4 T3 z+ \; W* ~, v! K scanf("%d",&amp;num_st);
$ {* q2 d1 n/ L0 i7 N for(i=1;i&lt;=num_st;i++) {
% G& M! M4 d; D, P. L5 a  printf("st.%d:\n",i);
  while(1) {
   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
   scanf("%d",&amp;re_st);
( e. r0 n, ?  c& B   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
   else break;
4 k+ M( T  R. l3 B! w  M) J* l( g  }
  printf("输入技术系数：\n");
+ v4 O- l' j$ D; g  \! G  for(j=1;j&lt;=num_x;j++) {  
: M* y# \- E/ x1 `6 x   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
  }
( A- o" E, w* l0 D6 m1 S  printf("输入资源拥有量：\nb%d=",i);
- M$ [# @3 z" `& ]$ B( L1 b( U  scanf("%lf",&amp;b);
5 d2 J1 m, D/ b6 n. u  
  printf("st.%i:\n",i);
  printf("%lg x[%d]",a[1],1);
  for(j=2;j&lt;=num_x;j++) {
( Q$ ^# M0 g3 F4 S: _$ o& G   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
   else printf(" %lg x[%d]",a[j],j);
- P( i: p3 o0 Y5 A( A* P9 ], l. E3 ?  }
. J: ~' K3 f  I/ s$ x" @( c  switch(re_st) {
   case 1: printf(" == %lg",b); break;
   case 2: printf(" &gt;= %lg",b); break;
+ y% e+ g8 f2 p  p) k) C' [   case 3: printf(" &lt;= %lg",b); break;
  }</P>
<>  while(1) {
: B& j' Q5 v" Y- R- \7 [4 S   printf("\n正确吗？(y/n)");
   getchar();
6 H, e$ U0 L; H   confirm=getchar();
7 A1 j) U3 m2 A5 o2 z4 E. T   if (confirm=='y') break;
0 c8 D, d3 r/ A8 U- g$ Z) F# W   else if(confirm=='n') {i-=1; break;}
# t& N; H. D; k, s$ U7 L- @  @  }
}</P>
<>//显示输入的方程组
- ~/ u9 k0 k+ _' x- {: y/ E printf("\n原问题为：\n\n");
+ z( i8 D1 _$ f% \  g& K8 G if(ma_mi==1) printf("max z= ");
else printf("min z= ");
1 u. z2 o# D" A: U printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
  I: R7 l: L1 D6 f3 M) H, e9 Z% Y9 @  if(c==1) printf(" + x[%d]",i);
1 P/ ?3 {9 K# M8 ^  else if(c==-1) printf(" - x[%d]",i);
  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
" M$ A8 L4 O; c7 c% U# n: |  B+ I5 s  else printf(" %lg x[%d]",c,i);
# k" B8 I9 K+ ? }</P>
- E& G7 i# d+ p) F% X# V- `( j<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
9 a% e3 n' ]+ |) C% x' s6 B- b  k=0;
0 c# W" }; Z& a4 v& |& J  for(j=1;j&lt;=num_x;j++) {
: z! o& O1 W1 l2 g: I   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);
, O1 t& }8 \& Y5 O' w1 N$ S% d    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);
% v0 a! f6 R8 \    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;
9 t( P. O: C! s. N3 _& c' X3 X: k3 j   }
5 y. w2 p( U. G2 v+ x! E  }
' N2 H& L4 o  b) o$ Q+ N. w  switch(re_st) {
. }+ D- T3 C# o8 H$ ]) m, O8 F   case 1:
4 b8 D6 h& N2 d% C( w5 M    printf(" == %lg\n",b);
6 R# J5 Q: l; [" m, s; l0 y. h    break;
' z$ k' z8 L% S4 C) R8 p: t! c+ q   case 2:
+ v# K! `0 k6 R; P6 R' ~    printf(" &gt;= %lg\n",b);
    break;
   case 3:
    printf(" &lt;= %lg\n",b);
    break;
, B. W* ]# T$ {: B  }
4 I" ]: x: T, R2 V, R2 [7 E }
printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
<> tnum_x=num_x;
for(i=1;i&lt;=num_st;i++) {
5 R& V  }6 Z6 Y$ c- b+ r& L: z* I  switch(re_st) {
7 o/ p" j8 T/ u6 j  case 1:
0 Q! C& J3 r, W9 T* Z, }  case 3:
, S3 K& s6 I' P* y! U   num_x+=1;
% k/ n& F; t! u# m   break;
  case 2:
5 M+ E( h+ K3 c: I6 [6 v5 T   num_x+=2;
0 L4 ^; W" c; S6 d1 D1 U   break;
% d% d# \7 N7 `; v( `, V$ z+ {3 s7 m  }
$ A7 @1 J3 h7 K8 r) P3 u }</P>
" ], D) K3 y% ~7 d3 k<>//化为标准形式
- ^2 f" s# @$ g5 @7 i: i) o if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
" x* M" u: Q0 r2 e  } for(i=1;i&lt;=num_st;i++) {
/ B/ C+ S/ l- }5 G5 Q$ B  switch(re_st) {
   case 1:
    num_addv++;
; }) F! n9 V$ F& Y0 l% o- `. y6 G    num_ba++;
    num_ar++;
    c[tnum_x+num_addv]=0;
    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
- Z# y' D& g$ N1 ]( T8 i    for(j=tnum_x+1;j&lt;=num_x;j++)
& S% n+ ~% _$ c/ Y# x. Z, f     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
' r, D7 j* l) p3 }" q/ {- l3 f    break;
   case 2:
    num_addv++;
: ^. y2 ?" z% y9 X8 A    c[tnum_x+num_addv]=0;
    num_addv++;
: E0 B5 m. x# D( e  a/ G    num_ba++;
4 a" y0 Z4 a2 S& I8 _! @" Y: O, a3 o    num_ar++;
    c[tnum_x+num_addv]=0;
% H! P. B8 \6 Q& I. C    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
    for(j=tnum_x+1;j&lt;=num_x;j++)
4 y$ |2 v0 ]; k3 L     if(j==tnum_x+num_addv-1) a[tnum_x+num_addv-1]=-1;
2 X. B) t2 o  _* c& }2 b* U     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
! p: K1 z- M9 @& G* g0 b4 n     else a[j]=0;
9 m, G2 D& v% o, E" @3 v: }/ D5 v5 J! f    break;
* s* {0 {. K3 s4 b4 D& j   case 3:
    num_addv++;
    num_ba++;
    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;
     else a[j]=0;
    break;
0 E8 o' W6 R; T4 a  }//switch
: a. T; }' y5 C' Q0 B4 a1 U" P }//增加松弛变量、剩余变量、人工变量、确定基变量</P>
9 }6 P% o" t8 e4 Z4 ~<>//显示标准化后的方程组
, q7 q& @- X  @9 G8 n printf("\n化为标准形式后：\n\n");
if(ma_mi==1) printf("max z= ");
! q. R4 l2 C" }6 j3 j. G, c9 U else printf("max z'= ");
printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
5 S- Z- m: Z0 |7 \3 C. F  k=0;
4 l3 u* S) H, g# D) }9 Y1 f; J) r  for(j=1;j&lt;=num_ar;j++)
: j. g/ Z( }# u% k5 ~  if(i==arti[j]) k=1;
5 k  ~8 c% X5 J! {: U; s" s  if(k==1) printf(" -M x[%d]",i);
+ R% h2 ?0 m8 l  else if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
0 m9 R- Z, L8 N  ?4 i, X  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=0;
% u- G8 ?) k9 d4 r  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
2 U+ S: z' ~% x( L1 r    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);
# f$ }$ u3 s& ^    else printf(" %lg x[%d]",a[j],j);
    k=1;
   }
' ?& T, ]0 w) v* g! U  }
" w* i& d) e& w: O, Q# M, }/ f  printf(" == %lg\n",b);
}
printf(" x[1]~x[%d]&gt;=0",num_x);
. p6 \. J/ ~/ ?4 s  V% d% U}</P>
<>void iterative() {
/ O9 P' m" ^0 u1 t. b- f 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;
( F* ?2 N0 I3 {5 a2 z double sigma[MAX],temp;
& R9 _' J( m" W. y9 M0 _ double value_be; //高斯消元里保存主元素值</P>
<> printf("\n\n第%d次迭代：\n\n",stop);
) D' N7 m& c/ U  ]# s for(i=1;i&lt;=num_st;i++) {
  printf("c%d=%lg\t",base,c[base]);
! e, z% M  H; _! V% k5 \2 `  printf("b%d=%lg\t",i,b);</P>
<>  switch(step) {
   case 1:
/ I, l+ t; s. q$ Y    for(j=1;j&lt;=num_x;j++)
9 R! F& H& A! B- X9 q    {
     printf("a[%d][%d]=%lg\t",i,j,a[j]);
  R0 K8 \  `  C" E    }
+ [; B4 t& |2 |1 `( R    printf("\n");
* o7 `7 i$ [- Y    break;
7 d3 M2 v& y! m; ~; G* y  A   case 2:
    for(j=1;j&lt;=num_x;j++) {
9 |0 f! f" g2 [% p1 f" I, T! {5 ^& J     k_a=0;
) B3 a. E3 C- H. C1 z& o     for(l=1;l&lt;=num_ar;l++) if(j==arti[l])k_a=1;
2 Q! U7 T" c8 J  M( j0 L2 G6 G     if(k_a!=1) printf("a[%d][%d]=%lg\t",i,j,a[j]);
9 P9 e5 R, f3 L' ]' j4 L2 B6 N    }
    printf("\n");
, w& c0 B5 s+ ^, s" ?/ T% a    break;
( ^; Z$ a1 L5 z( f$ w  `3 f' I  }
}
//求检验数sigma
1 E3 B, y1 ?- } for(i=1;i&lt;=num_x;i++) {
  sigma=c;
  for(j=1;j&lt;=num_st;j++) sigma-=c[base[j]]*a[j];
' t  w* h, Q2 U5 U7 W. x7 Z  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
& R% K+ j  @  M1 f  switch(step) {
6 ], `$ l3 {4 f   case 1:
    printf("sigma[%d]=%lg\t",i,sigma);
% f+ ^  n& m- C+ s3 r' E    break;
& i' u* {2 I" X   case 2:
    k_a=0;
    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
+ n$ R6 {# S9 v( I$ S( @+ V    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
( S# y6 r. L) r9 ?' L9 d- v    break;
* `' r( Q5 A' Q1 c, G" T  }
4 `: [/ ^# _- d; O  k+ h8 w2 U }
putchar('\n');
//检验检验数sigma是否全小于等于0
3 E- Q: M  a' u k=0;
for(i=1;i&lt;=num_x;i++) {
  if(sigma&gt;0)
  k=1;
3 L) Y; Q3 }) J: O+ U% Y }
3 g$ e( p; o6 @$ i if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
  for(i=1;i&lt;=num_x;i++) {
   k_f=k_a=0;
   for(j=1;j&lt;=num_ar;j++)
   if(i==arti[j]) k_a=1;
   if(sigma==0&amp;&amp;k_a!=1) {
- G5 N, z) i/ x/ v: S6 [2 B    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
   if(k_f==0) {status=-1; return;}
: d5 S4 K9 N4 C: b# ~7 Z+ S0 v. p   }
6 Y9 v6 B" i: b3 d" N- V  }
" L5 v1 f( r- E) X7 D. N8 ^  status=1;
  return;
: L7 `4 B! l+ d: w$ `# I6 Q8 A }
& v! P8 l' j! k. v1 P8 w) c0 ?: k//检查是否为无界解
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;
   if(k_f!=1) {status=0; return;}
  }
}</P>
<>//确定换入变量
for(i=1;i&lt;=num_x;i++) {
* B, Q& }5 Q% O# W' e5 c4 `  k=0;
+ ~7 `6 T" I7 u# A  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;
% c! J* b/ g8 j# m }//temp赋初值
+ d+ A3 c: l/ Y$ D! d/ a# W for(i=1;i&lt;=num_x;i++) {
  k=0;
3 e8 S- w# [6 y" O$ o  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  if(k==0)
  L' @9 n8 S7 `/ S. R+ V; Y   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
    base_in=i;
, l9 y  p9 n' U: r    temp=sigma;
   }
}</P>
<>//确定换出变量
2 g) _% ]6 c( W. H for(i=1;i&lt;=num_st;i++)
  if(a[base_in]&gt;0) {
   temp=b/a[base_in]+1;
   break;
  }//temp赋初值
for(i=1;i&lt;=num_st;i++) {
  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
0 g; U2 H! }2 z9 h+ U- X/ r3 g" X   for(j=1;j&lt;=num_ar;j++)
    if(base==arti[j]) {
     base_out=base;
2 c( H1 |8 D! f) _' K0 V6 h     base_elem=i;
     temp=b/a[base_in];
     break;
# k7 v: R- x4 I9 r    }
  }//人工变量优先换出
  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
- E0 |0 ]  e6 a9 t4 c8 o8 T# Q3 x   base_out=base;
' y( M" ]# b1 y1 A( L0 A: J   base_elem=i;
! M2 a1 G% q, ?) L4 d   temp=b/a[base_in];
  }
}</P>
" ?# H$ @) M. H  z$ }<> printf(" 基变量：");
# k- W  S8 i7 q& K8 @- v for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
, `5 U+ e1 M# p5 w, B  e! z9 @: e printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
//基变量变换，进行新方程初始化后迭代
9 S# [% f" ?) }% e8 f- K/ Y8 f  l" M for(i=1;i&lt;=num_st;i++) {
$ E$ A) Q! d2 P+ n+ I  v  I- a3 b. D  if(base==base_out) base=base_in;
}
//初始化主元素行系数
% M5 `6 H% f) t( h4 z) P+ l- f value_be=a[base_elem][base_in];
1 X& F3 A7 S% |2 H" X2 {5 `5 Z b[base_elem]/=value_be;
8 l- n* G4 f" e0 W' k* u% a5 O for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
<> for(i=1;i&lt;=num_st;i++) {
; H, v& P  l& x" W) ]  if(i!=base_elem) {
, e: B( I9 t, i( |# H0 @   b-=b[base_elem]*a[base_in];
: i5 Q( y3 K, H- Y. A1 I- x# t   value_be=a[base_in];
   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  }
}
  W+ K9 M* A7 A* U0 k+ u% H stop++;
if(stop&gt;STP) {status=-2; return;}
iterative();
}</P>
' N- T+ }9 G7 n4 e+ X- V" E<>void output() {
int i,j;
" c! @: z5 v, D" x% X double X[MAX];
printf("\n结果如下：\n");
printf("\nX=(");
0 M1 @( x  F$ \$ C& H for(i=1;i&lt;=num_x;i++) {
  for(j=1;j&lt;=num_st;j++)
$ Z$ Q4 L" D$ Z" K2 K  if(i==base[j]) {X=b[j];break;}
- h' g8 k" A6 h# [  else X=0;
  printf("%lg ",X);
/ |/ a- v( p& }/ Q/ r9 F8 d }
0 i7 B* S2 i# A7 ~. `1 G& ^ 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);
* x4 f# u% S& g}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了