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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
) A) q. S8 r4 i- [#include&lt;math.h&gt;
) j# Z: ^/ V, x5 U#define X 5
) c: K' G$ A2 C6 u/ B$ P#define Y 7
; r6 |  a2 C0 d) m+ U+ Lvoid xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
{
  int j;
  *c=0;
  for(j=1;j&lt;=*is;j++)
    {
      if((*a)[*ir][j]-*c&gt;0)
9 R+ B: }' ~: U  O( T/ D$ w0 w    {
      *c=(*a)[*ir][j];
$ r$ u* f: u, g. d8 p! E      *j0=j;
4 U) s" T6 k; @" l* v! V    }
7 x: K2 @: F2 @" K) \% N' u( \     }
}</P><>      /***************** 参数说明 **********************/
( ~! S: C8 q7 _$ f: q8 r      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
+ _. g4 w7 P( A( Q3 i      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
; K7 `. [; g+ j9 |  H5 M7 U      /********* <a href="http://happyyangxu.home.sunbo.net*******/" target="_blank" >http://happyyangxu.home.sunbo.net*******/</A>
; v- `* z: S" D* y% P2 Xint xi_sm(int m,int n,int m2,int mn1,int l1,float (*a)[X][Y],
2 a, P6 M" @2 L8 z) O2 }0 _8 @       int(*k)[],float(*x)[])
& H, }: X2 P6 y) T2 r{
  int i,m1,mn,j0,i0,j;
  float c,g;
  m1=m+1;
3 w7 _: ~% Y1 X1 o$ L& b  mn=m+n;
/ h# X( R. G: s, w  for(i=1;i&lt;=m;i++)
* W4 f) l7 ?9 G, i, e! ?1 i    (*k)=n+i;
) j# _6 o! `+ k5 Nleap1:
  if(l1-1==0)
% b5 ]( ]; W, _! N) R4 g( L# F     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  else
8 W& ?* k2 f8 x1 n: B     {
% v7 ?  z7 B: |9 {0 y+ Zleap2:   if(l1-50==0)
" k1 @! n' @- P# `  q0 \% p      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
) d+ I$ W% D7 z% u       else
      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
     }
! m) {2 k# f9 s1 z7 s   c-=1e-8;
+ W$ t& e$ c  Z4 a: [; l+ A: M   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
      {
' c5 G" {5 _9 z* J    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
9 _9 z  ]: T# p4 i; w    return(-1);
* z# M  m$ ?9 ]& _% c      }
   if((c&lt;0)&amp;&amp;(l1==1))
       {
     l1=50;
: _7 H3 v; @: V     goto leap2;
       }
    if(c&lt;=0)
5 n& G3 f2 J- S  }       {
9 u: v8 }% x' b! s* J      for(i=1;i&lt;=m;i++)
         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
/ r' _$ z  O- z      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
2 p  H/ E7 l: W! I      printf("\nF=%f\n",(*a)[m2][mn1]);
1 y& d! a% V3 Q6 x      return 0;
    }
     c=1e8;
; L: D7 S- n. P% J5 s     for(i=1;i&lt;=m;i++)
: ^" c! M* c' b- e& k, K* {       {
     if((*a)[j0]&gt;1e-8)
        {
; y( H+ G) ]( h% }          g=(*a)[mn1]/(*a)[j0];
) t$ B# M* s! D# @' i# J8 z! g          if(g-c&lt;0);
         {
4 R1 q, C3 y  U  {+ G7 j           c=g;
           i0=i;
7 g) I, v6 a: P& Z# t          }
# v/ ]3 P. m, d) b" r         }
! b2 ]( O  ]+ I3 ^% u4 J    }
      if(c==1e+8)
. A" m. }% u2 Q* T    {
6 ]) q! c7 ^: {) u  ^& X* K+ C; f      printf("\n\t*******Lp no solution********\n");
      return -3;
     }
' A! G9 s" W  b, K       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
9 m6 M) g) T) J4 ]- F& f      {
" i+ e2 Q' w5 V1 ?        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
0 x# M0 Q. H  D( {           continue;
         g=(*a)[i0][j]/(*a)[i0][j0];
; F$ T- v$ U$ _+ s9 G' f) i           (*a)[i0][j]=g;
- g, N" v% d% l9 E         for(i=1;i&lt;=m2;i++)
! P+ D9 d# E* ]1 W0 {           {
: \9 \7 {% b, ~, M) E2 L) G0 q4 |         if((i==i0)||(!(l1==1)&amp;(i==m1)))
  G- t5 X" m" z" T% {            continue;
