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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
* z+ {# j, I# E8 q  E1 z: j#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])
2 F* `, w+ W: E- n3 I{
# s; ^8 a1 x1 C3 G  int j;
' Z, \- \% [$ |  *c=0;
  for(j=1;j&lt;=*is;j++)
    {
$ T0 w$ V! `% O3 Z1 a0 W9 ?      if((*a)[*ir][j]-*c&gt;0)
. ]5 @3 {* u) n' g' `    {
: i: |# @! k% ^" B+ e      *c=(*a)[*ir][j];
      *j0=j;
* O" f) H# ^- }+ A  V) d" }    }
8 B% w% b- k% ~: w9 H* v5 w     }
* s5 f/ n' ^3 o8 m}</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>
" b7 I7 }' j' {0 e0 `! r+ uint xi_sm(int m,int n,int m2,int mn1,int l1,float (*a)[X][Y],
" U. ^+ r' v, o! N* D7 h       int(*k)[],float(*x)[])
+ G: M# H9 F* c( Y0 ~{
  int i,m1,mn,j0,i0,j;
  float c,g;
  m1=m+1;
8 x' h- {9 V4 I- x5 a& K  mn=m+n;
  for(i=1;i&lt;=m;i++)
; u) M2 o( u; F    (*k)=n+i;
, G- b0 m- I3 \  H$ k/ y. T+ Fleap1:
  if(l1-1==0)
4 W1 i: L( Z, R  A. a" T* d     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  else
     {
$ y0 R: l  G9 F8 z# g2 y# c) g7 qleap2:   if(l1-50==0)
; Y1 T, M, O0 T* P$ Z. @, W      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
       else
4 x9 x; b! ]1 h* H      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
$ [+ Y$ o& g) g; y% e. B- K  P     }
   c-=1e-8;
   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
      {
/ p+ [6 K- F5 U( ^7 E" A    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
    return(-1);
      }
   if((c&lt;0)&amp;&amp;(l1==1))
       {
     l1=50;
& C8 p. ?" R5 ]$ Y% X     goto leap2;
       }
/ O3 E1 J' O: k* q' f& \8 E    if(c&lt;=0)
       {
$ B+ G8 Y) E4 {, d! g  c7 s: {0 y      for(i=1;i&lt;=m;i++)
         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
      printf("\nF=%f\n",(*a)[m2][mn1]);
0 K! |+ Q- s2 B      return 0;
- A2 r1 g: L0 J5 _. K2 g" R    }
  u2 N- s; f! Z+ q; {     c=1e8;
% b! }! I  w% o  P5 Z7 E) W1 t     for(i=1;i&lt;=m;i++)
       {
3 U/ J4 d" c; k" b9 y     if((*a)[j0]&gt;1e-8)
/ {/ h2 U9 d: ?/ }# ]) J$ ~5 n        {
          g=(*a)[mn1]/(*a)[j0];
          if(g-c&lt;0);
         {
           c=g;
. d  ?: w' \9 Q  A" G3 I' T           i0=i;
, l3 E8 t& ^8 q1 t6 ]          }
' C! ]& N4 b8 V! r+ x( q         }
    }
2 {( G3 T- A! F9 t# t      if(c==1e+8)
    {
      printf("\n\t*******Lp no solution********\n");
      return -3;
% n- F5 P* \$ h4 f6 H     }
8 n4 o# @" s1 O. a$ l3 X       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
      {
# C& P4 ]% M( _1 U2 S% r! ?8 K        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;
& y5 ?. s  }$ x. D! @) c         for(i=1;i&lt;=m2;i++)
( \5 p  p7 R* f( X/ v- j2 U           {
         if((i==i0)||(!(l1==1)&amp;(i==m1)))
            continue;
4 x( D0 ]; b" T- O' z# A# h0 I3 _         (*a)[j]=(*a)[j]-(*a)[j0]*g;
6 H2 X# E& e6 n( `! W           }
7 N2 t8 Z) N! S7 \9 L) {       }
- N$ T/ V( q. [+ \( l/ L       for(i=1;i&lt;=m2;i++)
( j9 s! [* ~2 K* S      (*a)[j0]=0;
       (*a)[i0][j0]=1;
/ A* t7 B/ R" G: A8 v% r  U       goto leap1;
+ \" [- U! W* o+ U, U# J* L}  
! R$ f  p* ~1 n2 x. H7 bmain()
: Q& J- t2 E  a{
  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},
" L" n; M4 F! E8 _+ v2 y3 i* d' N                 {0,0,0, 0,  0,0,0},
                 {0,69,0,144,0,0,300}};
  int k[3];
  float x[3];
- k1 ~8 V  Z5 `; ?. e! c$ S, i" P% ^  clrscr();
  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
}
</P>

plgatc

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

编程环境：VC++6.0     
方程组输入说明：
/ Z. f6 B9 e& b* R+ [% T% L/ R! { 变量非负，按提示输入相关参数。
  l# @/ D3 ~+ V7 I1 t *************************************************************************/
# R5 i# k; x4 i3 D( W#include &lt;stdio.h&gt;
#include &lt;stdlib.h&gt;
& E$ h0 c" b: _" H5 a! y#define MAX 100
3 w* P& T6 g" m+ Q1 c' }7 H1 m  |#define STP 100</P>
( V) b" j' G0 e3 P  N! V<>int stop=1; //迭代记数变量
. a* X5 s$ O6 Kint status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
int step=1; //目前阶段</P>
<>double a[MAX][MAX],b[MAX],c[MAX],temp_c[MAX],max=0; //方程组相关系数
int num_x; //变量个数
int num_st; //约束方程数
int num_ar=0; //人工变量个数
int arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
<>void create(); //建立方程组
void iterative(); //单纯型法迭代
void output(); //输出结果
void banner(); //打印程序标题
- b6 }( P. E/ j1 T% q8 D4 cvoid exchange(); //交换两阶段价值系数
1 T/ F8 y+ g8 R! Q( }void show(); //输出方程组</P>
, z  a; ?% n& o$ ]8 X; l9 B: e<>void main() {
int i,j,k;
( ?1 V: c+ S, g6 R$ k; I$ ^ banner();
create();
8 [1 Y  \9 D! Q( A' t //保存原价值系数，转换为第一阶段价值系数
for(i=1;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
& x( S6 O& j- [! X) M* w1 [( K& M  if(k==1) temp_c=-1;
  else temp_c=0;
}
0 C& N9 G1 x+ c+ @ exchange(c,temp_c);</P>
<> printf("\n\n第一阶段问题为：\n\n");
show();
) p2 R" I; K! d step++;
/ X+ O! [7 L' _, g) y% Z0 F printf("\n\n按回车开始第一阶段迭代");
# i7 K# u" Y2 x+ w' z" X4 M+ p: i getchar();
getchar();
iterative();
if(status==-2) {
  puts("迭代超过限制次数强行终止！\n");
/ u1 {) S& I& ]/ a7 g  puts("\n按回车结束");
  getchar();
1 S6 s0 O8 j, {( x  exit(0);
}
' P+ D- X" }+ P% ]6 o output();</P>
7 r! S- R: ?6 s9 T8 d  z<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
  puts("\n按回车结束");
  getchar();
$ d) T9 G8 @5 P1 |* k% M  exit(0);
}</P>
<> //转换为第二阶段价值系数
( H0 A# Q# F; Z exchange(c,temp_c);
//把人工变量列全设为0
for(i=1;i&lt;=num_ar;i++) {
  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
/ j  n' v% ~  `$ u# F/ w }</P>
0 l) A3 C! o& u) a( T1 N<> puts("\n\n第二阶段问题为：\n\n");
+ f4 c  K, y1 V9 S3 C1 R7 O: A+ j5 |, | show();
puts("\n\n按回车开始第二阶段迭代");
getchar();
iterative();
2 I4 q2 ?2 U$ f switch(status) {
case 1:
  output();
  puts("\n\n原问题有唯一最优解。\n");
6 t: R0 a5 J3 y) S  puts("\n按回车结束");
, `8 h8 @/ E6 l  getchar();
  exit(0);
5 G$ w* u) p1 Z/ Z case 0:
  puts("\n\n原问题为无界解。\n");
) _- u3 T! O8 G  puts("\n按回车结束");
  getchar();
: h) M0 p4 c! V% w5 R  exit(0);
case -1:
8 x1 j% l: e5 k* X9 m! m  output();
  puts("\n\n原问题有无穷多最优解。\n");
7 \& Q1 V/ X: ~, h7 T1 N; ^  puts("\n按回车结束");
  getchar();
  exit(0);
case -2:
. w1 \6 _. [9 g1 H# r0 d2 O$ N  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
% s0 z& e8 C" f% a, j+ z0 P8 S }//switch
}</P>
4 l+ G2 [5 R+ Y% x7 }<>void banner() {
$ {5 O4 z3 j$ f8 N, t1 s. O8 ^" l printf("\t\t****************************************\n");
2 v  o0 |0 p6 {- X printf("\t\t         单纯型法解线性规划问题\n");
printf("\t\t                         作者：Thunder\n");
printf("\t\t****************************************\n");
printf("\n");
}</P>
! k; W0 q7 t. m<>void show() {
( S3 k% v- g2 O: Q7 B//对方程组以自然的格式输出，系数为零的x不显示
//为1的不显示系数1，-1系数只显示负号
int i,j,k;
) _- n' I# b, D+ l$ t% ^0 u switch(step) {
case 1:
; u* R" C) b6 _; w+ _  printf("min z= ");
6 b% F5 y1 q& ]( u2 K$ f0 k  printf("x[%d]",arti[1]);
  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
: t+ _1 ?; \( y4 M" ^, C case 2:
& M  j9 q( |) L  G, F0 ?& G printf("max z= ");
  h6 U) e9 c" I$ ]) G$ s4 t printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
3 _1 p  M  o- X0 o/ u, m; W8 N  else if(c==-1) printf(" - x[%d]",i);
9 r( h$ s2 b& h4 X0 p7 y" @  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
+ X# g3 o2 W8 a; j+ \  else printf(" %lg x[%d]",c,i);
7 {- w% y! m& c4 R }
: J8 H8 s7 ]. H- x5 A4 e  l  c break;
}</P>
1 D' y- }0 I; P$ x* r9 l5 W4 y<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
6 ?& \& r/ L8 e$ n2 B; F9 U   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
8 v, K0 ^6 C: K+ g3 c    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
    else if(a[j]==-1) printf(" - x[%d]",j);
% Z- X% z" B. o7 j2 l    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);
) g" ^- t! X- M/ w- c    k=1;
   }
, i1 T2 A1 x7 e. t  }
  printf(" == %lg\n",b);
% }" I& @& R$ F' z: W }
printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
<>void exchange() {
2 K, g. n% u; s- \7 H int i;
double temp[MAX];
# g- |; e1 D( o$ D/ q3 X) A for(i=1;i&lt;=num_x;i++) {
" }+ z; Z2 ~# J* J  G' s$ s/ n  temp=temp_c;
  temp_c=c;
  c=temp;
}
( n8 Q, V; _* H/ S4 j. j4 f}</P>
<>void create() {
% H/ q0 `4 y% z5 _//输入方程组系数，每个方程输完后回显确认
7 Y1 q$ E, W" z" y int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
char confirm;
: h! o% B. z5 J# V  q/ M0 a5 f7 t
6 j- J" A# L5 |! o& K while(1) {
3 L2 d* e3 V6 T  printf("请选择：1、求最大值，2、求最小值：(1/2)");
% i; i' C4 T( d# O5 S! T  scanf("%d",&amp;ma_mi);
  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
  else break;
, x2 X$ B7 {, t! D }

while(1) {
2 J* _3 p6 j; D& W% G  printf("指定变量个数：");
) C5 r2 E2 z* i5 \, o  I2 ~  scanf("%d",&amp;num_x);
8 X6 w9 o0 I- h% w( W6 \, L  printf("输入价值系数c1-c%d:\n",num_x);
7 n6 d! u2 E$ G/ ^  I  for(i=1;i&lt;=num_x;i++) {   
   printf("c%d=",i);
2 j* n: {2 }) O# ]5 `2 k   scanf("%lf",&amp;c);
  }
) Y6 Z% L+ C0 a2 F  if(ma_mi==1) printf("max z= ");
. K9 d  ?4 j; t) g  u  {. c  else printf("min z= ");
  printf("%lg x[%d]",c[1],1);
  for(i=2;i&lt;=num_x;i++) {
# V; {  r! r" a6 K+ l   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
) Q7 O, a* t( {% x. a   else printf("  %lg x[%d]",c,i);
  }
+ {* y% l, Z7 R9 m# C4 D( t3 |' Y printf("\n正确吗？：(y/n)");
$ k. f% n, t+ |+ D" t. P1 Y getchar();
confirm=getchar();
if (confirm=='y') break;
6 ?5 W3 ?8 i9 A5 _1 n- i# U else if(confirm=='n') continue;
2 F+ ?9 V. l8 F# g, C( H* H: C }</P>
<> printf("输入约束方程组个数：");
4 F! i/ [) R0 O7 {* _ scanf("%d",&amp;num_st);
, d: |* p" L; k4 S) ~1 M) P7 v for(i=1;i&lt;=num_st;i++) {
/ D9 b+ J7 C' X" t# ^) y/ v  printf("st.%d:\n",i);
  while(1) {
   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
' w4 n, {0 E3 e- C/ m+ K& P& y& F   scanf("%d",&amp;re_st);
( l+ P* Y- }& J6 Z! k8 N   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
, U; {) K6 i/ ^   else break;
  }
9 N9 M) R0 S0 n8 W' n  printf("输入技术系数：\n");
7 A- {. t) b, @" B) b  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
  }
3 d8 \/ F; X, `* l  R  printf("输入资源拥有量：\nb%d=",i);
# d( s* I- O! M. S. d& O  scanf("%lf",&amp;b);
  
  printf("st.%i:\n",i);
" ?3 @' u7 n9 Q- 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);
   else printf(" %lg x[%d]",a[j],j);
% h8 A6 l7 x" f! O# T4 G- X  }
% s# c0 Z8 o5 m% ^5 ]  switch(re_st) {
   case 1: printf(" == %lg",b); break;
. c% v. A: _% Q# h   case 2: printf(" &gt;= %lg",b); break;
, T3 [% F3 B# m   case 3: printf(" &lt;= %lg",b); break;
  }</P>
. E3 ?6 J. _, J, }' ~<>  while(1) {
$ `* X6 w, E" _& V: f$ o   printf("\n正确吗？(y/n)");
   getchar();
) S1 [! V8 e/ z3 g: P   confirm=getchar();
+ _6 h% z1 U, y8 O2 i   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
0 x* _- A5 Q! |; G2 z; K+ F7 W  }
}</P>
- R( M& X: H2 g6 y1 @<>//显示输入的方程组
' {! u8 U% F* P9 w" G printf("\n原问题为：\n\n");
8 f* y5 i" z2 C- l% ^+ U: P9 B0 X if(ma_mi==1) printf("max z= ");
& k6 L" x$ ~1 \- |9 N else printf("min z= ");
printf("%lg x[%d]",c[1],1);
' N+ c/ B7 c+ o3 w: C/ c5 h6 @ for(i=2;i&lt;=num_x;i++) {
' {, j1 @0 V" f7 [" P5 q: k8 ~, ^  if(c==1) printf(" + x[%d]",i);
* p9 Z/ w0 K; q4 o1 v& P  else if(c==-1) printf(" - x[%d]",i);
7 u) `* S+ b) u0 W  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
; @0 G2 ?  l! j3 ]2 d0 p# w8 N  else printf(" %lg x[%d]",c,i);
, p5 |0 I2 N+ @9 a7 R }</P>
* B+ ]. s% d5 ^6 b# l* m3 t<> printf("\nst:\n");
2 W9 J3 D& _1 G for(i=1;i&lt;=num_st;i++) {
  k=0;
2 s2 h5 s6 @6 `4 `- W& {  for(j=1;j&lt;=num_x;j++) {
9 ~( V0 v3 [( H. |- F9 q   if(a[j]!=0) {
7 S1 c. q1 U4 Q, L( P1 b    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
- j" w- t" p0 V- x5 q. S# c    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
8 X- t! o7 t7 |0 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);
- `  E, W2 b2 a5 B; g    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
! w4 C$ a( Z) r# V! Y5 v    else printf(" %lg x[%d]",a[j],j);
1 r, K- X" H& _- S0 _3 M7 ]9 L) X    k=1;
   }
$ {  w1 g; L2 x# K# r& j. c  }
! \: c& ]6 h% A! u) O  switch(re_st) {
   case 1:
8 t8 l9 Q% W$ T. }# k! V6 b4 l6 R    printf(" == %lg\n",b);
    break;
! R7 S4 V3 G2 b: H' a: x- m# d   case 2:
" @+ c: B8 J( N9 P  v    printf(" &gt;= %lg\n",b);
    break;
3 i: y; R% e% g7 g- u' ^$ {   case 3:
9 Y- p! C+ l0 h/ {4 x$ ^; l( M+ p    printf(" &lt;= %lg\n",b);
2 ]# D5 k- `) g( s6 ?2 O, j4 V0 f    break;
  S# q4 C' j. s: d9 s  }
4 \4 g" w5 M$ X }
printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
<> tnum_x=num_x;
for(i=1;i&lt;=num_st;i++) {
& s0 {" k, z- k8 R2 i/ c  switch(re_st) {
  case 1:
, n! m( @, K& K/ V" T  case 3:
7 x6 u6 y6 l& d   num_x+=1;
   break;
  case 2:
   num_x+=2;
   break;
  }
/ X- ~& {8 o& u7 n }</P>
: i: r/ p& Z) ]; s' D<>//化为标准形式
) d+ G, f7 r: w, | if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
for(i=1;i&lt;=num_st;i++) {
3 }- y$ b3 |/ P* u$ M  F  switch(re_st) {
" T( |6 h" G0 y* A  s" V% I* C   case 1:
    num_addv++;
    num_ba++;
    num_ar++;
& k" `' p: m( {: k) x8 d    c[tnum_x+num_addv]=0;
( D* @, p7 z1 T3 j: b* n$ J    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
    for(j=tnum_x+1;j&lt;=num_x;j++)
: |  I% ~5 a* Z( P     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
% z) q" v* Q# u! z6 |: I     else a[j]=0;
    break;
# c+ O. l0 G( E7 @1 r6 n) y   case 2:
    num_addv++;
    c[tnum_x+num_addv]=0;
6 [: j. L: h6 t! ?/ {- }    num_addv++;
) V( ^( t- |+ [0 q, e" v    num_ba++;
+ f/ h% u% w( G2 {$ a    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++)
     if(j==tnum_x+num_addv-1) a[tnum_x+num_addv-1]=-1;
1 h7 q. y* U- i) `     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
$ y+ M8 {; m4 H0 W2 j# {     else a[j]=0;
    break;
   case 3:
    num_addv++;
    num_ba++;
    c[tnum_x+num_addv]=0;
    base[num_ba]=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;
) A* I- o. D  h; t     else a[j]=0;
    break;
  }//switch
0 F$ Y- d0 u: i1 v6 e }//增加松弛变量、剩余变量、人工变量、确定基变量</P>
1 e0 A7 s  J+ W6 d# k# Z2 u4 ]<>//显示标准化后的方程组
printf("\n化为标准形式后：\n\n");
3 F4 N+ u+ i4 K8 W* u% |; e if(ma_mi==1) printf("max z= ");
1 b+ K) m. s9 l' t; s else printf("max z'= ");
5 _9 M; }" H9 z printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
  k=0;
; L/ S+ i  Q0 {& n( r  for(j=1;j&lt;=num_ar;j++)
9 G$ y1 T: z' L  if(i==arti[j]) k=1;
9 E5 |4 ]3 {& K+ o% U7 ]6 G  if(k==1) printf(" -M x[%d]",i);
/ S0 P$ R# D5 o$ P  else if(c==1) printf(" + x[%d]",i);
: K0 H2 d+ y# ^: a  ]  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);
}</P>
<> printf("\nst:\n");
, v) B3 ]$ c  Z4 W& j! v0 Y/ x1 j for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
5 u* s1 G' K- @0 c- I/ \  Y1 h& d' x   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);
; m& t1 J. e5 b+ q1 K; P& P! r( @    else if(a[j]==-1) printf(" - x[%d]",j);
3 ~' j, m1 {4 \& w! e/ H4 N    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;
# h0 Z' M1 ~% c3 m2 G) W   }
  }
0 i' q) `1 P2 N8 h! `: r  printf(" == %lg\n",b);
}
1 m* R2 {  n$ o9 k2 W printf(" x[1]~x[%d]&gt;=0",num_x);
1 t0 K* E, b1 a5 L5 y: g; _}</P>
<>void iterative() {
( o, h9 v5 ]: n) A- ~* ?; Y int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
int base_elem;
* Q" _- ?- \1 u1 G& P int base_out,base_in;
double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
" h4 l! c2 k( ?4 Y/ Z( d  _<> printf("\n\n第%d次迭代：\n\n",stop);
. _* J8 U. }+ y5 m$ j4 V for(i=1;i&lt;=num_st;i++) {
  printf("c%d=%lg\t",base,c[base]);
  printf("b%d=%lg\t",i,b);</P>
2 v: V8 u9 ?9 R& X3 }$ O; ~- d<>  switch(step) {
   case 1:
2 w+ H7 i7 f& J! n2 D; m# J* J* W    for(j=1;j&lt;=num_x;j++)
    {
  \( Y# N5 A0 y; p" W9 Y$ J, |! B     printf("a[%d][%d]=%lg\t",i,j,a[j]);
    }
* _8 y' ~+ ]. Y9 e5 H  m& Y    printf("\n");
    break;
9 m5 B& G8 T8 w6 y8 _2 I   case 2:
    for(j=1;j&lt;=num_x;j++) {
     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]);
    }
    printf("\n");
    break;
  }
}
% b# F4 I& X  _//求检验数sigma
* x. q0 B& {) L. L6 v, \" a& x for(i=1;i&lt;=num_x;i++) {
  sigma=c;
: a* ]7 F5 G! I. z4 N  for(j=1;j&lt;=num_st;j++) sigma-=c[base[j]]*a[j];
3 B: N: t" s. g. z4 F  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
  switch(step) {
8 r" f' t# b' U1 w4 r   case 1:
    printf("sigma[%d]=%lg\t",i,sigma);
    break;
* j8 i* A; k) b! O2 n( K# B! E   case 2:
    k_a=0;
. R9 T* v% N5 }" a    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
$ v5 Z; E. {* {& v    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
    break;
  }
* ^3 S7 y% v; S% v3 U }
putchar('\n');
//检验检验数sigma是否全小于等于0
k=0;
for(i=1;i&lt;=num_x;i++) {
  if(sigma&gt;0)
8 r, K( r" V0 ^  k=1;
}
% Q6 _# }$ i6 ^! k8 F if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
  for(i=1;i&lt;=num_x;i++) {
& J- H9 l  F4 ^/ e8 V1 l0 o2 L   k_f=k_a=0;
3 l5 c$ N, S2 g. i! J: \   for(j=1;j&lt;=num_ar;j++)
   if(i==arti[j]) k_a=1;
   if(sigma==0&amp;&amp;k_a!=1) {
% L: L' C6 o6 h) H. o7 K9 y    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
   if(k_f==0) {status=-1; return;}
  M* g7 X3 x* K0 d) [+ Y   }
) z, {) Y/ D" T, V' D, @4 ~. C  }
3 d' g4 {3 Q' I3 s  status=1;
; n$ W* ~3 ^& x% P: K  return;
, t( a  `% U! r) k$ g  l }
//检查是否为无界解
9 a2 f  X9 {, ]- @ for(i=1;i&lt;=num_x;i++) {
  k_f=0;
  if(sigma&gt;0) {
+ m6 Z. Q0 N9 J: U   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
   if(k_f!=1) {status=0; return;}
# l  p* `4 Z, s) _  _  }
}</P>
% D/ ~& O0 }2 D& h9 E& s5 Z0 g<>//确定换入变量
for(i=1;i&lt;=num_x;i++) {
4 o, {+ |8 L- ]5 G- J" i9 \  k=0;
4 e- }1 h4 S; o1 T) F- ?! \) a  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
' _( Y8 u1 i. g  }/ f! z9 D8 N( t+ S  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
- ~, @7 T! i  X# Q }//temp赋初值
) y) L' l  `- n* f/ C! f4 V% T for(i=1;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
3 g, F. n  q. m3 v/ |  if(k==0)
- J2 [7 V' F9 L   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
    base_in=i;
    temp=sigma;
4 ]+ S+ C: j. \) h   }
}</P>
<>//确定换出变量
for(i=1;i&lt;=num_st;i++)
  if(a[base_in]&gt;0) {
5 X' X; @5 ?5 a1 O, C   temp=b/a[base_in]+1;
   break;
  }//temp赋初值
