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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
#include&lt;math.h&gt;
#define X 5
#define Y 7
: ~7 V# ^$ ~2 d( L6 K$ l) ^void xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
8 ~5 z3 ?' ~6 X9 L{
6 i5 N& k7 C2 }# x1 h% m  int j;
  *c=0;
  for(j=1;j&lt;=*is;j++)
6 C" _3 w) [9 L/ ~    {
      if((*a)[*ir][j]-*c&gt;0)
    {
. c# U3 q- `) p# J4 f      *c=(*a)[*ir][j];
      *j0=j;
3 G) x6 y- l. A/ e    }
     }
}</P><>      /***************** 参数说明 **********************/
      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
2 {$ S( k: e/ n0 p0 ?      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
      /********* <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],
7 p/ d) p$ M8 J       int(*k)[],float(*x)[])
' F" W5 v3 p) G. {3 T( H. r. a{
  int i,m1,mn,j0,i0,j;
/ l( U! N1 r- D# Y  float c,g;
  m1=m+1;
  mn=m+n;
  for(i=1;i&lt;=m;i++)
% G* n$ Q' \. _' f7 l) |$ H0 g    (*k)=n+i;
- z  w0 h: w! Aleap1:
2 K% q, q# G/ f4 A! g4 T) j  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  else
3 j2 M- X9 d, X4 i, w# a     {
# _- Y8 K& _& D# d8 c* U# Cleap2:   if(l1-50==0)
% I: s* Y/ @1 y      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
       else
( L( i  v% H- P# O+ U      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
     }
) D2 h1 U+ l6 \/ o, b  w   c-=1e-8;
; M' j4 U. j) o2 H   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
3 X5 v# R1 f* G0 f4 [9 [, V      {
2 x1 Q) B- o7 t  J5 [    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
" Y$ ~8 C) J  v) S9 q- ~& |3 S' I    return(-1);
* h% a% @. S' l/ e      }
   if((c&lt;0)&amp;&amp;(l1==1))
& Y7 ?8 j1 j# N       {
; M( _1 b' b2 m$ |     l1=50;
$ ]* X- \3 m3 e6 D7 H5 ~. E6 r1 R3 ^     goto leap2;
# [! U8 y" m4 i7 s7 b0 ?! M       }
    if(c&lt;=0)
$ a; {/ ]6 s2 [1 v6 z       {
      for(i=1;i&lt;=m;i++)
         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
3 G* |( @( }$ \$ [( Y      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
/ h( L( {9 e8 Q  c5 h1 z      printf("\nF=%f\n",(*a)[m2][mn1]);
      return 0;
    }
     c=1e8;
  h# I6 V& m! Z! G8 e$ D, ~  Y     for(i=1;i&lt;=m;i++)
! W! N4 \1 K3 D# }6 K2 O) F9 o       {
     if((*a)[j0]&gt;1e-8)
        {
/ X6 t1 a" o4 h+ {          g=(*a)[mn1]/(*a)[j0];
          if(g-c&lt;0);
         {
           c=g;
/ D  F4 S" c0 Q$ h           i0=i;
          }
         }
# e0 A0 }% s. G& `% ^6 g    }
      if(c==1e+8)
# Y6 j& k7 L' P8 j+ I3 E    {
" F8 x& f! U2 n: w: ^      printf("\n\t*******Lp no solution********\n");
0 r. I1 Y/ O6 Q; W8 ?6 {) ^  Y- M      return -3;
     }
       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
) Y. t. Z/ k9 _1 @( D      {
1 D9 _- N9 M6 U3 G2 M. P        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
           continue;
; `9 _: {6 q* y) [( c  }" {- x         g=(*a)[i0][j]/(*a)[i0][j0];
3 Q! @; T* ~+ [/ {           (*a)[i0][j]=g;
3 o' f9 [/ s7 }9 u         for(i=1;i&lt;=m2;i++)
; f/ {0 M4 ^  [2 J4 P           {
         if((i==i0)||(!(l1==1)&amp;(i==m1)))
            continue;
         (*a)[j]=(*a)[j]-(*a)[j0]*g;
2 D7 R/ u$ D/ i, Z0 M6 S# z           }
1 K' l' d7 P) P, O       }
       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
! J; F1 D( C) ^4 S/ ^2 j: l       (*a)[i0][j0]=1;
, {' W+ L( ?+ `" \; r: X# J0 ?       goto leap1;
}  
( @1 U" H: Q* _8 z, z6 P" [; B/ m& dmain()
; ~  v' F. ]# c' v{
+ z/ y! e8 m( a, Y# U) x  float a[5][7]={{0,0,0, 0,  0,0,0},
6 `( i+ a7 r7 M* }5 O" x) E                 {0,3,-4,3,  1,0,12},
" {0 m, ^* Q) a( j" a                 {0,3,0, 6,  0,1,12},
+ N4 ~; }4 @0 N1 w2 h                 {0,0,0, 0,  0,0,0},
                 {0,69,0,144,0,0,300}};
- u+ G! }2 l5 w: w  int k[3];
  float x[3];
  clrscr();
- j  d- n2 s0 l0 x7 g* v3 E  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
}
</P>

plgatc

<>/*************************************************************************
9 r" y; y0 {( y% K6 p* u' `! E                    单纯型法解线性规划问题（两阶段法）        
6 p7 d! D! M% e9 E
编程环境：VC++6.0     
方程组输入说明：
变量非负，按提示输入相关参数。
*************************************************************************/
#include &lt;stdio.h&gt;
& A9 J6 L/ {0 W; E, p#include &lt;stdlib.h&gt;
#define MAX 100
#define STP 100</P>
; T1 u: U$ p4 K<>int stop=1; //迭代记数变量
int status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
) V0 G% d* K" ]9 a7 m2 t* ^int step=1; //目前阶段</P>
- r: U# _( _9 e; }$ ]<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
int num_x; //变量个数
% j9 o5 y5 E8 s% Pint num_st; //约束方程数
int num_ar=0; //人工变量个数
; d$ r2 V# D) x, gint arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
<>void create(); //建立方程组
void iterative(); //单纯型法迭代
. D, x# t6 B0 t0 A$ Tvoid output(); //输出结果
void banner(); //打印程序标题
void exchange(); //交换两阶段价值系数
$ I) X/ w  h$ P' M3 M, qvoid show(); //输出方程组</P>
! T; a9 o5 r5 ]# c4 c<>void main() {
7 `! g! U" e* F2 | int i,j,k;
banner();
% O9 y, e9 f& b  }2 R" d; f0 { create();
) v4 d$ ~' j, j+ n5 N //保存原价值系数，转换为第一阶段价值系数
for(i=1;i&lt;=num_x;i++) {
/ W, V" v! \5 b# b! S0 L4 R  k=0;
2 k( h' j( s  O" [! B  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
  ~2 V$ k$ i) h$ Q/ y  if(k==1) temp_c=-1;
  else temp_c=0;
9 Z) r) l9 _2 z/ { }
exchange(c,temp_c);</P>
<> printf("\n\n第一阶段问题为：\n\n");
show();
step++;
printf("\n\n按回车开始第一阶段迭代");
getchar();
: Z0 S) k% ~9 O  k5 G getchar();
iterative();
if(status==-2) {
  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
  i' A# j5 `/ N3 `7 i  getchar();
  exit(0);
}
output();</P>
& M0 f) x% u& ~1 R<> if(max!=0) {
, O5 k8 a! _8 x9 l7 {# G& K% A  puts("\n\n原问题无可行解。\n");
& c9 |/ r( y7 |( V" J& I  puts("\n按回车结束");
  getchar();
  exit(0);
( \( }% Q+ J* y }</P>
<> //转换为第二阶段价值系数
exchange(c,temp_c);
//把人工变量列全设为0
for(i=1;i&lt;=num_ar;i++) {
( m5 f8 B9 L1 y  c[arti]=0;
  }, i0 F( v: e. p. `. I8 y4 x  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
}</P>
<> puts("\n\n第二阶段问题为：\n\n");
show();
puts("\n\n按回车开始第二阶段迭代");
getchar();
+ F/ s4 ^# X8 T  s iterative();
switch(status) {
& \! C( ^* M) y8 P0 Z case 1:
" o$ Q( {# G; g  output();
  puts("\n\n原问题有唯一最优解。\n");
" W8 N) R8 b+ j) ~4 c  puts("\n按回车结束");
1 I$ l* E% s% r: O& k  getchar();
, v0 E% b0 }0 n& C' e  exit(0);
( g, c5 y6 ~$ ^& X1 N case 0:
$ G: A8 u3 z  L6 \, o# \& [  puts("\n\n原问题为无界解。\n");
  puts("\n按回车结束");
0 f4 K( j* z( O- E  getchar();
  exit(0);
& i/ N) J- B& S, @  ?6 g6 t8 B case -1:
  m: b3 Z% c0 {8 m$ r0 f! V  output();
4 ]4 b  S! G$ o  j+ \  puts("\n\n原问题有无穷多最优解。\n");
, j: J# q  l9 z) n" Y. Q  puts("\n按回车结束");
  getchar();
  exit(0);
( q9 I3 i( |# y9 Z3 A case -2:
5 Q7 u7 o) j! O# n8 O9 G  puts("迭代超过限制次数强行终止！\n");
* t/ c' Y% ~$ e* l) A/ Y0 {# ?$ v- f( S  puts("\n按回车结束");
  getchar();
/ L1 t: [7 H! n% p. R  n3 x! }  exit(0);
}//switch
  P) T- ?; c. x, t( t% y}</P>
7 r1 b( \; `* w6 Z  ^3 ]<>void banner() {
7 [. \7 F" ?# F# x% g printf("\t\t****************************************\n");
printf("\t\t         单纯型法解线性规划问题\n");
printf("\t\t                         作者：Thunder\n");
printf("\t\t****************************************\n");
printf("\n");
}</P>
  }* l6 e$ y# `! W<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
//为1的不显示系数1，-1系数只显示负号
int i,j,k;
switch(step) {
8 |5 h& }/ ~! {# Q2 d0 E9 u case 1:
. ~8 H0 D0 g1 l: @; Y  printf("min z= ");
  printf("x[%d]",arti[1]);
  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
$ `$ j4 ~, b3 _/ Z$ F! V& ^ case 2:
9 ^: f9 H: ^4 l printf("max z= ");
& R7 x2 w/ f: k) n$ w3 _  ` printf("%lg x[%d]",c[1],1);
8 g# y; C2 C' C$ y4 q/ Y for(i=2;i&lt;=num_x;i++) {
. p% o4 m4 L) m4 X+ q4 S: K7 c  if(c==1) printf(" + x[%d]",i);
. a( \. W0 H$ V7 @  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);
2 }6 u# M0 r0 l5 ` }
break;
}</P>
2 f- k  ?; W& b+ U4 P/ x- N) V<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
1 {2 E/ d+ O5 k1 P  k=0;
( ~& a5 w3 ~1 w  for(j=1;j&lt;=num_x;j++) {
. ?2 `- @4 m& Z( `   if(a[j]!=0) {
5 \. K8 G$ M7 l0 m5 f/ ~8 p    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
. _% ?! N+ I, ~" G7 n    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
    else if(a[j]==-1) printf(" - x[%d]",j);
: R, n' p6 f1 q/ s. i7 t0 u    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
0 ~2 V/ n8 `1 `* H    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);
* a' F# T2 V* p9 p% a    k=1;
   }
' H- h$ Q  o- w# b( @" E  }
7 k  P) n9 E+ c  printf(" == %lg\n",b);
}
2 |- `0 G. ~* v' x% j printf(" x[1]~x[%d]&gt;=0",num_x);
6 W! u" l5 w! g7 z}</P>
<>void exchange() {
. Q# L8 }& B( u+ o& U int i;
% T( C! I3 k8 H7 d# e: r, ~ double temp[MAX];
for(i=1;i&lt;=num_x;i++) {
3 \2 j; I; q/ D6 ]( A  temp=temp_c;
7 u4 o8 }9 A+ i) [/ G8 @, |, K  temp_c=c;
# T( H8 P- r9 n: Y' J2 p  c=temp;
}
( n! U& m- g' U( Q4 ]2 v9 f}</P>
3 Z2 A  o$ V- K, |# m<>void create() {
% W8 e) h3 a! ~1 Z//输入方程组系数，每个方程输完后回显确认
8 r2 u7 F1 j# X. x* R( D9 y2 S int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
+ _* \1 e* E1 e3 ] char confirm;
2 x& n0 M) s% |6 L
; _/ A& b" Q2 M$ v- }1 R% D while(1) {
" J3 J( g+ J0 ~1 j! s% x  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  scanf("%d",&amp;ma_mi);
  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
  else break;
( r7 }4 Y; e! @1 I8 t( w }

while(1) {
8 L) F' k# i% ~: X( o$ V% ^  printf("指定变量个数：");
* R$ q' r0 j6 w% H0 I  scanf("%d",&amp;num_x);
  printf("输入价值系数c1-c%d:\n",num_x);
) F3 z: v5 O7 |% p) D  for(i=1;i&lt;=num_x;i++) {   
8 s9 q2 F1 q* U" P   printf("c%d=",i);
   scanf("%lf",&amp;c);
  }
: d" u3 @/ l7 K8 K' U  m0 j5 K  if(ma_mi==1) printf("max z= ");
  else printf("min z= ");
; @) k  \( K9 P0 ?  printf("%lg x[%d]",c[1],1);
. Q" Z' m% g2 y* `7 |  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);
+ X7 }& B& W% W" Z  }
printf("\n正确吗？：(y/n)");
% w& |  A' W$ H6 j: L: r3 Y! h getchar();
confirm=getchar();
if (confirm=='y') break;
  Z# r: S: R5 `9 m, A2 b else if(confirm=='n') continue;
}</P>
<> printf("输入约束方程组个数：");
scanf("%d",&amp;num_st);
" J4 g- e: c. r$ D" E for(i=1;i&lt;=num_st;i++) {
  printf("st.%d:\n",i);
* n% G+ Y/ l# f- ^8 \  while(1) {
2 p7 \" t* C; q! T/ o5 O: {   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
   scanf("%d",&amp;re_st);
9 {8 `3 O- E6 [1 A+ W( ~7 M* |   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
   else break;
  }
  printf("输入技术系数：\n");
  for(j=1;j&lt;=num_x;j++) {  
  Z5 r9 e& `; h) w   printf("a%d=",j);
3 _* z% P/ u- O/ [) I8 Q) B   scanf("%lf",&amp;a[j]);
/ ^6 F+ |8 v$ V+ a9 d5 I  }
  printf("输入资源拥有量：\nb%d=",i);
