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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
3 A, s$ }1 j6 j) T#include&lt;math.h&gt;
0 r1 j! r" H9 U6 `- r8 K- y6 s#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])
{
  int j;
  *c=0;
  for(j=1;j&lt;=*is;j++)
    {
      if((*a)[*ir][j]-*c&gt;0)
    {
9 a! W& v# ~# D! O0 S+ t9 v# L      *c=(*a)[*ir][j];
      *j0=j;
    }
$ o. N+ b# `$ y# G3 [, b( D     }
3 R" ^, ]) d- j/ i/ \}</P><>      /***************** 参数说明 **********************/
( G/ Y5 v0 r9 s* h% r      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
: a5 L" ?* _* [; r1 X+ A! H      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
# A; X) @" y& \/ A" \2 n- u      /********* <a href="http://happyyangxu.home.sunbo.net*******/" target="_blank" >http://happyyangxu.home.sunbo.net*******/</A>
, }. v! X' ?+ q- y# C* d/ Lint xi_sm(int m,int n,int m2,int mn1,int l1,float (*a)[X][Y],
) @3 H: q& h! {9 e       int(*k)[],float(*x)[])
9 _( c. W& Q  y. k: U4 V{
. [5 R: r* e. _0 I3 E  int i,m1,mn,j0,i0,j;
, V, c6 W; v6 o  {1 t% A  float c,g;
! ]/ E* H, j7 V7 Z. `1 A/ E; b" v  m1=m+1;
  mn=m+n;
  ]( U3 M4 X+ {) B1 U6 ]3 N" h/ g* m3 p  for(i=1;i&lt;=m;i++)
    (*k)=n+i;
leap1:
  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  N9 Y, H, Q4 [7 `" h% z  else
     {
leap2:   if(l1-50==0)
      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
       else
      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
     }
   c-=1e-8;
0 m1 `- a0 `! [/ j7 f7 r  x5 q   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
      {
' r1 [& J+ z& j6 W) \6 E    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
    return(-1);
% e% Q' ~$ l/ ~      }
   if((c&lt;0)&amp;&amp;(l1==1))
. Y' c" l0 u( S+ v$ ]       {
     l1=50;
+ o0 ]4 C9 p5 t' I     goto leap2;
7 I4 l' u1 b$ e! m       }
    if(c&lt;=0)
       {
6 e0 _" |+ l9 Q$ v. d      for(i=1;i&lt;=m;i++)
( n* D/ j. Y$ l7 j' P         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
  P$ t  r( M: t2 p9 E      for(i=1;i&lt;=m;i++)
* H1 b3 Q6 [3 S6 ~7 [# Q         printf("\ni=%d,x=%f\n",(*k),(*x));
6 w! h" a' i# z" f. \1 x      printf("\nF=%f\n",(*a)[m2][mn1]);
4 z) c5 |1 i8 }- f      return 0;
/ B7 w% ]+ Y% L$ G% F% o# ]+ D) V    }
! }$ f/ A; Z7 e  }& J9 @1 `     c=1e8;
     for(i=1;i&lt;=m;i++)
1 @. H$ X1 y5 z/ R% W: T& i1 m       {
     if((*a)[j0]&gt;1e-8)
+ \+ x: J' }+ R# ~' |  V        {
          g=(*a)[mn1]/(*a)[j0];
          if(g-c&lt;0);
; T$ d" J* q2 c         {
1 [4 w1 ?# U# z/ P           c=g;
! \( q$ G% m/ S9 a  Z           i0=i;
2 T) y. L) L. M9 A; [. w. [0 Z  b          }
# M- p/ d  P4 ]& d* Z         }
    }
      if(c==1e+8)
( Y& G# f1 m' }$ |    {
  y  B+ L5 c. H) V* u+ p  [1 I/ n0 V! S      printf("\n\t*******Lp no solution********\n");
      return -3;
9 N: r3 H) R4 W* ?+ j     }
       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
4 r) v% q* u# `+ t+ L5 ]1 W- C      {
% P7 v3 [2 C6 O7 X6 Z        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
           continue;
         g=(*a)[i0][j]/(*a)[i0][j0];
           (*a)[i0][j]=g;
1 [& g: e% c( h; v- n+ X+ L" F# U         for(i=1;i&lt;=m2;i++)
           {
) o" O" }: u  [' F# ~+ z         if((i==i0)||(!(l1==1)&amp;(i==m1)))
. F5 d( J$ H: X. t9 L            continue;
: P* s6 f. x) C* i         (*a)[j]=(*a)[j]-(*a)[j0]*g;
           }
       }
       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
  _. k, f( _5 g       (*a)[i0][j0]=1;
0 C0 p" d$ E3 T/ t- N       goto leap1;
  @) t. D. v6 r. H}  
# w& X; [0 Q: ?, V) ~' u; Imain()
{
$ B( C9 _3 u; S7 V3 R: l  float a[5][7]={{0,0,0, 0,  0,0,0},
, D8 ~0 L7 b7 v; [                 {0,3,-4,3,  1,0,12},
2 _1 a- p; `% v* x                 {0,3,0, 6,  0,1,12},
, P' n. q1 [% h8 T2 I  ~5 X                 {0,0,0, 0,  0,0,0},
                 {0,69,0,144,0,0,300}};
  int k[3];
7 |; k( A# T0 R2 i1 d! H; W7 V" b' A  float x[3];
  clrscr();
  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
}
</P>

plgatc

<>/*************************************************************************
                    单纯型法解线性规划问题（两阶段法）        
$ p/ p# g7 w  o3 Y" x8 A$ a
编程环境：VC++6.0     
方程组输入说明：
/ U, K* P( r% h* N9 e- i* B% N 变量非负，按提示输入相关参数。
*************************************************************************/
#include &lt;stdio.h&gt;
#include &lt;stdlib.h&gt;
#define MAX 100
4 j; _( z9 s  z4 Z3 @#define STP 100</P>
7 Z' t4 I, e8 H. j<>int stop=1; //迭代记数变量
int status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
int step=1; //目前阶段</P>
<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
% N' C+ u# K; y8 w0 k- k( Tint num_x; //变量个数
7 K  ?8 e* n7 P7 _& E) @( bint num_st; //约束方程数
int num_ar=0; //人工变量个数
7 q- s0 @8 P1 I. dint arti[MAX]; //人工变量下标
% `# ?: E" U9 W  S4 j$ xint base[MAX]; //基变量下标
8 c5 n& |5 B. f: ?) @int ma_mi; //1为求最大值，2为求最小值</P>
<>void create(); //建立方程组
+ L; ]/ I3 E8 b* T' [, s' S6 Q# wvoid iterative(); //单纯型法迭代
2 r2 Q2 y( Q- h6 I# Uvoid output(); //输出结果
* j! p, b9 Y) L" h, Mvoid banner(); //打印程序标题
void exchange(); //交换两阶段价值系数
& A1 B$ n/ f- F* [; X6 M1 hvoid show(); //输出方程组</P>
<>void main() {
int i,j,k;
) ~6 X2 n$ q- s4 G0 m' H' N banner();
6 d+ p7 [4 W! M) I) \- ^2 i create();
0 G, E5 c" f# \" Y! s0 s2 n7 Q //保存原价值系数，转换为第一阶段价值系数
; Q2 B( C" _7 ^' S for(i=1;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
1 h0 I: E6 w, {* `0 E3 {  if(k==1) temp_c=-1;
  else temp_c=0;
}
exchange(c,temp_c);</P>
<> printf("\n\n第一阶段问题为：\n\n");
& M/ ^7 P; o2 k- h7 M5 @8 q show();
step++;
printf("\n\n按回车开始第一阶段迭代");
% y( j% h) j2 m: T; M getchar();
0 o, @& T  j  W- i# b! H9 M getchar();
  C/ A4 w3 F2 X) I+ t1 I  F$ l iterative();
if(status==-2) {
  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
! I0 M1 V. R) K& D  getchar();
1 [2 E% x& ]3 d, X! n  exit(0);
}
output();</P>
<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
  puts("\n按回车结束");
  getchar();
% j( ~. U% q+ J" r, Y2 B/ a  exit(0);
( ^4 A* a7 Z; j6 r) J2 G+ j }</P>
$ B* D, }. A& N3 ]; R( i<> //转换为第二阶段价值系数
exchange(c,temp_c);
/ v+ y- h  t# ]1 _+ v# V //把人工变量列全设为0
for(i=1;i&lt;=num_ar;i++) {
3 `/ P, x( F. y+ [  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
9 ~9 ?0 Y: }& h) s4 V4 S* P2 p }</P>
<> puts("\n\n第二阶段问题为：\n\n");
7 h- u8 F+ \0 s! _! _4 K& Z show();
* a# ]! S- q$ D/ g puts("\n\n按回车开始第二阶段迭代");
getchar();
iterative();
. w0 F- O0 f! ]  c switch(status) {
case 1:
# M' u1 u( ?1 g  output();
1 [. ?, Q/ [6 [9 p: x4 e! ?4 w  puts("\n\n原问题有唯一最优解。\n");
  V7 C  e7 q8 E  puts("\n按回车结束");
9 S  U3 Y( {4 D4 D! E2 G  getchar();
' B8 I4 ^9 C9 T/ A  exit(0);
7 C/ o9 _5 Z' w case 0:
  puts("\n\n原问题为无界解。\n");
  puts("\n按回车结束");
  getchar();
: |7 d8 |5 o4 y# E3 A" r  exit(0);
case -1:
  output();
* q( ]8 u. t- P- v5 j  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
case -2:
  puts("迭代超过限制次数强行终止！\n");
$ \% y4 C# H1 a( T$ W( m, {  puts("\n按回车结束");
  getchar();
  exit(0);
}//switch
}</P>
<>void banner() {
printf("\t\t****************************************\n");
! ?) ~9 Q# \( d; y+ T3 }: L$ }5 s printf("\t\t         单纯型法解线性规划问题\n");
printf("\t\t                         作者：Thunder\n");
" A, v/ q& |8 G: d" P5 K9 _ printf("\t\t****************************************\n");
printf("\n");
" ^! F4 v- A$ T3 X9 S+ C}</P>
  N5 g5 f; l+ s  O+ \+ U$ w<>void show() {
. Y7 @6 c  h8 I* q- e4 C" c0 ?//对方程组以自然的格式输出，系数为零的x不显示
6 V3 I  P/ |* @; f//为1的不显示系数1，-1系数只显示负号
int i,j,k;
1 @4 A2 p4 }* f5 v7 U* J switch(step) {
case 1:
  printf("min z= ");
( a5 V. t- L8 N, n4 n  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++) {
1 W( p  a/ ^& T' s4 b2 n1 q' L  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);
4 R+ y1 x% E) p, a  v  s }
break;
4 `& u: `( k4 ^8 F9 ?: F }</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  k=0;
! N# M# u8 l, t7 h  for(j=1;j&lt;=num_x;j++) {
0 v3 F% s% f+ v/ U) \( G) ^   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
: U2 u" q; C) [% I    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
; T- Z/ s% F9 w+ F4 C$ u% M! E    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);
9 Z. ~1 a  A4 Z) L    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
" o- [9 u8 t2 h8 w* `    else printf(" %lg x[%d]",a[j],j);
4 n' I1 t) r/ J, _, U. h: |    k=1;
   }
  }
  printf(" == %lg\n",b);
  Q1 |; l$ i3 U; v& O }
printf(" x[1]~x[%d]&gt;=0",num_x);
- j+ [  f$ @& I5 |* X2 @}</P>
! I! V" K6 N7 y& Q+ S& K<>void exchange() {
int i;
" v8 c+ K' m7 V: Q2 e9 L double temp[MAX];
7 J7 Z1 v/ ~5 t4 h5 y; x for(i=1;i&lt;=num_x;i++) {
! D6 q% u: V- J& r; A9 B* k  temp=temp_c;
  temp_c=c;
  c=temp;
3 s5 G% |5 G  |. } }
' Q* N( ]( E4 b- p* A0 \}</P>
<>void create() {
' b* V5 S4 c, D; Y//输入方程组系数，每个方程输完后回显确认
& s. i8 R9 D8 L8 \1 ] int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
) R$ S7 L+ x/ k+ u. J* y* p& o7 j$ w7 ~! k char confirm;

while(1) {
  printf("请选择：1、求最大值，2、求最小值：(1/2)");
/ g# ^, X) S& r  scanf("%d",&amp;ma_mi);
  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
+ G3 T3 b. P4 j0 ~' |4 ?  else break;
}
2 F" u7 Q% \- j, p
while(1) {
" h! e. d- @0 A9 A9 x3 P  printf("指定变量个数：");
5 }5 g6 g4 c, l1 o  scanf("%d",&amp;num_x);
  printf("输入价值系数c1-c%d:\n",num_x);
' f/ u, [1 ~2 ~: e" Q& X5 D* q$ }& N% A  for(i=1;i&lt;=num_x;i++) {   
- m6 B6 a) ]* U- {8 A# I   printf("c%d=",i);
   scanf("%lf",&amp;c);
  }
- a3 X. t% D( Q) _  t3 h  if(ma_mi==1) printf("max z= ");
  else printf("min z= ");
) O: R' r* h3 y- i& O  printf("%lg x[%d]",c[1],1);
  for(i=2;i&lt;=num_x;i++) {
9 }- o4 \6 O% W   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
. X; E3 N' A3 ~: u+ c, t# ?4 D1 A   else printf("  %lg x[%d]",c,i);
6 d8 w" }* B  w  }
# J* `7 u$ Q. P. G printf("\n正确吗？：(y/n)");
* h8 b% R: I( F! H1 b getchar();
% I8 V6 A6 W; h  a0 | confirm=getchar();
if (confirm=='y') break;
  _, b) L& H9 |7 z9 H! N: Z else if(confirm=='n') continue;
}</P>
7 I' I3 j5 i& D<> printf("输入约束方程组个数：");
scanf("%d",&amp;num_st);
for(i=1;i&lt;=num_st;i++) {
' S0 G( c5 w# [. M  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("输入错误，请重新选择。");
   else break;
  }
