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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
# {, `4 r" g8 m2 S9 A2 `#include&lt;math.h&gt;
9 B' q: f+ |/ t#define X 5
5 \2 O" _: i7 s9 ^$ \1 v5 I% n#define Y 7
4 I) @9 f- O# T2 b- V, evoid xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
{
6 G0 J3 _: a7 _* i( @1 J  int j;
  N1 b1 P' B2 D& z, ?3 H  *c=0;
( ?+ J% l  [$ f7 H2 f  for(j=1;j&lt;=*is;j++)
3 N% H  x$ W7 W3 e/ H    {
" V6 c7 U7 B% F' ^      if((*a)[*ir][j]-*c&gt;0)
    {
      *c=(*a)[*ir][j];
      *j0=j;
    }
     }
}</P><>      /***************** 参数说明 **********************/
      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
. v* y! n  X6 D      /********* <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],
       int(*k)[],float(*x)[])
{
  int i,m1,mn,j0,i0,j;
  float c,g;
3 n3 p0 h$ f: i9 r9 D% J  m1=m+1;
  mn=m+n;
  for(i=1;i&lt;=m;i++)
    (*k)=n+i;
leap1:
  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
7 v0 a/ `8 p4 M' L' v7 @1 u  else
: S& ~+ T& j  u% s# r     {
leap2:   if(l1-50==0)
  c7 H) W3 }( G# Z, ]/ o      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
/ h8 g2 E( d3 l( D9 v       else
      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
! Q) j6 }/ u/ e" n# `6 k3 |# b1 d( P     }
# C  U% F# _9 R5 P1 C& Z   c-=1e-8;
; g9 c0 H, q) I& W& e4 L3 E9 Q   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
2 I$ C  |& s- t) p      {
    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
    return(-1);
      }
/ o. h  a6 l" J) @$ C9 S8 J! s   if((c&lt;0)&amp;&amp;(l1==1))
% V4 \! g+ C5 j! U9 m, R; h6 h7 l7 w       {
; c1 L) B$ U: s     l1=50;
  s+ i" \! x& R0 ]) _$ Z     goto leap2;
1 R2 n. z3 ]- J% u6 C       }
    if(c&lt;=0)
: a" i3 |" \5 D/ y( C4 v       {
, y0 j# V( S: |# ^1 N      for(i=1;i&lt;=m;i++)
         (*x)=(*a)[mn1];
5 a, f5 E, D. }* y( A         printf("\n\t********Optimal solution********\n");
      for(i=1;i&lt;=m;i++)
1 h; w6 ]% I: N" _# m. Q         printf("\ni=%d,x=%f\n",(*k),(*x));
! q# l2 \. f- g# d      printf("\nF=%f\n",(*a)[m2][mn1]);
: Q0 Q3 W! R: K  w! X      return 0;
! e6 X4 {- ~* V- T0 y! D    }
     c=1e8;
     for(i=1;i&lt;=m;i++)
1 E' Z4 m8 K9 }7 b) w' t       {
, v5 @4 F/ a" O/ i     if((*a)[j0]&gt;1e-8)
        {
0 h  G, N" I+ X* z4 y          g=(*a)[mn1]/(*a)[j0];
4 v' n# o5 q4 ~+ d6 a$ z7 i          if(g-c&lt;0);
         {
           c=g;
           i0=i;
5 H* ?/ F! @5 [  R. m% W          }
         }
    }
4 y- S* f* p: {4 B9 i0 ~      if(c==1e+8)
6 ~' t0 H# L; u/ d' y. S* z7 l6 R    {
/ _$ ~- C" n  z: A3 m# y2 L$ A2 C      printf("\n\t*******Lp no solution********\n");
' Z* B$ a6 U& _      return -3;
3 n( n- ~$ L+ c# L8 N- _9 f  y     }
       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
      {
, ^. Z9 }' E  Q! Q3 a        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
           continue;
         g=(*a)[i0][j]/(*a)[i0][j0];
' G: l, B+ x" N, o1 z5 ?* s           (*a)[i0][j]=g;
         for(i=1;i&lt;=m2;i++)
* B; v, }  P) ]  n# h+ D' K% a7 u           {
         if((i==i0)||(!(l1==1)&amp;(i==m1)))
* X& M- F, t3 s9 I) s0 [7 U            continue;
         (*a)[j]=(*a)[j]-(*a)[j0]*g;
/ g$ F2 B5 U  `$ [           }
$ J9 R+ n9 t. v" a5 R( u       }
       for(i=1;i&lt;=m2;i++)
. k$ q8 Q7 I: g0 U      (*a)[j0]=0;
       (*a)[i0][j0]=1;
% |. {0 n2 l. {! X# }       goto leap1;
: J' b* w% x# u1 D}  
main()
3 D) c% N& Q' B{
9 j  {5 P$ ]5 Z5 {% b  float a[5][7]={{0,0,0, 0,  0,0,0},
- e& V( Z9 y% E' F0 b( g5 I                 {0,3,-4,3,  1,0,12},
  Y9 I! W. D8 t9 G                 {0,3,0, 6,  0,1,12},
                 {0,0,0, 0,  0,0,0},
                 {0,69,0,144,0,0,300}};
  int k[3];
  float x[3];
  clrscr();
3 p9 X# s* c/ q$ i- y6 R  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
% G) i! g! p% y4 Y$ O0 O, [, W}
2 E: S$ W& c$ o- v</P>

plgatc

<>/*************************************************************************
                    单纯型法解线性规划问题（两阶段法）        
5 L$ q8 d3 r, M) t, s  ~$ v0 ^- Q
2 Y+ L( v  r- I 编程环境：VC++6.0     
$ g" T, j- ?3 A. W7 c/ j 方程组输入说明：
变量非负，按提示输入相关参数。
) `$ ~- j* e8 n- Q *************************************************************************/
! o5 b( W+ Y8 w2 d% k#include &lt;stdio.h&gt;
$ ?; F3 b4 p1 R% C#include &lt;stdlib.h&gt;
#define MAX 100
. D' o0 V1 y/ _, D/ a& ~#define STP 100</P>
& a  W* I& Z0 `% T<>int stop=1; //迭代记数变量
6 N3 ]: O  ~; O0 e. H7 k4 Vint status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
4 E- P! i- {1 v9 x) ]0 T  @7 }6 J, iint step=1; //目前阶段</P>
<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
" _) ]+ j" ?! r) G. v* ?int num_x; //变量个数
" V. E6 I1 [( i+ S& n' Vint num_st; //约束方程数
$ ]9 H# e. ?2 F& V  @int num_ar=0; //人工变量个数
+ {) I) i1 w, h- h) C$ hint arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
<>void create(); //建立方程组
void iterative(); //单纯型法迭代
1 |% Y7 X! {0 evoid output(); //输出结果
void banner(); //打印程序标题
void exchange(); //交换两阶段价值系数
3 e. u4 k/ H) W& Rvoid show(); //输出方程组</P>
( R# M% B0 q3 ?2 x' _<>void main() {
int i,j,k;
banner();
create();
//保存原价值系数，转换为第一阶段价值系数
for(i=1;i&lt;=num_x;i++) {
0 o) T( I+ |4 `" ]: m  k=0;
1 Q3 ]8 f; O- t  ^% Y1 S7 F  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
% g0 _% ^5 S, {  if(k==1) temp_c=-1;
  else temp_c=0;
}
/ N+ V3 ]0 E" C exchange(c,temp_c);</P>
<> printf("\n\n第一阶段问题为：\n\n");
" v4 M9 C3 i6 m; E/ V, V0 u show();
) D; m$ {5 Q( D step++;
9 T% t! r, C/ I& Y- } printf("\n\n按回车开始第一阶段迭代");
getchar();
getchar();
2 ?6 l, J2 C+ B# a7 v9 O iterative();
4 Q: e. \* s2 C2 A0 } if(status==-2) {
. r* L* S% ]0 O9 u# A+ \3 Y6 h  puts("迭代超过限制次数强行终止！\n");
3 I4 k$ n1 m+ O& C! K9 b' H4 A  puts("\n按回车结束");
  getchar();
3 v! _- @- y0 w8 R  exit(0);
}
- Y/ I4 N8 \) U$ A output();</P>
/ i5 T! Z$ A9 {0 p* X* n<> if(max!=0) {
7 [! F0 x& q" q  puts("\n\n原问题无可行解。\n");
  puts("\n按回车结束");
/ o8 c  f$ n. X& D2 C1 y+ e1 u  getchar();
  exit(0);
' [# P8 ]4 t4 e5 ]/ U* n }</P>
<> //转换为第二阶段价值系数
exchange(c,temp_c);
//把人工变量列全设为0
: ?& R9 R8 S3 C* ^ for(i=1;i&lt;=num_ar;i++) {
, J, G3 B, \& T3 t8 r  c[arti]=0;
  d8 s, \! f4 _( Y  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
}</P>
<> puts("\n\n第二阶段问题为：\n\n");
# Q& H* d* `' D9 Q1 M show();
" K7 e6 X( `3 N0 f puts("\n\n按回车开始第二阶段迭代");
getchar();
' P# v: K, D' A iterative();
* E1 G) J" Z$ i8 m' Y switch(status) {
case 1:
0 A0 m7 B1 U2 K9 |  output();
  puts("\n\n原问题有唯一最优解。\n");
3 V% n! B; A1 g+ S, H! |7 Q  puts("\n按回车结束");
  getchar();
  exit(0);
case 0:
  puts("\n\n原问题为无界解。\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
case -1:
& H( t# Z! i9 K& P& A5 u$ k  output();
% w. y+ L, B7 z4 ~- m! A  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
2 c; u! |5 [0 p+ I- c7 D' N  getchar();
2 a5 I9 d: V6 w  exit(0);
case -2:
  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
  getchar();
' P! P( {8 s' B: r. ^  exit(0);
0 T% |, A/ i1 c% Q6 j# ^ }//switch
}</P>
<>void banner() {
. }. Y! Z& \& F5 x printf("\t\t****************************************\n");
printf("\t\t         单纯型法解线性规划问题\n");
printf("\t\t                         作者：Thunder\n");
' U! L; o6 z% `/ k# a. [5 m0 E) a printf("\t\t****************************************\n");
+ e0 A8 b6 h; [ printf("\n");
- M! f1 s; r) N+ C0 F: m}</P>
<>void show() {
: N/ C5 z2 U5 a: K# v8 g//对方程组以自然的格式输出，系数为零的x不显示
//为1的不显示系数1，-1系数只显示负号
0 n/ h; Z2 ?/ n% [, q int i,j,k;
switch(step) {
case 1:
! _' U& V1 B  I" z& r8 p. {  printf("min z= ");
  printf("x[%d]",arti[1]);
8 ]. U1 H% G. r6 B  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
2 U. n4 B' g( p3 q3 ~  break;
case 2:
/ b# k7 w8 m2 g- _& ?0 p3 o2 t printf("max z= ");
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);
- b$ {' U$ [( C% I  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
1 S! B! l+ S3 |3 s8 s6 p  else printf(" %lg x[%d]",c,i);
/ ^, `  k' p3 V$ R% @3 K6 L }
% W. y, f" j" j! `. d1 b7 {0 ^ break;
: Z1 S" v% ]; n" G* p" L }</P>
<> printf("\nst:\n");
; c0 v4 N% l" C# ]% M9 Y6 h for(i=1;i&lt;=num_st;i++) {
6 r0 \% K$ I# `7 X  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);
    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
7 z3 V/ o; h6 ]; O( h5 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);
    else printf(" %lg x[%d]",a[j],j);
0 M% @; R9 D; G  v' H, a    k=1;
   }
# D% g) S; y2 M$ N& @2 W2 t  }
  printf(" == %lg\n",b);
}
! E: \0 n7 _. u! ]: K" E0 { printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
# W% b) f( _6 m  r" t2 I<>void exchange() {
" e  K3 N1 n& F% A) Q/ b int i;
double temp[MAX];
for(i=1;i&lt;=num_x;i++) {
; F* a& w# r: f4 ?0 ]6 v9 t  temp=temp_c;
  temp_c=c;
- P% T5 x( y* Y  V  c=temp;
7 \! I- C' d1 R% q }
}</P>
<>void create() {
//输入方程组系数，每个方程输完后回显确认
9 n( A" E8 L3 y, G int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
char confirm;

while(1) {
  printf("请选择：1、求最大值，2、求最小值：(1/2)");
9 T9 y- o2 U+ N5 G  scanf("%d",&amp;ma_mi);
  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
: ~3 ]) x6 k( x  Y  else break;
}

while(1) {
, p  i: }; b- W% R' R0 y# x  printf("指定变量个数：");
  scanf("%d",&amp;num_x);
4 f: K3 B# B/ b& Z) l" \. G/ `: G  printf("输入价值系数c1-c%d:\n",num_x);
! ^: L( V2 z6 k, r. N  for(i=1;i&lt;=num_x;i++) {   
3 z6 ^9 \5 K/ A% h   printf("c%d=",i);
( E9 K1 m) `) M: S- ^5 [' R   scanf("%lf",&amp;c);
  }
2 ?# m+ K+ v+ k9 y. d7 y, r  if(ma_mi==1) printf("max z= ");
  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);
   else printf("  %lg x[%d]",c,i);
  }
printf("\n正确吗？：(y/n)");
& Z, c  {6 M/ W: h" ^9 F2 q% j getchar();
8 B$ _) v" W+ b( k+ v9 V3 ` confirm=getchar();
if (confirm=='y') break;
else if(confirm=='n') continue;
7 ?& p$ S. M  v- p! X6 J# U0 P }</P>
<> printf("输入约束方程组个数：");
scanf("%d",&amp;num_st);
for(i=1;i&lt;=num_st;i++) {
  printf("st.%d:\n",i);
' s* }2 x) _! d" `  while(1) {
$ C: P5 Z' z( `* @& y. }3 |1 H   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
$ R3 l$ K0 t, T7 S6 S7 ~1 Y- A   scanf("%d",&amp;re_st);
, c1 m, B0 u! P3 Q: z   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
   else break;
: O4 x3 c6 A) W7 ~4 c  }
4 S& b  K9 `+ b! U1 ?7 i; e. w* u  printf("输入技术系数：\n");
8 k: h' ~% e# {! U& R7 W  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
  }
  printf("输入资源拥有量：\nb%d=",i);
4 k; O# J$ R$ P) V. J) z  scanf("%lf",&amp;b);
  
4 H6 {& n, E8 K7 P! y0 G" s- v* c  printf("st.%i:\n",i);
, X- w* K5 |/ u8 {; j! K& u5 k8 o  printf("%lg x[%d]",a[1],1);
$ o! V* Y4 n- e% ~$ K  H  k  for(j=2;j&lt;=num_x;j++) {
   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
' _, G) |, P; m8 S+ b   else printf(" %lg x[%d]",a[j],j);
, u9 a0 L0 g' z$ a. `5 Q  }
. x: m0 e: j  b  switch(re_st) {
1 @4 A. L1 o! ^) _" q   case 1: printf(" == %lg",b); break;
# T& q0 k% ]0 z  q0 E   case 2: printf(" &gt;= %lg",b); break;
   case 3: printf(" &lt;= %lg",b); break;
  }</P>
<>  while(1) {
   printf("\n正确吗？(y/n)");
   getchar();
3 {1 ^- g) A" z5 F7 \1 ]   confirm=getchar();
   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
  o6 G# ^) O! W! w4 N* {  }
}</P>
<>//显示输入的方程组
printf("\n原问题为：\n\n");
, U  U1 p0 G. {" s+ {) P! F if(ma_mi==1) printf("max z= ");
else printf("min z= ");
3 B2 k) N; E( L* Q printf("%lg x[%d]",c[1],1);
8 ?& ]: |" |, {' ]% Y: ], L3 q for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
9 [9 ?! P( v' a3 U  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
  else printf(" %lg x[%d]",c,i);
3 n. Y; o/ v+ N8 l1 [  ~ }</P>
<> printf("\nst:\n");
9 `/ d2 V# n7 ~  T: | for(i=1;i&lt;=num_st;i++) {
+ J# K4 D' V& h' v+ U0 _  k=0;
  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
& E2 U8 T' a% g3 C5 Y7 S$ U( s! 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);
, K- b( X8 @8 h$ n8 m    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);
* c1 z* ]" C# }+ R+ J5 J    else printf(" %lg x[%d]",a[j],j);
    k=1;
3 T  i" _& l% C% s3 v) n   }
0 l0 m# ~: |7 X  }
  switch(re_st) {
4 Y0 a. [% X3 g   case 1:
' q* N  l* e+ P% k6 W9 V    printf(" == %lg\n",b);
0 _% O# E' {. t! H" K+ B/ X" L- l    break;
+ s4 |3 M) S: E. i0 y   case 2:
    printf(" &gt;= %lg\n",b);
    break;
   case 3:
% w7 z5 T9 D% }6 S- j  X  [    printf(" &lt;= %lg\n",b);
% B! ]* U, D! ~, v- Y1 G    break;
0 ]; f+ I' s6 G  }
}
' t; d9 g6 ^# z# i printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
5 m- [% D8 T4 w- f0 F<> tnum_x=num_x;
( i# ?0 ]4 _( E: X5 [* ~2 A& w for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
  case 1:
  case 3:
   num_x+=1;
, K# f6 [* d& F; x   break;
  case 2:
   num_x+=2;
7 B4 g4 C, `7 e% O* C. c- |   break;
  }
; T; B3 s% _1 Q7 Y0 W }</P>
<>//化为标准形式
if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
8 B+ m6 i3 g7 o3 v( m8 } for(i=1;i&lt;=num_st;i++) {
9 y' H% l" W3 o" l8 ?; m2 `  switch(re_st) {
   case 1:
    num_addv++;
    num_ba++;
: r) J4 d: _1 F0 n    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++)
8 S" F4 y) x) l  L* f: F; Q$ W     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
    break;
   case 2:
    num_addv++;
    c[tnum_x+num_addv]=0;
. W4 R) b; B( A( U; {    num_addv++;
    num_ba++;
3 O+ I5 M; z  o# m1 ~# ]) z; d    num_ar++;
    c[tnum_x+num_addv]=0;
    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
. Q8 J( n5 w) q) t& E& Q9 h    for(j=tnum_x+1;j&lt;=num_x;j++)
) F/ Q$ G) n! n+ L/ W" C0 }3 V     if(j==tnum_x+num_addv-1) a[tnum_x+num_addv-1]=-1;
  L  Z2 l7 g3 E2 S* A     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
+ \, q2 n  {$ h- Q" C. D' T     else a[j]=0;
4 N7 Z" h8 \2 M/ m* e" u. Y* ^    break;
* [! q3 S6 N' r3 g+ z" I   case 3:
7 D& k" v4 \' n' \5 }    num_addv++;
    num_ba++;
    c[tnum_x+num_addv]=0;
3 s7 `  Q% s( j3 |    base[num_ba]=tnum_x+num_addv;
0 x5 C6 h; g# K) A# ]- ^    for(j=tnum_x+1;j&lt;=num_x;j++)
8 x' r) J& `7 Q3 M     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
1 K6 y! `: j$ W  G8 f2 B5 X    break;
6 E& n! W& l4 K+ C% r3 [! M  }//switch
}//增加松弛变量、剩余变量、人工变量、确定基变量</P>
<>//显示标准化后的方程组
! |, M) H" L) K& W/ k  v printf("\n化为标准形式后：\n\n");
# c+ ?  C' o  R if(ma_mi==1) printf("max z= ");
else printf("max z'= ");
printf("%lg x[%d]",c[1],1);
& R4 a, g8 s/ _ for(i=2;i&lt;=num_x;i++) {
  k=0;
7 {- `3 M/ B7 q1 S$ V" H# S+ I  for(j=1;j&lt;=num_ar;j++)
  if(i==arti[j]) k=1;
4 b7 c! s4 ~6 \+ I1 t  ?5 h  if(k==1) printf(" -M x[%d]",i);
6 z+ Z3 j; B4 d* \) T  else if(c==1) printf(" + x[%d]",i);
4 E- j  ~2 P  ~+ `% O2 N, ^$ f  else if(c==-1) printf(" - x[%d]",i);
/ @# I# {" X, O0 f  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
5 F2 H% g: X$ a3 U1 b- p" W  else printf(" %lg x[%d]",c,i);
}</P>
$ ?; Z0 i) ~6 k( u" U0 {<> printf("\nst:\n");
$ y2 v4 E1 [" v+ l$ n3 | for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
' c7 v& L& e5 H- m% T- n! I' ]% q6 T   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
) `5 x- |+ ]2 Q3 l; L7 y    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
( I& K, d: C- g8 J6 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);
    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;
   }
  }