+ e* T% }% C6 m3 u3 G; E         (*a)[j]=(*a)[j]-(*a)[j0]*g;
. q* m! S; g4 j  S- B$ V$ T9 r           }
. n- W( i; h2 G9 r$ M" \( Q& L       }
# b, N- z0 m+ W- x       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
+ s$ z; {, b5 Z& z1 W" J       (*a)[i0][j0]=1;
' }3 D: }' a7 p5 c9 m       goto leap1;
}  
main()
' p5 L3 c) f2 w* j; ]1 F2 s{
  float a[5][7]={{0,0,0, 0,  0,0,0},
                 {0,3,-4,3,  1,0,12},
; A" Q" ^! j% U, p  A1 {1 p1 }                 {0,3,0, 6,  0,1,12},
0 s6 B1 O5 m7 K3 T" R                 {0,0,0, 0,  0,0,0},
                 {0,69,0,144,0,0,300}};
4 K8 Y, J* z' i+ O) F+ R  int k[3];
- b1 \' v6 `: S- Z  ^" E2 L4 r# G  float x[3];
  g, P" H- u1 n  clrscr();
6 I8 Y/ v8 {. l  I0 O" u" `7 {  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
}
- q% M8 }6 s% o; l1 u; Q& _9 [</P>

plgatc

<>/*************************************************************************
) V. ]4 m6 a8 ^7 N" L1 F                    单纯型法解线性规划问题（两阶段法）        
2 @' O# o3 f+ r" i" B
5 i& R9 }1 Z) f' x 编程环境：VC++6.0     
方程组输入说明：
变量非负，按提示输入相关参数。
, v& B- M1 y& C, ]* ]  R0 K *************************************************************************/
#include &lt;stdio.h&gt;
#include &lt;stdlib.h&gt;
#define MAX 100
+ L7 `9 Z. ~/ w+ @- o#define STP 100</P>
+ O. q& E; R0 ]' e<>int stop=1; //迭代记数变量
/ Z+ J6 N( V) [2 lint status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
int step=1; //目前阶段</P>
$ f/ S. y# y+ d$ g6 _<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
0 a) Q& ~3 }) f- c7 U) f, |int num_x; //变量个数
! Y) Y% l- }0 ]8 Mint num_st; //约束方程数
int num_ar=0; //人工变量个数
3 u* @9 c) M' t/ F* e# g* jint arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
6 W( a$ N/ z* U8 [: U0 Tint ma_mi; //1为求最大值，2为求最小值</P>
<>void create(); //建立方程组
void iterative(); //单纯型法迭代
) G# x; a: g% w# ~4 `% i% }% Vvoid output(); //输出结果
9 L9 t# y# }! s+ e8 F( ?void banner(); //打印程序标题
/ R2 Y. y6 u0 T' Y7 W9 c. s6 qvoid exchange(); //交换两阶段价值系数
void show(); //输出方程组</P>
<>void main() {
int i,j,k;
banner();
* y* b: Z7 \7 g2 |. e+ f create();
//保存原价值系数，转换为第一阶段价值系数
! i' H' ]% {' w for(i=1;i&lt;=num_x;i++) {
6 @9 @7 s/ F1 Q& r! Z  k=0;
' o, `) @. M* K3 D5 P4 V: e9 f  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
  if(k==1) temp_c=-1;
$ K3 y% T$ s3 G* z, l  else temp_c=0;
}
  h. H# `! H& ?1 _ exchange(c,temp_c);</P>
<> printf("\n\n第一阶段问题为：\n\n");
; O# t- F4 a2 o5 C, j4 R2 g7 { show();
9 }# N7 }0 L% P' ^+ `+ y step++;
printf("\n\n按回车开始第一阶段迭代");
3 k1 \2 T1 N  N2 N  s getchar();
5 M  _  M6 e* E0 ~% N# o getchar();
; F8 e9 y, T- u+ f* `. R iterative();
9 f$ b0 C4 Q. N8 W" J if(status==-2) {
  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
+ G: ~7 X* b, H. G  getchar();
+ S" s4 ~% v2 p/ U8 i* K  exit(0);
}
8 _8 @$ ^5 S3 c' [/ D- a output();</P>
8 J: q/ b- M  n  I$ T  n<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
  puts("\n按回车结束");
1 h9 P% h7 b. n  s( E7 e  getchar();
" n/ `( D0 X. y2 h, }* _/ `. q  exit(0);
}</P>
<> //转换为第二阶段价值系数
" Y6 @/ W, _# l" q1 |8 O" ? exchange(c,temp_c);
//把人工变量列全设为0
( H9 S$ j" W% B  x$ o for(i=1;i&lt;=num_ar;i++) {
  c[arti]=0;
) \* T; h1 l" C- ^  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
3 x$ _. n; {( G }</P>
; N! ~- o) f/ W. d<> puts("\n\n第二阶段问题为：\n\n");
$ z( k5 Y6 i; I, \% |  p show();
, z8 q: k- K: }- j/ W puts("\n\n按回车开始第二阶段迭代");
& U+ c3 V( o5 K6 S' x# V getchar();
iterative();
. z" L0 T) w6 s( u switch(status) {
case 1:
! E8 ?) B2 P' L/ |. Y  output();
  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
( J  s6 q) y, f4 |3 M: k! e; R  getchar();
  exit(0);
% m3 @: P$ l2 x" r( w3 q case 0:
8 b, b* e' B7 @  puts("\n\n原问题为无界解。\n");
6 f/ H  N- z9 R- x: Y( Z  puts("\n按回车结束");
  getchar();
  exit(0);
% I- D, h& z% ^& _) o+ B case -1:
  A# B' w' v% k' M+ w  output();
  puts("\n\n原问题有无穷多最优解。\n");
, U( u) g- x/ x7 j' s  puts("\n按回车结束");
  getchar();
; f& C* p% u& V  D3 Z  exit(0);
case -2:
, j' X" O& v0 ^. W+ Y  puts("迭代超过限制次数强行终止！\n");
8 d4 f) [# T8 G* }0 n& O  puts("\n按回车结束");
" U) \, s# ]5 f% y  getchar();
  exit(0);
( T; Z: ?9 d% a }//switch
}</P>
9 m3 l9 P" s( J& e<>void banner() {
printf("\t\t****************************************\n");
9 B3 u  B; A% Z2 n; a printf("\t\t         单纯型法解线性规划问题\n");
printf("\t\t                         作者：Thunder\n");
printf("\t\t****************************************\n");
/ w8 e  y6 e+ P) J, h printf("\n");
}</P>
/ u' k4 P  x/ V9 w2 n" M+ `* b<>void show() {
& o* a/ g9 }9 y//对方程组以自然的格式输出，系数为零的x不显示
$ Y3 `! \5 l6 S+ t  }3 ~5 O+ \//为1的不显示系数1，-1系数只显示负号
/ n3 ^* s- J, U# m% w0 ?3 e int i,j,k;
/ C0 y: z$ Y7 ~3 X- M5 Z- V switch(step) {
case 1:
% S$ Z' J2 F4 u  printf("min z= ");
  printf("x[%d]",arti[1]);
  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
case 2:
printf("max z= ");
printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
4 C# X' y; u, u: N3 {2 r  if(c==1) printf(" + x[%d]",i);
; Y1 g7 ]- c7 u6 s- h  else if(c==-1) printf(" - x[%d]",i);
( l+ b: O& h7 L2 M) a0 `2 {  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
; g  }. e: K3 K  n& v# W2 A9 ]  else printf(" %lg x[%d]",c,i);
5 G3 }& R5 Y) T6 }( B6 N2 j }
break;
( _) U# I5 a4 L( A }</P>
<> printf("\nst:\n");
+ r& R1 i% R6 d& p2 O+ z4 i for(i=1;i&lt;=num_st;i++) {
2 S" w" S$ i; X  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);
    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);
    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