# N1 s3 L2 N% g: Y/ @" Z! F6 n8 c  scanf("%lf",&amp;b);
( n  I+ q* g3 b  
  printf("st.%i:\n",i);
  printf("%lg x[%d]",a[1],1);
  for(j=2;j&lt;=num_x;j++) {
   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
. p& j9 C* c9 R2 C1 S   else printf(" %lg x[%d]",a[j],j);
$ q8 I8 O& g% [! y  }
) @1 M0 q6 _0 n; L$ `* s5 {# r" z  switch(re_st) {
$ Y2 u. w) t2 q: |' q   case 1: printf(" == %lg",b); break;
7 q  C5 z  |9 c6 N   case 2: printf(" &gt;= %lg",b); break;
   case 3: printf(" &lt;= %lg",b); break;
  }</P>
& v6 H1 [+ z. u: b3 F: z<>  while(1) {
6 c7 ~  h  F5 j* M0 \   printf("\n正确吗？(y/n)");
: h% q" s8 ~6 }* W6 F& J6 A- s   getchar();
# @8 g, m3 N9 I2 q) k   confirm=getchar();
   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
  }
}</P>
<>//显示输入的方程组
printf("\n原问题为：\n\n");
4 I5 ^; X$ e5 N/ v  H: u& M if(ma_mi==1) printf("max z= ");
# K6 a3 `% s' [+ x& \' A; B else printf("min z= ");
printf("%lg x[%d]",c[1],1);
9 @9 }" G8 J! u! v& f9 } for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
* B8 j$ D3 p$ u- g& x+ @: R! {  else if(c==-1) printf(" - x[%d]",i);
+ E$ u. m# \) b% c$ h: K  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
9 b2 p3 |; _% [$ ~; y) L, a) |  else printf(" %lg x[%d]",c,i);
}</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
* X8 b9 h6 ~( {6 S5 E( c- a) r4 Z   if(a[j]!=0) {
  ~" n* q+ `' @% X; j/ ~  I/ H    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
2 `: Z6 ?6 _+ s$ G! j8 r$ I    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);
4 e( D0 N) ~1 {$ T, P% ]- K/ {    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
: F' m( |0 B3 O, l    else printf(" %lg x[%d]",a[j],j);
1 c1 }5 |$ O6 B; ~3 I    k=1;
   }
  }
  switch(re_st) {
   case 1:
/ s: v* R, Q) G7 [    printf(" == %lg\n",b);
0 {+ T" w, a' b- c    break;
   case 2:
    printf(" &gt;= %lg\n",b);
; q' @7 G/ ^: t5 H" k    break;
% f* S/ {& D9 M  T   case 3:
/ g. ]9 K9 O' ~    printf(" &lt;= %lg\n",b);
" m0 m9 f+ u  h    break;
  }
, ?6 _( I2 Y0 H' V; O }
printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
' T3 V2 G! b% ]7 y3 V<> tnum_x=num_x;
- _. X' L: ~! l, I7 P$ Q6 _# r for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
% _: l7 O2 {  i2 ~1 @  case 1:
  case 3:
5 w) r! F) s( [( k- L   num_x+=1;
   break;
  case 2:
   num_x+=2;
+ W) `# N1 W, A7 v( i   break;
  R9 k& H1 z; [" S  }
$ n( @0 J, d$ @) C. \8 D6 c& R* O }</P>
! `. M& L8 u0 Z# K8 j<>//化为标准形式
6 o3 N0 }( K6 _3 X. M) O' Y if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
8 K+ q9 c( W) o7 i5 l6 O* B- T2 X for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
. P1 |+ e. a0 r' g   case 1:
6 e. m3 {6 N% l4 s+ ]6 G    num_addv++;
" }: E, K% U' A! C3 `- {    num_ba++;
3 T5 x0 p, A& @, a1 v% n* K    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 _! c/ t2 U1 T+ j' @) L- t7 V$ L$ {     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
; j+ D; H9 x. d     else a[j]=0;
    break;
! V# x/ X6 b- n: u& o8 a- U   case 2:
    num_addv++;
: J) F1 q* m. M2 x4 O2 P    c[tnum_x+num_addv]=0;
; _7 }* X5 ~5 \4 q! ^    num_addv++;
    num_ba++;
8 K* G& T" n; ~7 v) y2 i& ?* Z+ k    num_ar++;
    c[tnum_x+num_addv]=0;
    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
3 r8 x7 M1 Y! R) \) h    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;
. N  O6 A$ y9 Z$ I0 U4 R     else a[j]=0;
6 x% |; x; c* {' S3 n    break;
* d. e) Z2 A& D2 U; a( E2 m# M   case 3:
    num_addv++;
7 w4 @# P9 Z0 B. K  B    num_ba++;
  t" H( |0 h" E/ `" c    c[tnum_x+num_addv]=0;
    base[num_ba]=tnum_x+num_addv;
. I2 {$ E1 g6 C8 H2 L    for(j=tnum_x+1;j&lt;=num_x;j++)
     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
* G1 Q' k7 p5 v4 \     else a[j]=0;
  M( C# v9 B& L    break;
  }//switch
  s3 a) W: X/ ^) u/ E9 a% H& E/ [ }//增加松弛变量、剩余变量、人工变量、确定基变量</P>
<>//显示标准化后的方程组
printf("\n化为标准形式后：\n\n");
if(ma_mi==1) printf("max z= ");
: ]/ z1 ]( H, g9 t) c' V/ W( D+ K else printf("max z'= ");
, x& n; Z% j1 V1 T+ W9 u4 O5 @ printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
; Y5 g& j% I6 G# _0 `. @  k=0;
  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);
  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);
