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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
5 p# g2 I3 ~3 O2 E" G9 d# J  e#include&lt;math.h&gt;
#define X 5
. L% T# S: F) h5 l#define Y 7
; }9 t7 |/ v4 tvoid xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
{
7 ~5 L1 h+ p/ c8 ?4 K3 Q  int j;
$ G$ S& R6 j2 S2 E9 o- E9 p* `  *c=0;
  for(j=1;j&lt;=*is;j++)
7 B: O  ~* x# n* C( A8 Q$ A    {
8 z9 o+ u% Q0 s      if((*a)[*ir][j]-*c&gt;0)
- G9 [3 G% w5 s* q5 L    {
      *c=(*a)[*ir][j];
/ i- O% ?. e6 ]; X7 d      *j0=j;
    }
     }
% J7 I5 J1 H9 }' |}</P><>      /***************** 参数说明 **********************/
5 S1 h: p2 m8 O7 U$ ?      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
, g( g% A. ?3 I+ j( ]: o# {      /*   (*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],
' l' T( P' z( a3 h# O       int(*k)[],float(*x)[])
9 F8 C7 H0 k  l" v3 p6 i{
  int i,m1,mn,j0,i0,j;
  float c,g;
; d; `: P! f9 ]9 i+ O1 \. O/ |  m1=m+1;
+ h8 A+ ~! q6 g# F; s* J/ K4 A  mn=m+n;
2 F* `. k4 x( S  for(i=1;i&lt;=m;i++)
  U3 Q, L7 P0 \; b/ l    (*k)=n+i;
0 q1 O& G' k# v/ @. Y2 ]0 _" |( Kleap1:
  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  else
* C0 ]3 e! Z, Z, {- f0 g     {
, G( S5 ^/ B) {0 |2 C. k' mleap2:   if(l1-50==0)
      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
3 I9 }$ f% F5 ]& n       else
      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
" H3 Q( C  ?$ G8 @! L6 G     }
" g: U! Q- c# S* w: x7 g   c-=1e-8;
# Y  @0 M3 A% W1 n- b5 V   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
; W  F; \$ g8 ?& r      {
2 b5 R, V: Y# ~2 C    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
    return(-1);
      }
   if((c&lt;0)&amp;&amp;(l1==1))
- E3 f  v8 g% t. w6 }       {
& b$ M  P/ `0 n, X1 e2 x/ T2 K/ d5 q     l1=50;
     goto leap2;
. ]& ^8 {1 H2 B( }       }
    if(c&lt;=0)
5 k8 U+ Q# |# M2 d0 J       {
      for(i=1;i&lt;=m;i++)
$ G; E4 Q8 E, O4 q9 }; q: o         (*x)=(*a)[mn1];
# x7 P7 E. F' B( ], B         printf("\n\t********Optimal solution********\n");
8 n& N7 s4 o$ C% `) H      for(i=1;i&lt;=m;i++)
; T, I0 S) D* ^3 S$ _         printf("\ni=%d,x=%f\n",(*k),(*x));
      printf("\nF=%f\n",(*a)[m2][mn1]);
      return 0;
! q4 q/ I' s3 L# F% C  g    }
     c=1e8;
" k+ ?( W  m6 o% s     for(i=1;i&lt;=m;i++)
: |! h0 _2 |7 T" j       {
     if((*a)[j0]&gt;1e-8)
        {
$ \: Q% z# `5 G  i$ ~. z2 E, ^          g=(*a)[mn1]/(*a)[j0];
4 L$ F7 s0 M: Z; k  S( q          if(g-c&lt;0);
' B; }8 v( A9 r; T" y3 F7 m7 i         {
           c=g;
           i0=i;
          }
- q: ?. X+ H. J9 i. A         }
    }
      if(c==1e+8)
7 a* L' h6 |# h, Z5 s2 l* f9 Y3 u    {
      printf("\n\t*******Lp no solution********\n");
      return -3;
