求助：单纯形法的程序(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])
/ A) D4 e! L( C{
  int j;
: z! F5 _) {) y4 _- Q2 ~  *c=0;
5 o( e9 _3 ~3 f/ p# {  for(j=1;j&lt;=*is;j++)
    {
" Z9 k) }6 s( P" Q! z- ?      if((*a)[*ir][j]-*c&gt;0)
( N2 t  R: T) o% x    {
      *c=(*a)[*ir][j];
3 a8 x# k2 l; x  U4 i6 E$ Y3 `; B6 g      *j0=j;
' b, P( V" [, \8 a' b9 F    }
     }
) N0 A1 I# Z5 A% {" d* e; g; O}</P><>      /***************** 参数说明 **********************/
      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
* @: h% M0 F( d( @# O0 o+ z      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
) V$ F) j1 l; I      /********* <a href="http://happyyangxu.home.sunbo.net*******/" target="_blank" >http://happyyangxu.home.sunbo.net*******/</A>
1 E9 V# F, `6 m% g6 \int xi_sm(int m,int n,int m2,int mn1,int l1,float (*a)[X][Y],
       int(*k)[],float(*x)[])
+ ^9 O" E  g" d{
  int i,m1,mn,j0,i0,j;
  float c,g;
5 P5 f3 l7 F- x/ |2 k" C# j  m1=m+1;
" A1 `. F+ a! P  mn=m+n;
+ V' A+ Q' s+ z- l' g  for(i=1;i&lt;=m;i++)
, l7 @% a  t* e& K9 o* d; p    (*k)=n+i;
/ q5 t9 D/ G! |leap1:
2 }9 F) |/ l3 {4 h  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
9 }8 ^6 z: y1 h2 w1 G2 G6 |( t8 V  else
* x) N+ U( ]( G& |1 l     {
leap2:   if(l1-50==0)
# q7 C  m! ?' H      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
       else
8 U& }6 P0 M6 `* t7 |# j3 Q      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))
/ ^" K1 C6 n8 D' n) Y      {
    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
    return(-1);
      }
% S- W  ?9 |' F* ]% L0 O- o9 S8 O   if((c&lt;0)&amp;&amp;(l1==1))
4 E& l! W) \: r: n. @4 j       {
) v% }' k% E: `% c     l1=50;
( \! E- G9 D: Q. G" I/ v# e     goto leap2;
: m8 @) d2 Q0 ?, |: ~       }
    if(c&lt;=0)
7 ]0 b; F& E" I0 a' Q+ f- T       {
* n# z: k2 ~! l; J$ _- i$ u7 ~      for(i=1;i&lt;=m;i++)
         (*x)=(*a)[mn1];
$ J$ `9 y4 u1 U! v3 [         printf("\n\t********Optimal solution********\n");
+ g* }6 i- C; H/ w# T- y1 q      for(i=1;i&lt;=m;i++)
: A2 f& W0 F6 F3 A8 M& p  ^         printf("\ni=%d,x=%f\n",(*k),(*x));
) h- L' V. h7 w/ c1 Q7 G, `/ X      printf("\nF=%f\n",(*a)[m2][mn1]);
      return 0;
    }
     c=1e8;