2 C2 M4 x0 M1 P  printf("输入技术系数：\n");
  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
  }
  printf("输入资源拥有量：\nb%d=",i);
$ M7 I" h* J( V  scanf("%lf",&amp;b);
, y/ |1 U) H4 X/ V( u: c: b+ ~3 X/ m5 r  
  printf("st.%i:\n",i);
  printf("%lg x[%d]",a[1],1);
  for(j=2;j&lt;=num_x;j++) {
' d( P$ Y7 t% l7 U6 K" r% `   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
   else printf(" %lg x[%d]",a[j],j);
  }
  switch(re_st) {
   case 1: printf(" == %lg",b); break;
/ {# f+ P, A' j: ?   case 2: printf(" &gt;= %lg",b); break;
* |' W! M5 L# Y" @5 r. K5 s   case 3: printf(" &lt;= %lg",b); break;
  }</P>
2 R7 t3 r; @# V- F) C<>  while(1) {
   printf("\n正确吗？(y/n)");
   getchar();
   confirm=getchar();
   if (confirm=='y') break;
! `* s2 M& z1 ?2 q+ S; N   else if(confirm=='n') {i-=1; break;}
8 s% d" U# c/ p0 B: _& r% I1 c  }
}</P>
) p' \8 U1 |. N4 `  c! O5 m<>//显示输入的方程组
printf("\n原问题为：\n\n");
if(ma_mi==1) printf("max z= ");
else printf("min z= ");
( n) S5 v4 N+ t* a# O 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);
' O8 G3 V/ x' F  G  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
5 Z7 D3 |$ K' ~$ x- @& y  else printf(" %lg x[%d]",c,i);
}</P>
# U6 ^" p$ O. C: |+ v<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
9 T# U, a, _9 S. w  k=0;
! `8 f5 ^1 P$ y& Q* i4 I  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);
- o- T. ?3 F/ W! `. W& N/ v    else if(a[j]==-1) printf(" - x[%d]",j);
! s, A' o) Z- n- ]8 c; Z    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
+ a" I" u( b/ q4 w! v1 f    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;
   }
  }
) ?6 B9 `; Z7 D) E  switch(re_st) {
   case 1:
    printf(" == %lg\n",b);
$ x, ^: V9 S7 c& a7 a2 ~4 }1 }    break;
   case 2:
' B' s/ \+ w- `; S# ~( e6 h) l    printf(" &gt;= %lg\n",b);
9 R: `$ R1 [- Q! }; P( F/ N7 @    break;
   case 3:
0 O( i# P* X4 }( {    printf(" &lt;= %lg\n",b);
    break;
8 _( r9 M+ N5 w' ]* g  }
. E& \% A! U) u5 b5 {8 \) r8 Q }
- W. W" S# J) Y( O0 S( ]# m9 r printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
) k/ v+ R/ w6 C6 n& k+ L8 `<> tnum_x=num_x;
& w$ X1 i3 C% A4 p, s. ?5 g4 b3 L for(i=1;i&lt;=num_st;i++) {
) b  D8 `+ ~" B- ]$ Q! P  switch(re_st) {
  case 1:
6 _& {+ E" V  P6 H$ N4 i) p/ c% a  case 3:
3 A6 u) a. z% t   num_x+=1;
/ Q/ Q9 Y* f& [: I   break;
  case 2:
5 K) u7 S) e) V: F   num_x+=2;
   break;
3 [) U) D! G' G, R$ u  }
}</P>
; A' R, X% O( c! ]$ G7 g2 F<>//化为标准形式
: Z7 Y, I- X2 t6 ]& ?6 A 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) {
   case 1:
    num_addv++;
    num_ba++;
    num_ar++;
    c[tnum_x+num_addv]=0;
7 \8 `; S: S0 ~' U    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
! a3 Y" w' A# u' O$ G    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;
) R: r+ ~; m, x- k" ?( z0 G( @( W( v    break;
   case 2:
) V& h- m( _4 B/ R  V    num_addv++;
    c[tnum_x+num_addv]=0;
* b7 ]" [9 L* f2 L    num_addv++;
    num_ba++;
