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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
#include&lt;math.h&gt;
* L, j6 V) m3 [3 P#define X 5
$ t6 k7 J) n$ ?6 E#define Y 7
8 ~" v- c8 N" P& o% p8 }void xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
{
, r( R: `/ U. L6 ^4 W+ e' J  int j;
  *c=0;
  for(j=1;j&lt;=*is;j++)
    {
      if((*a)[*ir][j]-*c&gt;0)
    {
      *c=(*a)[*ir][j];
      *j0=j;
    }
" c' n8 K" x" F& ]  Z( k3 t. C     }
1 F/ a6 i0 [9 [# w: |}</P><>      /***************** 参数说明 **********************/
      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
# m+ f+ |/ L0 V4 u1 U# Z! v      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
0 I1 `# v( J0 h3 k4 r% n      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
8 V+ V" t2 a0 j  c      /********* <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],
9 f1 }- i  I2 S. Q% u% F) S       int(*k)[],float(*x)[])
) a- r. X+ b; F7 l6 t/ w' M( M! G{
  int i,m1,mn,j0,i0,j;
  float c,g;
  m1=m+1;
# A& T" u1 ^7 T7 f  j" H  mn=m+n;
  for(i=1;i&lt;=m;i++)
1 }2 I; n& B+ H8 b% o# \    (*k)=n+i;
1 m  r2 W4 m3 I% T, {8 jleap1:
5 |" I/ V, _' G0 g  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  else
     {
leap2:   if(l1-50==0)
7 T1 M" W( b, U# A$ k6 p      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
# e9 b% [1 V4 b       else
      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
     }
   c-=1e-8;
   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");
+ w3 K4 h1 H, }. e2 ?3 i( F    return(-1);
      }
# @$ _+ v4 K' H0 [! D$ W7 N' C( Z   if((c&lt;0)&amp;&amp;(l1==1))
       {
     l1=50;
     goto leap2;
3 L; w8 x" u+ T; i* [% d* y! V       }
$ a, q* m9 {4 D% B+ O% x/ S' J    if(c&lt;=0)
       {
  @( w) R% T: M' Y+ z) }      for(i=1;i&lt;=m;i++)
2 W" p) Z" O9 c/ j+ i% v* o0 [         (*x)=(*a)[mn1];
: e  I& T; ^7 v- D! u9 K4 W! [1 q         printf("\n\t********Optimal solution********\n");
3 l+ [$ z* u4 c- Z! f      for(i=1;i&lt;=m;i++)
0 S# y: J9 D  H  r* a2 R         printf("\ni=%d,x=%f\n",(*k),(*x));
      printf("\nF=%f\n",(*a)[m2][mn1]);
      return 0;
% _' |  s- W" F/ F! {- ^. i    }
     c=1e8;
     for(i=1;i&lt;=m;i++)
       {
7 x( q; O4 g& K% X; U" t& b     if((*a)[j0]&gt;1e-8)
" k  V. _) A6 r, i+ e' V, [        {
          g=(*a)[mn1]/(*a)[j0];
          if(g-c&lt;0);
; a: ], Y6 I7 I0 y+ ~         {
$ S3 A5 {+ J* \7 F+ M4 ]; }7 g# {+ H/ q           c=g;
           i0=i;
% r* w% d7 R1 n$ V          }
. o! w) o5 A, p; Z/ Q4 j         }
5 B7 c7 n7 t5 y) K% e: T! h    }
" K1 M$ m! G4 Z, f  a0 g1 m+ v. z0 f      if(c==1e+8)
* V, I8 D5 P1 V    {
      printf("\n\t*******Lp no solution********\n");
) W, W/ P5 J+ D  q& o* o* o) i      return -3;
5 q; k4 o% X, b  ~     }
       (*k)[i0]=j0;
! Z6 u) k: L% @/ g! K9 ?" b1 O       for(j=1;j&lt;=mn1;j++)
      {
        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++)
# ~" t  h8 [6 I, K           {
         if((i==i0)||(!(l1==1)&amp;(i==m1)))
            continue;
