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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
#include&lt;math.h&gt;
#define X 5
#define Y 7
void xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
{
9 B& g$ g- [3 m! \  int j;
  *c=0;
. q& [  T- A5 j5 D2 i& w- [, ?. O' y  for(j=1;j&lt;=*is;j++)
    {
      if((*a)[*ir][j]-*c&gt;0)
    {
      *c=(*a)[*ir][j];
! s5 d: D5 Q% n/ }( p      *j0=j;
7 l3 l7 ~5 ]! w+ Z4 O    }
6 N8 l, f1 w0 q     }
}</P><>      /***************** 参数说明 **********************/
      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
0 F* @/ }. A6 M      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
, _5 X, J. v7 ^5 P( U1 k; w      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
      /********* <a href="http://happyyangxu.home.sunbo.net*******/" target="_blank" >http://happyyangxu.home.sunbo.net*******/</A>
/ E9 n& k  _7 [6 N/ j& K, eint xi_sm(int m,int n,int m2,int mn1,int l1,float (*a)[X][Y],
7 X# E- L6 h$ J       int(*k)[],float(*x)[])
' U! M( P  M7 J2 p' D. ^{
# N6 x# F$ K) y$ e7 ?* [  int i,m1,mn,j0,i0,j;
  float c,g;
, b  T7 t4 c$ G% l  m1=m+1;
7 E- o: I) }9 V  \: h  mn=m+n;
  for(i=1;i&lt;=m;i++)
8 x- p3 ^0 m0 X( e$ z6 A0 X3 B& U. z    (*k)=n+i;
leap1:
  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
: @( }& z, v, L* x  B# @) {: a# t  else
     {
4 t* H$ m. [$ _* I5 ]5 O( ?leap2:   if(l1-50==0)
      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
' m5 _7 r1 h1 b& ?/ B9 K1 l2 M& r       else
7 ^# r- o" ^* p6 r      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
     }
/ Z) i# S" |0 u/ P   c-=1e-8;
   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
      {
6 x+ B0 r5 {' G1 X2 d    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
    return(-1);
) \  c7 J/ [' X. n      }
. M# K% K+ P3 D6 ]" J/ e   if((c&lt;0)&amp;&amp;(l1==1))
4 n: e* \  @' n; M# |$ f       {
     l1=50;
* I5 J+ ^% u, i1 U4 w     goto leap2;
( ?+ x1 _+ B$ n4 {% L       }
  I" R  x) ]: o9 ^7 x1 P    if(c&lt;=0)
# d  J0 v5 g0 M* y4 B( f& D       {
      for(i=1;i&lt;=m;i++)
         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
9 @/ f/ E4 J6 _2 w$ ^1 W      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
) l; [. N/ {' X1 |1 X7 t1 p      printf("\nF=%f\n",(*a)[m2][mn1]);
( f, H( ^% k& |5 H      return 0;
  p2 b5 b# E6 K6 `, y8 D7 }+ O    }
     c=1e8;
     for(i=1;i&lt;=m;i++)
- U# J% V, X2 `6 P% X       {
     if((*a)[j0]&gt;1e-8)
        {
( e6 E  c; ?' g' A$ I& Y1 R          g=(*a)[mn1]/(*a)[j0];
& X$ T/ N/ \" j4 L          if(g-c&lt;0);
/ p, j% X7 D9 h' m         {
# @( Y) {# D5 n9 q: H5 b           c=g;
           i0=i;
          }
         }
    }
      if(c==1e+8)
    {
      printf("\n\t*******Lp no solution********\n");
      return -3;
& n5 ~, y# n9 a0 j) f7 L     }
% x& V. z/ J2 r, M' k# f3 Z       (*k)[i0]=j0;
: @- i6 `0 {. [+ ]       for(j=1;j&lt;=mn1;j++)
      {
) h3 ~' |7 x/ }6 e  a# s        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
           continue;
0 J. V% @3 V6 h) m, c' Y2 ~         g=(*a)[i0][j]/(*a)[i0][j0];
/ H! t' ?7 V7 y6 f           (*a)[i0][j]=g;
         for(i=1;i&lt;=m2;i++)
. [& p9 y" c2 V" Z: ~           {
* t# o# c7 M* L  l, }% t4 w         if((i==i0)||(!(l1==1)&amp;(i==m1)))
            continue;
$ Y( ^! E/ y7 C) a& I         (*a)[j]=(*a)[j]-(*a)[j0]*g;
# {/ j* Y' Y' D" U  ~8 T5 m5 U           }
       }
       for(i=1;i&lt;=m2;i++)
$ p. g# Q8 ~- o. ]3 t' f      (*a)[j0]=0;
' u( ?4 Z! W9 a/ L# G       (*a)[i0][j0]=1;
       goto leap1;
$ e# a; @: R# L}  
" A$ q! b1 l  b7 I3 P. amain()
* B6 ^3 b1 X; K& U{
5 m  O( |2 o# D/ M  float a[5][7]={{0,0,0, 0,  0,0,0},
                 {0,3,-4,3,  1,0,12},
! M5 h- O: D2 d6 @6 a                 {0,3,0, 6,  0,1,12},
                 {0,0,0, 0,  0,0,0},
0 G4 {; w, X% u" c  C                 {0,69,0,144,0,0,300}};
  int k[3];
  float x[3];
  clrscr();
2 I  O2 L' R" c1 `  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
: Q% c  o# v1 {}
</P>

plgatc

<>/*************************************************************************
2 F( }% v. ?0 [' m! ]/ K3 Q                    单纯型法解线性规划问题（两阶段法）        
, J/ A, w' `) z- |) a" v
编程环境：VC++6.0     
" O5 ~: I( o/ y3 r$ |" R 方程组输入说明：
0 x8 X; K1 w2 e# \5 D4 S 变量非负，按提示输入相关参数。
*************************************************************************/
, h: s" Q- U$ G$ s, B#include &lt;stdio.h&gt;
( U& Y& `1 v5 v  V#include &lt;stdlib.h&gt;
8 E9 X; ~4 c- @& D$ N) C#define MAX 100
- ?. b$ {& A5 y4 X5 S. g& E2 j#define STP 100</P>
! L' _3 M, o, K* ^/ x/ L! J8 P<>int stop=1; //迭代记数变量
int status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
1 A3 ^& Y  G0 dint step=1; //目前阶段</P>
<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
% M, `% N6 w! t+ L6 _, f2 Kint num_x; //变量个数
int num_st; //约束方程数
- s! [. f! S' [5 Y# jint num_ar=0; //人工变量个数
int arti[MAX]; //人工变量下标
9 u: E7 T+ }/ S) v+ a5 j3 Fint base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
<>void create(); //建立方程组
void iterative(); //单纯型法迭代
void output(); //输出结果
$ b3 t4 }: l7 k9 w5 kvoid banner(); //打印程序标题
% M' j7 n. d8 Z8 t0 ~- W0 r) bvoid exchange(); //交换两阶段价值系数
void show(); //输出方程组</P>
4 F. |: H7 k9 `<>void main() {
: A: `2 h  w+ m1 K/ Q int i,j,k;
" h* V7 M/ `4 I" S1 z8 T banner();
% Y: \) h3 o+ w% e' b( ] create();
//保存原价值系数，转换为第一阶段价值系数
' |) B4 X# w+ u( H& D: `: K2 q9 u# W for(i=1;i&lt;=num_x;i++) {
  k=0;
/ |5 U" L; c7 X3 f6 d  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
0 [/ I3 O) E* D( a: u5 j; r( c  if(k==1) temp_c=-1;
8 E7 J8 K* S) C- O; G3 k" B  else temp_c=0;
- _4 Y$ V/ Q9 t) V3 R }
exchange(c,temp_c);</P>
& Q* f  ?5 p$ o; P<> printf("\n\n第一阶段问题为：\n\n");
show();
step++;
0 E! R: x) E7 J  R+ l) j printf("\n\n按回车开始第一阶段迭代");
( ?1 P+ ?8 n! T, O. ~( A: V( E) Y% N7 R getchar();
getchar();
$ m( }" B9 r/ ^6 E( h7 u/ N iterative();
if(status==-2) {
  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
  getchar();
7 f: i2 x  @: l% B6 c% U  exit(0);
}
output();</P>
4 `# ^  Y1 x! h5 G. z  F<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
  puts("\n按回车结束");
  getchar();
4 u: s) M% l) [( e; I( @5 [8 V$ x! j  exit(0);
3 {- H+ g! C1 D4 z+ j }</P>
<> //转换为第二阶段价值系数
exchange(c,temp_c);
; g8 F: ^( \( p8 U: b0 e0 A! D //把人工变量列全设为0
: E2 _. E* v8 v  B( P for(i=1;i&lt;=num_ar;i++) {
2 |) P7 g7 Y% Q/ ~! I4 P% M3 n  D  c[arti]=0;
3 ?# p2 `' g) z. i. z: C  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
}</P>
<> puts("\n\n第二阶段问题为：\n\n");
show();
puts("\n\n按回车开始第二阶段迭代");
  Y% x, C8 G# @ getchar();
& e- F, Y8 j% C1 _* ], b5 b iterative();
$ g7 }) S6 a9 Z) D& t switch(status) {
1 q5 X( D; q8 @+ \7 ~! I case 1:
  output();
, F6 B6 j; G4 R  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
/ U$ _% l! @6 j9 g9 n7 e1 _  getchar();
0 B9 o( A: U3 O  exit(0);
" P% W: R. @) A/ y" Z% r4 i( G case 0:
  puts("\n\n原问题为无界解。\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
5 a/ E/ s( E1 y  m. Y( A: J case -1:
  output();
6 q8 j0 D) t3 L( z. N% m  puts("\n\n原问题有无穷多最优解。\n");
# m& c$ |. ~8 ]  puts("\n按回车结束");
* r* P( ]# K* M$ F2 i7 @  getchar();
  exit(0);
' x# A) i9 W  j+ ^% d8 n case -2:
/ C  V  I( ~0 \7 J  puts("迭代超过限制次数强行终止！\n");
9 Y  l" k5 @# ~! q! \  puts("\n按回车结束");
  getchar();
  exit(0);
6 \- a$ g) p5 v1 m) m& }( ?; V: `' S }//switch
; Y# G, d6 u# M}</P>
<>void banner() {
printf("\t\t****************************************\n");
printf("\t\t         单纯型法解线性规划问题\n");
1 }) r6 m' j6 M- k, a printf("\t\t                         作者：Thunder\n");
printf("\t\t****************************************\n");
printf("\n");
}</P>
<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
0 S8 M8 [5 `  J1 [//为1的不显示系数1，-1系数只显示负号
int i,j,k;
% E. G. Q: i, }- V8 }. S$ M# _4 V switch(step) {
case 1:
4 E4 B- c; E& h  |! j  printf("min z= ");
' _/ ~: l! n5 l& K  printf("x[%d]",arti[1]);
( n6 c+ A0 n) m+ M, m! J4 E1 m  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
case 2:
printf("max z= ");
0 l! k! A& d- p' q 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);
  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
  else printf(" %lg x[%d]",c,i);
}
break;
}</P>
1 [) M% m4 i- V4 G2 F3 b# D<> printf("\nst:\n");
6 Q6 y8 |* M1 q7 V6 ^/ D0 S for(i=1;i&lt;=num_st;i++) {
2 V; L4 X! W: p- n  k=0;
  for(j=1;j&lt;=num_x;j++) {
: k+ v# m( A9 u) l   if(a[j]!=0) {
- q4 l  V7 }  R7 E+ K* a    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
; s2 U. g  \# m3 ^- p    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);
; U& Q7 q* }# c: z6 D+ j8 y    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;
   }
0 J) p; E! i9 K; h. x+ @- }! a  }
+ z  J7 Q# e% _% f. n0 B; N  printf(" == %lg\n",b);
: W# {! e* `3 O! \$ _ }
( A% C/ p$ R: J! ?+ D0 U" T printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
) d5 N; ~, o' B' n3 Y) x<>void exchange() {
int i;
, X2 T8 I* {3 a" C* ~; D1 U double temp[MAX];
, K" I* [* k/ \7 G4 Z! l for(i=1;i&lt;=num_x;i++) {
  temp=temp_c;
7 z7 ]5 S! ^9 z  temp_c=c;
: G" b# P6 Y: Q  D) s$ w  c=temp;
}
- H7 \% O) x4 L+ }* }}</P>
; [0 z& i3 l+ A7 q<>void create() {
0 g- e. |: K  O9 Q; ]//输入方程组系数，每个方程输完后回显确认
, O5 k# U; Z; l) |0 c# W6 t" a int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
7 y* X+ D0 d1 E& }* I  J+ Y char confirm;

