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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
: F$ @, U* H5 R4 E#include&lt;math.h&gt;
$ F# K6 t1 G! L. @% U#define X 5
% b! |$ Y9 t- d6 I8 Q" t- c#define Y 7
' u0 ]* ?  C: M3 z% |( K9 svoid xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
{
  int j;
  *c=0;
: |! B/ J3 \" \1 ~  for(j=1;j&lt;=*is;j++)
$ F6 h( q3 K# ^* V! v    {
      if((*a)[*ir][j]-*c&gt;0)
$ e) X7 ?( u. ]+ [: f- X    {
4 c+ `4 ^9 {5 b3 v3 k- Z6 N& C      *c=(*a)[*ir][j];
7 I1 x; e& E! u$ N( _$ c% \      *j0=j;
* R. i/ o- U/ B* }    }
/ T/ i- n* V: D% \0 H     }
1 x% T7 |& @  m: [) {( Q! L& i( T3 u}</P><>      /***************** 参数说明 **********************/
      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
; z) b9 p* J: Z3 g: L      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
      /********* <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)[])
" Z! c" S" F$ R$ A- ]( g5 L{
/ A2 m& u2 n" s5 |! L1 j  int i,m1,mn,j0,i0,j;
  float c,g;
  m1=m+1;
- X: H- y0 e7 I  mn=m+n;
6 S3 W( `: s/ _$ e1 [+ P) Y  for(i=1;i&lt;=m;i++)
    (*k)=n+i;
leap1:
  if(l1-1==0)
6 V2 M& ]8 B& Z; B4 g1 j- y$ V5 @1 @     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
. i- g0 E0 o( Z/ O  else
     {
  g- m; E% i/ ^) W/ j* d* z3 y# mleap2:   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);
     }
6 C) R0 |8 h# Y9 X   c-=1e-8;
   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
4 p% \( G( ^& _9 z% C      {
+ K0 z) ?8 d1 V1 u    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
" v% W7 G7 `. i0 V: |- u1 `, v    return(-1);
      }
- \& T' g' ~- o' H   if((c&lt;0)&amp;&amp;(l1==1))
       {
     l1=50;
     goto leap2;
       }
    if(c&lt;=0)
* r$ F, l$ V. g3 B       {
# b6 U/ O2 v9 t) t( R% X      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));
) W9 I% q) I8 T7 }: i: `0 a% `% K      printf("\nF=%f\n",(*a)[m2][mn1]);
      return 0;
