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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
#include&lt;math.h&gt;
#define X 5
, l# {# U+ U! p( U$ M#define Y 7
4 F# W* f; s) B% j: J  j" M, ~void xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
{
+ H( Y! r+ H2 b  int j;
  *c=0;
  for(j=1;j&lt;=*is;j++)
    {
      if((*a)[*ir][j]-*c&gt;0)
+ c8 ^# h, Y. e1 O! S    {
8 a/ F( \/ k4 b, @( ^  k      *c=(*a)[*ir][j];
. S/ M1 @) W; _' Y( I  [1 O) y- a      *j0=j;
3 E/ W3 _9 z) Y/ Q, I" ?    }
     }
3 v7 w, q, a4 s; A9 D- n4 u1 j4 ?# T}</P><>      /***************** 参数说明 **********************/
: H/ o5 y" r, p% c      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
$ B# }6 U. L* A& W% U5 }      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
# Q7 N  X, k  f0 E5 F      /********* <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)[])
{
  h7 X. v0 V2 [% _7 A5 ^  int i,m1,mn,j0,i0,j;
5 J- f. }/ T$ c3 |' V  float c,g;
8 Q3 K& }% P& c# n  m1=m+1;
  mn=m+n;
  for(i=1;i&lt;=m;i++)
. x3 f1 y  u7 }    (*k)=n+i;
! U) T& ~; d+ n  \2 cleap1:
  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  else
+ h2 ?( r! ]+ T. f     {
( p; C9 d: U6 X* U: U5 s2 Ileap2:   if(l1-50==0)
      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
       else
2 r* k; l  k+ r! w) i$ e      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
     }
$ R  y& n5 m/ O# v* s  _   c-=1e-8;
1 q& j1 E1 q# b1 Y   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
, w4 j' {0 y# m& [      {
9 Q; x* ^- _5 p6 j4 M+ k% }5 h    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
* p1 }& n# D. y0 d& r* B1 R    return(-1);
      }
  m* I8 e# b3 d7 X! y# [' e   if((c&lt;0)&amp;&amp;(l1==1))