- B# t7 |  e; q6 V$ G3 s" L     for(i=1;i&lt;=m;i++)
       {
     if((*a)[j0]&gt;1e-8)
        {
          g=(*a)[mn1]/(*a)[j0];
- {* Q3 m* k. W- ?+ G- T          if(g-c&lt;0);
. {( J* m+ A8 H% [1 U% Z3 w1 ?& s         {
! G' k  e9 b! A& |3 P           c=g;
. ]" w0 ^2 n+ o9 I" |           i0=i;
          }
         }
) V; z$ J; s( R# y9 q    }
      if(c==1e+8)
    {
: r; E) P1 }/ f3 }      printf("\n\t*******Lp no solution********\n");
      return -3;
7 L  B: C9 d3 i     }
5 u+ K# X: G& A8 p       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
8 d  B: r# p4 C) F7 C7 o" r      {
        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
% C* l1 v6 ]6 L3 y2 t; d8 L* T           continue;
         g=(*a)[i0][j]/(*a)[i0][j0];
           (*a)[i0][j]=g;
         for(i=1;i&lt;=m2;i++)
           {
1 [# Y3 r, f6 B0 s1 a/ N, z         if((i==i0)||(!(l1==1)&amp;(i==m1)))
" D) u% C3 P# T) E" Z            continue;
! U, O  O3 L3 y2 T! t( h& s         (*a)[j]=(*a)[j]-(*a)[j0]*g;
  ^$ c% m( A& j/ v; w4 y9 B           }
       }
       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
# R8 E3 u7 k& W* A7 T- r% Y* q6 K       (*a)[i0][j0]=1;
  @8 e2 `  o  D7 D& i; @       goto leap1;
7 R; m3 z* v( _* Y0 u}  
main()
, R; F6 X% B( \; `, L# _$ R/ g{
3 F2 p  {7 x4 n- y2 @' ?  float a[5][7]={{0,0,0, 0,  0,0,0},
6 q7 y$ O' i3 l2 @8 Q3 t  `                 {0,3,-4,3,  1,0,12},
( B; l+ q) X+ y/ M1 r7 j9 q9 f5 d                 {0,3,0, 6,  0,1,12},
                 {0,0,0, 0,  0,0,0},
2 H: S1 s# Z# q, [6 C                 {0,69,0,144,0,0,300}};
  int k[3];
; b; D8 D/ p. B$ ^8 _6 S  float x[3];
: \& q3 S; J2 P0 v  clrscr();
9 B. U/ m) v! K: W( |. O  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
}
</P>

plgatc

<>/*************************************************************************
                    单纯型法解线性规划问题（两阶段法）        
, W4 N& L& }' i7 N
: `4 E0 w/ J( P( n% U 编程环境：VC++6.0     
方程组输入说明：
变量非负，按提示输入相关参数。
4 R4 a4 N) {. y. ~) |1 X *************************************************************************/
. ~' S% i8 k$ W#include &lt;stdio.h&gt;
2 a' ~) J# Y  p- t#include &lt;stdlib.h&gt;
% R) T' I- D' w+ ?, c* f$ X* j#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; //方程组相关系数
, x7 q4 s' r2 u# `/ xint num_x; //变量个数
: ?- a" ^. B  Q/ A$ h$ {% Qint num_st; //约束方程数
int num_ar=0; //人工变量个数
int arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
0 a# y% Z( ?. ?, k, }<>void create(); //建立方程组
void iterative(); //单纯型法迭代
void output(); //输出结果
void banner(); //打印程序标题
void exchange(); //交换两阶段价值系数
void show(); //输出方程组</P>
<>void main() {
int i,j,k;
banner();
create();
//保存原价值系数，转换为第一阶段价值系数
4 z( {1 o/ k- S' c1 r1 k for(i=1;i&lt;=num_x;i++) {
) X; v0 m1 y& z7 q; s3 f% X  k=0;
  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
9 R/ X  x3 O( j6 Y- M9 L3 t9 Q0 o  if(k==1) temp_c=-1;
1 a; V' d6 E5 X5 g! p& r- z  else temp_c=0;
}
' G) H7 s8 Y" m exchange(c,temp_c);</P>
" w  f. r( X4 U; l2 k1 {8 e<> printf("\n\n第一阶段问题为：\n\n");
4 A# d' Y: P9 A5 N& c  E4 m show();
1 o- h: y( u# `6 l- a7 ? step++;
printf("\n\n按回车开始第一阶段迭代");
5 L3 I" G' S+ I# Z+ E" g getchar();
1 q, S# a9 U( V( z2 C: g9 e1 M getchar();
iterative();
if(status==-2) {
( D$ b) p! X2 N4 v9 A0 Z# f  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
" Z8 W- J7 M4 t: G8 l9 Z2 h) I$ W* Z  getchar();
1 W. M+ M8 R& _3 Y  exit(0);
8 }% P* B( [* Z8 @# \ }
output();</P>
<> if(max!=0) {
! z* G7 `# \% C2 b% x/ i7 b9 e  puts("\n\n原问题无可行解。\n");
( V! S/ e" ?7 M7 }1 w4 B$ g  puts("\n按回车结束");
  getchar();
5 T8 |- t* F' [: X6 Z+ ^0 y  exit(0);
}</P>
<> //转换为第二阶段价值系数
exchange(c,temp_c);
2 o8 a4 W9 M: V' T/ Y //把人工变量列全设为0
for(i=1;i&lt;=num_ar;i++) {
  g: ^. |/ c& F/ Z  c[arti]=0;
  s2 @+ u% [$ i, _  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
}</P>
<> puts("\n\n第二阶段问题为：\n\n");
$ `  r3 g! v" z show();
puts("\n\n按回车开始第二阶段迭代");
' @6 S8 ^) e$ u2 W getchar();
+ u8 F6 U& m- l3 @ iterative();
switch(status) {
, c8 ?2 I5 u* { case 1:
  output();
  |7 o2 I$ `7 s  puts("\n\n原问题有唯一最优解。\n");
$ }3 H* U1 w" D, f  puts("\n按回车结束");
  getchar();
5 L( z' g$ n$ z  exit(0);
case 0:
$ h6 A# X- K1 C; W6 `" I4 u  puts("\n\n原问题为无界解。\n");
4 {& a# A% e. _" l8 v) E& y. t  puts("\n按回车结束");
9 d9 n/ O# U4 ^2 {7 m  getchar();
  exit(0);
/ d+ Y5 c$ Q. U% H: ]1 f! z2 E; o case -1:
0 x0 B: b5 N# `) B/ @2 C  output();
  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
% Q( M& j$ A6 M' {  getchar();
1 N$ U9 L" t' L8 i; v# J7 A) s; b  exit(0);
case -2:
  puts("迭代超过限制次数强行终止！\n");
; y, \0 c& e  U) ]. g% r& t  puts("\n按回车结束");
5 O, l/ Y0 ^' F/ @% K  getchar();
3 @; d% [. k3 G9 k8 a* q5 H# p  exit(0);
$ G9 A* w- t6 A& h! D5 E+ F6 r }//switch
3 D0 r. [; {. b: s9 J}</P>
<>void banner() {
printf("\t\t****************************************\n");
printf("\t\t         单纯型法解线性规划问题\n");
7 i0 s4 T/ a3 I printf("\t\t                         作者：Thunder\n");
" K- W! `  R: e; a+ ] printf("\t\t****************************************\n");
# a' O* I" |6 C# G printf("\n");
}</P>
4 I; o* j& A! q<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
4 }+ C: J( ^# v4 N//为1的不显示系数1，-1系数只显示负号
  q  s" r* r' G0 A% }' A int i,j,k;
switch(step) {
case 1:
+ }) p$ _: B+ M& n7 H, O  printf("min z= ");
  printf("x[%d]",arti[1]);
  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
+ @% b) e# d) F! Q. p! n  break;
case 2:
" @: g) N  D- z8 \6 w printf("max z= ");
9 P# j2 |" v. I printf("%lg x[%d]",c[1],1);
1 S; i; l9 ?. ?9 g$ i$ [: U for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
/ M$ \( a2 k1 n  |1 X7 h  D* E  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
, n& a7 E; i% q& ]% J  else printf(" %lg x[%d]",c,i);
}
break;
) W* d2 r3 i+ L& b }</P>
. v% O2 j3 q+ \7 z" X<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
+ h: J, n1 Q+ ?( q7 y% P  k=0;
  for(j=1;j&lt;=num_x;j++) {
6 g: _, S1 S, G, v   if(a[j]!=0) {
  `5 H/ T1 w7 ?/ H+ A- O    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
  B8 w0 D! R- x+ {: w9 r! M    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);
8 c' c& n' ?8 e7 C    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;
   }
  }
  printf(" == %lg\n",b);
9 M5 K8 C. [) b4 M9 H6 S1 e6 F( ?) r4 l }
printf(" x[1]~x[%d]&gt;=0",num_x);
! ]/ ?$ A* b" X. s}</P>
<>void exchange() {
" q. q) g( ?# `0 V! ]3 Z- t+ d: y int i;
double temp[MAX];
' \# R# r+ g/ |9 w8 Q/ I for(i=1;i&lt;=num_x;i++) {
% i% D* C6 x* s  Q2 u5 J3 {/ l; m# U  temp=temp_c;
  temp_c=c;
' {  s' W- i, E$ Z' ]" [  c=temp;
/ w6 j5 n4 s* i5 Q }
) s% }7 d# Y+ y) U2 n}</P>
<>void create() {
//输入方程组系数，每个方程输完后回显确认
int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
4 |) [( p; o4 }' T$ D char confirm;
! I4 _' {0 B; K
/ H( {/ y" Q4 \  \6 T, m2 L while(1) {
  printf("请选择：1、求最大值，2、求最小值：(1/2)");
1 `  }) j* ]* a& ^6 F  scanf("%d",&amp;ma_mi);
# }& r8 n% [6 r% m' L8 D: A  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
6 m0 o( f7 c" m  else break;
}

& K$ d$ V* a, S5 _5 u while(1) {
  printf("指定变量个数：");
6 y: ^, k4 e. I* m4 O6 J. e  scanf("%d",&amp;num_x);
  printf("输入价值系数c1-c%d:\n",num_x);
6 k8 J; b* m! N  for(i=1;i&lt;=num_x;i++) {   
% k/ S& n* g5 s4 e8 l   printf("c%d=",i);
! l) R, V$ S. u- Y& `1 D) D" X" `/ p( J* [   scanf("%lf",&amp;c);
  }
  if(ma_mi==1) printf("max z= ");
2 {: V! j- f2 T2 B% L) @6 N  else printf("min z= ");
3 [  K. \9 @5 E7 Y/ C, I# q  printf("%lg x[%d]",c[1],1);
3 J  A3 Q) ]; H7 R) o' ~& x  for(i=2;i&lt;=num_x;i++) {
+ ^3 A9 {% U* K. }   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
) k1 i4 E: S& K0 s5 d$ h( r7 Y) \   else printf("  %lg x[%d]",c,i);
  }
printf("\n正确吗？：(y/n)");
5 v/ o% H) V6 l getchar();
! ~; {9 Z1 J* W+ D5 y9 ~ confirm=getchar();
if (confirm=='y') break;
else if(confirm=='n') continue;
}</P>
<> printf("输入约束方程组个数：");
scanf("%d",&amp;num_st);
for(i=1;i&lt;=num_st;i++) {
/ s8 ?9 @0 M/ W1 w" F1 b  printf("st.%d:\n",i);
4 D; f* P0 T3 k! m$ ?2 i  while(1) {
( `( S. A: v3 T; W+ d   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
   scanf("%d",&amp;re_st);
   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
   else break;
  }
  printf("输入技术系数：\n");
9 |3 K; v- q1 D! `# q0 G# }  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
( ?3 W' a2 ~$ Q  }
8 n7 k0 u! b8 v. e" K  printf("输入资源拥有量：\nb%d=",i);
  scanf("%lf",&amp;b);
  
  printf("st.%i:\n",i);
  printf("%lg x[%d]",a[1],1);
2 X& A6 L) y: _4 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);
9 t  S; |/ B: Q% ]; R  }
  switch(re_st) {
6 }* d% R' U& Y+ |8 V   case 1: printf(" == %lg",b); break;
   case 2: printf(" &gt;= %lg",b); break;
) Z# L& O4 O" N" o9 a! @   case 3: printf(" &lt;= %lg",b); break;
  }</P>
) _% g* R$ U; T5 i6 l' ~8 Q+ Q<>  while(1) {
   printf("\n正确吗？(y/n)");
   getchar();
, t7 }  e2 I6 f, l0 R   confirm=getchar();
   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
& x2 Q: U8 w- J& K3 [  }
}</P>
/ ?# @/ e* N- g<>//显示输入的方程组
printf("\n原问题为：\n\n");
if(ma_mi==1) printf("max z= ");
else printf("min z= ");
  o7 [) F( j2 J6 \7 Q6 d( ]3 U printf("%lg x[%d]",c[1],1);
; U* @( K; f7 n' i( G" [0 B/ }3 t for(i=2;i&lt;=num_x;i++) {
  X% W( F% N( _0 g8 D! Y1 _4 ?6 @  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);
}</P>
<> printf("\nst:\n");
8 \: P9 t5 {3 X, u for(i=1;i&lt;=num_st;i++) {
* l8 Y! D& r: o- @$ o" H4 p( ^  k=0;
  for(j=1;j&lt;=num_x;j++) {
7 k8 ^; t) K6 u* `/ v' }' y7 ?+ q   if(a[j]!=0) {
- e8 i0 J0 w0 Q: E( l  F# h7 I0 S    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);
4 l- _- g* b4 p    else if(a[j]==-1) printf(" - x[%d]",j);
" i4 m8 n* A- z) ~' x    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
9 q8 A) q2 e+ I- ?    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;
. }4 V, k5 u$ ^* J& J   }
  }
. j9 u- w7 K. y3 A) x  switch(re_st) {
: ~" i( c1 f6 r- U" e   case 1:
- R  H4 Q8 z( s8 M8 z    printf(" == %lg\n",b);
: n+ U5 q3 r. e3 s    break;
0 d, \! R: |0 o   case 2:
) K1 w& D4 M  n  |3 t& R" m( c    printf(" &gt;= %lg\n",b);
    break;
   case 3:
    printf(" &lt;= %lg\n",b);
" x$ ^; X8 T' h- U0 s8 P    break;
* T$ N+ v0 @' i1 a# ^6 E+ ]- k  }
) n7 l# e+ X) Z- c+ a9 X }
7 R' o+ u8 p0 B/ } printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
<> tnum_x=num_x;
; T" V4 l3 C' j2 |8 l; s2 u% w% m5 \ for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
  case 1:
. s7 j% q0 q/ I9 r& Z  case 3:
3 G* P( _( c" X( V& m9 s: c   num_x+=1;
; e1 M' b! H( G0 k   break;
  case 2:
' A% |( m( a. H/ m9 ~   num_x+=2;
4 y2 P3 O- J  t3 R0 `* s7 K5 E9 S! l   break;
) e; b( X  a7 R. k! H+ M) W7 W7 Q& p  }
}</P>
- `5 y0 D: B: G<>//化为标准形式
* r# }! ^7 U5 |. C# Q/ J2 O' Q0 ` if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
for(i=1;i&lt;=num_st;i++) {
0 ?( k5 G" }  E8 p* T0 q) i  switch(re_st) {
   case 1:
9 v) f5 `# `( j# F8 Y! g    num_addv++;
    num_ba++;
, n1 H& Z1 {1 \7 u6 S% }    num_ar++;
    c[tnum_x+num_addv]=0;
9 W9 g( ^# w( {6 s    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
" J( a. O6 r7 @" m+ M; X    for(j=tnum_x+1;j&lt;=num_x;j++)
     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
4 V7 _) B* L; |* v" Y     else a[j]=0;
    break;
& S  [/ b5 t& R8 I   case 2:
5 v4 w4 G7 `! M. H# o. ?    num_addv++;
& M; P$ x& `- A9 w+ w; f7 E    c[tnum_x+num_addv]=0;
    num_addv++;
& E& j! B. {( d% \' S    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++)
8 d. o2 W+ i7 m: r  R) {     if(j==tnum_x+num_addv-1) a[tnum_x+num_addv-1]=-1;
/ |9 K( k# N! V5 T2 o     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
2 q9 B5 g- c0 B- x7 ]     else a[j]=0;
3 x+ x6 s" N, Q7 I    break;
+ @: c- l5 M0 F& o0 f   case 3:
    num_addv++;
    num_ba++;
0 W. I# x' `2 g$ o- Q    c[tnum_x+num_addv]=0;
) @$ ?) `% \1 o- m' k    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;
9 V5 \$ Y% t7 n9 u; y  y* W5 e     else a[j]=0;
3 `  g6 Z( m4 \- q2 }( ^    break;
$ [1 r/ Q, O7 x8 r) Q# h  }//switch
}//增加松弛变量、剩余变量、人工变量、确定基变量</P>
1 ^* e2 @: g$ ?7 I/ X: \<>//显示标准化后的方程组
" Z0 W5 Q. f! d4 o printf("\n化为标准形式后：\n\n");
if(ma_mi==1) printf("max z= ");
; v: i% N) J2 b$ i7 f1 O5 U! D else printf("max z'= ");
8 D# Y* X' G) a6 |6 | printf("%lg x[%d]",c[1],1);
. I4 l8 D4 p% d! e for(i=2;i&lt;=num_x;i++) {
  k=0;
* {" W8 w0 Z/ v2 h9 f  I( q9 Y+ S  for(j=1;j&lt;=num_ar;j++)
  if(i==arti[j]) k=1;
" }5 j3 S0 c# F: m) `( g+ Z  if(k==1) printf(" -M x[%d]",i);
& P, G$ A6 ]: r) F  Q/ N- Y( r  else 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);
}</P>
<> printf("\nst:\n");
2 T0 v, }( X+ O3 T% `/ F9 @9 {- l/ P 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);
+ q% v: F: \0 Z    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
    else if(a[j]==-1) printf(" - x[%d]",j);
- Q7 A# a1 B- y1 ^6 @% a, G    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
. s( K4 @  K8 L, c8 g% V    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
$ j3 }& S7 R) [, h    else printf(" %lg x[%d]",a[j],j);
4 U6 W1 G0 z  \* I    k=1;
   }
% ~! E3 ]0 L9 n  }
  printf(" == %lg\n",b);
}
7 ~2 ~8 H( f; [$ q: z printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
<>void iterative() {
5 X' Z9 Y5 J3 } int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
: ?( `8 Q0 q% I, _- l  {3 ?9 W int base_elem;
int base_out,base_in;
! i) A' |5 z0 M: D. W8 \ double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
6 K7 x) S* z: w/ `% D5 @+ c<> printf("\n\n第%d次迭代：\n\n",stop);
3 j1 V7 R4 w7 p6 E0 e for(i=1;i&lt;=num_st;i++) {
  printf("c%d=%lg\t",base,c[base]);
  printf("b%d=%lg\t",i,b);</P>
, Y7 K% g% p' Y; M<>  switch(step) {
2 w! J* _6 ~! G6 y0 o   case 1:
    for(j=1;j&lt;=num_x;j++)
. s/ U2 {; k% `5 N2 g# P( i' c' K    {
# u2 G$ S  S$ o8 F) f     printf("a[%d][%d]=%lg\t",i,j,a[j]);
    }
    printf("\n");
$ N8 h( E; k8 v+ y; ~' _" i    break;
   case 2:
    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;
' d' _- l# A& O* M$ r     if(k_a!=1) printf("a[%d][%d]=%lg\t",i,j,a[j]);
    }
8 J# m8 O8 c' L% G( q    printf("\n");
    break;
, T5 u" d3 |# Y8 d% i  }
}
& G3 z: z& I! }% C//求检验数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];
9 R& m& f0 y0 m1 R9 s* q8 P  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
  switch(step) {
' y1 U+ F* ^* C  F; c4 R) }   case 1:
$ I: l& k, w1 F: i. m& o9 Q. O    printf("sigma[%d]=%lg\t",i,sigma);
    break;
   case 2:
- D. j5 @; e2 W2 ^5 l    k_a=0;
    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
, A* e5 i: x$ N' D, D2 p. L    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
) E) j0 x& |0 o1 d9 s    break;
% v$ x4 d- w' }! X- G: r  }
+ @, b: X( T- f; J- G' c1 e }
putchar('\n');
//检验检验数sigma是否全小于等于0
k=0;
# X: p4 M# X1 d3 c1 X5 E for(i=1;i&lt;=num_x;i++) {
  ~- P: P% Y' o4 s$ ^  n# a  if(sigma&gt;0)
  k=1;
}
if(k==0) {
6 L1 p$ n0 v8 C: x2 C4 g  g) n6 i' A  //sigma是全小于等于0时，检查是否为无穷多最优解
  f7 N6 j4 x. ]. L  I0 x4 N$ ]+ `! v  for(i=1;i&lt;=num_x;i++) {
   k_f=k_a=0;
   for(j=1;j&lt;=num_ar;j++)
$ _& n- g  p; Z   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;
, J3 a* i2 }# b" `* ]  E   if(k_f==0) {status=-1; return;}
   }
  }
  status=1;