; D  x" v7 B! r    else printf(" %lg x[%d]",a[j],j);
    k=1;
   }
9 o( d- p. ^8 J7 ~6 T- U: c, i5 ?( f  }
% }( x+ G" S7 p1 Z0 ]/ ^* G  printf(" == %lg\n",b);
}
printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
% `- D$ ~8 ~* n5 ?! f* |3 A3 g<>void exchange() {
3 l/ {( C# s* g$ j  f int i;
& Q- S5 J4 m6 L0 s. |4 J4 y* j double temp[MAX];
for(i=1;i&lt;=num_x;i++) {
3 k  j; o5 A; w( ^, K8 m  temp=temp_c;
  temp_c=c;
  c=temp;
( P0 [. m$ P: L }
}</P>
! v2 w- T% Q% m4 J  r; S<>void create() {
//输入方程组系数，每个方程输完后回显确认
int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
char confirm;

while(1) {
- x7 n/ c0 ~  e  ?: k  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  scanf("%d",&amp;ma_mi);
% o* e& V# c0 Z: z, a8 B  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
/ t( R1 k0 a/ `  else break;
}

2 ?% H$ w  N, e! Q  ^ while(1) {
  printf("指定变量个数：");
% x3 C2 ^1 l, p- T* l  scanf("%d",&amp;num_x);
. E% d( U* L) @: F  K- `! ]  ?3 O  printf("输入价值系数c1-c%d:\n",num_x);
+ Y" y5 p3 T2 O3 Y& a8 r1 T* n- h  for(i=1;i&lt;=num_x;i++) {   
/ ~, f) k7 ~6 g( o   printf("c%d=",i);
5 }0 Z5 c3 _8 Y" L6 E1 z% H   scanf("%lf",&amp;c);
  }
4 E) `) r% C6 E+ P  if(ma_mi==1) printf("max z= ");
# ~, |3 r/ y) U% K1 F) K1 m3 ~  else printf("min z= ");
. E4 A: N6 {" a' b+ P, H: _  printf("%lg x[%d]",c[1],1);
; z9 ?: V0 z1 I+ E% U. j6 g, o  for(i=2;i&lt;=num_x;i++) {
   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
5 ~6 K: k. ?4 y- @) t2 Y9 p( z   else printf("  %lg x[%d]",c,i);
  }
printf("\n正确吗？：(y/n)");
3 r# p* o3 m( c. ?' C* W" Z* B getchar();
4 K+ J$ v) Q; V+ i1 h confirm=getchar();
; [+ q7 ]/ ^/ G6 h. N& x! d5 \ if (confirm=='y') break;
: B" [, W8 {& C+ y9 H5 v/ N else if(confirm=='n') continue;
}</P>
0 Z: B0 a4 r2 v# B' V9 @) ^: p# M7 T<> printf("输入约束方程组个数：");
scanf("%d",&amp;num_st);
for(i=1;i&lt;=num_st;i++) {
+ R4 R; b& p# u+ y# |) j  printf("st.%d:\n",i);
  while(1) {
" Z' n) U$ Y, p: m: C0 L: `   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
: n/ |( T2 i/ w8 M7 B   scanf("%d",&amp;re_st);
   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
   else break;
  }
2 ]5 _5 H, ]5 h; \% o5 h* k  printf("输入技术系数：\n");
  for(j=1;j&lt;=num_x;j++) {  
5 d& W9 V6 t/ g1 o% p   printf("a%d=",j);
0 c5 _" H4 ]8 r; ?; l   scanf("%lf",&amp;a[j]);
  }