/ A  `% S4 z, S    }
     c=1e8;
! {2 B8 z0 ]# I2 }/ Z     for(i=1;i&lt;=m;i++)
       {
8 J) i0 M' A* h( L: w" v# z     if((*a)[j0]&gt;1e-8)
- ^+ v; h( [8 e  u2 E! a        {
- g1 K9 W3 a1 z( S' x6 X! P) b          g=(*a)[mn1]/(*a)[j0];
          if(g-c&lt;0);
         {
4 p% R8 s, _4 I( A+ z# p           c=g;
+ s3 T6 ~& ]; G/ }1 A7 o9 f           i0=i;
' D5 W- j2 }* S' ?9 ]( [# N+ e' }          }
         }
; i* a9 S: ?9 X+ I    }
9 r, ?. j! c  E, k      if(c==1e+8)
' a; M$ h9 Y1 @# A* B, b    {
      printf("\n\t*******Lp no solution********\n");
      return -3;
     }
       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
      {
        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
, M0 D8 {; ?: l% }* ~           continue;
         g=(*a)[i0][j]/(*a)[i0][j0];
8 v7 y% S( Y0 W" ^% o* s+ P8 \! w           (*a)[i0][j]=g;
         for(i=1;i&lt;=m2;i++)
5 \3 [- u+ {, @5 R           {
         if((i==i0)||(!(l1==1)&amp;(i==m1)))
8 w& O, A! u2 ^' g* U' ]+ B" T) X0 w7 x            continue;
         (*a)[j]=(*a)[j]-(*a)[j0]*g;
           }
       }
       for(i=1;i&lt;=m2;i++)
) T7 w& `2 e, W% t4 U      (*a)[j0]=0;
1 w5 o" f/ |7 V       (*a)[i0][j0]=1;
. z  n8 \& v, m! B       goto leap1;
}  
main()
{
  float a[5][7]={{0,0,0, 0,  0,0,0},
  s* O' c2 `5 A# i& n                 {0,3,-4,3,  1,0,12},
                 {0,3,0, 6,  0,1,12},
                 {0,0,0, 0,  0,0,0},
9 C( _  ]; q$ @, s1 L                 {0,69,0,144,0,0,300}};
* }9 i+ T! Z! n5 {2 \  int k[3];
  float x[3];
  clrscr();
  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
( Y+ M) P$ S( r7 n5 O! o}
</P>

plgatc

<>/*************************************************************************
$ ]0 `/ Q8 m/ R7 w# S# L                    单纯型法解线性规划问题（两阶段法）        
3 i, z# t9 `7 }2 T& y
编程环境：VC++6.0     
方程组输入说明：
变量非负，按提示输入相关参数。
7 `4 S5 X0 A) T7 X6 l$ D" g7 g *************************************************************************/
7 p0 J& Q% Z; D( z2 v#include &lt;stdio.h&gt;
+ k4 C9 w9 t7 ]" K2 e/ K( {1 {$ E#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>
<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
int num_x; //变量个数
int num_st; //约束方程数
) z7 t) J# ?5 G  H$ q$ `3 t2 i& }int num_ar=0; //人工变量个数
* A/ x2 |6 [# I4 u6 v! ^1 W/ Yint arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
<>void create(); //建立方程组
void iterative(); //单纯型法迭代
' |, l9 K4 v6 [4 `void output(); //输出结果
- G6 F; a4 j- Q3 O5 Ivoid banner(); //打印程序标题
9 f$ N; q: d' N/ O$ P. Nvoid exchange(); //交换两阶段价值系数
: b; y0 I3 T* {4 Gvoid show(); //输出方程组</P>
) `5 w1 o) q5 h% q% C# m<>void main() {
int i,j,k;
. M" i' H8 X) r+ x banner();
+ K5 f$ M4 @0 ?+ L; o! k% A create();
8 \! F1 U( O4 ^+ @8 w+ p //保存原价值系数，转换为第一阶段价值系数
! H4 x& o0 @, w6 j9 H# `! S. `3 D for(i=1;i&lt;=num_x;i++) {
- l+ t4 ^, K: t# V" |& J" I  k=0;
' a7 c0 N' `( }  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
  if(k==1) temp_c=-1;
  else temp_c=0;
}
exchange(c,temp_c);</P>
<> printf("\n\n第一阶段问题为：\n\n");
: \& S' p/ W* @3 _( j show();
+ \' W  b' `0 `: v step++;
5 f, x& V) O7 J printf("\n\n按回车开始第一阶段迭代");
; h; y6 y1 L" \; X$ {( [4 D& U getchar();
getchar();
' K- V/ r2 T' X6 \ iterative();
+ y. k7 J* E& `  p: N- t! r if(status==-2) {
% f. A" c' d- ~# X7 A7 ]  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
! l/ g  r: T+ e$ Q# T  getchar();
  exit(0);
}
+ f( Z- h5 k. @6 ?% x( t/ }9 u output();</P>
<> if(max!=0) {
  t' R6 E8 N* k4 c1 @0 s7 U1 e9 r; v  puts("\n\n原问题无可行解。\n");
: ~, i5 T; c8 r" d4 g  puts("\n按回车结束");
3 X2 Y/ @6 n+ N/ Q, @; }  getchar();
( x- ?- V1 t6 n7 L2 V. X; }  exit(0);
}</P>
<> //转换为第二阶段价值系数
exchange(c,temp_c);
//把人工变量列全设为0
! f8 S5 M# S4 Q& j' R# m for(i=1;i&lt;=num_ar;i++) {
  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
+ [* d; w( }, C1 {4 g3 B' u& u  n }</P>
! S2 j0 i2 q) {9 u  \& r! _<> puts("\n\n第二阶段问题为：\n\n");
6 e* ^0 n4 e' q) O4 _ show();
- }/ n1 d5 f$ x' G puts("\n\n按回车开始第二阶段迭代");
getchar();
5 E' ?+ e; }, \9 a" X iterative();
switch(status) {
case 1:
  output();
, m  j2 R' E! I' E. @  puts("\n\n原问题有唯一最优解。\n");
' @, `9 e* I! L  puts("\n按回车结束");
9 s2 s) ]* f) i6 s3 Y  getchar();
  exit(0);