* L, V# H% o+ E& D2 s- @3 f. w9 z while(1) {
5 L! |7 _3 P. `5 o  S; A- z  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  scanf("%d",&amp;ma_mi);
2 ]+ i# C: c' d4 K" S2 H% G) S% {% H  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
2 o: u- X3 c. l  else break;
}

7 \- Q: V: {# M0 O while(1) {
  printf("指定变量个数：");
6 ?1 h8 C8 ]8 r8 J$ W0 O  scanf("%d",&amp;num_x);
  printf("输入价值系数c1-c%d:\n",num_x);
  for(i=1;i&lt;=num_x;i++) {   
   printf("c%d=",i);
   scanf("%lf",&amp;c);
5 i# V" v: Q+ o4 e* w9 [# M& \  }
! a: w. d0 {2 D6 I  if(ma_mi==1) printf("max z= ");
2 f6 d1 ?# @* R$ \3 ~  else printf("min z= ");
; D; d# w, ^$ j  printf("%lg x[%d]",c[1],1);
( c  t. T% `  l, v/ t% F" C6 p  for(i=2;i&lt;=num_x;i++) {
   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
3 F! N0 o. |; _4 k   else printf("  %lg x[%d]",c,i);
. |/ |. L. i: \* I, ^" U1 ^  }
printf("\n正确吗？：(y/n)");
) X) D* e3 w  [/ B: M6 |& C getchar();
% d( h2 Q& I. @, x- O7 a) n1 j confirm=getchar();
0 t% R' a6 O, [% V! r if (confirm=='y') break;
else if(confirm=='n') continue;
}</P>
<> printf("输入约束方程组个数：");
scanf("%d",&amp;num_st);
for(i=1;i&lt;=num_st;i++) {
$ `7 m; u: r8 G1 f% M  printf("st.%d:\n",i);
, h) Y, ^8 f2 E) |* M" A  while(1) {
: X- R3 o8 {1 z- n3 O: N) Y- E; }   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
9 `/ ]* p( e! U' I9 D   scanf("%d",&amp;re_st);
3 V6 N. B. y( _. _0 N( |3 a   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
: A/ d. u: Y6 J% l   else break;
7 R( G7 l/ T' q# @  }
* O- s: d( Q5 t) ]* `/ B  printf("输入技术系数：\n");
, F- K  H6 Q# Z4 e$ o& ^  for(j=1;j&lt;=num_x;j++) {  
- L7 ]( d, K& Z9 n! v   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
  }
/ R; P$ V' X! h7 Z# @1 O- v  printf("输入资源拥有量：\nb%d=",i);
  scanf("%lf",&amp;b);
  
  printf("st.%i:\n",i);
" Q; g8 N8 R2 V2 U: o  printf("%lg x[%d]",a[1],1);
2 E" I! L5 s- \$ R  for(j=2;j&lt;=num_x;j++) {
4 |" V$ p" {: o( t2 ?  p   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
, h. R7 Q, J( f7 i   else printf(" %lg x[%d]",a[j],j);
  }
# k% \  J4 N/ v4 P4 a  switch(re_st) {
* Y, w3 p! V* g8 C! g. ~! H   case 1: printf(" == %lg",b); break;
   case 2: printf(" &gt;= %lg",b); break;
   case 3: printf(" &lt;= %lg",b); break;
  }</P>
<>  while(1) {
   printf("\n正确吗？(y/n)");
, v6 X: I" N: @( Z   getchar();
( |, o5 D9 R! D) P- u' _2 N: Z" B   confirm=getchar();
( H" ^) |1 c* x0 s* t0 E   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
; K4 l# q0 X* e8 H8 d  }
}</P>
+ m9 x+ X& P. h& y' T) x, j: ~1 `<>//显示输入的方程组
0 T2 ~5 b5 o6 ?/ L printf("\n原问题为：\n\n");
if(ma_mi==1) printf("max z= ");
# ~$ A/ u& v# o7 `1 T else printf("min z= ");
printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
9 o7 Y% S8 b  y* m% U  if(c==1) printf(" + x[%d]",i);
% c* B3 G/ X6 W% }( m7 C( E- Y  else if(c==-1) printf(" - x[%d]",i);
- U# {+ a0 G9 ~5 a7 e  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
: B- G3 Y$ n# S  else printf(" %lg x[%d]",c,i);
}</P>
3 N; s8 Q. u5 V' A( T* `, Z<> printf("\nst:\n");
: \( M% F5 f' \% y# P) S for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
! S# s) n6 \3 \" v   if(a[j]!=0) {
$ t1 H0 M( i* |2 X    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
1 w- L4 M. f4 M; M/ L3 H    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
9 e" E1 n$ r* c# z0 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);
( I' V* p8 n  k# 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);
: d/ P1 }8 {; I; w" h: D    k=1;
2 E3 }& I. c7 {$ L' _1 p8 v: ^   }
  }
5 z1 L% }3 ]; ~* ~4 t/ _; i  E3 ^9 y  switch(re_st) {
   case 1:
    printf(" == %lg\n",b);
% I" E: D$ j; J5 ^, s    break;
/ W% j# O) G. A9 ?   case 2:
, y. Z- a$ z/ J    printf(" &gt;= %lg\n",b);
    break;
   case 3:
% Q% ^5 u5 W/ p( i7 h    printf(" &lt;= %lg\n",b);
    break;
  }
}
printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
9 J, X! ^3 x0 @; G# e# a<> tnum_x=num_x;
for(i=1;i&lt;=num_st;i++) {
2 X$ V. T, c/ J9 M  T- Y- X# Q. M  switch(re_st) {
. ?+ F  ?$ L1 {4 B( _- b  case 1:
3 L- i* J0 d8 {9 Y  case 3:
   num_x+=1;
; i' O* u4 u( I+ e5 H- c( |+ N; I   break;
& X% N: i' d$ E9 q' y1 ]" n  case 2:
   num_x+=2;
   break;
  }
}</P>
<>//化为标准形式
- {; H: C  r$ b- p+ h) H. u! I if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
for(i=1;i&lt;=num_st;i++) {
' N+ s+ Y6 H4 {5 U6 O8 A2 J8 `  _  switch(re_st) {
   case 1:
    num_addv++;
    num_ba++;
6 Q; R& p: A( q+ g# K' T/ V: ?+ J    num_ar++;
) \) h# R) }- E* J8 r    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++)
! ?# j0 A: m$ R     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
' s0 V( H- C0 l0 b+ p    break;
! {2 @! ]) g6 t9 ]   case 2:
% I. z! n4 m% `% c; ~# |4 Q9 J    num_addv++;
) s/ ~* d, @5 f0 _/ v" g; W: x    c[tnum_x+num_addv]=0;
    num_addv++;
    num_ba++;
# s. {4 I0 H. ^9 b$ R" e" K8 `9 \  O    num_ar++;
5 E2 \- \! ^3 E1 G& [5 h    c[tnum_x+num_addv]=0;
+ ?3 D9 O# W7 O. f* y: u. `    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
3 x$ O9 R5 ]0 ~+ P. ?4 K& 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;
5 ]0 x. t3 _: F& r& \0 X" I$ {# j% `     else a[j]=0;
    break;
) U6 u7 m1 l5 \   case 3:
    num_addv++;
    num_ba++;
) T% r: [- s/ }4 e7 ~/ v- [    c[tnum_x+num_addv]=0;
, b6 `7 K2 H" k. F+ N  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;
/ E- r0 r: Q. b+ ]3 H" S     else a[j]=0;
; }" X/ }$ _6 e: F" s    break;
  }//switch
6 w8 w5 a4 G3 a( V1 u+ V- J6 s }//增加松弛变量、剩余变量、人工变量、确定基变量</P>
7 z) F5 s4 L) m. D- Z5 Z! L<>//显示标准化后的方程组
printf("\n化为标准形式后：\n\n");
if(ma_mi==1) printf("max z= ");
$ F" C! z* w* J: ?# o) J5 a else printf("max z'= ");
printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
1 p) A( N/ S% e" M: n9 W% O/ {  k=0;
  W, `/ R: v) P+ g" B  for(j=1;j&lt;=num_ar;j++)
