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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
+ V' m8 g# O: x3 {% y6 M) K#include&lt;math.h&gt;
  C' v, S5 J) S#define X 5
! L- `5 @; z2 m* F#define Y 7
void xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
% S* J/ _% e! M6 k5 b) P4 \/ O{
& h% U6 C$ w5 _  int j;
  *c=0;
  for(j=1;j&lt;=*is;j++)
    {
      if((*a)[*ir][j]-*c&gt;0)
: J# W$ X. h" L: S& [    {
5 }; h0 Q! X! N0 }# O: m. F      *c=(*a)[*ir][j];
      *j0=j;
    }
     }
}</P><>      /***************** 参数说明 **********************/
4 Y( v- R1 S' P$ O+ X: V5 A0 H      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
0 n; [0 w" X% i/ t( @      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
0 A" {" s% r+ X8 i9 Z& z$ O7 D      /********* <a href="http://happyyangxu.home.sunbo.net*******/" target="_blank" >http://happyyangxu.home.sunbo.net*******/</A>
2 B/ U2 }/ c! ^% f+ T7 Zint xi_sm(int m,int n,int m2,int mn1,int l1,float (*a)[X][Y],
       int(*k)[],float(*x)[])
4 N. T8 e) `8 j% k6 z! V  @{
9 e0 Q6 J+ k( Z" H; `4 x! |8 V' }  int i,m1,mn,j0,i0,j;
  float c,g;
' f" ^; ?; V4 K  a  m1=m+1;
  mn=m+n;
  for(i=1;i&lt;=m;i++)
8 ]/ Z* s6 G, W$ q0 q  g: r    (*k)=n+i;
leap1:
4 R, a1 h6 y+ ^) }2 m7 n  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)
' N2 q4 o% i' \; A( V2 O" B      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);
5 P2 m4 f- @! \: o+ O7 |- Q( }; f     }
   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");
    return(-1);
" }3 k5 I# u. e' r& p3 U7 }      }
   if((c&lt;0)&amp;&amp;(l1==1))
% n2 L# _* [# s7 n2 L4 J3 k' z       {
     l1=50;
     goto leap2;
       }
    if(c&lt;=0)
       {
' e. v! \: H! @! \- Y9 G: Q# c      for(i=1;i&lt;=m;i++)
         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
! r- q5 j% N/ @0 }/ Y      printf("\nF=%f\n",(*a)[m2][mn1]);
      return 0;
& ]) U, o5 ^3 R# n* U4 @( z    }
7 ]% I  ~9 Q- r- w; ~     c=1e8;
! t, q; G# k+ T5 {' F8 S+ j" k) L     for(i=1;i&lt;=m;i++)
9 t/ s* O. X; s1 ]6 E- B       {
     if((*a)[j0]&gt;1e-8)
        {
          g=(*a)[mn1]/(*a)[j0];
          if(g-c&lt;0);
         {
  I+ V8 U$ u+ P! G) K% O% o4 _! o           c=g;
, |+ {9 ~) ?6 u  U: x+ H           i0=i;
          }
1 f/ L! m% e1 {9 e; o         }
    }
8 D! `/ B! ~2 s: R  ^, m( S" \) Q      if(c==1e+8)
    {
# ^) {3 C" D( v+ ~      printf("\n\t*******Lp no solution********\n");
8 X9 k' u# O( z: h3 o      return -3;
* J; a$ H2 F. g     }
- b! Z: A6 C( E       (*k)[i0]=j0;
- q. U+ e  K5 G: R$ V% c       for(j=1;j&lt;=mn1;j++)
: Y; v' ]' W$ }- H, e4 |      {
6 H! `! @, b. O        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
           continue;
         g=(*a)[i0][j]/(*a)[i0][j0];