case 0:
; h& b: J' L; a  X  puts("\n\n原问题为无界解。\n");
: ?3 d: n" [5 ?* d) E& N  puts("\n按回车结束");
  getchar();
- J+ t! a- ^2 T" l  exit(0);
( @2 U* \2 c0 j" h& P0 t case -1:
& D* J2 {+ A) W; y( Q( H4 @& {  output();
  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
  getchar();
' w9 U' T# h( D6 D  exit(0);
case -2:
$ T4 H& h8 S$ x& w  p  k6 a  puts("迭代超过限制次数强行终止！\n");
7 P) `  m) a( x6 p! M7 U  puts("\n按回车结束");
5 |$ k7 z6 B+ F5 |; X  getchar();
6 P6 q; z: K5 ?) I+ [4 S  exit(0);
7 I) a0 y% N3 ~6 l }//switch
+ e' h7 U5 F3 y8 m9 Q. X}</P>
: ^: `' H: v; M0 H9 H<>void banner() {
9 M5 x! D9 U" A! b6 ]3 ^, Q printf("\t\t****************************************\n");
printf("\t\t         单纯型法解线性规划问题\n");
printf("\t\t                         作者：Thunder\n");
printf("\t\t****************************************\n");
printf("\n");
! O, w8 h7 M8 u: y6 d0 }" L}</P>
9 m0 l2 M7 |. o6 H( j4 a<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
//为1的不显示系数1，-1系数只显示负号
2 r  |- C. h0 z7 P$ [6 S. u int i,j,k;
% Z! q3 G& z1 M6 U& B switch(step) {
case 1:
4 c  ?7 T" c6 |3 S$ D+ n. M- `  printf("min z= ");
. O9 R3 C$ r4 n% K" ^  printf("x[%d]",arti[1]);
7 D! K8 _0 }; n5 a5 D  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
, {* ?# ~# C. a, O% u- C7 n case 2:
printf("max z= ");
8 {/ x* i8 d$ Q8 T printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
6 g  A! g; S( K! K  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);
}
& {1 T' T, k+ t+ z; y! R break;
$ c7 a: A1 L+ g7 G/ F6 R }</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
9 v$ H& I4 L' u5 _  Q  k=0;
9 z+ D* ?8 ]& J5 q2 u  |  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
+ F: G  b8 {( o9 D2 X$ Y    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);
, O0 P$ M5 U4 U! 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);
8 f4 K3 H5 B0 F' x0 D    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
3 A- F) R9 z5 O8 M' ~1 h5 M    else printf(" %lg x[%d]",a[j],j);
/ G# q1 [' `: G( a6 |8 U    k=1;
   }
  }
  printf(" == %lg\n",b);
' v1 ~$ |5 g$ x# A  i- t }
+ a, Q  j6 T7 \% M# l2 X printf(" x[1]~x[%d]&gt;=0",num_x);
( }  u& r6 s# p! h}</P>
# [9 G/ f1 [0 t8 \+ Y+ ?( o<>void exchange() {
" u9 h# k" Y8 e0 m- U int i;
/ [: M7 a- ?3 h' l  ~! |7 L  c$ O6 R double temp[MAX];
9 k1 D+ Y$ |" L+ O: V for(i=1;i&lt;=num_x;i++) {
  temp=temp_c;
  temp_c=c;
  c=temp;
7 V2 F& ?7 {5 I5 ^6 p- D* V }
7 ]* I" W0 E3 l- [2 G}</P>
<>void create() {
- p/ Z$ f1 U: h; T' i0 U% x, h# [//输入方程组系数，每个方程输完后回显确认
, n; X! u# n+ H+ M int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
7 |% C5 z& Q3 h1 K char confirm;
" k' T' J. Q9 ?; s, a
, o; M% O, s* B0 M while(1) {
  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  scanf("%d",&amp;ma_mi);
  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
  else break;
" R, I3 A: Q5 ]1 {0 c0 J }

5 D4 x4 A- s' i while(1) {
6 z% p, e2 ]  S$ d5 P4 V  printf("指定变量个数：");
1 v& }& a1 I" _, p; |+ {, o  scanf("%d",&amp;num_x);
" ~# w8 k! W; V7 c0 ~4 t$ O& q  printf("输入价值系数c1-c%d:\n",num_x);
1 S* W% v6 N, Q# k3 o  for(i=1;i&lt;=num_x;i++) {   
" Y5 _5 ?5 _: w; ~, n6 ~) `   printf("c%d=",i);
   scanf("%lf",&amp;c);
  }
  if(ma_mi==1) printf("max z= ");
  else printf("min z= ");
  printf("%lg x[%d]",c[1],1);