. ?' ~5 E8 J  K; _# ]! H         (*a)[j]=(*a)[j]-(*a)[j0]*g;
           }
) G% Z( s" M$ V9 k' E% O3 D/ v       }
0 y1 v2 W' S8 Z$ ?2 x       for(i=1;i&lt;=m2;i++)
+ _) k1 _0 o8 y6 i% P7 u8 P& S      (*a)[j0]=0;
       (*a)[i0][j0]=1;
' }( _6 w  u8 H) G. ]8 t8 o       goto leap1;
- c. ~6 {/ M' x- M- x2 e1 S6 Y$ x}  
( I/ |$ X, v8 ~% w7 x3 q) z, y- Imain()
{
  float a[5][7]={{0,0,0, 0,  0,0,0},
                 {0,3,-4,3,  1,0,12},
                 {0,3,0, 6,  0,1,12},
                 {0,0,0, 0,  0,0,0},
: G4 i2 B- c, w* v$ g6 z                 {0,69,0,144,0,0,300}};
0 q3 p, S4 r; v  int k[3];
  float x[3];
  clrscr();
; V, K# p* u) f% p! c/ }. d9 o  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
6 r1 n2 r2 ^8 I( E8 J- t}
</P>

plgatc

<>/*************************************************************************
+ \! I/ m% r( `, \                    单纯型法解线性规划问题（两阶段法）        
# C% Y/ l) {" r/ R7 X, n
! f1 s  {- H  w  n2 c% I 编程环境：VC++6.0     
# N) [( Z  i* T9 j) s5 M. L 方程组输入说明：
变量非负，按提示输入相关参数。
9 d$ z+ \* M; I9 ^  J *************************************************************************/
( h; Z0 o4 G" u5 |$ |1 k- g3 Q1 z; ]#include &lt;stdio.h&gt;
5 M) j/ K4 A) q8 `: {9 G#include &lt;stdlib.h&gt;
#define MAX 100
#define STP 100</P>
" {  @! J3 f$ F# ]! D& I$ F<>int stop=1; //迭代记数变量
* m! T% ~1 e/ M; ?int status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
* @1 B5 Y. o  I6 b( bint step=1; //目前阶段</P>
<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
/ E( g2 j/ {' s& N  ^1 Xint num_x; //变量个数
int num_st; //约束方程数
int num_ar=0; //人工变量个数
) _1 @4 q" m6 p+ }1 Mint arti[MAX]; //人工变量下标
% c9 d7 `; P) _. t) F7 z  rint base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
0 n: ?- e4 F4 l: b. ]- o4 G<>void create(); //建立方程组
6 o4 c( x; K- i7 W1 S( \! w8 gvoid iterative(); //单纯型法迭代
* H3 m- o; R2 g* R, a. q& dvoid output(); //输出结果
void banner(); //打印程序标题
0 F% d3 x" p* @1 Ivoid exchange(); //交换两阶段价值系数
void show(); //输出方程组</P>
<>void main() {
int i,j,k;
banner();
3 A3 ?& d9 W/ B& i* h create();
//保存原价值系数，转换为第一阶段价值系数
for(i=1;i&lt;=num_x;i++) {
  k=0;
( S" X* F8 t, b* E' B( `- q  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
# N) Q& Q& L3 t4 ?4 ?/ e  j  if(k==1) temp_c=-1;
4 \. \; Z8 f0 \0 n7 Z  else temp_c=0;
}
0 P" Y7 d6 X! q; X8 ^- G% Q( _ exchange(c,temp_c);</P>
! S# r5 D0 {- _! D  H1 Z1 m: A<> printf("\n\n第一阶段问题为：\n\n");
show();
* L; M* n/ A% i+ d8 N step++;
: N. I' {2 h8 [  v: w printf("\n\n按回车开始第一阶段迭代");
getchar();
getchar();
iterative();
if(status==-2) {
  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
; z7 R0 s! C' r" F: }" s  getchar();
  exit(0);
}
output();</P>
6 T5 |9 {- b1 N, n7 |" v<> if(max!=0) {
+ ]/ I2 i3 E9 x9 t! z) `4 Q  puts("\n\n原问题无可行解。\n");
  puts("\n按回车结束");
  getchar();
, r  F6 d* i! v8 X  exit(0);
}</P>
<> //转换为第二阶段价值系数
" S" T( q& G- m exchange(c,temp_c);
//把人工变量列全设为0
& o- c' U9 Y( x/ p0 L for(i=1;i&lt;=num_ar;i++) {
# \+ w# ]3 i: h4 ^$ I! }  c[arti]=0;
& X" `) b+ g. s& z4 k9 h2 l( d  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
}</P>
<> puts("\n\n第二阶段问题为：\n\n");
show();
puts("\n\n按回车开始第二阶段迭代");
getchar();
3 M. P7 e/ y# x% f" T# k" w$ E iterative();
switch(status) {
case 1:
  output();
; e# N$ S% T% k0 ^4 X' s4 v  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
# \9 p2 ~8 l4 W  getchar();
  exit(0);
+ j' F$ ?3 ~' k- R. m2 r, J case 0:
  puts("\n\n原问题为无界解。\n");
: ]; S1 Z9 J8 t# H, T/ i, E  puts("\n按回车结束");
9 B. C0 k2 R! Q! \: U6 h  getchar();
  exit(0);
( c7 C4 c5 ]+ y! z case -1:
  output();
  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
  getchar();
- I. J% g' G/ V5 G5 d$ ~  exit(0);
3 S3 b4 N" Z  c; a) ^1 C/ j case -2:
5 U9 R, R8 N. j/ o& D% J. M  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
) P( h: g. j. B/ ^+ B- _, ?  getchar();
  exit(0);
- i2 V* _# K* n }//switch
( F6 [; k+ V0 q; p6 L; i9 }3 b: S' s}</P>
<>void banner() {
% q4 ?8 I8 E# d# _; g  Y printf("\t\t****************************************\n");
7 H5 M% h5 h" i* e printf("\t\t         单纯型法解线性规划问题\n");
printf("\t\t                         作者：Thunder\n");
8 B! j9 A1 Y9 u7 w printf("\t\t****************************************\n");
printf("\n");
; x5 w. k" \1 Z, O8 ]$ D}</P>
<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
3 i4 q, g! }2 _  D1 k: ]8 S% ]//为1的不显示系数1，-1系数只显示负号
int i,j,k;
switch(step) {
4 D' r# P* e. f case 1:
% l0 E/ K7 [; J9 T  printf("min z= ");
: n+ n6 V# O" m* m/ Q" P6 `  printf("x[%d]",arti[1]);
4 j& p3 N. |$ u0 s# @0 }. Z: H- H* Y  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
case 2:
: f, U! O+ I3 T# c6 B* u% {9 c printf("max z= ");
- x( Q, q4 O. X8 d0 I printf("%lg x[%d]",c[1],1);
# W1 F, Y6 r3 `- ~3 K for(i=2;i&lt;=num_x;i++) {
& R. T" s$ f9 ~" V  if(c==1) printf(" + x[%d]",i);
1 s: J4 B- |" D% p. G' y2 L' p  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);
% c7 b+ P( Y  {! N8 i: V" d% { }
break;
}</P>
<> printf("\nst:\n");
7 a) G+ |7 H$ }5 s. i5 J% X for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
; J/ ~1 b3 x0 z: i4 g/ j0 u   if(a[j]!=0) {
4 C3 y- K$ t) p# x9 {$ r. v) B    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);
/ M1 V  w' l1 Q  w" a    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);
1 a0 ]4 W* M; O" x( t0 D( _    else printf(" %lg x[%d]",a[j],j);
    k=1;
   }
4 N7 e# c9 j1 y( x  }
7 p  L; a, F: ~$ V1 Q  printf(" == %lg\n",b);
}
printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
<>void exchange() {
5 `+ u# N( f6 f- R" J% w int i;
double temp[MAX];
' [1 T, b' G3 G5 R: y# u for(i=1;i&lt;=num_x;i++) {
  temp=temp_c;
& f: Z+ A0 b9 h3 Z5 f" k  temp_c=c;
0 }7 H3 n9 \- c& i4 e; k) E  Z  c=temp;
}
. |3 a! U9 T; |+ _/ s  U6 b1 A" n+ Q}</P>
, l- g# ]+ o+ f& Y4 _) K<>void create() {
//输入方程组系数，每个方程输完后回显确认
/ x: p% G  Z: S* W  P7 u; a int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
8 o* j  }/ g5 b0 F; X' i char confirm;
; W9 A  |, X) Q, u7 ?
while(1) {
- p% R4 K" u! @8 w. y5 k6 }; [  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  scanf("%d",&amp;ma_mi);
- S8 C- h0 b) N; T( Z  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
  else break;
1 V# j2 a$ O" q$ w5 F2 ? }
9 P3 F* ~9 r, J1 K$ t% H" g- j
while(1) {
$ j8 N$ X( j' D  C9 y/ k: z  L  printf("指定变量个数：");
/ f! R2 s$ q, `8 Z  scanf("%d",&amp;num_x);
8 ~$ d$ n' @+ B2 i1 ~1 O  printf("输入价值系数c1-c%d:\n",num_x);
  for(i=1;i&lt;=num_x;i++) {   
   printf("c%d=",i);
, F9 a) t: H3 T9 }$ L' q   scanf("%lf",&amp;c);
  }