* y% W: _; ]0 |; N* o- t6 u  return;
" p6 H7 Y0 B% f, o }
//检查是否为无界解
5 ^* d4 ?! Q9 V$ P3 M8 w for(i=1;i&lt;=num_x;i++) {
) z& P' z3 V$ R. y/ e4 I  k_f=0;
) l& _2 x1 E. {( ^# t- ~0 L  if(sigma&gt;0) {
5 S  ^' p7 v3 a# x/ J# |6 M   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
# e  u% t: T( y3 G   if(k_f!=1) {status=0; return;}
; H- V9 W" u' j* w, w6 h  }
) h% u3 I/ e, M' e7 a }</P>
<>//确定换入变量
for(i=1;i&lt;=num_x;i++) {
, Q9 s# _  \! a& m  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;
+ l: h' D; M: f }//temp赋初值
for(i=1;i&lt;=num_x;i++) {
  k=0;
6 X+ d$ n8 L$ c2 i. t' I: i' c  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  if(k==0)
* I6 {2 p1 X. A7 {   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
    base_in=i;
# V8 }% }. W, B! d  Z& x# P6 g    temp=sigma;
   }
}</P>
<>//确定换出变量
for(i=1;i&lt;=num_st;i++)
  if(a[base_in]&gt;0) {
  E0 W3 [/ d  `! S$ l# L0 ^( t   temp=b/a[base_in]+1;
   break;
/ j: L" ~; u; ?3 E- j  }//temp赋初值
* z1 G  M5 h% @" l1 }- n$ { for(i=1;i&lt;=num_st;i++) {
  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
   for(j=1;j&lt;=num_ar;j++)
    if(base==arti[j]) {
     base_out=base;
     base_elem=i;
     temp=b/a[base_in];
: U& r8 l( g+ b  x; X( y     break;
    }
  }//人工变量优先换出
6 ~+ P& K; t% r' y4 G  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
   base_out=base;
% @- g; M2 w) V! q   base_elem=i;
   temp=b/a[base_in];
5 p+ B& P; e4 o- {  }
6 d/ {, V* R: X3 ?' i. i, ` }</P>
; ~$ J8 A5 ~% J8 k9 ?" E8 L: F<> printf(" 基变量：");
for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
//基变量变换，进行新方程初始化后迭代
9 P# u2 ^: L( i0 a% D& l3 f; x for(i=1;i&lt;=num_st;i++) {
+ g6 p$ l! [4 V& i& I  if(base==base_out) base=base_in;
}
//初始化主元素行系数
value_be=a[base_elem][base_in];
8 A8 S3 g( e- C2 h, P+ l+ K b[base_elem]/=value_be;
2 h" P3 D" r& V. ^ for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
<> for(i=1;i&lt;=num_st;i++) {
, F& q8 Z+ b4 w1 S  if(i!=base_elem) {
% A  D1 h# `1 Q+ w6 Y2 _   b-=b[base_elem]*a[base_in];
6 }0 m8 ?; @5 C3 q7 l   value_be=a[base_in];
. C" J; }0 `, q- a1 Q   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
9 I& O& U3 N. h. \  }
' u) G6 k; j5 d9 ~7 m) w }
1 a6 h; U% l1 m& N$ ]/ m$ \5 B stop++;
if(stop&gt;STP) {status=-2; return;}
" c7 j9 T& c1 h3 U iterative();
}</P>
<>void output() {
1 b- E- y: N8 B5 m int i,j;
% n" R* r  T6 f7 X/ O double X[MAX];
' F  a- K2 f6 G printf("\n结果如下：\n");
printf("\nX=(");
for(i=1;i&lt;=num_x;i++) {
  for(j=1;j&lt;=num_st;j++)
  if(i==base[j]) {X=b[j];break;}
  else X=0;
. Y- b/ v; e/ \, n5 |  printf("%lg ",X);
}
printf(")");
8 M# p1 v$ X. ?7 ? for(i=1;i&lt;=num_x;i++) max+=c*X;
: B2 k6 F8 T6 u8 i" j0 X; X% T if(ma_mi==1) printf("\nMax z= %lf\n",max);
; b6 m1 V1 w/ @' j/ B) e& ~9 }6 j* ~ else printf("\nMin z= %lf\n",-max);
}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了