1 c: N( [! Z1 p1 Y% k" h: q' m1 [7 u  printf("输入资源拥有量：\nb%d=",i);
  scanf("%lf",&amp;b);
  
  printf("st.%i:\n",i);
  printf("%lg x[%d]",a[1],1);
) E- Q' t+ B; T' ~+ Y3 `  for(j=2;j&lt;=num_x;j++) {
   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
- Q1 B, y9 @4 J1 E5 M) f" O   else printf(" %lg x[%d]",a[j],j);
6 [7 X2 C" h# \- H* g! l+ A7 L  }
  switch(re_st) {
   case 1: printf(" == %lg",b); break;
3 A6 l+ N' j  E. y1 t3 L; s   case 2: printf(" &gt;= %lg",b); break;
   case 3: printf(" &lt;= %lg",b); break;
6 o% g  k0 F( {3 n, z  }</P>
, D5 R' w) D' S  k  W* y$ D/ x<>  while(1) {
& G' J* Q" {5 J( {   printf("\n正确吗？(y/n)");
$ [! z/ n+ P4 y% u   getchar();
( i; O  I5 ~: w) F$ K% ?/ U   confirm=getchar();
, O3 w5 f9 Z& z   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
+ q, H* [* v' q3 y) L- l  }
& M% m# C, i5 s$ x: ^- l }</P>
<>//显示输入的方程组
8 ?/ J( P% ~6 P9 V printf("\n原问题为：\n\n");
if(ma_mi==1) printf("max z= ");
else printf("min z= ");
printf("%lg x[%d]",c[1],1);
$ t8 u8 \6 g7 D' u. K7 H for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
2 x1 ?# j% i% W9 `  U  else if(c==-1) printf(" - x[%d]",i);
8 ]( K' W) p! d5 ^  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
  else printf(" %lg x[%d]",c,i);
}</P>
/ k% f* j$ W; V. d. q<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
( b. J. L% e' [) `  i  k=0;
  for(j=1;j&lt;=num_x;j++) {
0 k: [0 k4 ]; _- i1 s   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);
    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);
    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;
   }