: k8 e7 C, f* P  if(ma_mi==1) printf("max z= ");
  else printf("min z= ");
  printf("%lg x[%d]",c[1],1);
. k9 z) |6 l+ ]: p4 [. h5 [  for(i=2;i&lt;=num_x;i++) {
   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
; e; L1 _9 T" B% H8 W   else printf("  %lg x[%d]",c,i);
  }
printf("\n正确吗？：(y/n)");
getchar();
7 O+ Q- ?5 y2 M9 P confirm=getchar();
5 |0 k- E5 a6 G: w* w if (confirm=='y') break;
else if(confirm=='n') continue;
2 G7 o8 w/ |) h5 _- w4 N }</P>
<> printf("输入约束方程组个数：");
scanf("%d",&amp;num_st);
for(i=1;i&lt;=num_st;i++) {
  printf("st.%d:\n",i);
  while(1) {
   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
, H* o( {+ Q' ]& o   scanf("%d",&amp;re_st);
   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
   else break;
! s' F5 r% A6 t  }
  printf("输入技术系数：\n");
+ [7 i% M$ _; n8 o  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
( T# m4 G7 l/ o5 i' p' F# e   scanf("%lf",&amp;a[j]);
  }
! L) E- |- C% L2 d1 F, H  printf("输入资源拥有量：\nb%d=",i);
0 J2 u% g+ ]9 v4 R7 t" W% a  scanf("%lf",&amp;b);
8 p/ _5 `) c6 O& k4 P7 q# A  
2 y3 w) R/ S+ N" C+ z( R  printf("st.%i:\n",i);
8 H+ {5 @( n% ]0 Y5 X  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);
   else printf(" %lg x[%d]",a[j],j);
  }
8 y" p0 {5 k! b# l/ ?  switch(re_st) {
   case 1: printf(" == %lg",b); break;
   case 2: printf(" &gt;= %lg",b); break;
   case 3: printf(" &lt;= %lg",b); break;
+ u/ ?# {* e3 h$ b* k. @0 U  }</P>
2 y4 N3 X* N) i- T* k, r- M9 E<>  while(1) {
1 z# c! b6 u; W   printf("\n正确吗？(y/n)");
   getchar();
% v3 ]  j$ t3 M" _   confirm=getchar();
   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
  }
# `, B$ v2 [; j+ e" g }</P>
<>//显示输入的方程组
6 F3 u. u! i1 C7 B& h  ` printf("\n原问题为：\n\n");
' N: Z, v9 Q; g, Q2 |- M4 T if(ma_mi==1) printf("max z= ");
6 @  ]% m0 Z1 V  u: R5 n* ]7 v else printf("min z= ");
  d1 c; A8 H0 f$ k  l% |) R/ ^ printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
1 O7 f- a9 c' p5 h* G% d" F  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! Q- N0 V# H7 y1 q5 x3 P0 r }</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  k=0;
- r; ]- X7 s: @  n  for(j=1;j&lt;=num_x;j++) {
5 x* O8 D4 f9 _6 ^, D   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
4 w" D! h& i5 {4 |! r4 q3 }; c    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
    else if(a[j]==-1) printf(" - x[%d]",j);
  z  v9 [4 p2 U: {- v+ z' ?) R    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);
# w; m6 z9 x- p  x3 S: {0 z    k=1;
   }
& |' z1 J  z9 B/ j  B3 M  }
! ^7 o! K* ?8 T% t4 ^  switch(re_st) {
7 \  ]0 ^. }3 X3 H( B+ T   case 1:
    printf(" == %lg\n",b);
! f4 G! T/ ~* t, G4 |8 A$ X' z    break;
8 q5 i; l4 V% @   case 2:
( R4 k1 h2 k* `! ~    printf(" &gt;= %lg\n",b);
    break;
   case 3:
/ X5 {1 m7 R$ N4 a7 C  r) n    printf(" &lt;= %lg\n",b);
* A  |2 i" G# T. \9 }    break;
  }
; U3 c+ H% ^! r% n: A9 U8 A) w }
+ g# B: F6 M- p$ l5 ^& J printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
<> tnum_x=num_x;
for(i=1;i&lt;=num_st;i++) {
7 C; Q9 ?2 A- n  switch(re_st) {
+ `  J8 K/ h/ o  case 1:
& h. T6 u' |: q2 N6 o9 p4 s+ p  case 3:
   num_x+=1;
   break;
  case 2:
3 `/ a8 I+ _/ d% E* D   num_x+=2;
& Q' _4 ^' J& w; Z1 `   break;
% [- E, A' c4 c$ }; V  I  }
" F+ O) H5 K7 s! ~- \; Z" \ }</P>
<>//化为标准形式
0 J5 D1 y9 ]+ x: b5 n: | if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
for(i=1;i&lt;=num_st;i++) {
3 g' T  ?& y8 @, I* E  switch(re_st) {
; {" _+ o& x3 v8 U  n7 U' Q* Y   case 1:
! B% l& J& Q$ D3 P5 b, Q" E+ s9 Q/ M    num_addv++;
# q0 \8 t" O4 x( C    num_ba++;
    num_ar++;
7 r: T8 H  P( m* r    c[tnum_x+num_addv]=0;
6 l5 D9 |' D& k+ y, I& O    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
    for(j=tnum_x+1;j&lt;=num_x;j++)
" b9 v& \7 E: N. z8 B6 h+ K, m     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
    break;
/ y# o" o  r$ P9 p  _; Q* }   case 2:
' v3 ?8 N$ Y2 n7 O) \    num_addv++;
    c[tnum_x+num_addv]=0;
    num_addv++;
    num_ba++;
/ P& e. q/ ~1 D/ A8 c    num_ar++;
8 E$ A1 C0 i3 ~% W9 X* h0 @    c[tnum_x+num_addv]=0;
9 z! E) L" c% R6 N    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;
8 k  Y$ B7 V) C# ?8 ]     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
7 C, W' y  K( t2 r; L' d, D    break;
' k9 D8 ^; R' a' Q4 ], `   case 3:
    num_addv++;
    num_ba++;
    c[tnum_x+num_addv]=0;
/ o1 c4 q% Q0 m2 w    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;
( p* t$ p- Y7 f! C7 o1 \" L! H     else a[j]=0;
# y' h2 m8 @. @: i( M! ?4 l    break;
  }//switch
}//增加松弛变量、剩余变量、人工变量、确定基变量</P>
<>//显示标准化后的方程组
. g4 u( d9 c9 {4 m! c& Z printf("\n化为标准形式后：\n\n");
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++) {
  k=0;
  for(j=1;j&lt;=num_ar;j++)
3 F* k0 x, o5 R# d  if(i==arti[j]) k=1;
  if(k==1) printf(" -M x[%d]",i);
  else if(c==1) printf(" + x[%d]",i);
- W" v+ @4 N: r% Y2 s+ h  else if(c==-1) printf(" - x[%d]",i);
  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
; x0 F7 V8 ^/ H, T  else printf(" %lg x[%d]",c,i);
5 W7 }; |+ a3 Z, q }</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
1 K4 |$ O0 a/ W- _7 |1 z8 O  k=0;
  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
, q. Z3 m5 ]: ]8 S/ X. V    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);
' {; w5 `6 b! z: _: `7 I    else if(a[j]==-1) printf(" - x[%d]",j);
$ Y4 Y4 S& U  o0 S3 p3 Y$ f) R% a4 k    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);
# @, m1 Z  b! V3 B& b    else printf(" %lg x[%d]",a[j],j);
    k=1;
   }