7 H7 G; `; W* c: K( W, S( ]  for(i=2;i&lt;=num_x;i++) {
) u3 y! h7 O1 U* ^   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
   else printf("  %lg x[%d]",c,i);
$ P5 Y$ S: k# V" x7 n7 m  }
8 }' g  A- \/ x, v  H( N printf("\n正确吗？：(y/n)");
getchar();
confirm=getchar();
" f  ?" J3 G7 `* D' }& D- m! T if (confirm=='y') break;
5 D1 S. n7 {! n: t else if(confirm=='n') continue;
9 x) R' ?( Z* d, S) Q" o" V }</P>
<> printf("输入约束方程组个数：");
scanf("%d",&amp;num_st);
) ?) k6 s3 @3 l; c. j3 ? for(i=1;i&lt;=num_st;i++) {
" Y: g. c& r8 h4 J: o4 V  printf("st.%d:\n",i);
  while(1) {
- c  A. j2 Q& t& m9 g8 h   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
   scanf("%d",&amp;re_st);
( W% o5 K3 o# Q& {+ d8 }   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
3 u( `) f) d1 h0 j   else break;
  }
  printf("输入技术系数：\n");
  for(j=1;j&lt;=num_x;j++) {  
( J. @4 v) d, b0 W; P6 o8 h   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
7 z+ f, x6 }4 D( G& F1 P  }
& Q; I% o8 ]. }' C) L4 `  printf("输入资源拥有量：\nb%d=",i);
  scanf("%lf",&amp;b);
  
  printf("st.%i:\n",i);
  printf("%lg x[%d]",a[1],1);
* _5 P, j. q+ `+ o0 X# T, i  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);
0 J3 Z9 I1 q2 T6 o3 Y  }
  switch(re_st) {
+ f7 m/ H- H9 f- W6 a/ D% ]   case 1: printf(" == %lg",b); break;
   case 2: printf(" &gt;= %lg",b); break;
9 p, c- g( V& c+ e   case 3: printf(" &lt;= %lg",b); break;
  }</P>
<>  while(1) {
   printf("\n正确吗？(y/n)");
+ y! J# [/ m" A) G+ C   getchar();
; b" @# E8 h% f1 V   confirm=getchar();
   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
  }
}</P>
<>//显示输入的方程组
printf("\n原问题为：\n\n");
if(ma_mi==1) printf("max z= ");
else printf("min z= ");
printf("%lg x[%d]",c[1],1);
: Z& X7 q9 B. @9 m9 \& l for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
3 O. O' g9 o/ k7 L3 j, Y  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);
# c" Q' t6 m) i7 w) e- n }</P>
<> printf("\nst:\n");
: N0 J0 S: v* z3 N' L6 O for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
7 B# p$ v' Q( o( ~( X% q+ n   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
* K/ r; H/ u$ V% [2 h# p    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
/ S9 ]# p5 G6 }9 a3 v    else if(a[j]==-1) printf(" - x[%d]",j);
( u+ o% h- H8 z% m4 c) N    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);
! O( y4 p- x8 _. [* |: \5 Q0 o9 g9 m    k=1;
' k( r' a! y& Y9 d( P' w% \   }
) H/ L% \4 w$ M- ?  }
  switch(re_st) {
2 {3 g4 P3 ~# v7 V* L1 k   case 1:
. T- U$ a& q; [    printf(" == %lg\n",b);
3 w% c6 S$ l1 ^1 l4 H/ P    break;
   case 2:
    printf(" &gt;= %lg\n",b);
# z  y* [+ C+ [- c/ L* z0 d' S7 c    break;
   case 3:
! |/ o9 K8 o: p1 O% [+ P! p    printf(" &lt;= %lg\n",b);
    break;
7 C3 J" E# c$ ]" U  }
% F$ d% N  [1 n. Q }
printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
# H- Q8 b8 k, Y6 Q/ _% W* n<> tnum_x=num_x;
for(i=1;i&lt;=num_st;i++) {
3 q/ f5 y- M2 T' N& p; v1 Q+ s/ {  switch(re_st) {
& }( ~% ~$ A+ |/ w; J" b/ v. z( [  case 1:
3 B% N: `6 g% i5 H  case 3:
   num_x+=1;
   break;
" U$ f. i& w' k9 E5 I7 q$ A4 s  case 2:
   num_x+=2;
+ d1 f/ u8 G1 q. a   break;
  }
8 J  H$ A! n. U' y/ d4 i7 `7 H; K }</P>
0 m1 E3 v9 D9 g5 r  ]% k3 A<>//化为标准形式
if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
' i% G0 W5 G8 x' E& o; Q8 u. \! } for(i=1;i&lt;=num_st;i++) {
! c9 @$ R3 r6 r2 }9 T- b8 z  switch(re_st) {
. Z* t6 V* c# C! v) r7 ^: w& J   case 1:
    num_addv++;
  K6 Y3 F/ @4 U4 D9 w    num_ba++;
    num_ar++;
    c[tnum_x+num_addv]=0;
, I1 z& U/ x* F, Y7 ?& O! Z  z, Y    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
    for(j=tnum_x+1;j&lt;=num_x;j++)
+ G& e4 V$ }( D& Y1 @& D+ E' i7 _     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
    break;
, M% M4 Q/ P; J" i( S  u3 h+ C% W! u   case 2:
( H. U' a8 U7 G7 A& `" E    num_addv++;
) `% `" y/ ^% q: W    c[tnum_x+num_addv]=0;
4 N' P3 B% q9 n4 V+ F* C    num_addv++;
    num_ba++;
/ y4 {4 |  I( M$ W; O& g    num_ar++;
6 U( w8 ]* G! c8 t8 t    c[tnum_x+num_addv]=0;
% g* }. O# f& X  Z$ @( Q    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
8 B" h  e3 w8 r  C) o$ }- c9 C    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;
+ A; w: f0 D0 d; v    break;
   case 3:
    num_addv++;
" q0 @3 w+ R* v6 C( \- ~9 c    num_ba++;
3 ]+ M( Z+ J8 t% O, f( x  U/ B    c[tnum_x+num_addv]=0;
    base[num_ba]=tnum_x+num_addv;
% n6 c# X# _/ k  k& ?% N# I    for(j=tnum_x+1;j&lt;=num_x;j++)
% [1 E+ X; `) D" T* k     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
, G7 x7 C# g9 b! H- l     else a[j]=0;
: i+ v8 h4 D8 A% I    break;
3 _# ]+ V& Z6 o& i  }//switch
) T+ g# {& @  N }//增加松弛变量、剩余变量、人工变量、确定基变量</P>
; w4 |( ]- F. i0 I! G2 i<>//显示标准化后的方程组
* W# d$ K) b" ?$ N printf("\n化为标准形式后：\n\n");
if(ma_mi==1) printf("max z= ");
else printf("max z'= ");
) q: w6 J6 B+ e; Q& O' K( i8 K printf("%lg x[%d]",c[1],1);
7 X9 K5 ^6 ~' g& ]+ Y for(i=2;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_ar;j++)
  if(i==arti[j]) k=1;
+ d! M" E8 R2 A4 B  if(k==1) printf(" -M x[%d]",i);
  else if(c==1) printf(" + x[%d]",i);
, B2 c4 Z) Q+ P" u: u+ l/ c  else if(c==-1) printf(" - x[%d]",i);
  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
8 y" `! K# T. a7 G4 o% }) K' B  c# y  else printf(" %lg x[%d]",c,i);
$ Q+ b' p. q! k3 @" m }</P>
<> printf("\nst:\n");
6 N8 c; E  t) s% `5 s& H, ~ for(i=1;i&lt;=num_st;i++) {
+ Y  Y2 S0 d' F; y5 M; ~! a  k=0;
2 I! J) H5 Z3 p0 s  for(j=1;j&lt;=num_x;j++) {
2 |3 r. N4 [  B* m   if(a[j]!=0) {
/ m6 J: s1 M; x. W$ }; H    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
1 n) F+ ]3 b2 v5 q8 G8 g( X    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
" g+ |5 s: V" l! W    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);
: f0 {. I7 p- i0 @& ^% p2 Z1 f    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$ }/ T) g+ m8 S    k=1;
( f- P& H1 c) g5 T( @/ ~   }
: n: h, T2 u' G/ p* R5 {  }
% G4 M$ r' z; R, {, \- y  printf(" == %lg\n",b);
$ k+ F  e; X+ Y* w9 d  a! A) s. }+ s }
printf(" x[1]~x[%d]&gt;=0",num_x);
. Q! |3 w, \. ^: L# x! h}</P>
$ k3 V: `2 ]+ [% I& I+ V& d<>void iterative() {
int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
int base_elem;
2 g/ b8 u" c7 s% n int base_out,base_in;
, G, X. J* t# E# W double sigma[MAX],temp;
, c$ [1 U0 l, z; e" j' S double value_be; //高斯消元里保存主元素值</P>
! T4 Y0 w- G. I% Z8 E; Y6 `<> printf("\n\n第%d次迭代：\n\n",stop);
for(i=1;i&lt;=num_st;i++) {
) p: E7 l3 Y5 ~  printf("c%d=%lg\t",base,c[base]);
  printf("b%d=%lg\t",i,b);</P>
- U  H9 q' ]( s$ N/ U( B<>  switch(step) {
. ?' n9 ]; b; E& D! P) u   case 1:
    for(j=1;j&lt;=num_x;j++)
    {
     printf("a[%d][%d]=%lg\t",i,j,a[j]);
    }
    printf("\n");
    break;
   case 2:
' |( G% v0 I+ S, I    for(j=1;j&lt;=num_x;j++) {
     k_a=0;
     for(l=1;l&lt;=num_ar;l++) if(j==arti[l])k_a=1;
& @  `8 k+ F! }/ Z0 {* M. T     if(k_a!=1) printf("a[%d][%d]=%lg\t",i,j,a[j]);
0 q. M% a4 F6 E    }
5 c; D  @0 B4 y$ q! y3 ~) V    printf("\n");
: a$ c: _6 m$ i& K8 }% T# R8 t; H' U2 ]    break;
+ p/ D$ q& u" E, l  }
' h8 A- i8 A, E6 r% g7 J }
//求检验数sigma
for(i=1;i&lt;=num_x;i++) {
  sigma=c;
& N& s2 y3 Z# E& z3 ~  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) {
+ y$ H( G$ ]% u. o* i   case 1:
7 o2 X0 L" q0 g    printf("sigma[%d]=%lg\t",i,sigma);
    break;
   case 2:
  K% h9 {3 b0 K5 W5 H    k_a=0;
( q1 N& e3 Q- i5 K    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
9 [6 F. \; Y+ K9 r) r5 p  r    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
& n+ H& v' i2 C& i    break;
# C' s) ?9 O. P& ?( a  }
7 W1 A2 u$ b% _* N4 m# O }
. t/ F5 H3 b& D4 W. Q putchar('\n');
//检验检验数sigma是否全小于等于0
k=0;
$ R) r- C1 B9 _: q$ N for(i=1;i&lt;=num_x;i++) {
; [9 N. n1 H' C5 U$ I6 \  if(sigma&gt;0)
# o! M4 Y# C1 P: U  k=1;
/ Q3 X. N( a' F2 w/ D  q, I }
: ~0 |4 k  N, s if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
) b3 ]8 }: C  B, A0 d8 C  for(i=1;i&lt;=num_x;i++) {
& g5 J+ s3 x# w0 m9 [0 m   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) {
- B( D& K. G/ f3 Y    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
   if(k_f==0) {status=-1; return;}
   }
  }
  status=1;
4 e: ]  B6 L0 D9 F, b0 {2 v% g  return;
# ]! j# `3 S, D }
5 p! g7 d+ I! b//检查是否为无界解
for(i=1;i&lt;=num_x;i++) {
3 Z0 B) n0 o: y" Z# U0 S  k_f=0;
% o* t+ F7 d1 Q1 q0 ~! g9 W4 J  if(sigma&gt;0) {
% c4 B; @: q: X8 m8 H9 h8 }   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
   if(k_f!=1) {status=0; return;}
  }
) u( D1 }& i3 J3 C }</P>
0 N! N+ x# b9 l5 ^& \5 z6 `<>//确定换入变量
2 G5 V( J5 G" n! s3 N; P1 Y6 T0 T for(i=1;i&lt;=num_x;i++) {
: P6 o9 g( c0 N5 e7 g3 u  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
2 Z+ n& l; r! y6 F# ~& K  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
}//temp赋初值
; Q/ _+ w( A4 D; r, i/ ~1 B. n for(i=1;i&lt;=num_x;i++) {
1 h% w" A1 B$ D) y  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  if(k==0)
   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
    base_in=i;
- Y2 {3 D+ e# j& i% V    temp=sigma;
% N+ ]8 X9 G: U2 e, x9 y0 f   }
}</P>
<>//确定换出变量
$ [$ |* P: K6 [3 U% l$ z$ A for(i=1;i&lt;=num_st;i++)
  if(a[base_in]&gt;0) {
& D6 _- |' B& U( U. Y% b   temp=b/a[base_in]+1;
! }5 B" p. u9 e8 @* d   break;
$ @* L& P; G  ^8 k  Q9 T# n  }//temp赋初值
for(i=1;i&lt;=num_st;i++) {
* O9 o9 Y# y* F5 l  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
  u9 a( z( ^  @7 G" P0 J   for(j=1;j&lt;=num_ar;j++)
    if(base==arti[j]) {
' }0 k0 m, a* i     base_out=base;
4 ^/ p# t4 r" d6 b8 H     base_elem=i;
     temp=b/a[base_in];
     break;
    }
  }//人工变量优先换出
  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
   base_out=base;
   base_elem=i;
3 k5 Y6 c# {+ t# d2 g   temp=b/a[base_in];
% ^4 I1 o8 Z, b, [$ o  }
}</P>
<> printf(" 基变量：");
3 C+ |0 F* ?  c, J  p( ~+ i& L for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
" _5 l0 H/ ?9 {' h//基变量变换，进行新方程初始化后迭代
for(i=1;i&lt;=num_st;i++) {
  if(base==base_out) base=base_in;
}
7 }# T: d& S$ u3 e$ G# e//初始化主元素行系数
value_be=a[base_elem][base_in];
b[base_elem]/=value_be;
# g* k) S0 q5 U/ z; P  d for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
<> for(i=1;i&lt;=num_st;i++) {
  a+ [  a8 O1 W4 M  if(i!=base_elem) {
   b-=b[base_elem]*a[base_in];
   value_be=a[base_in];
   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  E* }2 Z% S9 a0 P2 W) H6 n# ]9 ?  }
}
stop++;
$ r% @# }' Y  \: q: f9 ] if(stop&gt;STP) {status=-2; return;}
iterative();
' K+ P6 _+ A. Q) {}</P>
<>void output() {
int i,j;
double X[MAX];
printf("\n结果如下：\n");
printf("\nX=(");
for(i=1;i&lt;=num_x;i++) {
9 R" M  g! i6 i9 d  for(j=1;j&lt;=num_st;j++)
  if(i==base[j]) {X=b[j];break;}
  else X=0;
  printf("%lg ",X);
" v4 c1 ?: t  b- m }
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？帅呆了