' N7 A  _( C, |  printf(" == %lg\n",b);
}
printf(" x[1]~x[%d]&gt;=0",num_x);
8 j, b' z* F- `4 f# n: ^1 g2 M, m}</P>
8 i, Q) P& n3 u<>void iterative() {
int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
5 e, h+ c/ |( w! W4 w, o7 T int base_elem;
int base_out,base_in;
double sigma[MAX],temp;
8 n% y; R: M( d9 x( s double value_be; //高斯消元里保存主元素值</P>
: r9 A, V/ G6 b<> printf("\n\n第%d次迭代：\n\n",stop);
, _: t7 z3 g( ?3 N) f, ]1 X9 C for(i=1;i&lt;=num_st;i++) {
  _2 Z5 y) n, q3 d2 k  printf("c%d=%lg\t",base,c[base]);
* d0 |$ K4 D* |! V% s  printf("b%d=%lg\t",i,b);</P>
<>  switch(step) {
1 |! \4 F! }. \6 k0 q" @6 W! E   case 1:
$ C1 o! x7 w* i; C( _4 P+ I0 n    for(j=1;j&lt;=num_x;j++)
    {
# R  G) M8 i8 P' a) O, ~( K+ U6 b     printf("a[%d][%d]=%lg\t",i,j,a[j]);
; x2 O4 _& P4 [9 V5 N+ O    }
    printf("\n");
: L; v& ]  w) t    break;
" ~4 h  M# R8 S0 Y- ^' @# v* Y* P   case 2:
+ r  f( B. }- J4 a* F  A% ^    for(j=1;j&lt;=num_x;j++) {
5 d4 ~2 G5 S$ D9 R     k_a=0;
     for(l=1;l&lt;=num_ar;l++) if(j==arti[l])k_a=1;
4 j  v/ I! B' W* s. ^( v5 o: U     if(k_a!=1) printf("a[%d][%d]=%lg\t",i,j,a[j]);
: V$ r: L5 B7 z& Q9 V3 j  T( R    }
0 T" s6 {6 c/ e0 [    printf("\n");
/ K) ^: @7 Z/ ?" l+ a' o/ i, _% _    break;
* i2 A6 i3 x( ^1 ~  }
) n+ }7 S, ]% j, w+ i' g+ ^" w8 E }
//求检验数sigma
$ @% v( k  F& E# }( s for(i=1;i&lt;=num_x;i++) {
/ a. V8 k% V- d1 W2 h  sigma=c;
  for(j=1;j&lt;=num_st;j++) sigma-=c[base[j]]*a[j];
2 Q% l. y- Y- d9 U1 h& T; }9 a3 u  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
  switch(step) {
+ k% {6 u9 b6 ]1 d( o/ G6 R  B   case 1:
    printf("sigma[%d]=%lg\t",i,sigma);
    break;
   case 2:
6 }& C  G+ ?' X% G    k_a=0;
    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
# l$ k( @. u9 y; J9 n6 a    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
# F; C5 _( P7 i4 I7 P6 Q" l    break;
( b; v9 L. g5 f% D$ P4 c  }
}
9 v& a1 J8 D5 C8 e putchar('\n');
8 F9 L- N, s1 ]//检验检验数sigma是否全小于等于0
$ `; o: `% z# D4 U* e$ l& O5 R- \ k=0;
for(i=1;i&lt;=num_x;i++) {
  if(sigma&gt;0)
) ^1 ~" x" Q7 i! ?" B7 W2 H  k=1;
}
7 a. |' P4 F1 z+ a) V: q if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
& [0 u0 \- Y3 ^6 l7 `* _  for(i=1;i&lt;=num_x;i++) {
   k_f=k_a=0;
+ \- i$ A# G: t: g/ B' ~  X$ s   for(j=1;j&lt;=num_ar;j++)
3 G* B. U* @6 ?' q) `+ p9 T( u   if(i==arti[j]) k_a=1;
   if(sigma==0&amp;&amp;k_a!=1) {
! a! A* A& r( V& a6 M8 z7 f! Y& _    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
. S  L( Y" J3 D6 G7 i$ h2 w   if(k_f==0) {status=-1; return;}
   }
  Z7 [+ i" X3 Y  }
% c4 v8 m. q4 g: ?' t- ?! y  status=1;
  return;
}
//检查是否为无界解
for(i=1;i&lt;=num_x;i++) {
4 J% G& a- @8 H& z2 z  k_f=0;
; a, _( N. U2 s2 Y, T$ c) t  if(sigma&gt;0) {
7 Y$ Y. \/ {0 S6 a$ p# d0 _   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
   if(k_f!=1) {status=0; return;}
! q) d: ?! Q& P% F) R  }
}</P>
& Q; U% R7 X% W<>//确定换入变量
, z7 @  W8 U1 U' o for(i=1;i&lt;=num_x;i++) {
- s9 i+ m' `1 o) d3 @& e3 A* 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;
* o$ d2 C9 T5 }5 O }//temp赋初值
0 L0 [0 F: e# U; W& H! w+ L for(i=1;i&lt;=num_x;i++) {
  Z* _  _# ~2 g  C# ?) x  k=0;
" p! w1 @; Y. }  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  if(k==0)
  P: ~: G9 O4 Z/ [   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
% ^% ?0 Z4 ^2 B$ i3 x6 R    base_in=i;
' \! G" }0 B2 K0 z9 R    temp=sigma;
   }
}</P>
<>//确定换出变量
0 c7 N9 E: K; k! H/ p* ~ for(i=1;i&lt;=num_st;i++)
6 `1 g; _- c2 R4 R  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) {
$ q9 b! |8 v. V  ?$ G1 N% f& e   for(j=1;j&lt;=num_ar;j++)
6 C' q" M0 p: N5 B: }/ c5 G, _    if(base==arti[j]) {
     base_out=base;
     base_elem=i;
     temp=b/a[base_in];
     break;
    }
  }//人工变量优先换出
  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
0 I2 l/ n: T3 _   base_out=base;
# R! L( V3 `/ M) ~( Q/ w. {9 k! i; p' n   base_elem=i;
   temp=b/a[base_in];
' ^9 t; a- {% P7 j  }
1 ^: G5 q- `4 V0 h: v) M }</P>
<> printf(" 基变量：");
for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
# @. C* f' L8 z* v5 G% T//基变量变换，进行新方程初始化后迭代
. ?% }* R, U# c) W2 M3 L for(i=1;i&lt;=num_st;i++) {
  if(base==base_out) base=base_in;
}
+ u: C. g3 c5 ~' v/ S/ j2 ?//初始化主元素行系数
value_be=a[base_elem][base_in];
& Y" U3 b+ Y; }; `# f# E 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++) {
9 a) G! i! d& c' H3 m  if(i!=base_elem) {
   b-=b[base_elem]*a[base_in];
   value_be=a[base_in];
/ G8 a. O7 _* _  M   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  }
4 ~6 _- i2 Z$ s& r }
8 E# S& |, i) T5 U. e stop++;
if(stop&gt;STP) {status=-2; return;}
iterative();
}</P>
<>void output() {
4 \2 \5 `8 S/ v6 V6 J2 d- o- N int i,j;
0 V' ^* h0 t7 I4 h; O- N8 q! M9 j9 V2 I double X[MAX];
8 O; z/ |4 e5 i* d: ?. i9 _ printf("\n结果如下：\n");
, W. Q& f0 g, U" L7 g6 R6 F* v printf("\nX=(");
4 y6 Q! Y/ }5 d9 i3 T6 k for(i=1;i&lt;=num_x;i++) {
  for(j=1;j&lt;=num_st;j++)
8 t5 N3 o, e9 z) h( P$ U9 D  if(i==base[j]) {X=b[j];break;}
2 O) g; k1 r5 K5 q: {% S  B  else X=0;
  printf("%lg ",X);
2 P1 d$ p" d: c" p }
7 H8 s9 z( E1 v; Z3 J' g4 A5 D* } 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);
}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了