% j5 d  {! _* r6 y  if(i==arti[j]) k=1;
  if(k==1) printf(" -M x[%d]",i);
  else if(c==1) printf(" + x[%d]",i);
1 f( W9 C/ p) e+ `" A  else if(c==-1) printf(" - x[%d]",i);
  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
$ p9 U' R* R& s8 v2 P  else printf(" %lg x[%d]",c,i);
}</P>
: z4 y& c2 ~. \: o8 W# s$ `! V2 i<> printf("\nst:\n");
, i2 r, }' z, m4 z! t for(i=1;i&lt;=num_st;i++) {
  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);
! w9 E# @  q3 W$ ^9 u0 Z    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
    else if(a[j]==-1) printf(" - x[%d]",j);
3 C+ u+ i# U8 U; t3 d    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
1 I" ^/ j- V- l' R, s$ 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);
; y  {6 B$ f& g$ Q' s! I; [' b0 ~    k=1;
   }
; s& O# h) `( M4 o) g# H( M  }
  printf(" == %lg\n",b);
& T* `3 `. v) v- o$ O+ W# \- K }
printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
: b4 l: r- f; D<>void iterative() {
  }3 {( ], j* I; Q int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
& Q! n) B: z* \, {% r7 C1 I int base_elem;
int base_out,base_in;
2 e3 z* k! u! v- B. b; A5 }* X double sigma[MAX],temp;
, n1 n9 m- a% ^ double value_be; //高斯消元里保存主元素值</P>
, I  K. G8 Y2 m( ?6 I. F<> printf("\n\n第%d次迭代：\n\n",stop);
for(i=1;i&lt;=num_st;i++) {
; B. M) I0 D8 P2 r7 @  printf("c%d=%lg\t",base,c[base]);
8 D$ I. j/ V& ]* X; ]6 ~  printf("b%d=%lg\t",i,b);</P>
<>  switch(step) {
   case 1:
    for(j=1;j&lt;=num_x;j++)
2 w( L! k' ?" @; Q2 w, r    {
6 c; `2 l5 s1 h* P7 d     printf("a[%d][%d]=%lg\t",i,j,a[j]);
" g0 B. q9 L1 R1 K6 N1 E& a    }
% }9 r1 q! c. s  G, ^    printf("\n");
8 i5 U: q& Y; `% j4 J" k    break;
   case 2:
  u/ j9 S6 n8 T; X& v& l/ k' Z% n    for(j=1;j&lt;=num_x;j++) {
2 C0 p- a, c3 \% K: B( A, @+ ?% b     k_a=0;
     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]);
/ ?' N: b: I) O    }
    printf("\n");
    break;
  }
}
, H0 r/ y% h* T) U+ J0 ?7 M//求检验数sigma
8 y, S7 o; x, W; { for(i=1;i&lt;=num_x;i++) {
5 E: j* K, d2 h( ?1 b2 d( n; Z( F  sigma=c;
8 j9 _1 i) T2 G# X" s  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;
, Z! j; d3 \  n  switch(step) {
+ `9 E" \2 `. B9 i+ U   case 1:
    printf("sigma[%d]=%lg\t",i,sigma);
    break;
   case 2:
    k_a=0;
    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
5 y& N: O' R& d1 q    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
3 W# o* m( `+ b" ^' @( e' e    break;
  }
# B" M! e* I+ n5 v: t/ O+ d0 W% d }
+ m5 k: o4 Y; h9 e! p0 l, f putchar('\n');
* I, O. k3 Q8 E4 c3 u//检验检验数sigma是否全小于等于0
) ~: b& M, W9 b3 J9 V% L) i8 X k=0;
for(i=1;i&lt;=num_x;i++) {
5 c/ S; I& w# S0 T  if(sigma&gt;0)
3 M( _: ^, o( s8 M8 E4 n+ i  k=1;
}
if(k==0) {
- ^; ?  u. ~. q: w9 F& x7 _) |+ C0 \  //sigma是全小于等于0时，检查是否为无穷多最优解
  for(i=1;i&lt;=num_x;i++) {
   k_f=k_a=0;
   for(j=1;j&lt;=num_ar;j++)
% ^- g+ B' X( Z5 y& D   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;
" x+ N6 k3 X7 _' j/ r. \   if(k_f==0) {status=-1; return;}
4 Z  O; @- s+ S2 m   }
( L" T% f- [3 U0 f  }
& A( }. V# A; T/ b3 J. c  status=1;
! T. S' x& j+ \6 H  r4 L5 \  return;
}
//检查是否为无界解
for(i=1;i&lt;=num_x;i++) {
  k_f=0;
- F8 B* F) _5 A5 n4 u  if(sigma&gt;0) {
; e1 h6 ]/ _% U! g0 g* `2 u   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
* [* o: C" L* J( @5 b. d# }   if(k_f!=1) {status=0; return;}
) |% I# P: m+ K  m4 P9 `  }
( F: ]# O& D, V4 o: G- u3 U& J }</P>
7 ]! p- D) j& O$ F" K<>//确定换入变量
; w* ]% f* s( `5 k; b for(i=1;i&lt;=num_x;i++) {
: E, \" ^/ H$ @  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
2 A7 n$ V9 S: I; N0 E5 j5 o+ u  I  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
}//temp赋初值
8 r3 n9 R! S8 {- K# _  s for(i=1;i&lt;=num_x;i++) {
5 C% q6 Q  m3 P! S) d% I% p  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  if(k==0)
0 U- ~3 J, M/ J: [   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
  I/ R- m; F3 o8 X; {    base_in=i;
2 B: n# E/ G; ]3 _( [; X    temp=sigma;
   }
}</P>
* D5 k* Y* ]/ |  Y<>//确定换出变量
for(i=1;i&lt;=num_st;i++)
, i$ s: z5 t8 j8 j1 V8 g  if(a[base_in]&gt;0) {
   temp=b/a[base_in]+1;
4 y8 H* {; V7 X  f; P  B   break;
  i8 j8 f9 S, j) N/ z8 h9 `  }//temp赋初值