+ g8 n/ \" m9 j% c       {
+ Z& S# u+ y  {: @, S     l1=50;
; Q4 |- A( J, s+ H( X0 _7 H% T     goto leap2;
       }
) C) J  x$ C+ M3 L8 D    if(c&lt;=0)
       {
      for(i=1;i&lt;=m;i++)
& x9 i5 q2 C2 `( E/ L, n2 G         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
9 u7 P. }0 R8 s9 d      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
      printf("\nF=%f\n",(*a)[m2][mn1]);
; @9 ^, v: l  x3 H8 b      return 0;
    }
     c=1e8;
4 j8 P9 r# ?, {7 x, H8 e8 C     for(i=1;i&lt;=m;i++)
/ N  F8 k1 u& t" n2 j       {
  ^8 d0 R$ o) Y$ j$ z) H     if((*a)[j0]&gt;1e-8)
& u3 u1 U$ H9 o2 ^. q- o! s- y6 S        {
          g=(*a)[mn1]/(*a)[j0];
          if(g-c&lt;0);
" X5 x  I1 v3 U5 A; b3 h         {
$ ?( G; K8 I- ~0 u9 x) Q           c=g;
9 o* u1 ?4 C( D" `2 A% `           i0=i;
          }
: s8 B! L8 j( ~5 w' J% e1 p         }
    }
+ ^; ?) O# k0 t& _7 |9 J. s; ^      if(c==1e+8)
    {
      printf("\n\t*******Lp no solution********\n");
      return -3;
     }
! T% h1 w" H$ l" Q2 }+ @& u0 u       (*k)[i0]=j0;
( [4 s4 d! [' m; u9 x, V       for(j=1;j&lt;=mn1;j++)
8 F/ o3 M3 d: ^8 p* G3 a9 ^      {
        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
2 P- [1 @. O# E  i0 l1 X           continue;
. g' b( f& K: ~8 y         g=(*a)[i0][j]/(*a)[i0][j0];
3 c# p8 L8 r& H( b( q- W* ?% R, {           (*a)[i0][j]=g;
         for(i=1;i&lt;=m2;i++)
           {
         if((i==i0)||(!(l1==1)&amp;(i==m1)))
3 i2 i+ A$ l5 H% Z4 F2 k  m            continue;
         (*a)[j]=(*a)[j]-(*a)[j0]*g;
1 d3 [2 @  B$ [( I8 p% g           }
       }
7 O/ [7 p6 Z. N  ~8 Z       for(i=1;i&lt;=m2;i++)
; M+ W, `" N# ^0 a      (*a)[j0]=0;
       (*a)[i0][j0]=1;
2 H) t' w- \5 l: P# o4 N' D' A% A       goto leap1;
}  
main()
" T& S$ Y5 ^1 `& V4 G{
  float a[5][7]={{0,0,0, 0,  0,0,0},
8 x7 [. L# D  J                 {0,3,-4,3,  1,0,12},
                 {0,3,0, 6,  0,1,12},
% g6 H2 Q2 j2 _1 m! W                 {0,0,0, 0,  0,0,0},
                 {0,69,0,144,0,0,300}};
  int k[3];
+ }9 A0 t, M9 I- x! L- J, c7 J  float x[3];
  clrscr();
  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
  z1 [7 R  B* C. F' Y, K9 t1 C}
# v' @+ k4 v) A" k" t  U! X. p</P>

plgatc

<>/*************************************************************************
                    单纯型法解线性规划问题（两阶段法）        
8 L1 w  ?2 s' \8 Q) {7 n' I
编程环境：VC++6.0     
方程组输入说明：
变量非负，按提示输入相关参数。
0 T* |  l4 p; A! y( p7 e *************************************************************************/
- z5 B5 B3 k" a& L+ N! b; h#include &lt;stdio.h&gt;
9 |, e- [, W, a, T- W  \1 {: }#include &lt;stdlib.h&gt;
; h4 I2 [1 U& T1 F( C. V#define MAX 100
; ?  [  `& l) @#define STP 100</P>
<>int stop=1; //迭代记数变量
4 @: w" ]0 W0 `8 F+ A9 `' G; k8 eint status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
$ H0 d8 \" b$ S. Aint step=1; //目前阶段</P>
5 {7 R( e% D4 a- N1 f; q<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
int num_x; //变量个数
int num_st; //约束方程数
; q$ Z/ w# B# r! xint num_ar=0; //人工变量个数
# [0 K3 H# D# M# v! J5 xint arti[MAX]; //人工变量下标
, K1 L- U# Y) z5 tint base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
% _9 y: K, _" a$ w, E6 l<>void create(); //建立方程组
void iterative(); //单纯型法迭代
void output(); //输出结果
6 Q1 S" w( a: X# J) G! C5 w& [1 q( mvoid banner(); //打印程序标题
void exchange(); //交换两阶段价值系数
void show(); //输出方程组</P>
<>void main() {
3 b' f4 d! C3 b" [! ~9 v& \ int i,j,k;
/ K# r9 _- i1 e( a0 d/ ^. l banner();
8 Q. N1 g( y3 C3 m, P create();
' ^0 d) X& S0 b/ j  s8 B* w //保存原价值系数，转换为第一阶段价值系数
: A0 k  n) x: j" c' N+ ?8 U) C) ^ 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;
, t% B1 g% l& b' l3 S }
exchange(c,temp_c);</P>
<> printf("\n\n第一阶段问题为：\n\n");
1 y) l- O: e/ C  t show();
+ o: a0 ?4 V- t5 p# f3 O step++;
printf("\n\n按回车开始第一阶段迭代");
getchar();
3 g6 |; q( m4 q/ S# y6 I getchar();
iterative();
if(status==-2) {
  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
  getchar();
$ n' c8 H- d3 R  P! A' S5 l: \  exit(0);
* Q% }( @( R& S  z* R }
output();</P>
<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
" w' X7 ^( l/ U( {  puts("\n按回车结束");
5 N( k$ j& C; f/ T  getchar();
  exit(0);
}</P>
<> //转换为第二阶段价值系数
exchange(c,temp_c);
" q6 Q& m% @, a8 C1 n& j //把人工变量列全设为0
for(i=1;i&lt;=num_ar;i++) {
/ Z% [! z. P0 |& f/ A* W0 `3 q3 I# q  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
+ f8 x7 ]" o; O$ B }</P>
<> puts("\n\n第二阶段问题为：\n\n");
show();
puts("\n\n按回车开始第二阶段迭代");
' O5 J% V& {7 N. S) Q" V2 w getchar();
" z* K, C% ?! {# ?2 s" s+ X iterative();
switch(status) {
case 1:
  output();
  puts("\n\n原问题有唯一最优解。\n");
4 B& a; `/ V5 @1 i  puts("\n按回车结束");
  getchar();
  d, }. {/ O% R" ~1 @" r7 x  exit(0);
case 0:
5 d% F- n) d$ p+ U+ g  puts("\n\n原问题为无界解。\n");
  puts("\n按回车结束");
  getchar();
( I4 g6 `* a  K6 f) C. T  exit(0);
case -1:
  output();
  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
4 r4 B4 S" C8 Z' N) Q5 W  getchar();
  exit(0);
5 o" @: e5 H" q" Q case -2:
. z  x+ r' e6 \6 E# V0 Q2 D  puts("迭代超过限制次数强行终止！\n");
/ A% F! v; w1 ?0 K- T1 W  puts("\n按回车结束");
  getchar();
  exit(0);
}//switch
8 b0 O* y' w9 [: [}</P>
! x' G9 P9 V' b9 s# K7 o<>void banner() {
7 ]3 C+ \& d2 c2 b" C% t  J$ b; n! k$ c printf("\t\t****************************************\n");
% ?2 z7 S# v% S) r' y$ Z printf("\t\t         单纯型法解线性规划问题\n");
printf("\t\t                         作者：Thunder\n");
' g: i+ D+ s( ?) @$ g# s printf("\t\t****************************************\n");
printf("\n");
5 ?; `: O" V3 U; F- D}</P>
/ n9 y% f- C, O: z! p+ N: m; b<>void show() {
1 v2 X( C; L* _1 B//对方程组以自然的格式输出，系数为零的x不显示
$ h: u* x( M5 |. L* X: d! U//为1的不显示系数1，-1系数只显示负号
9 b  r0 Q  j, k& u+ m int i,j,k;
switch(step) {
0 L) I6 y) b4 D; {8 C case 1:
  printf("min z= ");
  printf("x[%d]",arti[1]);
/ t+ M& k' [& }  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
case 2:
- q  V' ]3 l; X+ l printf("max z= ");
4 b# J0 M' u% o printf("%lg x[%d]",c[1],1);
8 K9 f5 k0 [/ B for(i=2;i&lt;=num_x;i++) {
  x" `" B1 t' B. P. i: K  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# E6 z! Y: H/ u( H }
break;
}</P>
" ?/ I! N  M9 ?<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
& @% ^3 H0 S5 |# w+ l: d5 m7 j  k=0;
1 n$ n8 b. g7 `5 u; R0 ~  for(j=1;j&lt;=num_x;j++) {
0 O8 U, e: U! @/ c& ^9 O& c1 e5 \4 I   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);
' f$ V: Z+ F& Z! ]2 {: f' @1 ^    else if(a[j]==-1) printf(" - x[%d]",j);
& Z! y/ p6 S" Z4 q5 f9 o3 R# H    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);
+ M6 d( x3 f  D# |0 C& U    else printf(" %lg x[%d]",a[j],j);
    k=1;
, L! H2 N" p8 e   }
& I+ E$ j( w/ G: ]  u6 T3 @3 k  }
  B& H# |- k2 w: O# U% D2 D4 j2 Q9 t  printf(" == %lg\n",b);
}
0 p6 E" P+ X; f! |1 M( |) ] printf(" x[1]~x[%d]&gt;=0",num_x);
9 [0 _. }: e, i- P4 `. X}</P>
' ^$ A, ^2 h+ v( M' G<>void exchange() {
int i;
double temp[MAX];
for(i=1;i&lt;=num_x;i++) {
  temp=temp_c;
/ X9 W4 W8 _& H+ m2 X  temp_c=c;
, z) Y8 M: {5 W/ b6 c  c=temp;
: `! E; T& O5 [; Z }
9 j9 t! |& H0 A- t9 R' B( ]}</P>
; V6 B" ^) N* Y& K( [  |<>void create() {
//输入方程组系数，每个方程输完后回显确认
int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
" t; `5 N, t, b( z( p$ l char confirm;

# v- |, D& M7 ~- C8 d while(1) {
  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  scanf("%d",&amp;ma_mi);
  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
' Y6 |& [% y; G- @. [$ U  else break;
}

while(1) {
4 ?* K9 M; L6 p6 n$ P  printf("指定变量个数：");
4 G( i- `2 S1 {! M7 x% U# `  scanf("%d",&amp;num_x);
  printf("输入价值系数c1-c%d:\n",num_x);
, z" s3 G  \& y. Y& D  for(i=1;i&lt;=num_x;i++) {   
   printf("c%d=",i);
3 f* F1 b* i3 r8 K0 i( c; j+ r   scanf("%lf",&amp;c);
( {% H6 e/ }% I  g& V, B# U0 b  }
$ n2 {5 T& T% N7 e6 l" o0 M5 f0 Y  if(ma_mi==1) printf("max z= ");
  else printf("min z= ");
* T6 J9 r& w" ^( J) d* k& l  printf("%lg x[%d]",c[1],1);
1 z; W  ?5 D4 p5 _+ Y; x  for(i=2;i&lt;=num_x;i++) {
   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
; |- E, R4 H' U! U/ ]% C   else printf("  %lg x[%d]",c,i);
2 c6 J) B$ l8 a# y" v  }
printf("\n正确吗？：(y/n)");
  u( @" o3 C+ T0 w0 F  F. [ getchar();
, l1 n# u( S' t confirm=getchar();
if (confirm=='y') break;
. p8 G  t2 a: j9 g2 k1 M& F  z else if(confirm=='n') continue;
' c2 v6 X, ?/ O. ~ }</P>
<> printf("输入约束方程组个数：");
( t, Q! a  l8 u, u( K scanf("%d",&amp;num_st);
for(i=1;i&lt;=num_st;i++) {
  printf("st.%d:\n",i);
+ J$ u; W8 ?1 r1 k6 s  while(1) {
   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
   scanf("%d",&amp;re_st);
9 C1 y3 d" d6 W: N   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
   else break;
9 i1 i% E- w* u  }
# u- h  b7 U% l: k1 v: k1 C/ }8 |  printf("输入技术系数：\n");
  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
  }
  printf("输入资源拥有量：\nb%d=",i);
  scanf("%lf",&amp;b);
  
+ y5 J8 d1 Q# g( c/ v7 j4 [  printf("st.%i:\n",i);
3 u6 x8 I; [# |1 F  printf("%lg x[%d]",a[1],1);
, S( Q$ J$ G6 R+ \8 W  for(j=2;j&lt;=num_x;j++) {
   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
7 E0 r% P3 x4 a6 E   else printf(" %lg x[%d]",a[j],j);
4 p2 b) S* A9 K; V# c9 L  }
  i$ M- V) n. M  switch(re_st) {
- f" N/ i( e% X* x6 @7 r   case 1: printf(" == %lg",b); break;
2 b3 E& h" k' Z5 f: Z9 q5 N   case 2: printf(" &gt;= %lg",b); break;
+ X1 ~7 M2 V1 y   case 3: printf(" &lt;= %lg",b); break;
8 u& ?! z6 H- U& ]& k  }</P>
7 A. {5 b& r* p, D2 d<>  while(1) {
   printf("\n正确吗？(y/n)");
2 n- ^5 s/ M! d$ O$ ?   getchar();
   confirm=getchar();
   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
- n+ l2 W) P1 N8 D  }
9 O  n9 \$ n8 t7 f' ^/ v }</P>
<>//显示输入的方程组
8 {9 `* l. m& o0 r) H, m, P0 j printf("\n原问题为：\n\n");
) _8 |( s. b( c" f- U% v if(ma_mi==1) printf("max z= ");
0 l. o/ H6 g' \# B2 O else printf("min z= ");
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);
! [0 k1 F1 r) s% q1 X. d4 ^& h' S  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
8 W1 b1 ~: N+ j2 t# H' Z  else printf(" %lg x[%d]",c,i);
! r9 V- k2 A, k$ _3 s. X* b8 A }</P>
/ \3 y: Z8 v$ Y) A<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  k=0;
: O' t1 I" E7 Q, Q! F0 D  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
  S1 I, @& ~0 {    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
9 `7 r$ t6 i8 g    else if(a[j]==-1) printf(" - x[%d]",j);
; @, C% u4 U7 c( U# x( `1 U2 A    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);
0 C, N% v1 z. \6 ?    k=1;
   }
  }
  switch(re_st) {
   case 1:
    printf(" == %lg\n",b);
/ d. @( u" ^- G5 }    break;
& X' ^2 ~/ t8 J; ?- F   case 2:
    printf(" &gt;= %lg\n",b);
    break;
& y/ o" x2 O7 I5 ~4 r   case 3:
  a) L4 z  r+ \& ~9 N. ~    printf(" &lt;= %lg\n",b);
    break;
9 m, ]4 _5 w9 n3 g  }
, A$ m* E1 g( Q. y# ]  ^ }
printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
/ q% T  m1 I, O0 |( k2 d<> tnum_x=num_x;
/ u5 x& Q/ N' p8 T3 Q' h3 h for(i=1;i&lt;=num_st;i++) {
0 K/ U3 W7 P4 x9 J* F$ V  switch(re_st) {
, g( d3 R9 M# k2 i* Q6 t  case 1:
' G! }: u8 A- D0 l8 V: D  case 3:
   num_x+=1;
4 h4 Y4 r$ F3 g* c* U   break;
  case 2:
   num_x+=2;
   break;
  }
4 f# T- N8 B. S+ P; `7 K. d }</P>
7 M. U) ?0 n7 K; ^+ A<>//化为标准形式
- C  |6 u- M/ F: X6 l9 [7 u2 G6 Y: _* _ if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
for(i=1;i&lt;=num_st;i++) {
: o6 `% }, v0 L4 _* [9 o- b% l& h  switch(re_st) {
2 ]  M' [. d" W; i/ t* y% K8 S   case 1:
  D+ E8 O" A9 q3 I* G    num_addv++;
- R1 c  y+ F) i& X, p    num_ba++;
% F0 H2 m8 H* F8 k2 q) g0 u6 D    num_ar++;
    c[tnum_x+num_addv]=0;
4 h$ y8 V: s' j* A    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) a[tnum_x+num_addv]=1;
. r1 W! e; E% L& ]: B     else a[j]=0;
" V, q+ b0 F9 C6 Q: F    break;
   case 2:
    num_addv++;
) Y- J! e* K, V5 j2 r8 i$ J4 o7 G3 ]    c[tnum_x+num_addv]=0;
/ ^* W6 {2 e. |" B* z) r3 Y6 j- Q; I    num_addv++;
    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;
% D/ R! V8 D8 j6 M8 P+ ~, b; \2 A* q) Y     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
    break;
   case 3:
7 w1 e. `6 x5 N7 i# c% e    num_addv++;
: F. Q$ A+ w$ u/ G& z, Y    num_ba++;
    c[tnum_x+num_addv]=0;
    base[num_ba]=tnum_x+num_addv;
& u' v% ]% {6 t, n1 c! v2 b( q    for(j=tnum_x+1;j&lt;=num_x;j++)
' ]+ k7 c0 |( K9 a3 ~% Y' C! w     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
    break;
  }//switch
. f4 o2 Q; O# \2 Y7 d1 ~ }//增加松弛变量、剩余变量、人工变量、确定基变量</P>
<>//显示标准化后的方程组
0 s; j5 j4 }% S4 T9 Q# j* k1 r) C" o printf("\n化为标准形式后：\n\n");
& E" I. O& ^% ^4 t. o9 i if(ma_mi==1) printf("max z= ");
; E% d: _* B" Y1 k else printf("max z'= ");
printf("%lg x[%d]",c[1],1);
/ X( s! R/ {8 z7 |# w$ i for(i=2;i&lt;=num_x;i++) {
  k=0;
4 ^/ U3 [3 q$ Z) J  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);
9 T  j* D. Z% r+ o+ l  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>
) o+ x& K1 P7 q5 q  u+ Q<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  k=0;
8 W5 \# t: C* e+ ]5 s  l  for(j=1;j&lt;=num_x;j++) {
% L- v' F' |; C   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
; Q3 l' X2 S+ s5 a# S3 M    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
' b% Y) ~4 e& x) O: l% ~, r    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 |% x$ l0 F0 H7 _! 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);
- l/ `- L0 K$ B& M6 Z& E# M, N! p    k=1;
/ [7 h% q( j) K   }
  }
2 c. M1 L* D$ i# w  printf(" == %lg\n",b);
1 ^8 M5 z$ ]/ W0 } }
3 Z( u. {) W2 P  l) ~5 Y printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
! A4 N/ ?  U% l) C5 }4 v<>void iterative() {
int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
0 O- B( d4 d2 a int base_elem;
int base_out,base_in;
double sigma[MAX],temp;
! B& ?# C- C+ k( k5 W double value_be; //高斯消元里保存主元素值</P>
3 Z0 g9 f8 A( Y  j5 d# K<> printf("\n\n第%d次迭代：\n\n",stop);
( V9 w3 m1 q6 e3 y4 N% c: n for(i=1;i&lt;=num_st;i++) {
  printf("c%d=%lg\t",base,c[base]);
1 G% v: q3 v4 ^/ ^" T( c  printf("b%d=%lg\t",i,b);</P>
<>  switch(step) {
   case 1:
    for(j=1;j&lt;=num_x;j++)
    {
$ F, U. `+ F8 ]     printf("a[%d][%d]=%lg\t",i,j,a[j]);
    }
    printf("\n");
    break;
( D& W/ Q. I5 W9 [0 {; d( r   case 2:
, |# D4 r9 k: \    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;
     if(k_a!=1) printf("a[%d][%d]=%lg\t",i,j,a[j]);
" |) M1 L) E: E# J6 a+ ~    }
; J3 t) C& T$ A! K- n" D( b; O    printf("\n");
    break;
  }
}
//求检验数sigma
% e! _$ W8 m% @. O for(i=1;i&lt;=num_x;i++) {
6 }8 M- Z# r1 Q, d* Z4 @; v% T  sigma=c;
  for(j=1;j&lt;=num_st;j++) sigma-=c[base[j]]*a[j];
, Z9 F. P" q/ I* j. w  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
3 @- j1 A+ z6 x  S# K' B  switch(step) {
' p$ z2 {7 K7 |7 f   case 1:
    printf("sigma[%d]=%lg\t",i,sigma);
& q/ l* L' o. }' @2 S" r6 G* U    break;
/ f. Q) d" c0 b- ~; F0 N; }   case 2:
+ V7 q7 j; F6 q. \6 W8 y( g    k_a=0;
# B6 g% N1 Z; E( a) E5 [: z# t    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
! B- G1 n* a& D/ p, ]* U    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
. o* w& K, y/ i# b8 l1 T+ V: L    break;
  }
' Z4 X% U/ l- A  R }
putchar('\n');
//检验检验数sigma是否全小于等于0
k=0;
for(i=1;i&lt;=num_x;i++) {
- K0 ~* q0 ^) U* f$ F; W  if(sigma&gt;0)
- T, }; O2 i* H& V  k=1;
! \( s/ y+ z0 z5 X+ ?0 q. j }
" W7 l" E( c5 }9 z9 k. x: ]+ X' p* c* J if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
7 k  j- |; B0 `  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;
. Y* f) B  L7 b   if(k_f==0) {status=-1; return;}
   }
% h- V! d6 ?' `  }
; M" j& L0 `& [! k% c' p" Q$ h4 y  status=1;
  return;
4 Z4 I6 B: b* U/ \7 U' @ }
3 L# P" H7 }$ K5 R//检查是否为无界解
" b( ~, q( J2 J) E+ B3 i9 k% d for(i=1;i&lt;=num_x;i++) {
  k_f=0;
  if(sigma&gt;0) {
) H/ D- U$ E: q   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
) T5 A' I9 y* C4 O- ^   if(k_f!=1) {status=0; return;}
) o9 x, J6 p. n1 T  }
1 y) T. z7 x; p# M }</P>
<>//确定换入变量
for(i=1;i&lt;=num_x;i++) {
  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;
: s: I" F# J7 C0 x }//temp赋初值
2 y& ]. q8 V- L4 Q) n for(i=1;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  if(k==0)
  a: P) l! \9 g" R+ Y# a; v8 ]   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
/ g& r* Q% _$ p    base_in=i;
3 s( F( h9 y2 [" z7 _    temp=sigma;
   }
. g- c2 R( @6 b" k. j8 E1 u }</P>
7 R7 j9 Q* a% j/ ^<>//确定换出变量
for(i=1;i&lt;=num_st;i++)
; @( m  q0 c% b  if(a[base_in]&gt;0) {
* F; v, p; S# v9 A  O. c* r   temp=b/a[base_in]+1;
8 s" r, p& C3 t* e6 O$ r   break;
  }//temp赋初值
$ k7 ^3 t0 O% A; ]! `8 U0 Z% a8 T for(i=1;i&lt;=num_st;i++) {
  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
$ b% c* K+ q3 R& V/ t   for(j=1;j&lt;=num_ar;j++)
    if(base==arti[j]) {
     base_out=base;
+ s- q% {/ e- d  E6 Q, P; [; `& \     base_elem=i;
$ P% i: G/ A4 Z7 Z  z- R     temp=b/a[base_in];
7 L8 D( ]' [3 q2 ~3 `, r) D' Y     break;
& L2 F2 |/ p& S4 B2 q8 [9 A1 W    }
  }//人工变量优先换出
$ v! U: c/ H: A  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
   base_out=base;
" K* }! _& |5 j2 d6 `& U   base_elem=i;
   temp=b/a[base_in];
  }
}</P>
% F& q' r2 W' X9 e8 c  {<> printf(" 基变量：");
2 B% S1 _+ m8 f* U" c3 i for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
: r% O4 \" _: G3 [! S8 F/ R# r printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
+ m4 N+ A: d9 r, s# z) z$ w# Y1 f! v//基变量变换，进行新方程初始化后迭代
for(i=1;i&lt;=num_st;i++) {
8 P( G4 i; g' R0 A: o/ p  if(base==base_out) base=base_in;
( {7 t  j: N  |2 l" D% d+ D/ r5 Z; X }
//初始化主元素行系数
8 g$ h& R' W8 N& N$ n; x3 w value_be=a[base_elem][base_in];
5 S$ ?0 M1 p) y8 _ b[base_elem]/=value_be;
7 s' o% u. f' y" _5 v for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
<> for(i=1;i&lt;=num_st;i++) {
) `* c" P4 H0 q1 h, o  V0 g  if(i!=base_elem) {
0 H' ^4 M0 k( e- U, f) {' Z   b-=b[base_elem]*a[base_in];
! J; t% l" L- C) d8 w; u0 K   value_be=a[base_in];
( t4 d9 b9 `5 E8 d7 R/ D1 g7 K0 W   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
. |! {5 W1 I! E2 v! D& T/ |  }
. w( h' K. _' y }
stop++;
if(stop&gt;STP) {status=-2; return;}
iterative();
! A9 Z. X2 M2 g8 P; o}</P>
1 z+ y) N. h: d' i7 A; K5 K6 i<>void output() {
int i,j;
" T; w8 g! ^& [! T double X[MAX];
+ \& ]; Z% _) C: S( ]1 A% N4 m printf("\n结果如下：\n");
printf("\nX=(");
3 q- f7 Z) t9 O) N, v for(i=1;i&lt;=num_x;i++) {
# j' q  Z$ R6 _7 M7 a$ v, C1 N  for(j=1;j&lt;=num_st;j++)
, j% d; ]+ L% I( X9 m' ^  if(i==base[j]) {X=b[j];break;}
  else X=0;
  printf("%lg ",X);
}
' @8 E6 B( I3 F- j6 \; J printf(")");
. j! d4 H. S% \+ @4 y! F& m for(i=1;i&lt;=num_x;i++) max+=c*X;
# f; E) Y. m1 P4 y if(ma_mi==1) printf("\nMax z= %lf\n",max);
$ n2 k  x- H6 p+ t1 @ else printf("\nMin z= %lf\n",-max);
! p( {. }7 n8 o1 O; b  P2 G}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了