: H. p4 X" u$ ?2 y. x% t" k) D }</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
2 }. S( H, U/ D4 g5 c0 F  k=0;
  for(j=1;j&lt;=num_x;j++) {
8 C% P. d) I+ _: |) _6 G   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);
" V' R+ i6 z  f, W: J' o, H    else printf(" %lg x[%d]",a[j],j);
' U8 m% ^8 s; L6 I+ ^    k=1;
   }
  }
+ d- d  l( N. |/ n6 G5 Z  printf(" == %lg\n",b);
}
( g6 V7 L+ i/ ?" L, k printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
<>void iterative() {
int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
/ P; }; D- L  X int base_elem;
int base_out,base_in;
# S. q9 y# y# v+ K7 O" x& B double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
# r( I0 Y$ u& V& E9 X+ u' c" L) V3 I<> printf("\n\n第%d次迭代：\n\n",stop);
2 c1 Y# L0 o3 l5 ?5 g! x  ] 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) {
: r% M5 x$ J8 K# z/ E- Y! F   case 1:
, x6 X- W8 m* d; \/ Z7 y9 {5 C    for(j=1;j&lt;=num_x;j++)
* h" v7 i1 [% \5 T    {
  z& k% U$ Z% G7 F- G3 J     printf("a[%d][%d]=%lg\t",i,j,a[j]);
9 f) r7 t1 ]6 W4 ]$ h    }
    printf("\n");