% c4 x" W. Z, F% m  }
  switch(re_st) {
   case 1:
    printf(" == %lg\n",b);
    break;
/ {8 w; I/ k6 g# a/ }   case 2:
    printf(" &gt;= %lg\n",b);
    break;
8 S8 E/ G% v' r; o   case 3:
* r" L5 m$ g4 @    printf(" &lt;= %lg\n",b);
5 a7 w. H) ]: n8 K    break;
  }
}
& u9 Z0 n' u$ O3 i printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
<> tnum_x=num_x;
for(i=1;i&lt;=num_st;i++) {
7 I3 k; c; S% ^. t; }1 s9 I  switch(re_st) {
  case 1:
  case 3:
   num_x+=1;
   break;
  case 2:
   num_x+=2;
- s& N3 i' l, j- l   break;
  }
8 B) O% `  K1 `3 x, \1 ^ }</P>
! a4 T& ^; N9 d' Z8 A* U<>//化为标准形式
, i$ W' A* M; ]8 Z3 h* w" J if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
9 \) J5 @, `- l$ o" `   case 1:
$ D5 A- ?  z6 o# }. }) X* c/ I    num_addv++;
    num_ba++;
  e  j4 f! L- N& ~0 g) D    num_ar++;
6 o) @; [# ~, ^' I/ I2 f    c[tnum_x+num_addv]=0;
- _. M7 s5 _- _  L. B4 E) @% l    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
  @5 c' R, g) d8 `8 `4 Z2 \    for(j=tnum_x+1;j&lt;=num_x;j++)
$ A- W7 |; D7 R, S9 S8 U, X8 }     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
    break;