; X3 w4 |5 q; I5 _, U( k! [4 e. `    num_ar++;
( }4 z, o" o7 U2 o+ E    c[tnum_x+num_addv]=0;
  n  `0 i; u, W4 X' r' k    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;
     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
8 W0 h, R+ m- ~1 [     else a[j]=0;
    break;
   case 3:
9 ~; @* Y0 _& C" E  F2 P& p    num_addv++;
    num_ba++;
    c[tnum_x+num_addv]=0;
    base[num_ba]=tnum_x+num_addv;
$ h( H4 K% X2 t3 l$ s% q, s    for(j=tnum_x+1;j&lt;=num_x;j++)
" e2 G- U6 W/ Y" Y- d3 \     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
0 D* ^7 Y7 n2 b! _4 l  N     else a[j]=0;
    break;
# C+ q. b! m# L" s% v  }//switch
}//增加松弛变量、剩余变量、人工变量、确定基变量</P>
<>//显示标准化后的方程组
printf("\n化为标准形式后：\n\n");
; K% q% z/ b' J; Q' g if(ma_mi==1) printf("max z= ");
else printf("max z'= ");
printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_ar;j++)
  if(i==arti[j]) k=1;
  if(k==1) printf(" -M x[%d]",i);
. Y* B8 S8 c9 o  j8 F5 f  else if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
$ A$ x$ H. X' s! C6 F7 a" K  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
  else printf(" %lg x[%d]",c,i);
& N! H3 e0 g( h9 w1 @ }</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  j2 M( f# p% l5 |1 ]  k=0;
  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
2 p  N# m( p6 \/ ~) }    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);
% l, {' D4 R  S* r6 T) C0 {1 `    else if(a[j]==-1) printf(" - x[%d]",j);
* ~/ l, |: l4 O* F    else if(a[j]&gt;=0&amp;&amp;k!=0) printf(" +%lg x[%d]",a[j],j);
* B' |3 u! F% y9 m* u    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
& g. P5 Q4 a3 `* k. B$ h    else printf(" %lg x[%d]",a[j],j);
$ |; `. \" E; s. Y0 _    k=1;
   }
  }
  printf(" == %lg\n",b);
" \% s+ l1 F, K  ~+ U1 _& c }
" p; Z& M6 K) B2 @4 r0 j printf(" x[1]~x[%d]&gt;=0",num_x);
9 z6 W+ h" {# K2 z8 W- `6 [}</P>
<>void iterative() {
int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
3 h1 v) V( N/ Q4 R int base_elem;
int base_out,base_in;
double sigma[MAX],temp;
! g6 M3 i6 U! ?0 M0 w0 X9 X5 f6 X double value_be; //高斯消元里保存主元素值</P>
% F1 q, j! R7 p( j2 D7 `<> 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>
. K7 U/ m: F" a- P3 S1 t9 S<>  switch(step) {
: f7 L  R/ u1 C   case 1:
6 ]5 @6 _) u- ~/ d7 k# _    for(j=1;j&lt;=num_x;j++)
8 J* O, r" c+ c; I5 Z    {
' S) K; P# F2 d3 B: j     printf("a[%d][%d]=%lg\t",i,j,a[j]);
! I+ x/ k0 v; o9 h! v' e( p    }
    printf("\n");
    break;
   case 2:
    for(j=1;j&lt;=num_x;j++) {
' b. V- U4 |1 D3 c     k_a=0;
     for(l=1;l&lt;=num_ar;l++) if(j==arti[l])k_a=1;
- [- Q0 L& ~8 {7 B* \     if(k_a!=1) printf("a[%d][%d]=%lg\t",i,j,a[j]);
    }
    printf("\n");
9 c* b1 E6 ~. T    break;
- R" S- ~- i( F- D* b0 p. e" U, N  }
8 |% J& _# u( V }
//求检验数sigma
. g, f+ X# S4 G: [; C  E+ ` for(i=1;i&lt;=num_x;i++) {
  sigma=c;
  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 L) s5 ~' [. R9 m  switch(step) {
' d0 F! x+ s* H/ H   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;
    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
4 W$ R8 a0 @3 H5 [    break;
0 V, T! X. [8 c$ O" \' ^) b5 ~  }
}
8 R* {3 t& @7 {# Q; D2 X, Y putchar('\n');
% q) d/ n' O* s5 C1 W8 A* A//检验检验数sigma是否全小于等于0
k=0;
  C) s7 b7 h" m' D for(i=1;i&lt;=num_x;i++) {
  if(sigma&gt;0)
$ ?" t6 B5 @# e6 p- a1 ~# c  k=1;
4 u  x$ Z9 S) H! x7 A9 S  h }
& r  j- H' V, d+ ?! z if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
  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) {
0 B  H9 t% v; g5 N$ h    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
, Y5 ~/ F  f& x- s   if(k_f==0) {status=-1; return;}
2 A, y- ]$ O% A; h( J   }
  }
  status=1;
+ N" n+ T' o5 y5 B  return;
/ V- P) Z- r' o$ s" L }
" w2 P: y9 c2 t! y  M//检查是否为无界解
" q! w6 T2 Y1 v9 O  E for(i=1;i&lt;=num_x;i++) {
# W5 f& m$ W& t% ]+ N) X1 U% U3 h  k_f=0;
) F; K* J8 ?# U  if(sigma&gt;0) {
   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
