求助：单纯形法的程序(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])
{
/ h+ d3 R8 B4 F  M, E% L  int j;
! c# k+ w& t7 t0 x# o$ g* Y  *c=0;
, F! T( v* W7 M, A. D- w  for(j=1;j&lt;=*is;j++)
  h1 r, @. E' r6 ]6 k    {
      if((*a)[*ir][j]-*c&gt;0)
8 s6 l8 ~+ G4 N: b' s( p) p* f* X    {
! X9 k0 Z+ A: u) k6 T: E      *c=(*a)[*ir][j];
5 i3 j5 a% T; j+ D      *j0=j;
    }
' c# C, x0 e% ~8 y/ k     }
}</P><>      /***************** 参数说明 **********************/
      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
      /********* <a href="http://happyyangxu.home.sunbo.net*******/" target="_blank" >http://happyyangxu.home.sunbo.net*******/</A>
$ s6 t- S& C6 {int xi_sm(int m,int n,int m2,int mn1,int l1,float (*a)[X][Y],
: _( a7 u$ s" P$ N# }, y       int(*k)[],float(*x)[])
{
* |2 m2 V' \2 B7 L) T, A" H  int i,m1,mn,j0,i0,j;
. ~7 w% D" y9 Q* O% H# v+ P  float c,g;
! m) ~6 P  H; Z' D, F' `  m1=m+1;
% \; J7 j" {9 f6 T$ ~  mn=m+n;
% R- ^! y% ?( t' e5 S  for(i=1;i&lt;=m;i++)
' V. H  l! T9 R. v: }4 h+ [7 A    (*k)=n+i;
+ ^0 G/ K3 v: v+ [+ zleap1:
$ R$ M( ^" I7 C, U8 a3 F  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  else
     {
2 N$ G( b4 {% p# g6 ileap2:   if(l1-50==0)
      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
) J3 v( s: u/ b( T, g9 M2 G       else
      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
* c' [8 U4 t5 w$ M% e' v* A     }
+ `8 j2 p, P* [% Y9 b; w) x   c-=1e-8;
, O. X- v+ i; i( h5 @1 p   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
      {
    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
    return(-1);
      }
   if((c&lt;0)&amp;&amp;(l1==1))
       {
     l1=50;
     goto leap2;
' B2 W# D6 j8 ^       }
- o. u5 v9 h+ @& p    if(c&lt;=0)
       {
( l. G1 ~. S7 P. ]" ~4 ^      for(i=1;i&lt;=m;i++)
( E" y( W) u# u3 W7 F4 {         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
      for(i=1;i&lt;=m;i++)
1 d0 v8 {5 s8 B4 m% d         printf("\ni=%d,x=%f\n",(*k),(*x));
      printf("\nF=%f\n",(*a)[m2][mn1]);
% s# Y% o/ j$ ]5 `% x5 G" Q" }      return 0;
' Q" t1 F5 i; H/ @' @    }
     c=1e8;
* H5 a$ @; \, Z0 x( S0 m* W4 ]8 U     for(i=1;i&lt;=m;i++)
; G; i! O$ O# k0 ^' {5 F       {
5 K, }6 O# n! P. ^) M% ?     if((*a)[j0]&gt;1e-8)
& H$ X7 f# P! ?) z        {
          g=(*a)[mn1]/(*a)[j0];
          if(g-c&lt;0);
3 d. P# Z" x' Q8 }0 e         {
           c=g;
           i0=i;
          }
         }
    }
7 y0 P5 C: ^1 y7 o' s1 G! X& {      if(c==1e+8)
6 T5 B& }4 l9 s* z' {    {
      printf("\n\t*******Lp no solution********\n");
      return -3;
     }
+ @3 K% \( j7 r+ D9 I       (*k)[i0]=j0;
4 s5 c/ ]2 O. s. q       for(j=1;j&lt;=mn1;j++)
      {
- D$ F% N6 [7 Z. e6 k$ ~! P        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
           continue;
# {$ O" b$ }) o# r         g=(*a)[i0][j]/(*a)[i0][j0];
7 l: Y0 }5 z0 _/ r8 O/ e# [           (*a)[i0][j]=g;
' |4 m  s& r5 n8 S, |( _- h) J         for(i=1;i&lt;=m2;i++)
           {
) j. G4 q  x  t# e         if((i==i0)||(!(l1==1)&amp;(i==m1)))
+ {: g, r" G% @            continue;
9 }# u" K/ y- k) R: s+ X6 Y) n         (*a)[j]=(*a)[j]-(*a)[j0]*g;
% c( b' O3 t  ?  d& d% k5 o           }
       }
       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
1 L0 ]: f% A, _; P5 h, f0 O       (*a)[i0][j0]=1;
       goto leap1;
}  
main()
% l) d: V# C! ]8 m, z5 Y{
5 N, v, t- L3 v4 ?9 l  float a[5][7]={{0,0,0, 0,  0,0,0},
2 {3 J+ }6 x# _                 {0,3,-4,3,  1,0,12},
                 {0,3,0, 6,  0,1,12},
                 {0,0,0, 0,  0,0,0},
: B+ {8 T* V% b' f* ?% S  w                 {0,69,0,144,0,0,300}};
  int k[3];
) H7 i6 g& Z6 }/ [/ H2 v( _+ D8 M  float x[3];
  clrscr();
  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
& b2 l7 f: W6 e! y}
0 m. L2 a" F; x- F# b8 P0 G</P>

plgatc

<>/*************************************************************************
                    单纯型法解线性规划问题（两阶段法）        
9 J0 P. }/ U4 a, P/ B1 p0 ?) ~! H
! g. C$ i. z# D" H3 X. r/ a5 p9 ` 编程环境：VC++6.0     
# m- `9 d4 ^* s  z2 m( l 方程组输入说明：
变量非负，按提示输入相关参数。
% g* W. n- K5 U' B# b) b *************************************************************************/
4 T$ ~+ L* Q, h6 }$ F$ f#include &lt;stdio.h&gt;
  B: P) i: q! ]  x) y#include &lt;stdlib.h&gt;