+ ]' \3 ?8 ~: z! c$ a1 y2 I for(i=1;i&lt;=num_st;i++) {
  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
% J  l3 l" ~3 e, z  b- {6 H# \   for(j=1;j&lt;=num_ar;j++)
    if(base==arti[j]) {
     base_out=base;
     base_elem=i;
$ ^* [! |" s, N5 P     temp=b/a[base_in];
     break;
! w2 t) N$ c5 ]' `    }
  }//人工变量优先换出
  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
3 p) n& I- _* f6 R2 \$ \   base_out=base;
   base_elem=i;
" I# |, A$ O# F9 p. J9 M   temp=b/a[base_in];
( V4 j9 {& C+ P/ E3 U  g  }
  D: L) C% J4 R4 t+ J9 R7 p3 l# ?) ~ }</P>
7 [( \! r! H* x7 S* h<> printf(" 基变量：");
9 }' \. V5 V% s, U" n for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
  `) Q5 g: Z% ?% W7 k1 i! Y printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
//基变量变换，进行新方程初始化后迭代
for(i=1;i&lt;=num_st;i++) {
4 s* ]: L& z) D0 B! _8 Q& {& F; f  if(base==base_out) base=base_in;
}
//初始化主元素行系数
+ k, B$ E& H5 X: Y) Q" } value_be=a[base_elem][base_in];
b[base_elem]/=value_be;
9 ~, t. f5 ^, U7 Z2 g for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
<> for(i=1;i&lt;=num_st;i++) {
  if(i!=base_elem) {
   b-=b[base_elem]*a[base_in];
& r5 X6 Q  o1 A' |- V   value_be=a[base_in];
   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
) p3 H+ V, W, C# B4 q  }
}
) A3 T! u* s& I4 \0 Q7 _# n stop++;
if(stop&gt;STP) {status=-2; return;}
* j; y0 s7 ]$ L8 Z, |8 T% E! @  R iterative();
}</P>
: F, v$ k4 @2 a: f4 [1 Z<>void output() {
9 T; A. w1 P# A; _0 _. g4 H int i,j;
; ?  A* S- [7 `# V double X[MAX];
4 A. N8 r, g! `. B* |* S6 \+ ~0 X printf("\n结果如下：\n");
6 O/ R. S: N: b+ A printf("\nX=(");
5 t) K% o8 e2 W. `/ z- f for(i=1;i&lt;=num_x;i++) {
  for(j=1;j&lt;=num_st;j++)
3 a  n8 [5 ?2 q+ Y8 Q+ R! ~  if(i==base[j]) {X=b[j];break;}
  else X=0;
  printf("%lg ",X);
1 D5 y3 `6 C  D/ O" n  N4 u }
& r( ~$ X5 ~' O6 [7 N; I7 k printf(")");
6 O) v3 @% W6 S1 M# O2 l" o4 J1 o for(i=1;i&lt;=num_x;i++) max+=c*X;
  c9 W' ~9 P7 Q: t+ }' A 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？帅呆了