+ @) d5 n% Q+ [+ [8 g     }
       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
      {
        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
+ l( |; o! f0 O& s3 s6 p  U: J           continue;
         g=(*a)[i0][j]/(*a)[i0][j0];
% k1 A3 v2 W- x9 j( W           (*a)[i0][j]=g;
: ?" l8 L5 y# U# E' O% J+ M) t. ?& g9 }         for(i=1;i&lt;=m2;i++)
           {
         if((i==i0)||(!(l1==1)&amp;(i==m1)))
            continue;
9 {& w) p: ^8 ?: @8 d' _         (*a)[j]=(*a)[j]-(*a)[j0]*g;
           }
       }
       for(i=1;i&lt;=m2;i++)
. ?5 G# h+ v0 |% |* q& s      (*a)[j0]=0;
       (*a)[i0][j0]=1;
( ~. x: V) a$ b       goto leap1;
3 Y, D. s4 y( [" W}  
1 W' |% t7 z+ _1 W1 A4 J2 Hmain()
  ^2 c$ C* u2 o( C! p# M{
  float a[5][7]={{0,0,0, 0,  0,0,0},
                 {0,3,-4,3,  1,0,12},
                 {0,3,0, 6,  0,1,12},
- x7 @& u$ s" A, m2 S                 {0,0,0, 0,  0,0,0},
' z9 K) I' [* @2 j                 {0,69,0,144,0,0,300}};
  int k[3];
. U3 y/ t# b% V5 F/ X2 `  float x[3];
  clrscr();
) q. W+ C6 U6 m1 V5 W# l  z  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
}
</P>

plgatc

<>/*************************************************************************
                    单纯型法解线性规划问题（两阶段法）        

' a+ }6 X$ |! K; Y3 }6 [ 编程环境：VC++6.0     
方程组输入说明：
, l# P( Y0 Y. U- r 变量非负，按提示输入相关参数。
*************************************************************************/
% e4 ~+ ]: q2 x1 K5 R2 y5 a% F#include &lt;stdio.h&gt;
: {2 c# j3 \( B" b% `#include &lt;stdlib.h&gt;
2 ?1 S: v" s$ l  [# b9 @4 K0 L#define MAX 100
#define STP 100</P>
<>int stop=1; //迭代记数变量
int status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
! c- R8 b  M  ~) ^, j+ k% ^int step=1; //目前阶段</P>
<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
" R2 H; ?4 h) \* D% x% _1 R; G6 @8 Eint num_x; //变量个数
int num_st; //约束方程数
int num_ar=0; //人工变量个数
, p4 [- u2 F! a4 _3 I, C5 ^9 b' Iint arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
4 w$ o: \! X( Jint ma_mi; //1为求最大值，2为求最小值</P>
$ s/ U& `3 n/ |8 v. f<>void create(); //建立方程组
void iterative(); //单纯型法迭代
void output(); //输出结果
void banner(); //打印程序标题
" `0 g7 \" h7 h+ Hvoid exchange(); //交换两阶段价值系数
void show(); //输出方程组</P>
<>void main() {
; m, d' U2 ], _8 W! D5 p* a int i,j,k;
% }  c3 O5 c, [" a: F, a3 e banner();
- I7 d1 ~! V9 {# \+ M. y$ h" m+ [. [ create();
: I& c* k& R5 j' J* `' I$ [4 E2 J# D //保存原价值系数，转换为第一阶段价值系数
for(i=1;i&lt;=num_x;i++) {
  k=0;
# M& o' J, v  [- c1 a! |  d1 t  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
  if(k==1) temp_c=-1;
  else temp_c=0;
}
% \2 R; w7 k7 y( D, {, _* N exchange(c,temp_c);</P>
) \( m9 _7 l6 Y1 A+ [3 O<> printf("\n\n第一阶段问题为：\n\n");
2 n7 y- y  g$ c% C; i8 t show();
step++;
printf("\n\n按回车开始第一阶段迭代");
getchar();
getchar();
iterative();
/ X# Y) Z( q% z0 E: a& x if(status==-2) {
5 N1 I5 ~) ~+ |* b3 p  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
  getchar();
& {* \/ j: x% d  exit(0);
}
output();</P>
# {, u6 H2 v0 A& e" Z<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
  puts("\n按回车结束");
0 n2 _. i. W$ E& q) M  getchar();
  exit(0);
}</P>
<> //转换为第二阶段价值系数
; Y2 C( c- a; u4 a; ` exchange(c,temp_c);
% p$ @8 ^( B$ q3 B/ Y //把人工变量列全设为0
for(i=1;i&lt;=num_ar;i++) {
% @  y0 |" d& n! h+ s) ~  c[arti]=0;
1 N5 c( ?9 u) s4 r  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
}</P>
$ B, ?; O9 c! C5 l" \0 R<> puts("\n\n第二阶段问题为：\n\n");
show();
puts("\n\n按回车开始第二阶段迭代");
) k  o& o, ^5 ]" g% _! T9 e& L getchar();
  [' T! ]: t: m! d( L, { iterative();
switch(status) {
case 1:
  output();
  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
7 U  }) p0 W- ^) c2 z  getchar();
  exit(0);
: W1 m# G' C  T0 C case 0:
+ a; G5 a( K* Z) ?, W6 a3 x  puts("\n\n原问题为无界解。\n");
  puts("\n按回车结束");
: l+ V* \8 ?5 O4 r) q! _, y4 Z  getchar();
  |6 K8 u6 Q0 p; w! n  exit(0);
case -1:
  output();
9 X# v5 Z' g2 H4 T. I" I9 I, A* ]  puts("\n\n原问题有无穷多最优解。\n");
& C; K5 g6 U" m2 _# y9 U- v: h  puts("\n按回车结束");
* Z# m2 a! S7 E- \/ R: G  getchar();
  exit(0);
4 b' n- X7 n7 f! A/ V4 W case -2:
  puts("迭代超过限制次数强行终止！\n");
1 i) I* P4 B; Y' c. X2 E3 q) o; j' R  puts("\n按回车结束");
  getchar();
4 L5 S# g% |/ j9 ?  exit(0);
# o% F& P8 ]; J! q* |- ^ }//switch
  b. H9 `2 e" z& j3 q}</P>
% `3 d" P0 X7 v3 a- }<>void banner() {
' z5 F3 m2 Q- ` printf("\t\t****************************************\n");
printf("\t\t         单纯型法解线性规划问题\n");
0 l2 K5 ]. m" o* h printf("\t\t                         作者：Thunder\n");
printf("\t\t****************************************\n");
printf("\n");
0 C( [+ w2 F1 N/ O}</P>
$ {3 M3 r. Z; k3 r. ]<>void show() {
8 Q9 M8 ^8 k* O2 \9 a//对方程组以自然的格式输出，系数为零的x不显示
% T, H- A$ |: W//为1的不显示系数1，-1系数只显示负号
# z: X# W9 P2 y. n; ]0 a& V% k int i,j,k;
switch(step) {
6 J, L# g: O$ [  `6 A# B case 1:
& p, Y1 I- f1 q6 A, D  printf("min z= ");
  printf("x[%d]",arti[1]);
; R$ f& t3 F/ }1 M2 J% V5 X  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
: N/ c6 r; W" T2 t' E8 g1 T7 M* b; { case 2:
+ b2 V: Q  o) a 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);
# c8 \+ \- {0 g* Q- Y8 {  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
# [' E: R& u: N  else printf(" %lg x[%d]",c,i);
+ o2 r; u: H! ~6 j2 i# ] }
' j5 q0 Q5 y7 q/ J- [* }' w break;
}</P>
<> printf("\nst:\n");
: \! m) P) K* ]; W2 O3 u$ ?3 |/ \ for(i=1;i&lt;=num_st;i++) {
/ ^* L8 ?" z, n+ V9 v4 Q  k=0;
, @% S" @; s( R! ~  for(j=1;j&lt;=num_x;j++) {
2 l* w3 Q1 e% P, \) ~2 r6 D   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);
2 |6 K4 Y# E8 Z1 n+ j5 C  X4 ]    k=1;
   }
9 F+ A  i$ y; G1 |  B& m; X' s  }
4 F9 K8 j+ z- P: u$ y  printf(" == %lg\n",b);
9 {4 L, [) z3 h }
$ _9 K" J+ L6 S& b; q% M$ x  H printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
3 d% A+ ]3 ?# K) X<>void exchange() {
: R! V5 e, t" [  X6 U int i;
double temp[MAX];
; s5 W" U: Q+ n) `7 ]6 F4 q for(i=1;i&lt;=num_x;i++) {
  temp=temp_c;
: M" ]$ `0 a% \; V: n0 k; a" \  temp_c=c;
3 b: s% C' q& p! W% K( U6 v+ S+ e% G9 @8 }  c=temp;
" J3 `7 [2 h# b4 h5 E  l/ o" X }
: m( D' O, d' l/ m) c: W; V}</P>
<>void create() {
2 O1 ?4 t, l' w  m//输入方程组系数，每个方程输完后回显确认
int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
char confirm;
# i3 b5 h0 U( G  l5 P7 j
while(1) {
7 i# O3 w) p0 x" O- \: V& H! g  printf("请选择：1、求最大值，2、求最小值：(1/2)");
5 N# j) C0 y* J6 j  scanf("%d",&amp;ma_mi);
# ?" F& M5 B8 `* P+ l. c  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
! R! l; K. D* f9 }5 R9 i5 R& I  else break;
}