# X4 S8 V5 T% `) V: a    break;
; z" g6 b4 r4 b( _   case 2:
/ M2 l( v, Y3 {2 f6 a    for(j=1;j&lt;=num_x;j++) {
$ r( j) F5 A9 C$ F# q% B0 ]* X     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]);
    }
- p2 \2 b! w, G3 r2 W. Q    printf("\n");
: K+ k9 N$ e4 }    break;
' L2 [+ t4 V- @- Y) D  }
. `9 f* Y8 g6 ^1 W3 n }
) F- D4 }. v1 T//求检验数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];
. I, O: E& g  b" p  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
  switch(step) {
   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;
2 E5 E5 ?% J8 i5 M" ?0 u    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
    break;
  }
  q/ a5 {# \5 S3 \/ N }
putchar('\n');
# @0 _' c/ `/ s- [3 L//检验检验数sigma是否全小于等于0
6 D3 R+ ]' M! R) ^! {7 A9 k k=0;
# ^- L2 |0 X+ L! E for(i=1;i&lt;=num_x;i++) {
( ^+ k. h+ f, L/ T7 C  if(sigma&gt;0)
7 N  `& S* w+ C( h- c2 r  k=1;
. |' X6 |: F* y. m- o& d! s+ y }
if(k==0) {
: R1 h. w' x7 |' a# ?& a  //sigma是全小于等于0时，检查是否为无穷多最优解
  for(i=1;i&lt;=num_x;i++) {
& K# J$ d2 p6 p, g$ R# ^: t4 b   k_f=k_a=0;
9 ^" S3 p: p4 h1 ?2 b   for(j=1;j&lt;=num_ar;j++)
. n4 F) i2 j. e! v7 b( ^   if(i==arti[j]) k_a=1;
   if(sigma==0&amp;&amp;k_a!=1) {
' s' d. _' p, A' v5 {; t5 [  N    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
& X9 l; q, g* M' j8 E   if(k_f==0) {status=-1; return;}
( R7 e. Z* B. j   }
& ?8 _% w3 \3 _9 x$ ^! c  }
  status=1;
  return;
}
//检查是否为无界解
7 G( z& T4 f+ } for(i=1;i&lt;=num_x;i++) {
  k_f=0;
  if(sigma&gt;0) {
   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
   if(k_f!=1) {status=0; return;}
- \% z& d9 g5 o  }
}</P>
+ f& h$ a" c: m& i% d$ m5 U<>//确定换入变量
for(i=1;i&lt;=num_x;i++) {
  k=0;
! U2 u1 W3 u% q9 @  d  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;
}//temp赋初值
! m! T( g9 L6 H$ V  `, z for(i=1;i&lt;=num_x;i++) {
  G2 G- O4 e0 t. k  k=0;
' A& x  E# v2 S$ P' u0 e  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  if(k==0)
  @* V$ I+ `: h: t4 d4 |! v   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
    base_in=i;
    temp=sigma;
   }
$ a) ~' k" n: {9 O0 Z- ~ }</P>
+ l$ G  I' l7 M: v+ T5 t# w<>//确定换出变量
for(i=1;i&lt;=num_st;i++)
  if(a[base_in]&gt;0) {
1 ~5 p; F8 ?9 u! I, v1 `$ ]   temp=b/a[base_in]+1;
: p' U! ~$ Y* O" G4 a# {   break;
  }//temp赋初值
* o8 H# J. T0 G9 Z- _- ?0 N for(i=1;i&lt;=num_st;i++) {
+ e& X! q2 b1 h1 R1 R  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
. \6 d2 ]; b  i$ C1 e   for(j=1;j&lt;=num_ar;j++)
3 N8 a3 L* L$ p  S7 {    if(base==arti[j]) {
1 c: v8 F: U. n9 a. _     base_out=base;
     base_elem=i;
: M3 s: P6 w, J& a3 T: p     temp=b/a[base_in];
- x) t# ~; i. \3 c/ l, R* g     break;
    }
  }//人工变量优先换出
  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
   base_out=base;