4 H+ I+ ]; a6 j0 \5 N( U6 B   case 2:
( t3 I9 `7 w8 l5 r$ T    num_addv++;
) T1 w2 A, m* i* d  \    c[tnum_x+num_addv]=0;
  E$ G! J" P& |7 o1 h    num_addv++;
    num_ba++;
    num_ar++;
' X/ q8 e: Y8 ]% @: H$ E    c[tnum_x+num_addv]=0;
    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
2 ~- R) V3 |7 o& f5 Y" ~, H    for(j=tnum_x+1;j&lt;=num_x;j++)
, v. O" O4 S4 ~3 `% r' p: {     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;
) P1 {8 [1 r5 l6 p" w. L     else a[j]=0;
    break;
& ~4 q9 c  I# A1 p) f   case 3:
4 p9 g' \7 n7 N1 Q- U* F7 e2 @  f    num_addv++;
1 N, M; L( X# D2 J6 {2 |    num_ba++;
    c[tnum_x+num_addv]=0;
& {( y& T1 j( s9 Z$ v    base[num_ba]=tnum_x+num_addv;
    for(j=tnum_x+1;j&lt;=num_x;j++)
  E+ W9 i. e; k5 P  n7 t4 h     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
" b% Z. {0 m: u     else a[j]=0;
7 h) l' X+ ^8 P% Y$ v6 A8 w    break;
9 n; [' `+ E  ~' W: ?& ?  s  }//switch
}//增加松弛变量、剩余变量、人工变量、确定基变量</P>
<>//显示标准化后的方程组
printf("\n化为标准形式后：\n\n");
7 M" y. U5 B0 T5 O& C5 s& S5 V if(ma_mi==1) printf("max z= ");
& r+ J1 S: ]: { else printf("max z'= ");
' X8 W0 ~+ T2 Q6 [ printf("%lg x[%d]",c[1],1);
3 k8 Y2 L: [, C: u- B7 \ for(i=2;i&lt;=num_x;i++) {
) f  K  ^/ z6 `- V) _5 a  k=0;
& o' J4 v- T, c9 [' M  for(j=1;j&lt;=num_ar;j++)
# ~; U' V1 c: z  if(i==arti[j]) k=1;
  if(k==1) printf(" -M x[%d]",i);
/ ^5 V; u1 I1 Q' Q2 @( `  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);
! ?9 z) Z: X) n  else printf(" %lg x[%d]",c,i);
}</P>
<> printf("\nst:\n");
, p$ e7 o; @; x6 n for(i=1;i&lt;=num_st;i++) {
* W6 w$ g0 U' P, c) U. [* u. J  k=0;
( J' k( k, }, p% S' X/ K! X  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
# H: ^5 V  @2 D; Z    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
2 v% r, j5 R/ |4 e6 C' Q    else if(a[j]==-1) printf(" - x[%d]",j);
. p7 m5 a( I2 X! O. o5 y2 l    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
4 D4 t( O8 L; v- P    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);
. N/ |& }9 Z% O    k=1;
  K2 m- p% S, B- ]   }
  }
  printf(" == %lg\n",b);
}
1 B) Z* B* d7 d) n printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
# Y7 D8 C- a8 |6 i5 V% S! {<>void iterative() {
int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
* Q% A* y9 B$ Z6 u: z# B int base_elem;
+ w' U8 G& C, L% g% H$ b int base_out,base_in;
double sigma[MAX],temp;
% t' u. T" B% h: i# ^- t double value_be; //高斯消元里保存主元素值</P>
<> printf("\n\n第%d次迭代：\n\n",stop);
for(i=1;i&lt;=num_st;i++) {
  printf("c%d=%lg\t",base,c[base]);
  printf("b%d=%lg\t",i,b);</P>
<>  switch(step) {
2 p1 P0 g1 B% O. c   case 1:
: t- [# Z% o" P/ h6 D7 ~) ?    for(j=1;j&lt;=num_x;j++)
4 {4 s& o$ a3 L0 o1 U    {
     printf("a[%d][%d]=%lg\t",i,j,a[j]);
0 i/ ~3 c! [& Y7 m# z1 L% J    }
5 N* M- B6 a  K- ?    printf("\n");
6 s, M# G* r. A' {% b    break;
# `/ {  ]6 X- v. B/ ]% e' @   case 2:
    for(j=1;j&lt;=num_x;j++) {
4 S9 O" V* U5 @     k_a=0;
     for(l=1;l&lt;=num_ar;l++) if(j==arti[l])k_a=1;
3 V! Y! I" W3 _* e     if(k_a!=1) printf("a[%d][%d]=%lg\t",i,j,a[j]);
$ [. U" T5 R0 x! x3 o5 |    }
    printf("\n");
  Y$ f" D" v- F: n- @: Q    break;
/ v! \% e8 y2 N# v: D' ], j  }
* E6 ?" F: a  j% p  e8 [ }
//求检验数sigma
0 y0 Y% q% A1 h+ P7 ^6 ]2 J0 e7 L for(i=1;i&lt;=num_x;i++) {
  sigma=c;
1 |6 W3 `- ^: l$ W2 ?+ ]  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;
6 p; m' Q, z4 D/ [' j2 B  switch(step) {
   case 1:
    printf("sigma[%d]=%lg\t",i,sigma);
3 ^9 I; V% N# Z) [    break;
+ S2 C% p; l" |' u# R! F& X   case 2:
' R7 x: n# a+ O2 \; |    k_a=0;
    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
8 K9 f* S1 {8 `% r& Q2 [, q    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
    break;
  }
6 L' I1 ^! T  Y, y2 g, C( } }
putchar('\n');
; B5 {% I: t0 a' I. O//检验检验数sigma是否全小于等于0
k=0;
for(i=1;i&lt;=num_x;i++) {
% \# u5 S+ ?+ t8 H$ `2 j% v  if(sigma&gt;0)
$ P4 H% T/ O& {/ ^$ b' M  k=1;
}
  H, U- R& O. @# o& D* r if(k==0) {
9 }" R$ l. b* B7 \# X% I3 C; p  //sigma是全小于等于0时，检查是否为无穷多最优解
, D. H9 s) v5 a& r8 U  for(i=1;i&lt;=num_x;i++) {
   k_f=k_a=0;
   for(j=1;j&lt;=num_ar;j++)
5 U6 Z/ V- l, a% s3 i- O  @+ T   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;}
- s& T( s+ Q, g: T5 M: H. W' \   }
  }
  status=1;
  return;
}
$ S2 \+ e  w; X8 a" o4 O0 @: C//检查是否为无界解
& d- {. _* D+ P1 d for(i=1;i&lt;=num_x;i++) {
, B* ?! G" }6 p9 ~8 M8 \0 C  k_f=0;
/ f1 f+ p! `+ q' ~7 `6 M  if(sigma&gt;0) {
4 ?- s8 Q" U/ o' T- y5 `   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
   if(k_f!=1) {status=0; return;}
+ `3 `4 `9 a3 U0 |0 B, s  }
}</P>
* U: ~6 `8 @0 n. r( k, h<>//确定换入变量
for(i=1;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
9 x& \" H# f! D1 ]- h/ Y2 h6 Y+ [1 g) D  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
1 F6 p  f9 U0 u8 Z& N2 b }//temp赋初值
! m1 Q6 M/ y7 ?6 t for(i=1;i&lt;=num_x;i++) {
$ O, G2 p: g# t# H4 |! b  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  if(k==0)
  D2 s4 N5 q9 o3 h; f+ M0 O( b   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
* u3 g, P; G  Z7 c# l7 W4 W    base_in=i;
    temp=sigma;
   }
5 A; d- `3 Y/ P, e }</P>
<>//确定换出变量
( c' c  g& s0 b! q. \ for(i=1;i&lt;=num_st;i++)
9 q8 c- @$ m& }2 @+ k( t  C% a  if(a[base_in]&gt;0) {
   temp=b/a[base_in]+1;
   break;
  }//temp赋初值
% I' b6 b# g4 r# X' b' Q% S% V for(i=1;i&lt;=num_st;i++) {
  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
  o7 Q& F0 K! F/ r   for(j=1;j&lt;=num_ar;j++)
* \/ c; i; ^; y0 l    if(base==arti[j]) {
. Q  |: H. Q  `- u" t     base_out=base;
     base_elem=i;
     temp=b/a[base_in];
     break;
    }
  }//人工变量优先换出
  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
# }+ J% K) d* J7 b6 m% S   base_out=base;
9 y' v3 U) m( X( R  n  N+ s   base_elem=i;
( G2 F3 h# \7 _* s6 _1 C8 n; w   temp=b/a[base_in];
  }
}</P>
<> printf(" 基变量：");
( y8 @+ Z+ y, y% F8 ], k5 k for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
' t4 O* o% b0 C5 e! { printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
//基变量变换，进行新方程初始化后迭代
for(i=1;i&lt;=num_st;i++) {
& M6 M7 ~/ F6 _+ P# o' P  if(base==base_out) base=base_in;
  d0 @) y. _+ m& ]/ N, a6 o }
//初始化主元素行系数
value_be=a[base_elem][base_in];
* M  I) V/ e/ y: T6 d b[base_elem]/=value_be;
" |6 W" o3 ~; N4 T( y for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
<> for(i=1;i&lt;=num_st;i++) {
' G5 \4 R; a& U6 A  if(i!=base_elem) {
   b-=b[base_elem]*a[base_in];
   value_be=a[base_in];
8 a: B# ?# `0 J, E/ i( |   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
( c( K/ Q0 B, |9 \8 Y  }
2 G# }. N/ m9 v- Z/ V4 `- ]+ l }
stop++;
' b2 |( K4 U' c5 \ if(stop&gt;STP) {status=-2; return;}
0 }" P& o. A: \% r* Z$ \6 y- `6 y" l iterative();
9 S" G$ N! B0 G# C2 v. D}</P>
, Q/ x4 `. Q6 W( H" u  P; \2 ^<>void output() {
int i,j;
double X[MAX];
) m7 v& E  r9 D printf("\n结果如下：\n");
2 B$ ^4 s! Z; j$ i1 s printf("\nX=(");
$ ^# \6 K+ n  q% m) B for(i=1;i&lt;=num_x;i++) {
  for(j=1;j&lt;=num_st;j++)
  if(i==base[j]) {X=b[j];break;}
$ {( T$ P( Q8 k& O( V$ q5 X& V  else X=0;
  printf("%lg ",X);
, Y) z8 O) R1 i+ ~8 n& I2 |) v1 i& G }
- N* y( p( o: x$ x+ B' \& C printf(")");
for(i=1;i&lt;=num_x;i++) max+=c*X;
if(ma_mi==1) printf("\nMax z= %lf\n",max);
: h& P1 ]7 r9 H9 a0 C) @! R5 }1 H else printf("\nMin z= %lf\n",-max);
! ^' B0 h# G8 y& J2 U  K/ }3 c/ w}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了