4 \' L: R$ Y. m2 B# R, {   if(k_f!=1) {status=0; return;}
  }
}</P>
<>//确定换入变量
3 s0 C) C- m) ~! Z$ t for(i=1;i&lt;=num_x;i++) {
, r7 x; o2 N9 S/ k  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
. L" J/ O  o" `* i3 i# B* [  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
+ }- b" w: g8 m1 w- G0 P- ]* d }//temp赋初值
. b1 N/ A, L: F1 M: I 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)
! E! ^% x! K$ \( _   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
    base_in=i;
    temp=sigma;
   }
}</P>
- c# G# g: p0 C4 W9 b% u6 r<>//确定换出变量
* b) g1 m( L! j' b. B6 E; g for(i=1;i&lt;=num_st;i++)
  if(a[base_in]&gt;0) {
, W. F- G2 f' E: s$ V$ Y( o4 \   temp=b/a[base_in]+1;
   break;
  }//temp赋初值
for(i=1;i&lt;=num_st;i++) {
* ^  Z7 s: z* o  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
- D( W4 j: v1 I& ?: K7 T& P   for(j=1;j&lt;=num_ar;j++)
3 P; r* K# a) g2 K: B* v+ i    if(base==arti[j]) {
5 N; @8 N7 _0 b; B     base_out=base;
8 @' {0 p: ]/ o& R6 @8 x& [     base_elem=i;
, Y' n% Z6 Z5 x. ^/ f! r     temp=b/a[base_in];
. y* }- h: g; e     break;
* x# j+ e; C5 k6 n6 Q' t3 y4 l6 O  t7 E1 Y    }
7 O- j( Z1 b; v: v6 O  }//人工变量优先换出
  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
. Z9 p2 q2 o3 R$ v  B   base_out=base;
   base_elem=i;
# Y% D2 r" |4 w3 Z6 r- x   temp=b/a[base_in];
% K4 \0 G0 z& |/ s2 ~  }
}</P>
, u6 s$ O: l. r' t<> printf(" 基变量：");
1 X, M% j. b  }! P* \ for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
8 X/ u( i. `6 ?4 W" s0 f! B//基变量变换，进行新方程初始化后迭代
- n/ L% W9 T" Z( k" c for(i=1;i&lt;=num_st;i++) {
" G7 w; U/ C3 T- m  d4 n) B: x  if(base==base_out) base=base_in;
}
9 F4 Z2 n; \9 \9 \9 o8 [//初始化主元素行系数
) h5 G# i+ \' c& u1 X value_be=a[base_elem][base_in];
b[base_elem]/=value_be;
4 ^) I+ q& ^9 l- E, l- y; @ 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) {
% i0 d# P" T& R  Q# b. U   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;
  }
}
5 _7 F: i6 @9 d4 y! j  Y4 U  Y# T stop++;
5 u' D; f$ ]! b* ^ if(stop&gt;STP) {status=-2; return;}
5 G; `9 T8 i: f# K4 }$ t9 d8 s iterative();
}</P>
<>void output() {
int i,j;
double X[MAX];
printf("\n结果如下：\n");
printf("\nX=(");
' @5 w8 _% u5 }! m" F, w3 _ for(i=1;i&lt;=num_x;i++) {
! L0 b% F6 b1 D2 Z. W" ^+ M3 Z. m! H- L  for(j=1;j&lt;=num_st;j++)
! P/ u0 j# k2 m1 m8 Q1 o  if(i==base[j]) {X=b[j];break;}
, V, q' q. g5 v( h) g  else X=0;
  p0 J' K4 t5 Z# c0 E& ]  printf("%lg ",X);
! p/ n" B. |! _/ u& H3 I }
printf(")");
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);
) a& o. k$ G2 m, V* D# s}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了