for(i=1;i&lt;=num_st;i++) {
  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
   for(j=1;j&lt;=num_ar;j++)
    if(base==arti[j]) {
     base_out=base;
     base_elem=i;
     temp=b/a[base_in];
! N0 R/ c! J6 v: A. q7 G     break;
    }
  }//人工变量优先换出
0 R% _* l: G) o! u% f; W% Q  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
, q( i+ }, E# B! \   base_out=base;
: C2 O# I" s; s; O0 d  A" Z8 F   base_elem=i;
   temp=b/a[base_in];
  }
}</P>
<> printf(" 基变量：");
for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
//基变量变换，进行新方程初始化后迭代
, L5 i% w* P4 S; b for(i=1;i&lt;=num_st;i++) {
  if(base==base_out) base=base_in;
}
//初始化主元素行系数
- c2 b, F5 H1 k6 { value_be=a[base_elem][base_in];
& M' c5 T: i2 D# M' i b[base_elem]/=value_be;
3 i. x% s: H6 D6 l( x for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
( U8 R* r5 [' M: a" E; G<> for(i=1;i&lt;=num_st;i++) {
$ {/ H. K: {# c  if(i!=base_elem) {
  `7 m+ {5 `+ \" r3 p1 q   b-=b[base_elem]*a[base_in];
   value_be=a[base_in];
- h& K+ h& J3 `9 W   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  }
: u% F9 p# d& S7 F }
" l1 ]) \3 G- W stop++;
if(stop&gt;STP) {status=-2; return;}
- r& c7 y$ ~) O7 z, [5 g/ L iterative();
}</P>
) a% x9 V$ F- `; ^<>void output() {
8 f  c6 x0 X! J2 g; d int i,j;
double X[MAX];
' Z- f' r, t  ~9 q printf("\n结果如下：\n");
+ h7 H- s( Z7 A' k& o printf("\nX=(");
  q  P* {7 }: o+ i1 y9 t2 C for(i=1;i&lt;=num_x;i++) {
  for(j=1;j&lt;=num_st;j++)
8 Z; j5 L& o  k! d  m  if(i==base[j]) {X=b[j];break;}
: R1 W5 u* v3 @. v* N  else X=0;
$ F: ?. ^. T& O6 P2 M0 n  printf("%lg ",X);
4 x% u6 }' K3 d" q& h; ] }
# L; Z3 Y+ D4 {$ J+ E  G9 R) J 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);
}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了