5 M  t1 x) a' y9 ]           (*a)[i0][j]=g;
3 x8 r, U/ R% M6 q         for(i=1;i&lt;=m2;i++)
           {
% V7 R6 d% @0 l2 K- y( Y, p  Z         if((i==i0)||(!(l1==1)&amp;(i==m1)))
% _: V1 R7 B1 b4 k            continue;
         (*a)[j]=(*a)[j]-(*a)[j0]*g;
0 `& u- ?. `* ?9 C2 C( ]  d( @           }
# \2 b/ p$ |, i8 X. F       }
       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
       (*a)[i0][j0]=1;
! K, `, ^  [. l4 p5 ^       goto leap1;
}  
- c# y' H, x! 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},
5 Y% g! i" m( v* b% u1 {0 b2 D                 {0,0,0, 0,  0,0,0},
5 j7 S/ M; U$ U( c9 I% j! p: }                 {0,69,0,144,0,0,300}};
5 p/ Y' M% k$ r# \  int k[3];
+ t6 u  K% h% g! ~; S  float x[3];
  clrscr();
, A3 Z8 F3 S+ r7 g  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
  Q1 p0 O7 B8 P2 O}
</P>

plgatc

<>/*************************************************************************
                    单纯型法解线性规划问题（两阶段法）        
0 y* ]# k: q% p9 m7 e- o
) }0 m% T( `* ~- ~ 编程环境：VC++6.0     
方程组输入说明：
, o9 ~" S/ u! y' c: m  ]0 \* L 变量非负，按提示输入相关参数。
*************************************************************************/
' S) k3 m$ c. K( H5 D  r; H: A' r#include &lt;stdio.h&gt;
#include &lt;stdlib.h&gt;
#define MAX 100
#define STP 100</P>
<>int stop=1; //迭代记数变量
int status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
int step=1; //目前阶段</P>
. y" F. m1 Y2 M) y# t7 V8 p<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
int num_x; //变量个数
int num_st; //约束方程数
int num_ar=0; //人工变量个数
+ W6 b8 T( ?4 y3 e) o( ~int arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
/ [7 H! a$ D- u8 I* o4 C' Aint ma_mi; //1为求最大值，2为求最小值</P>
<>void create(); //建立方程组
void iterative(); //单纯型法迭代
void output(); //输出结果
+ L9 W+ g5 b2 U' a& ]3 g4 l% }void banner(); //打印程序标题
/ C2 v9 |/ b7 B$ ovoid exchange(); //交换两阶段价值系数
( k5 f6 o( @4 a  y$ |6 D4 Bvoid show(); //输出方程组</P>
<>void main() {
int i,j,k;
banner();
create();
, R( N, \7 u' U+ S8 e/ \$ A, V //保存原价值系数，转换为第一阶段价值系数
' g& r1 g, V) P( F4 s for(i=1;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
  if(k==1) temp_c=-1;
  else temp_c=0;
6 ~0 c0 S0 p& t, i! | }
; v" ^0 [# R' k8 f0 v; J9 W0 `$ t exchange(c,temp_c);</P>
8 Z# ?2 F. K! K7 k9 i% d<> printf("\n\n第一阶段问题为：\n\n");
( [. G# j- F( m* a3 W. \: h show();
; H: r9 _) z+ r! s* z step++;
( S& o; a  Q# W5 J! e/ D printf("\n\n按回车开始第一阶段迭代");
3 S1 x& G/ j$ Q9 x getchar();
getchar();
7 r$ q$ Q# P# ~$ Z iterative();
if(status==-2) {
  l# X& X2 K* a1 g9 ]# c  puts("迭代超过限制次数强行终止！\n");
9 d" u3 O& _9 |* j7 }) q1 Z4 A/ {7 G  puts("\n按回车结束");
  getchar();
7 `/ w/ b6 W% w" ~0 G7 b  exit(0);
}
output();</P>
<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
}</P>
, @" z1 T5 s1 v; A1 ]: Q* T<> //转换为第二阶段价值系数
- f) x% Y/ m  X0 I$ t) D exchange(c,temp_c);
//把人工变量列全设为0
for(i=1;i&lt;=num_ar;i++) {
$ b! m8 g! J+ w- {  c[arti]=0;
1 I/ A, x( t0 z8 S  `: p) }4 g  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
}</P>
<> puts("\n\n第二阶段问题为：\n\n");
show();
; F" _$ b2 \/ @- P; C/ p6 b puts("\n\n按回车开始第二阶段迭代");
2 |! u! I8 {1 H  c0 G getchar();
iterative();
switch(status) {
case 1:
  output();
  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
" b. c9 `" _+ t2 R  getchar();
  exit(0);
1 N: m: C8 z1 P case 0:
  puts("\n\n原问题为无界解。\n");
% `$ G5 {# n: T2 Y( T5 B' z9 V( _  puts("\n按回车结束");
( x: I# `, O4 z  H. D+ f  getchar();
8 D: i+ a" R9 H' L1 H& W+ P  exit(0);
  Y7 Y5 A) [1 J% Q' K case -1:
  output();
  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
* I8 H+ @0 i/ Z4 u: E- {  getchar();
  exit(0);
; }1 ~4 h# B+ p( j: w& E case -2:
2 Y; }% t6 q' d1 n  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
, f# b6 w, c' R  getchar();
0 F+ c( v/ g' q2 z3 V- M  exit(0);
}//switch
}</P>
7 g4 x8 v9 K/ P1 Z<>void banner() {
printf("\t\t****************************************\n");
printf("\t\t         单纯型法解线性规划问题\n");
printf("\t\t                         作者：Thunder\n");
printf("\t\t****************************************\n");
printf("\n");
}</P>
<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
7 w) Y3 ]) \5 k% v//为1的不显示系数1，-1系数只显示负号
! ]) V6 w+ I/ }8 j; ? int i,j,k;
switch(step) {
case 1:
% k5 T) m2 e' x5 X* ~! G# F4 k  printf("min z= ");
- ]- X0 J5 f) {/ J, N5 P( s3 Z  printf("x[%d]",arti[1]);
  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
/ x2 k, q! Y8 h; c4 S$ k' i5 Q  break;
% Z4 Y* c- D2 y; ~# w+ B9 H; e case 2:
printf("max z= ");
printf("%lg x[%d]",c[1],1);
% `- ?5 I% f! s, D6 l) l) k for(i=2;i&lt;=num_x;i++) {
& k& V9 q( T8 S! @$ q$ [+ g3 M& U# g! t  if(c==1) printf(" + x[%d]",i);
+ L8 i' o0 w/ r3 n+ e/ M  else if(c==-1) printf(" - x[%d]",i);
- P0 W  }9 q2 D  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
6 @3 J: ^# D3 d  else printf(" %lg x[%d]",c,i);
- J; B$ H7 o4 S6 E( ^4 i3 u }
+ c! E, f( |) d) x) W5 q break;
4 @$ U5 ~: f2 I8 \ }</P>
1 x& ]! I0 m  {3 L2 I, \8 z0 ?<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  k=0;
7 t) q. N$ v) W% p, n* ?  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);
$ P" W/ e1 n% j  T0 W" E3 i    else if(a[j]==-1) printf(" - x[%d]",j);
  i2 g- e& A* N) ~# M3 @( b! o    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;
   }
# p9 b  t8 P' m; o$ w# r+ C1 n  }
  printf(" == %lg\n",b);
4 V9 `; ~, A+ o- y, Z }
printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
* U6 ^. _9 f1 g1 L, Z$ q<>void exchange() {
. M1 w2 U& n, _5 ^2 e int i;
double temp[MAX];
8 x2 i, M3 M; J, |8 ~' s1 Q for(i=1;i&lt;=num_x;i++) {
& S! o; }( t9 H6 Z; Z& _" Z3 }  temp=temp_c;
( O/ i1 |1 \% a1 l4 P, Y  temp_c=c;
5 p# \& W. g  u3 E  c=temp;
" j- H- s% T; Z9 U$ R }
' p* W- Z+ ~) l3 e* c8 x}</P>
6 \2 r6 f8 C( I+ N<>void create() {
//输入方程组系数，每个方程输完后回显确认
8 S9 c5 m7 Z7 x- w 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("输入错误，重新选择。");
6 w' ^( ~( u7 _( j" H- A$ m+ @) }  else break;
% E7 l6 {9 b8 ` }
5 ?$ P) h. o2 g' |9 ^
2 t. ]( i0 A% J while(1) {
  printf("指定变量个数：");
- H$ ^8 I' w1 r0 n2 q( l  scanf("%d",&amp;num_x);
1 T( ]" Y- ^# v7 U  printf("输入价值系数c1-c%d:\n",num_x);
  for(i=1;i&lt;=num_x;i++) {   
- Y0 k$ K! r# x7 d! p   printf("c%d=",i);
  l4 t; P8 z; P   scanf("%lf",&amp;c);
  }
% x4 h2 d: w( Z7 _( Z, P* h  if(ma_mi==1) printf("max z= ");
  else printf("min z= ");
0 ]" ^5 g" c7 m8 q1 H  printf("%lg x[%d]",c[1],1);
4 U7 [, f% l1 {( b! f  for(i=2;i&lt;=num_x;i++) {
6 |( |: [# E# s* g+ H. i   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
! Y# {6 ^5 R  Q   else printf("  %lg x[%d]",c,i);
  }
& i8 v. a/ p' _" M! N/ z. D9 J printf("\n正确吗？：(y/n)");
) M/ A+ D6 E* u+ C8 F7 i$ L getchar();
; b% ?# f" s$ R  i7 W confirm=getchar();
8 V- ]! b) D5 E8 X& h6 O& @ if (confirm=='y') break;
else if(confirm=='n') continue;
}</P>
# @! }$ o4 J: a3 s$ Z" y<> printf("输入约束方程组个数：");
scanf("%d",&amp;num_st);
. c+ K) r" z3 n, R) k3 q for(i=1;i&lt;=num_st;i++) {
3 k& f) L" @# D! y  printf("st.%d:\n",i);
  while(1) {
1 ?/ y! y# A  g, r) @& }   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
# T  N3 {3 y5 f( t; R! E   scanf("%d",&amp;re_st);
! @# ^* J* f% a2 H1 z; e1 ^/ R   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
   else break;
  }
- v; t; s0 N8 Z/ I8 H. ]2 M  printf("输入技术系数：\n");
6 p. \+ {/ P& y' s3 x  for(j=1;j&lt;=num_x;j++) {  
) |; @; b) _& {# i2 n* Q6 j- v/ {/ d   printf("a%d=",j);
, ]- U# X  d; u8 d$ G" J   scanf("%lf",&amp;a[j]);
+ X9 |( c! c6 t( z8 @' |  }
6 ?0 I! [- k1 |  V& [3 l  printf("输入资源拥有量：\nb%d=",i);
: o8 L  Q1 }7 z- @9 \* k! x$ \  scanf("%lf",&amp;b);
  
  printf("st.%i:\n",i);
: z* i3 N' P6 \9 Q1 a  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);
8 i! {, Q: p7 P/ `8 r1 _   else printf(" %lg x[%d]",a[j],j);
& g: \4 }; X* E  }
  switch(re_st) {
% C& E2 p( w# f  n3 y$ a5 S! |# t% ^   case 1: printf(" == %lg",b); break;
   case 2: printf(" &gt;= %lg",b); break;
4 n( r/ O& |8 D% k9 j" o* x   case 3: printf(" &lt;= %lg",b); break;
  }</P>
8 A3 W9 s- L! W<>  while(1) {
1 [4 `9 s( P" g, X3 e% z. v   printf("\n正确吗？(y/n)");
- D6 J) @( b$ a) h: C   getchar();
   confirm=getchar();
" M  f* ?/ m7 `   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
  }
  m/ c3 _3 L5 R  U/ u2 U# V }</P>
2 K! l. h# d2 p* s<>//显示输入的方程组
2 @3 I* T6 T3 j& f( s, X& k printf("\n原问题为：\n\n");
if(ma_mi==1) printf("max z= ");
else printf("min z= ");
& d6 b6 e/ E& W1 W7 k- y printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
2 T+ ]8 f$ W0 F9 d8 {  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>
& U9 E/ |" h/ t5 X<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  _4 N  B2 b/ _- T  k=0;
2 |* Y7 E; Q: N) P0 f/ [& K' {  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
, w9 O: o, `# F4 R6 f( 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);
! c7 U) ~7 @, u: h! B    else if(a[j]==-1) printf(" - x[%d]",j);
/ H, x% z  X. u4 ~" l- S    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;
6 j% X9 ]/ D  f   }
  }
( i: E) p6 X6 H3 e- O  switch(re_st) {
   case 1:
    printf(" == %lg\n",b);
0 R, g2 E+ E+ e# z: \- t    break;
   case 2:
$ E: m: a  H# E    printf(" &gt;= %lg\n",b);
: h& N5 m# ~  l0 ~+ C) p) e    break;
   case 3:
# L( C0 |# K8 l  S    printf(" &lt;= %lg\n",b);
    break;
  }
1 H! h7 L9 J! e. i( D, f }
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) {
. S% d# U7 x1 F: P  case 1:
& E* u0 A, p6 L( P6 d+ D  case 3:
   num_x+=1;
   break;
# x, U. Z1 V' d. F  case 2:
' ~3 l  T) H3 F9 Y! o# X: t   num_x+=2;
% w  S3 e& N# A" Y5 t4 P) k, A6 x   break;
( A0 h: V! k4 y# \  }
}</P>
<>//化为标准形式
: H+ S6 D4 v. Z6 P; q( i$ E3 g if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
6 [; t8 W5 R5 O. _8 E for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
) S* r- A& {) E3 g+ C+ E   case 1:
    num_addv++;
    num_ba++;
% C( \4 V6 c. Z6 P1 w, I    num_ar++;
    c[tnum_x+num_addv]=0;
$ N) _/ Q: U. c8 _9 M    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
* k  L% P: x4 ?. ]  ^, Q9 k    for(j=tnum_x+1;j&lt;=num_x;j++)
     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
$ C, \) T5 A: Z7 q# ^1 J9 V     else a[j]=0;
    break;
9 H* M! w3 P) C6 p/ ]" l& q   case 2:
/ o, q- n4 ^# z  @# f- V! |* f- q    num_addv++;
    c[tnum_x+num_addv]=0;
! X6 {, O$ R+ \: `$ K+ c2 r    num_addv++;
, u& p+ P1 A7 a    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-1) a[tnum_x+num_addv-1]=-1;
9 @& K! |8 p) t     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
8 T1 A! C; M- |+ f/ a( M  L    break;
   case 3:
: J" A- A' C* x0 h; Y( s& n    num_addv++;
    num_ba++;
, P, U9 U$ E) `/ P; r    c[tnum_x+num_addv]=0;
  i$ D/ Q1 q$ R& @    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;
/ [! h9 }/ r7 z: F' G, I     else a[j]=0;
* O0 E2 w4 Y- S* U0 m$ k    break;
  }//switch
}//增加松弛变量、剩余变量、人工变量、确定基变量</P>
1 d. y+ z1 M6 H* M<>//显示标准化后的方程组
* K2 b0 Y' I4 j( M6 c' |# B printf("\n化为标准形式后：\n\n");
# }3 Q! j, e; k! p; ?+ d( d# w9 s if(ma_mi==1) printf("max z= ");
. }! i: N9 \% s* [ else printf("max z'= ");
printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
1 r3 g! Y7 t; Y% e; D; K' U  k=0;
  for(j=1;j&lt;=num_ar;j++)
  if(i==arti[j]) k=1;
  if(k==1) printf(" -M x[%d]",i);
  else if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
4 ~( g2 y* S7 l  K( n9 K: c5 y  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
, I2 O/ ?, H5 g/ l% Y# }4 B  else printf(" %lg x[%d]",c,i);
}</P>
<> printf("\nst:\n");
" ?8 m: _! e: U1 h6 |9 s for(i=1;i&lt;=num_st;i++) {
  k=0;
3 @( [5 q. D" N- f/ N0 q; T  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
3 `% n9 M6 ?/ {+ d8 ]4 L# s/ N/ d    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);
. L6 i! Q* n$ H: I    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);
+ o* P* S, s+ d5 E" Q0 b    else printf(" %lg x[%d]",a[j],j);
  S. ~* l7 p! T- o% }    k=1;
' P; t( x9 Y* ^" E7 j   }
  }
+ y8 a/ f/ w9 p5 n! e8 i" C  printf(" == %lg\n",b);
2 {! l7 n* k9 J" g: v* ~ }
, [) F& r( T0 M- d printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
6 Y0 g+ Z* x% u, D  w* r8 |) [<>void iterative() {
+ c5 L9 l( _2 B 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;
$ R  \3 k1 Q4 a  h4 V double sigma[MAX],temp;
/ ^: x4 h/ r" a1 m9 I7 x double value_be; //高斯消元里保存主元素值</P>
# ?9 R0 W. ?3 T' ]4 |6 K<> printf("\n\n第%d次迭代：\n\n",stop);
for(i=1;i&lt;=num_st;i++) {
* _. g$ g. G' [" A  printf("c%d=%lg\t",base,c[base]);
  printf("b%d=%lg\t",i,b);</P>
<>  switch(step) {
   case 1:
; y6 [3 ?5 R0 f5 p- q# Y! L    for(j=1;j&lt;=num_x;j++)
    {
     printf("a[%d][%d]=%lg\t",i,j,a[j]);
9 K3 H1 }6 i, _. m7 x/ E( l1 u    }
    printf("\n");
    break;
8 \5 C. `7 N- p" d, f, w7 b   case 2:
2 @8 {% G8 v8 t3 y6 ^    for(j=1;j&lt;=num_x;j++) {
8 [: c3 S4 z5 D     k_a=0;
! [7 t5 F" }2 u4 t7 _: o1 {     for(l=1;l&lt;=num_ar;l++) if(j==arti[l])k_a=1;
     if(k_a!=1) printf("a[%d][%d]=%lg\t",i,j,a[j]);
    }
    printf("\n");
    break;
6 Y5 |4 F9 u% R3 z! {+ a  }
6 C7 _7 Q6 G  H! ~3 L" V4 J/ `/ K6 ~ }
# V6 Z6 b1 U9 X9 N+ D) F//求检验数sigma
for(i=1;i&lt;=num_x;i++) {
  sigma=c;
  for(j=1;j&lt;=num_st;j++) sigma-=c[base[j]]*a[j];
; B0 v. B$ d+ u( k& u; R2 A  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
+ Z$ z$ g  _" \9 p# x5 j8 d/ ?& ^  |  switch(step) {
7 ~4 N+ H( M( f. Z( `* H: ?% D   case 1:
    printf("sigma[%d]=%lg\t",i,sigma);
; A7 [; N+ d- t( D0 ]: n& B: P    break;
0 J  O9 G) f- F0 R4 w   case 2:
' C* z- |6 q& p& w    k_a=0;
+ @6 l& h2 S( ?) Q$ D* t    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;
  }
6 a' s( ?8 S& @3 d6 ` }
: V6 L: a9 M( Z- i putchar('\n');
. D1 r5 O! Y: t. I) M//检验检验数sigma是否全小于等于0
k=0;
! Z- `- V6 v  e8 q/ D) M for(i=1;i&lt;=num_x;i++) {
: V0 U4 R+ \  d7 W  if(sigma&gt;0)
  k=1;
( o5 q2 d2 @/ P& w8 A! g }
if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
  for(i=1;i&lt;=num_x;i++) {
   k_f=k_a=0;
: o2 w+ u$ I/ t   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;
   if(k_f==0) {status=-1; return;}
   }
  }
  status=1;
. w4 m( l+ N  _  return;
& _# R5 ^+ i9 y; `( c }
9 L& q% ^4 a+ r3 A: Z7 V//检查是否为无界解
for(i=1;i&lt;=num_x;i++) {
) Q8 C% }' _" l; K3 z, L& X  k_f=0;
( }" g6 O& \  R  if(sigma&gt;0) {
   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
4 ?4 `$ z2 Z& T% n   if(k_f!=1) {status=0; return;}
  }
8 ~$ Y6 g4 i& w6 p  @ }</P>
<>//确定换入变量
$ i' V) H5 @" _9 \8 E for(i=1;i&lt;=num_x;i++) {
. f6 Q( `& I! U+ f3 _  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
8 n' P) q# e- B3 W$ j  |  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
}//temp赋初值
& w4 z* W# O- d% ]+ w% { for(i=1;i&lt;=num_x;i++) {
  k=0;
7 v. W  \5 \2 k7 z' t  K+ ~0 d; t  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
: A6 v6 b) t% _( B0 A3 k  if(k==0)
   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
    base_in=i;
    temp=sigma;
   }
}</P>
( u) L7 Z  \$ w6 a. r/ B* x<>//确定换出变量
& g5 m/ z/ t" Z9 f/ U- I% C$ x% U for(i=1;i&lt;=num_st;i++)
  if(a[base_in]&gt;0) {
   temp=b/a[base_in]+1;
   break;
  ~3 X+ X5 g' I. h5 n& q, p; S& ~  }//temp赋初值
for(i=1;i&lt;=num_st;i++) {
5 y3 I" e  a6 @% `  z7 r  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
   for(j=1;j&lt;=num_ar;j++)
- a1 g! }; v3 n8 P0 K    if(base==arti[j]) {
     base_out=base;
     base_elem=i;
5 P# P' A( w' j  T' B5 G     temp=b/a[base_in];
: q) [. |1 Q" {     break;
' J$ `5 r- k  c6 q% m0 y3 O4 e    }
9 Q# b# Z5 O( E) @' m  }//人工变量优先换出
- G9 n4 I8 B' V' X3 z  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
" r. e1 n  R2 L( l& P, |3 X   base_out=base;
- i/ O, O; y6 v   base_elem=i;
   temp=b/a[base_in];
4 @, J7 ]3 S" @$ q* N3 H7 a# T  }
}</P>
<> printf(" 基变量：");
8 d+ f+ w8 i' H* r for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
; K! o! [8 j7 ^2 M' b* k! y4 n printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
//基变量变换，进行新方程初始化后迭代
for(i=1;i&lt;=num_st;i++) {
  if(base==base_out) base=base_in;
& B" z. b2 W" A/ C4 s }
: L* g' i  B  v$ l& U4 k. i2 [//初始化主元素行系数
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>
7 v1 s6 j* ]+ M5 ^: Q, W& r+ T<> for(i=1;i&lt;=num_st;i++) {
1 D( b% N4 ^$ E% s  if(i!=base_elem) {
   b-=b[base_elem]*a[base_in];
   value_be=a[base_in];
, g4 a/ \7 P' O. ]$ \" J+ ~   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  }
0 I/ W; T! V$ O/ Q, R: h4 b7 F }
stop++;
& Z0 N' N- T* k# g; R$ R4 b$ C if(stop&gt;STP) {status=-2; return;}
iterative();
! e( o. a) n9 `3 L5 n4 e}</P>
<>void output() {
int i,j;
  x1 j: M/ q/ g9 r" f double X[MAX];
) N. U7 E. r$ C8 n printf("\n结果如下：\n");
printf("\nX=(");
; T/ u! x7 V& Q+ v& v2 l% m for(i=1;i&lt;=num_x;i++) {
) {  [% B9 _4 A  f  for(j=1;j&lt;=num_st;j++)
! E0 x! |1 G  A9 ~/ ]  if(i==base[j]) {X=b[j];break;}
* [! p9 W( T- ^; ~6 `  else X=0;
' y# e1 K9 C' S- X9 z) [# m  printf("%lg ",X);
: Z. Z3 L% o6 a  ~4 y  } }
# Z4 ?; y. M* Y: ]" v- i& O 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？帅呆了