while(1) {
  printf("指定变量个数：");
  scanf("%d",&amp;num_x);
# t9 E" Y( i. b6 Y! g  printf("输入价值系数c1-c%d:\n",num_x);
  for(i=1;i&lt;=num_x;i++) {   
' v$ e7 _# L! h3 l5 L0 B$ }$ }- d   printf("c%d=",i);
' [8 G5 z; h& J3 ?4 h   scanf("%lf",&amp;c);
  }
8 `5 S7 V5 {& F0 j0 X$ E  if(ma_mi==1) printf("max z= ");
  else printf("min z= ");
( H  |6 S5 G/ a  printf("%lg x[%d]",c[1],1);
# b; I3 d3 f6 c$ }) l  for(i=2;i&lt;=num_x;i++) {
6 `( s: g: X5 [# S8 O9 v% U   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
   else printf("  %lg x[%d]",c,i);
$ d/ Y2 K( A- X0 ]0 b7 w1 n  }
6 v) t0 l/ b) z9 Y  @. u printf("\n正确吗？：(y/n)");
2 B, h& u) y) _$ e9 U' g getchar();
confirm=getchar();
: E  ^$ G0 ^7 l  w4 S0 ~ if (confirm=='y') break;
else if(confirm=='n') continue;
+ z* W. Y1 r3 F6 T( W' \2 u }</P>
<> printf("输入约束方程组个数：");
scanf("%d",&amp;num_st);
for(i=1;i&lt;=num_st;i++) {
  printf("st.%d:\n",i);
  while(1) {
   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("输入错误，请重新选择。");
# Q3 q* F) N! E: }   else break;
- I% P) Q- n* `% W! u! w: v  }
  printf("输入技术系数：\n");
- D. N. i; h2 \4 M1 v6 a1 F: `  for(j=1;j&lt;=num_x;j++) {  
4 B1 j+ v1 o0 X% C7 w; t   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
  }
  printf("输入资源拥有量：\nb%d=",i);
  scanf("%lf",&amp;b);
0 r2 K+ Z. B  m% r& l% i; X2 p4 j  
4 S6 T1 l5 v( _& e3 v( R, N( S. m  printf("st.%i:\n",i);
  printf("%lg x[%d]",a[1],1);
  for(j=2;j&lt;=num_x;j++) {
/ X5 f8 Z0 S" b   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
   else printf(" %lg x[%d]",a[j],j);
  }
  switch(re_st) {
4 L& S2 R9 T% W* s% q$ r   case 1: printf(" == %lg",b); break;
# O% W! A7 z/ W   case 2: printf(" &gt;= %lg",b); break;
8 O1 x& t; X* N: B% M   case 3: printf(" &lt;= %lg",b); break;
; h$ I/ s1 \/ x, Q7 m6 n7 ]  }</P>
<>  while(1) {
# d! q4 C: W) w5 H   printf("\n正确吗？(y/n)");
   getchar();
   confirm=getchar();
9 C: v5 k3 _! b! \- {" ~   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
/ D. J, y- R- p: ?! `! G. D  }
2 g, s! ]% j! M }</P>
3 ^( A8 }& R, [<>//显示输入的方程组
+ A; d6 V$ j3 b printf("\n原问题为：\n\n");
if(ma_mi==1) printf("max z= ");
else printf("min z= ");
: s) l2 I7 ?& T4 e) b4 o" o printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
1 K; R' o8 i/ y  O/ _  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);
' F* o' _' D8 @$ o& m  else printf(" %lg x[%d]",c,i);
}</P>
<> printf("\nst:\n");
4 i$ H/ w, v9 J5 L% X. q for(i=1;i&lt;=num_st;i++) {
# y$ Z3 Z' e/ f  g5 B' u  k=0;
! i" ]0 h% K: l  for(j=1;j&lt;=num_x;j++) {
- Y' n3 M: \/ l& h, F, m# b   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);
. O$ @5 c; ]: t/ U# C    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;
   }
" a8 l2 ~6 \, x. o  }
+ v' V3 q6 B( A  switch(re_st) {
   case 1:
    printf(" == %lg\n",b);
    break;
1 ^5 R7 c- `& J9 R6 m3 l: q$ K0 f   case 2:
    printf(" &gt;= %lg\n",b);
    break;
   case 3:
: _6 J4 D: [* u! X    printf(" &lt;= %lg\n",b);
0 V4 o! D( u1 X    break;
  }
6 P9 {0 \2 `, b. M5 e* H& s }
printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
<> tnum_x=num_x;
: o$ k" R' J* G  E  ^ for(i=1;i&lt;=num_st;i++) {
( ]7 X$ A- T7 D/ @( m3 J. C+ {  switch(re_st) {
  case 1:
  case 3:
2 i( L5 {9 L* }. f6 y( A9 M, b   num_x+=1;
/ Y( U6 C9 G0 Z! L4 I6 s   break;
  [. a4 e, \9 @  case 2:
( H: H7 I  Z4 m! j6 e7 _   num_x+=2;
( T- x% u" l# B* j; |, D   break;
  }
" f9 n* d+ k8 o2 S- T: A& s: Y) g1 L }</P>
# N: d& O4 s( Y" \<>//化为标准形式
if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
+ X0 ?- W0 N# k: a for(i=1;i&lt;=num_st;i++) {
9 C9 h8 c1 t9 f( Z; v$ n6 S9 z  switch(re_st) {
2 p6 `( W! t. A- j$ T   case 1:
    num_addv++;
; a  r$ P, E& b  O, }+ \- ?    num_ba++;
# X0 s, W/ ?) E: d% T* C    num_ar++;
; ~4 [) ^! G5 I$ V5 M    c[tnum_x+num_addv]=0;
( I+ {" q' N. ^' n5 B    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;
    break;
2 E" J, p8 J' a* f: \   case 2:
  x( v0 g2 a! T9 `8 v    num_addv++;
5 w5 j  i3 P% J9 j% C, U7 a& A5 v- o    c[tnum_x+num_addv]=0;
5 r; s' Z$ B8 O! i    num_addv++;
    num_ba++;
    num_ar++;
8 C: D+ `+ K) F! ]3 K: c- s    c[tnum_x+num_addv]=0;
    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
6 c# G0 Y1 A2 @! `1 ^    for(j=tnum_x+1;j&lt;=num_x;j++)
! z9 I, `* N4 w+ x     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;
    break;
   case 3:
! e) b; T" h5 z+ O4 g; G    num_addv++;
    num_ba++;
) J1 E" }8 T; A7 V$ L6 u7 K    c[tnum_x+num_addv]=0;
    base[num_ba]=tnum_x+num_addv;
    for(j=tnum_x+1;j&lt;=num_x;j++)
3 ]% w# k) H! C1 i     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
! H4 Y6 j& f: F. x( W% i3 j    break;
( j% O$ Z- w# \6 t1 Z3 f  }//switch
4 F% v  K8 w" z6 k }//增加松弛变量、剩余变量、人工变量、确定基变量</P>
& f/ ]; W8 w% @' T1 X<>//显示标准化后的方程组
: q4 q% M' n0 l/ p4 n0 v9 U printf("\n化为标准形式后：\n\n");
if(ma_mi==1) printf("max z= ");
# i. ?4 i* E& s/ ?" |; e1 T else printf("max z'= ");
printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
' ]/ N  b0 v3 j6 ~: l, Z" c. ?  k=0;
  for(j=1;j&lt;=num_ar;j++)
  if(i==arti[j]) k=1;
/ l' |/ b* H! ?; d  L  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);
# @. u6 |1 ?5 o7 U# Z- e; K& I  else printf(" %lg x[%d]",c,i);
}</P>
3 l$ U1 s2 F+ X; `) ]* o8 X<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
! J0 `7 v, g3 P! I# Q' A/ A* h1 q  k=0;
+ ^2 |9 ]0 j) K7 k. `% ~8 z  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);
$ `7 w- {- D0 q* `% ~1 ^( o! {    else if(a[j]==-1) printf(" - x[%d]",j);
; K1 E* ]) k6 R0 }* 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);
$ Q0 I3 `  m4 W1 p$ P    k=1;
  q0 h: |0 ]( O- g4 V; J' b   }
  }
  printf(" == %lg\n",b);
' D' W! v8 d/ E/ V- J& l }
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[]里的搜索结果
8 m2 W4 z2 N6 ^4 P4 Q9 W+ a7 G4 O int base_elem;
0 \: F" C  F' b" J int base_out,base_in;
double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
<> printf("\n\n第%d次迭代：\n\n",stop);
2 v0 ~' b; Q( y 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) {
& b& r. }7 P( H/ k# Q2 x" w   case 1:
    for(j=1;j&lt;=num_x;j++)
2 h" E, S9 Y# D% P, F% ^    {
( @* j1 d+ b) m9 C* Y     printf("a[%d][%d]=%lg\t",i,j,a[j]);
0 ~: g$ g# y" Z! P    }
    printf("\n");
5 G" `, f: `& B. B% Q    break;
   case 2:
9 b8 E% j  ]2 ~6 }) k& ^5 j% g    for(j=1;j&lt;=num_x;j++) {
5 w& O' |6 \1 Z* Y0 G8 Z% l3 Z" i. h     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]);
    }
8 e3 l- D1 ^" Y0 s) _2 p    printf("\n");
5 S  Z2 T+ b) r: z    break;
. m/ @0 r" X! S0 R/ r% {6 q/ w  }
1 X" [5 N/ k; l* ]" P3 B# v9 z& Q$ w }
7 j, Z9 ~; X, c+ O6 \: V3 T//求检验数sigma
for(i=1;i&lt;=num_x;i++) {
, ]3 h3 B) g2 D5 ]' w, T  sigma=c;
- v( T  z/ j: ]" L* j+ ^3 Y1 e  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;
  switch(step) {
" @1 L" a: R; f8 A! k/ _: U   case 1:
7 U8 M2 A, m( e. _" J$ k( Y1 s8 }4 ?    printf("sigma[%d]=%lg\t",i,sigma);
7 r% e+ R( i  M$ \    break;
! Y0 h" Z' H% _" t   case 2:
" a/ b4 h. _7 d$ w    k_a=0;
    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
3 r1 q' T0 m2 f    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
" Y7 U4 ^  H- Z    break;
3 N' r8 F5 s1 M' t" A. t  }
! ~+ p  S. z0 e+ J0 f# M2 B }
putchar('\n');
$ G, C- q8 ^& |0 h0 [  G7 e//检验检验数sigma是否全小于等于0
k=0;
for(i=1;i&lt;=num_x;i++) {
  if(sigma&gt;0)
  k=1;
! W  v8 e, @: S2 Q7 u }
if(k==0) {
6 B0 M6 g" S  l; Z  r5 ^8 ]# U  //sigma是全小于等于0时，检查是否为无穷多最优解
  for(i=1;i&lt;=num_x;i++) {
$ I7 s2 |( P( c7 Y& o' \8 m   k_f=k_a=0;
% U3 r3 X; s  ~: L& Z- c7 S# N   for(j=1;j&lt;=num_ar;j++)
" U' f- }7 j+ Y  J6 P   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;
, o, O$ D; C9 B. M) Q! Y# e   if(k_f==0) {status=-1; return;}
, A: e7 N8 F2 U! x2 x( b0 h   }
  }
  status=1;
1 m3 `) S9 v$ e( j  return;
}
//检查是否为无界解
for(i=1;i&lt;=num_x;i++) {
  k_f=0;
, d+ `. _' _+ m8 l7 B9 ]# U7 F  if(sigma&gt;0) {
" p! o6 c' ?8 h- L9 x   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
   if(k_f!=1) {status=0; return;}
, i2 u2 S" [, N/ \0 T! `  }
: Z- {; u6 B* n- e" ?0 m }</P>
: ^: H9 K6 L  V  m% z<>//确定换入变量
for(i=1;i&lt;=num_x;i++) {
  k=0;
% m; `8 p1 ]6 X1 g  d3 e; W6 t  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
( _: y6 W# P# a& Y1 s( X, ~  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
}//temp赋初值
for(i=1;i&lt;=num_x;i++) {
0 Z7 _2 B1 b, Z5 I5 v* `7 ^  k=0;
" y! f$ [# ]0 S; |  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  if(k==0)
   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
2 I( `. T) R% P4 X    base_in=i;
    temp=sigma;
" g' J  f; x3 b4 J3 M( T# ?   }
}</P>
<>//确定换出变量
) [# l- m/ s  P% | for(i=1;i&lt;=num_st;i++)
  if(a[base_in]&gt;0) {
   temp=b/a[base_in]+1;
4 D& |& T) F( w1 I' r   break;
4 g# {' v. Q1 a  }//temp赋初值
) d# @! K4 t4 Y8 ^6 I for(i=1;i&lt;=num_st;i++) {
: [: l7 e9 F% `) _2 [1 Z  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
   for(j=1;j&lt;=num_ar;j++)
3 F' C( m9 h2 b- \2 s9 x; v    if(base==arti[j]) {
7 o4 y+ _2 R$ Y7 J+ `, z     base_out=base;
     base_elem=i;
     temp=b/a[base_in];
     break;
    }
  }//人工变量优先换出
: H. w  P2 X2 s* v; h5 n  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
   base_out=base;
   base_elem=i;
   temp=b/a[base_in];
& V6 t# T4 y, ]6 o- d; ~. s' n  }
}</P>
<> printf(" 基变量：");
! r' C) t5 X3 C! A( p/ A7 R for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
" e9 ^3 B6 B; @$ A. e printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
//基变量变换，进行新方程初始化后迭代
for(i=1;i&lt;=num_st;i++) {
) a0 N' P3 Y, _' `6 Y- d/ k5 l+ ?  if(base==base_out) base=base_in;
! p& a) w% Y0 W3 m5 O }
//初始化主元素行系数
9 P& g8 q1 Z5 k8 ^' o- ~ value_be=a[base_elem][base_in];
b[base_elem]/=value_be;
- t  c" b5 D  B for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
<> for(i=1;i&lt;=num_st;i++) {
  if(i!=base_elem) {
  P0 x1 F. R' I8 P   b-=b[base_elem]*a[base_in];
/ G+ h) ?& W( R, ?& Y1 P- Q   value_be=a[base_in];
2 b( Z& I" A4 M7 U   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  f! N7 [7 N- T* H. J% ?3 K) L7 P  }
}
stop++;
: r9 g1 S7 |. o. }3 Z, ]" ~! R0 | if(stop&gt;STP) {status=-2; return;}
iterative();
}</P>
- b9 N0 g5 @# d  Y% z<>void output() {
, |$ C/ n, X" y- j4 e7 ~ int i,j;
double X[MAX];
. e" C- b& {& ~6 l& c printf("\n结果如下：\n");
printf("\nX=(");
. C, o3 @6 t& W' a; r 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;
" ]( n1 f- U5 ^  printf("%lg ",X);
( Q/ i- h5 T2 _: ^8 y$ m" t- L5 o# d, m }
printf(")");
) T# u1 V* i$ m6 h for(i=1;i&lt;=num_x;i++) max+=c*X;
if(ma_mi==1) printf("\nMax z= %lf\n",max);
$ S. q) c7 }( z" n8 J0 C7 g else printf("\nMin z= %lf\n",-max);
}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了