9 `0 k5 j5 M4 `" c5 n1 q. n   base_elem=i;
( C# D1 x# T9 i1 X! g   temp=b/a[base_in];
; P0 K. E$ A$ r6 s( B  }
}</P>
<> printf(" 基变量：");
& H0 [, A5 z* ?: z for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
1 H! n+ V: Y6 _ printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
4 U6 X* t5 {- ]# C# k//基变量变换，进行新方程初始化后迭代
5 h9 R4 z" _6 F for(i=1;i&lt;=num_st;i++) {
' H* K) u9 ]8 J5 D' x0 D+ C  if(base==base_out) base=base_in;
}
+ j9 g+ C* P% U/ I. W# H+ l: Z//初始化主元素行系数
value_be=a[base_elem][base_in];
: u# }) F; T4 Q* ^. F+ ~  |3 D/ } b[base_elem]/=value_be;
! ?) T9 V5 G0 l7 G5 L% g3 Z for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
<> for(i=1;i&lt;=num_st;i++) {
; e0 B" H( @8 i: k, K/ s1 m  if(i!=base_elem) {
   b-=b[base_elem]*a[base_in];
   value_be=a[base_in];
   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
) o* K6 I$ D2 H) L4 t  }
$ z' _2 t% f7 y/ S' L1 p }
stop++;
if(stop&gt;STP) {status=-2; return;}
iterative();
4 |* E; N9 F* ]. _}</P>
<>void output() {
$ h: o& M+ g$ S int i,j;
5 m' h  w2 ?; d9 r$ o- _4 R double X[MAX];
% |, R- I2 s/ o- q: g printf("\n结果如下：\n");
printf("\nX=(");
! H; p. s. v5 G' h 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;
- v  g/ _; {" ~+ E  printf("%lg ",X);
}
printf(")");
9 J6 M# ^) [: X for(i=1;i&lt;=num_x;i++) max+=c*X;
if(ma_mi==1) printf("\nMax z= %lf\n",max);
else printf("\nMin z= %lf\n",-max);
. j6 f) J2 F+ d, F7 D}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了