' t; u+ M7 M# q1 z4 O4 A  }
  printf(" == %lg\n",b);
$ L& D1 t+ ^5 l, S( e: W7 c }
printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
! _, P5 C& w$ y. R6 k<>void iterative() {
% |  t- R+ [$ S2 A( 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;
double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
8 i2 z$ K, j- k4 r. l0 g+ W- Z<> printf("\n\n第%d次迭代：\n\n",stop);
for(i=1;i&lt;=num_st;i++) {
6 }4 m7 z5 c4 L" m0 U  printf("c%d=%lg\t",base,c[base]);
  printf("b%d=%lg\t",i,b);</P>
<>  switch(step) {
   case 1:
    for(j=1;j&lt;=num_x;j++)
/ S( `+ M* ^/ E9 O+ J3 B    {
- |  X+ t' X9 l/ K7 v. m2 z* X% f" v     printf("a[%d][%d]=%lg\t",i,j,a[j]);
; g, K3 f8 @* |! @" \. o) d    }
9 x& Z! I' h. X( g    printf("\n");
% ?* K2 b! z! ~0 V. X    break;
1 |6 Z) d/ U2 T$ m5 p   case 2:
5 @5 b; I& h$ l: g' U# |    for(j=1;j&lt;=num_x;j++) {
+ c. E9 N! J5 L8 n' F     k_a=0;
     for(l=1;l&lt;=num_ar;l++) if(j==arti[l])k_a=1;
  z& V- _0 \9 q; c2 q     if(k_a!=1) printf("a[%d][%d]=%lg\t",i,j,a[j]);
. M1 T6 k( V4 k3 p    }
    printf("\n");
    break;
0 K+ W) S6 H, S! Y4 b' A: p  }
7 ], P$ W$ j, P6 d: x  `9 |$ V }
//求检验数sigma
# ~8 H5 T5 y8 S% \6 d# ]: T for(i=1;i&lt;=num_x;i++) {
  sigma=c;
  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:
    printf("sigma[%d]=%lg\t",i,sigma);
    break;
) _! t# Y7 S4 S) d/ }   case 2:
% S  w" K+ v, P; E$ ]# k' w$ U    k_a=0;
) S, F; I, N/ f  c    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);
    break;
+ Q6 ]/ ?$ X% ]  }
}
putchar('\n');
% i- U7 \- ^' o1 ~3 L5 N4 i  q//检验检验数sigma是否全小于等于0
% H! k# s! A$ B k=0;
) I+ m4 g! t+ M) Y- \6 C* s! h0 H for(i=1;i&lt;=num_x;i++) {
  if(sigma&gt;0)
  k=1;
}
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) {
    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
$ C9 N$ g5 m8 j: v. K   if(k_f==0) {status=-1; return;}
# p( u: G  O8 G1 L  v   }
& h8 I- r: M; ]6 B  }
" d- p# f5 a& q  c7 j* _5 _  status=1;
; n5 |- e% {3 U( g3 j  return;
}
//检查是否为无界解
2 l" X$ w" K- E9 L/ R2 r' R% m for(i=1;i&lt;=num_x;i++) {
; }& }! I( d4 ?; @- v( |( @$ W( ~  k_f=0;
: c* |; G2 k" e; m  if(sigma&gt;0) {
   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
! F+ i1 G8 p/ l, g& R   if(k_f!=1) {status=0; return;}
  }
1 V, {9 c+ ^2 H4 _: J! Z) E }</P>
  ]/ p5 Y, s0 a: J0 x, y4 p<>//确定换入变量
/ r& w: a7 J! g* C7 A for(i=1;i&lt;=num_x;i++) {
  k=0;
. `6 Q; I7 U3 r  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
- _3 n/ A3 `6 V+ p; t6 H  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
}//temp赋初值
* Q# F  b; R; K0 e$ _ for(i=1;i&lt;=num_x;i++) {
  k=0;
+ {9 Z  t1 y1 `$ ]" ]  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
# f, \% S$ P9 @: e; u  if(k==0)
   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
    base_in=i;
: ~; n2 Y& N0 A5 V" ~$ Z4 H    temp=sigma;
& Z* H2 J7 c$ b' y$ n' @   }
}</P>
. P) D4 R4 H4 w& Q" F8 s4 G: ~<>//确定换出变量
for(i=1;i&lt;=num_st;i++)
7 v7 s5 }! Q  ?: x' ~  if(a[base_in]&gt;0) {
   temp=b/a[base_in]+1;
; v/ v( n) o/ ]6 D, C% e   break;
  }//temp赋初值
4 Z% ^; _% c6 S7 z for(i=1;i&lt;=num_st;i++) {
9 K6 x( e9 Q; Q2 ]5 k8 ?  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
& B% J2 X0 E. Q0 b) |4 a# D. w0 c   for(j=1;j&lt;=num_ar;j++)
& Q. a  ?8 x! G4 Y    if(base==arti[j]) {
     base_out=base;
, U# V* i2 S/ m/ D! X+ y1 s     base_elem=i;
     temp=b/a[base_in];
     break;
) e5 H# _' m$ z4 s    }
  }//人工变量优先换出
. V! c, n2 A1 q! L! |3 L: y' m' s7 \/ r  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
, r, @  Y" b" R. \8 X9 `# f   base_out=base;
   base_elem=i;
   temp=b/a[base_in];
  ?5 I# k, g4 I9 C  }
}</P>
( J1 v8 {& r- O0 {<> printf(" 基变量：");
# i+ q5 r; K0 A1 Z  K# [ for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
8 C6 A! u) ?5 h7 Q' R1 e  x printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
//基变量变换，进行新方程初始化后迭代
for(i=1;i&lt;=num_st;i++) {
  if(base==base_out) base=base_in;
}
* y8 G0 J5 ^# x//初始化主元素行系数
value_be=a[base_elem][base_in];
b[base_elem]/=value_be;
7 \. r! W0 k3 `' @9 \9 [- Z5 C for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
4 s; D( p7 l3 w% K& l  C<> for(i=1;i&lt;=num_st;i++) {
8 Q- m7 X% N( U; r  c  if(i!=base_elem) {
1 R3 I: |: t+ E1 O# k$ U  |5 `   b-=b[base_elem]*a[base_in];
   value_be=a[base_in];
: K6 J" X* V. z" P* T2 |   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  }
}
stop++;
0 L* w6 N9 X: z* @ if(stop&gt;STP) {status=-2; return;}
1 p0 n) r8 Z1 _3 G) \: G& q iterative();
}</P>
9 ?8 ?2 @2 R, I7 @8 h7 ]1 V<>void output() {
int i,j;
double X[MAX];
printf("\n结果如下：\n");
printf("\nX=(");
for(i=1;i&lt;=num_x;i++) {
0 k- f0 B8 I7 h6 @, @  for(j=1;j&lt;=num_st;j++)
  if(i==base[j]) {X=b[j];break;}
  else X=0;
  printf("%lg ",X);
) ^5 P/ f3 R  V" A. q }
2 N! @' J# G% _) E' q printf(")");
3 T+ S/ H$ f. B- _# a6 x$ _ for(i=1;i&lt;=num_x;i++) max+=c*X;
1 S# O. p7 m. a if(ma_mi==1) printf("\nMax z= %lf\n",max);
else printf("\nMin z= %lf\n",-max);
" {! R8 [1 v  j4 K* T}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了