( z* ]) {6 X  U& U' Y, G0 |) G#define MAX 100
#define STP 100</P>
' ]0 q+ ~- ]1 c+ _9 M* o" j. w<>int stop=1; //迭代记数变量
" j7 j# {% K) [) W: Y+ o9 xint status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
int step=1; //目前阶段</P>
  O" @, s# ^7 {5 x5 ~<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
$ t9 q, ?& j; p! w' I5 s; {& h+ Vint num_x; //变量个数
int num_st; //约束方程数
! S+ z6 B# U/ G  z) Kint num_ar=0; //人工变量个数
int arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
$ b  E% L1 c+ C1 o1 O# ]$ r, T0 ?+ P<>void create(); //建立方程组
void iterative(); //单纯型法迭代
$ L) z: i5 b# c. w# U" wvoid output(); //输出结果
void banner(); //打印程序标题
void exchange(); //交换两阶段价值系数
void show(); //输出方程组</P>
4 F0 H! J/ @* r0 M1 r7 h<>void main() {
1 y2 d8 T9 x; [3 T# O; V: G int i,j,k;
# W8 {: @: z% A) T. F4 q( v banner();
* L) g' @! a/ ^+ y* ^6 R( q create();
//保存原价值系数，转换为第一阶段价值系数
0 v6 F  \0 L. ] for(i=1;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
  ?* v4 h9 E2 n7 T  if(k==1) temp_c=-1;
! B' Z2 Q5 y0 B" u' `  else temp_c=0;
}
. ?1 B! k3 ]# @: ~" o$ |2 P exchange(c,temp_c);</P>
% \0 J1 B: ?9 O% q, W<> printf("\n\n第一阶段问题为：\n\n");
show();
& d  N& ^4 m3 U( P( d- W  l) w8 Q step++;
printf("\n\n按回车开始第一阶段迭代");
getchar();
3 q- }5 a9 w$ _, n" J% F" v getchar();
" P# O4 e4 H1 m; W% y' _0 \& z iterative();
% `' @/ [: k' R5 E5 | if(status==-2) {
9 j) v( }6 z/ n1 [6 L  puts("迭代超过限制次数强行终止！\n");
8 P: s/ V7 q* E+ {6 i) V% ]  puts("\n按回车结束");
  getchar();
6 U; I  w5 o9 f1 v) d+ i  exit(0);
# y$ }" t. f$ @) u8 u }
output();</P>
<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
  g% `$ p* u) M  puts("\n按回车结束");
  getchar();
  exit(0);
}</P>
<> //转换为第二阶段价值系数
' n% L+ B1 f, I! G4 h) a2 B4 a0 i exchange(c,temp_c);
, @9 |" T) w! G+ \; I //把人工变量列全设为0
for(i=1;i&lt;=num_ar;i++) {
9 X. k4 u& z6 I1 L, _/ \  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
8 r* N  W- m) c& [& f4 g }</P>
$ p% [, N8 P0 \2 \<> puts("\n\n第二阶段问题为：\n\n");
show();
puts("\n\n按回车开始第二阶段迭代");
getchar();
0 P$ w' b; z1 s  e, n: R iterative();
switch(status) {
case 1:
! t1 \, Z3 a5 K, ?# t9 s) G& \, i  output();
( i5 u! k# ^' j( U0 W9 C. W  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
5 n( p# R# e7 W+ {" ^" L- s  getchar();
" {6 r3 j# u& F9 E7 Z  exit(0);
case 0:
! j0 ^* T/ V$ v  y  puts("\n\n原问题为无界解。\n");
/ J; w( b8 B$ m3 E  puts("\n按回车结束");
  getchar();
7 q! Z$ b" s" z! U5 ?% r. c  exit(0);
case -1:
  output();
% r4 `& c" {* @3 C" x  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
! i0 P3 Y6 k# P  getchar();
' O& c+ c. d: Z* M5 W3 b  exit(0);
case -2:
: k: @; H. {1 w2 q' b  `  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
( C: b7 X( a  [* k7 h! v% z  getchar();
7 Y0 v+ y% ?% ]& `( w; r  exit(0);
}//switch
}</P>
! _. o. E6 I5 H9 i9 N& m<>void banner() {
  o% d) n" z" p; ^+ v printf("\t\t****************************************\n");
$ o* s  q2 J2 z/ o" S printf("\t\t         单纯型法解线性规划问题\n");
printf("\t\t                         作者：Thunder\n");
printf("\t\t****************************************\n");
8 I* R7 I5 i( f( f" t: o printf("\n");
}</P>
<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
2 f) U9 B6 z  O2 J//为1的不显示系数1，-1系数只显示负号
8 t0 O- e" I( l" m/ x" ^$ X int i,j,k;
' T  B8 e) f2 h5 G+ {; P1 p% c! r: ] switch(step) {
case 1:
  printf("min z= ");
  printf("x[%d]",arti[1]);
  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
# M! u3 P& T- @  break;
case 2:
5 D- d' g1 {7 K0 ~ printf("max 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);
- f5 s; k3 u5 n$ J0 [, `, e  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
- [( ^/ |4 p0 U# b6 r4 l9 [: _7 _  else printf(" %lg x[%d]",c,i);
( {; f1 P0 P$ V% h- B; J }
6 A4 i" v& }* u7 E; H8 O break;
}</P>
<> printf("\nst:\n");
3 b) i; {* y& r7 m; N. m for(i=1;i&lt;=num_st;i++) {
  k=0;
9 I9 [: ?5 D! f! S  ?, [' k  for(j=1;j&lt;=num_x;j++) {
% u# h" ~9 W% m) W$ \+ t   if(a[j]!=0) {
; y4 U1 u- O; W- F0 V    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
0 e. N( o# h( y, o% U6 _5 @: n    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);
4 r4 i9 s; \5 T9 p/ b" u, }    else printf(" %lg x[%d]",a[j],j);
    k=1;
1 h* l- c# b. e+ l+ {   }
$ U6 {. v8 D1 z4 d# t0 P: Y  }
  printf(" == %lg\n",b);
}
" c/ B* v! v) O+ n& n6 u5 @/ x+ i* _ printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
8 \/ ^( J$ \5 J<>void exchange() {
int i;
double temp[MAX];
- N( r- l! Z4 K4 I/ u for(i=1;i&lt;=num_x;i++) {
4 c7 p( q) ]# ^1 O- q4 ^5 y; I( F  temp=temp_c;
/ l! x6 ?/ N1 b8 C  temp_c=c;
* J8 l2 q. q/ o8 H( V" `  c=temp;
4 U5 m6 q3 w6 [ }
& B  J0 D9 z( E}</P>
<>void create() {
9 x# S6 X5 {2 A4 u& T0 G: G# [//输入方程组系数，每个方程输完后回显确认
! p9 n& a7 M3 C0 X9 W int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
, V! z) s5 g# `5 ? char confirm;

1 ^" J6 H+ R; x while(1) {
3 |. ^' A3 Z/ o# x# r" I! D6 `  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  scanf("%d",&amp;ma_mi);
$ z" X- ~4 T1 o2 O  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
  W  m. g8 v# P; E  else break;
}

# r: B2 \: l3 y5 ^ while(1) {
  printf("指定变量个数：");
" t/ z8 O5 r1 ?! r! p) N  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);
. x" G/ k9 }0 Q& Z6 f  }
) A$ @- Z# r) d- e  if(ma_mi==1) printf("max z= ");
2 x$ L& n8 p  x, a% T+ X  else printf("min z= ");
  printf("%lg x[%d]",c[1],1);
* U6 s+ [* B2 W& X0 v9 T  for(i=2;i&lt;=num_x;i++) {
   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
   else printf("  %lg x[%d]",c,i);
  }
. o* e- `4 f. U8 I printf("\n正确吗？：(y/n)");
* t2 A) A# g. |# W getchar();
/ R4 y! Q) y7 f% P+ Z0 C; | confirm=getchar();
if (confirm=='y') break;
else if(confirm=='n') continue;
- d8 ~- z! I) W! q0 W& @4 _" j# S }</P>
  Y$ w+ {& z/ V9 k$ C" ?( q# U4 e<> printf("输入约束方程组个数：");
& I( d& i1 k9 Q& s, A  d4 x scanf("%d",&amp;num_st);
& V) I- v3 b+ c. M7 b5 q for(i=1;i&lt;=num_st;i++) {
  printf("st.%d:\n",i);
: u8 W- `* K/ x1 Q+ Z: U  while(1) {
# j: Q" w/ I. t% t- ^7 {4 e   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
6 B- j. U; e' {! d# m% b+ z   scanf("%d",&amp;re_st);
   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
: @7 x( N* n4 S/ U5 G8 g   else break;
  }
  printf("输入技术系数：\n");
% _6 H+ `: v+ D3 \2 G! b  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
; Q: ^- u9 G5 f+ m) K& ?  }
5 P) N1 a2 h% e, O1 m* c  printf("输入资源拥有量：\nb%d=",i);
  scanf("%lf",&amp;b);
' f% q! N$ E2 ~# |# d: G# Y5 y3 C" [  
  printf("st.%i:\n",i);
  printf("%lg x[%d]",a[1],1);
  for(j=2;j&lt;=num_x;j++) {
7 n& a; h3 s& [. T  F   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
   else printf(" %lg x[%d]",a[j],j);
  }
8 `7 c' Y: k- q! b2 _4 w  J  switch(re_st) {
9 M& ?" M% ^. Q6 G5 s4 W   case 1: printf(" == %lg",b); break;
8 Q' [: T# \# J: ^/ y4 }$ q   case 2: printf(" &gt;= %lg",b); break;
+ h9 [6 q- M1 g4 u: @, s: d# ]0 C: H   case 3: printf(" &lt;= %lg",b); break;
  }</P>
<>  while(1) {
   printf("\n正确吗？(y/n)");
   getchar();
7 p; W& D# i" g, _3 t! @   confirm=getchar();
   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
  }
: h. t3 S6 f# ~- d. h/ R }</P>
6 f, b4 J9 Q2 Z, S* a<>//显示输入的方程组
2 P8 d9 P& j( i2 g printf("\n原问题为：\n\n");
if(ma_mi==1) printf("max z= ");
else printf("min z= ");
printf("%lg x[%d]",c[1],1);
4 E$ l/ [4 ]  Y4 s- f. d0 J for(i=2;i&lt;=num_x;i++) {
5 o! B" D" j' @  if(c==1) printf(" + x[%d]",i);
' A( Z" p- [3 Z/ a2 g9 }  else if(c==-1) printf(" - x[%d]",i);
. ^' F" [4 ?$ r1 W/ I0 ]  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
+ R, s! _6 f8 c2 {! ?1 {, W3 j  else printf(" %lg x[%d]",c,i);
6 f) o/ z+ |( Y) ~ }</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
& |& }' [7 Z% I. i/ z2 t' Z  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);
    else printf(" %lg x[%d]",a[j],j);
    k=1;
   }
  }
$ |7 T& X, w' x% N1 b  switch(re_st) {
   case 1:
    printf(" == %lg\n",b);
  x+ [0 q6 ]1 G  c9 o; _    break;
1 t) ~% l0 X; x   case 2:
4 y7 l0 R7 V; P' v+ R6 D( E    printf(" &gt;= %lg\n",b);
: X; k- p+ Z* M# u" z    break;
# O1 g- b, F3 C4 N$ r   case 3:
    printf(" &lt;= %lg\n",b);
4 B& ?# d3 w& T! m+ f+ G2 ]    break;
2 k/ K1 a# K% l: b) n0 {  }
  J; n, R, c+ H }
1 K" Y- A1 v& b# F, D. [$ n* X printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
* ^0 f& F! s% O0 i* c/ L<> tnum_x=num_x;
for(i=1;i&lt;=num_st;i++) {
4 H& P  F" l! ]+ [' j  switch(re_st) {
  case 1:
( x2 v: H8 [! r5 _8 r& F- ^  case 3:
' Q, u( `" \! V   num_x+=1;
   break;
  case 2:
   num_x+=2;
   break;
8 |4 r. r. U2 p5 ^/ S  }
$ G8 o" q1 b: B2 W( Y. w0 h* [ }</P>
" t2 u# i) K4 ~  L  K1 L! ^% b<>//化为标准形式
if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
for(i=1;i&lt;=num_st;i++) {
0 r! M+ d; y, M. V  switch(re_st) {
   case 1:
' d/ {0 r+ V: k% X" i0 y9 l+ B    num_addv++;
0 E6 C# k/ x4 q    num_ba++;
- A, x( A1 F  _2 D/ z    num_ar++;
2 c9 S3 p+ j* u3 A& C4 a    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) a[tnum_x+num_addv]=1;
     else a[j]=0;
+ |% h8 u! {8 D5 v. b    break;
0 D* i" |8 Q* ?) Z   case 2:
; n$ I. T8 g& E5 J0 z, E  J2 Y0 @5 E    num_addv++;
4 [( Z0 g/ ~8 k5 U2 p/ }5 w3 W    c[tnum_x+num_addv]=0;
    num_addv++;
    num_ba++;
6 H) U; J5 f. K8 M. h    num_ar++;
5 V( [' s; x5 |8 p8 L    c[tnum_x+num_addv]=0;
    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
* y( {; I, ?+ T! P    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;
  v3 @& [9 l: A$ G+ J+ I    break;
   case 3:
% Q, P8 x# W$ {    num_addv++;
: E2 @. W" w$ m% g# q4 G: N    num_ba++;
2 b2 O2 R, D5 K+ b9 d, y8 l    c[tnum_x+num_addv]=0;
2 e( G. i) A! B    base[num_ba]=tnum_x+num_addv;
' j3 q! ^+ p& @+ T    for(j=tnum_x+1;j&lt;=num_x;j++)
     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
/ G7 X4 C7 z2 r7 R9 ^: T  A0 T/ S' [$ b$ ]     else a[j]=0;
) U3 z/ B& r$ f7 A6 |1 G* C    break;
  }//switch
}//增加松弛变量、剩余变量、人工变量、确定基变量</P>
<>//显示标准化后的方程组
5 \5 v4 A( m' N. I; E& a printf("\n化为标准形式后：\n\n");
if(ma_mi==1) printf("max z= ");
( n  s" k3 u/ e2 d8 E else printf("max z'= ");
5 ^6 {& K5 G8 A4 k printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
8 R& s; U) F1 @) W- @7 V  k=0;
  for(j=1;j&lt;=num_ar;j++)
5 t! w! _+ V0 e2 p( @  if(i==arti[j]) k=1;
  if(k==1) printf(" -M x[%d]",i);
  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);
. a, R. t# o+ X2 Y5 F; A! h2 J  else printf(" %lg x[%d]",c,i);
}</P>
<> printf("\nst:\n");
6 a4 l- H" D7 k7 {8 |5 j+ a4 { for(i=1;i&lt;=num_st;i++) {
+ J8 `% H/ }4 k% Y  _  z  k=0;
' Y0 R0 w$ O1 g6 K# Y6 k  for(j=1;j&lt;=num_x;j++) {
0 \, B2 L" c5 \+ D   if(a[j]!=0) {
+ P. s: t2 U2 i9 T9 J& _4 A$ E1 d    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
( i% `1 u8 a7 Y. J9 w) b; P    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
    else if(a[j]==-1) printf(" - x[%d]",j);
- \* V. b+ t$ X  l& G    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;
8 A. ?% x3 J5 ^3 v" s5 y0 C- G0 U   }
  }
! c6 g- a% F" X; i/ \  printf(" == %lg\n",b);
: ^$ r1 n7 i1 B3 i+ c& D; [ }
! s3 m; y1 y  j printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
" H- I; n) l( O. B<>void iterative() {
int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
5 D" J% u- d' b; g4 w8 q; r- |) C2 T int base_elem;
int base_out,base_in;
double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
<> printf("\n\n第%d次迭代：\n\n",stop);
for(i=1;i&lt;=num_st;i++) {
% u; l! X' N, l/ z! @  printf("c%d=%lg\t",base,c[base]);
  printf("b%d=%lg\t",i,b);</P>
! g- N& B' \2 X  j+ z  c! I<>  switch(step) {
9 l4 n/ W) {4 T: _' ~   case 1:
    for(j=1;j&lt;=num_x;j++)
" ~3 j$ \  G: X3 o- Z! g    {
     printf("a[%d][%d]=%lg\t",i,j,a[j]);
    }
    printf("\n");
    break;
5 `5 V" @/ @. E3 ~2 u   case 2:
1 z* K  x( Y1 s+ _/ b% I+ T0 Y8 B    for(j=1;j&lt;=num_x;j++) {
" _- V- Y& x8 L5 w! A1 N7 R6 W     k_a=0;
     for(l=1;l&lt;=num_ar;l++) if(j==arti[l])k_a=1;
8 e& z7 z3 V1 ]8 p* _     if(k_a!=1) printf("a[%d][%d]=%lg\t",i,j,a[j]);
    }
  C  q: s, W7 N* f    printf("\n");
    break;
1 G: Z2 M) g/ d6 G  }
1 Z! G' i/ j3 o; U4 n8 ^ }
3 l$ n7 n0 z2 w2 o5 g//求检验数sigma
" x7 q- v8 t4 G9 c9 ^% ]! r& N7 f& | for(i=1;i&lt;=num_x;i++) {
4 e3 x2 b& L/ G# n& `( i8 v  sigma=c;
, K4 ^) F$ _, Z6 e  for(j=1;j&lt;=num_st;j++) sigma-=c[base[j]]*a[j];
- _7 ?4 v( N& q3 i) O  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
  switch(step) {
" x$ P7 i) u; l1 V( ^" c9 j) m/ O   case 1:
    printf("sigma[%d]=%lg\t",i,sigma);
7 V  ^3 {* V( I; r7 c/ i3 n    break;
   case 2:
; u% k! l+ \0 I- @# r. A    k_a=0;
5 V! V4 r. [7 H3 g. H( }    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
0 X6 F, w7 Z# e  j1 H) M1 U) o' P    break;
6 A" ]% {! K* \8 z+ `, c+ b. M+ i  }
}
putchar('\n');
/ p0 B6 G) \1 ]% _& I4 Z2 H2 q//检验检验数sigma是否全小于等于0
* b2 {. @6 c/ |; E6 u k=0;
; }- ?8 e: p9 r" a" v% A% N- w% I! U for(i=1;i&lt;=num_x;i++) {
  if(sigma&gt;0)
; T  a/ c% L$ j6 c! O  k=1;
}
0 b0 H, u7 y: a if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
  for(i=1;i&lt;=num_x;i++) {
   k_f=k_a=0;
' Z. A8 m% p5 V, ]0 U0 A  D/ m% Q   for(j=1;j&lt;=num_ar;j++)
   if(i==arti[j]) k_a=1;
/ f4 t  x% S2 l1 w0 R  e   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;}
# h) o- n: {- a# B# g8 d" ~  Y   }
4 Z! e) U0 u$ ^- Q/ m2 l6 |  }
  status=1;
  return;
}
//检查是否为无界解
% w2 }5 w# s  }3 m$ n& f6 `9 f for(i=1;i&lt;=num_x;i++) {
6 ?9 ~/ @6 _9 F" p- R4 F5 }  k_f=0;
  if(sigma&gt;0) {
   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
; B. p" n& C% \  Q. w$ P- o3 O   if(k_f!=1) {status=0; return;}
+ v- f+ I9 r' k3 P) L$ D' t  }
8 t4 |' L  T7 t9 L, v# M }</P>
9 j# W: v' A8 F6 {' s; C1 t" r<>//确定换入变量
for(i=1;i&lt;=num_x;i++) {
2 u) X( |3 Q+ |: a! W/ g) u  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;
9 Y: h' h  Q- U6 Y }//temp赋初值
for(i=1;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
) |$ b( ]5 x: n! h3 X" Y8 o$ X  if(k==0)
) I9 }! ^+ O9 d- t) w7 T   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
$ a2 [, `; d* P0 |6 O' o    base_in=i;
. ~3 n0 ]: X! p    temp=sigma;
) v1 J7 \2 [& v1 l   }
0 c5 a2 \; I5 R3 v }</P>
<>//确定换出变量
for(i=1;i&lt;=num_st;i++)
. x. |  u: x, u0 Y  if(a[base_in]&gt;0) {
% c4 ^! a+ n* j: y2 s   temp=b/a[base_in]+1;
8 [7 _) s- b, J) I9 b   break;
  }//temp赋初值
for(i=1;i&lt;=num_st;i++) {
& g( g" ]6 J7 r* J4 W  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
   for(j=1;j&lt;=num_ar;j++)
" B# r3 k+ C# e    if(base==arti[j]) {
, y) V4 N2 v) [* T1 T     base_out=base;
- l# b" V7 i' B" Z     base_elem=i;
     temp=b/a[base_in];
     break;
    }
  }//人工变量优先换出
) H, _) Q6 I' ?- F2 r. Q6 c* N  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
  F' G7 _& o5 A$ G- n# e* w: l   base_out=base;
   base_elem=i;
$ o: Z8 Z, f( l7 @3 D0 R/ j4 o6 C$ ]0 K   temp=b/a[base_in];
  }
}</P>
5 }. f/ _6 l5 f$ U* K<> printf(" 基变量：");
for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
0 C3 [6 H0 r4 A# i$ R0 _* q1 q printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
2 I3 m1 g& g! Y" o//基变量变换，进行新方程初始化后迭代
for(i=1;i&lt;=num_st;i++) {
  if(base==base_out) base=base_in;
}
0 u: B- ?/ S  a; K7 C) _4 V0 f//初始化主元素行系数
! G" P, a: M  h- u" x! P  E value_be=a[base_elem][base_in];
b[base_elem]/=value_be;
for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
) t6 F* h7 o8 Z7 A6 s$ ~. R# @<> for(i=1;i&lt;=num_st;i++) {
  if(i!=base_elem) {
   b-=b[base_elem]*a[base_in];
" s- t. H, Q  k6 E! a6 S8 X6 }   value_be=a[base_in];
   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
- ~- X, K& r$ U/ }  }
: b" G$ `! C) k. O, w }
stop++;
! z2 \" \3 v- \: C/ x if(stop&gt;STP) {status=-2; return;}
2 v' N% I3 \% C: Y* O. c8 n0 d iterative();
}</P>
3 s6 T. h0 ?# f4 f<>void output() {
8 u1 n+ i' y) X. y& a" N) f int i,j;
' D9 w1 r2 h" g! Q! |# q$ T double X[MAX];
9 R; B/ `8 R9 w( f printf("\n结果如下：\n");
, }1 z% l( O, G2 [& z, i printf("\nX=(");
for(i=1;i&lt;=num_x;i++) {
' J; x3 |% A! l  p7 ~& Q  for(j=1;j&lt;=num_st;j++)
  if(i==base[j]) {X=b[j];break;}
  else X=0;
  printf("%lg ",X);
) f% S% [2 H. O }
, B: s0 E5 ?$ ~; v printf(")");
# o7 \* s5 }/ r' Z5 ~ for(i=1;i&lt;=num_x;i++) max+=c*X;
) q' j) t( ?- p5 O% A' ~& B& G if(ma_mi==1) printf("\nMax z= %lf\n",max);
: w; k6 V$ D/ Z else printf("\nMin z= %lf\n",-max);
3 C. V) T+ Z8 W; c8 b% ?% d4 Q' \5 G}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了