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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
3 x, y) V2 g- n7 q#include&lt;math.h&gt;
) I1 m* v4 C& ]5 y: S1 g! M8 B#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])
) ]9 y$ u6 p4 g1 @9 |$ g- @3 B$ t{
  int j;
+ J* j+ U6 R9 d+ E# _, K  *c=0;
4 d0 B# X0 s1 K7 S  for(j=1;j&lt;=*is;j++)
' P: W2 E$ w  |) Q    {
( l3 C1 p6 p) m9 Z4 C6 \* M      if((*a)[*ir][j]-*c&gt;0)
    {
; J( e* Z0 [% \! J9 o. @      *c=(*a)[*ir][j];
      *j0=j;
9 H! U- K  h+ ]) ]  K    }
+ n5 y6 W& J6 ?/ b; ~. J" D3 N     }
6 t  j- z$ ^& G# u- x}</P><>      /***************** 参数说明 **********************/
" x  s7 u4 A8 Q& e      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
! q2 L! G. P' q6 \, P      /*   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],
       int(*k)[],float(*x)[])
{
! f" ^; A$ Z; @! o5 }  int i,m1,mn,j0,i0,j;
1 Z" v( C- U; l# X- }  float c,g;
) g+ b. W+ P9 R$ P0 b  m1=m+1;
6 E" }- E" K  o% v0 @  mn=m+n;
' n! i8 Q1 y9 Y  s$ W( D8 S2 N2 Y7 ?  for(i=1;i&lt;=m;i++)
2 k) I  y1 I) G) M    (*k)=n+i;
leap1:
" [9 p5 Q; l9 E& n1 H! u" W  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  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);
8 z6 u0 C( v1 I# X, h5 A     }
   c-=1e-8;
   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
) }) [: t7 }8 X8 ]0 ]# A      {
    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
% ^% z7 P0 [; Y" {& x    return(-1);
      }
: X$ b9 e  a* n. L8 @   if((c&lt;0)&amp;&amp;(l1==1))
# I8 P" o0 z0 a. A( `9 ~       {
     l1=50;
/ L, @9 q  x- h! D     goto leap2;
) L7 ]: {% Z" I0 a# ?       }
    if(c&lt;=0)
3 b/ L; a8 u4 z9 v% f7 B0 R       {
      for(i=1;i&lt;=m;i++)
         (*x)=(*a)[mn1];
- ~, P3 H' A: e1 u- I. H4 A  [3 g         printf("\n\t********Optimal solution********\n");
& D* l  p! I  Z3 v4 A9 B7 H3 R" d      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
      printf("\nF=%f\n",(*a)[m2][mn1]);
9 [$ d+ e0 ~8 f( @2 N2 h      return 0;
    }
     c=1e8;
     for(i=1;i&lt;=m;i++)
( a* A; w; q; I9 {# B6 `" ?) p3 J       {
3 M8 q# G/ @- S3 _5 d; x3 O- N     if((*a)[j0]&gt;1e-8)
        {
          g=(*a)[mn1]/(*a)[j0];
7 ^- M8 \' u: u          if(g-c&lt;0);
         {
           c=g;
' `% W3 U7 Q$ f: P0 ]           i0=i;
/ t$ t( f0 Y3 t          }
" }& K2 ]3 \3 n7 q+ J7 y         }
    }
      if(c==1e+8)
1 p' H& Y( Y' A& y3 t7 t7 L    {
, B9 V: N! [( |6 U      printf("\n\t*******Lp no solution********\n");
3 o9 T7 x7 E5 }3 t4 {7 k- F      return -3;
     }
       (*k)[i0]=j0;
       for(j=1;j&lt;=mn1;j++)
* _( u; V* i3 r5 y! B+ W      {
        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;
         for(i=1;i&lt;=m2;i++)
           {
         if((i==i0)||(!(l1==1)&amp;(i==m1)))
            continue;
         (*a)[j]=(*a)[j]-(*a)[j0]*g;
           }
       }
       for(i=1;i&lt;=m2;i++)
0 {8 a' b1 r3 j) f+ E      (*a)[j0]=0;
* j8 S8 x% q" }       (*a)[i0][j0]=1;
       goto leap1;
2 |$ p% M5 u5 Y5 d( y9 j& {' A) P}  
main()
5 n% l- ~7 Q; j* s7 W{
  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},
                 {0,0,0, 0,  0,0,0},
, {& ^( C4 L$ f: l& d  R3 B                 {0,69,0,144,0,0,300}};
' }8 F3 h& J& G4 o7 d; K. H$ t  int k[3];
  float x[3];
  clrscr();
  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
7 ]; M+ d+ N. l1 ?, m1 m- g}
8 C$ X2 g' f, O</P>

plgatc

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

编程环境：VC++6.0     
* H: x5 Q. r+ Z- i& Q$ ? 方程组输入说明：
% g% ^" E* M' P7 e& L7 }: C3 b+ w' f5 [ 变量非负，按提示输入相关参数。
*************************************************************************/
#include &lt;stdio.h&gt;
#include &lt;stdlib.h&gt;
2 E% d( N' n& x0 u1 s. L#define MAX 100
#define STP 100</P>
<>int stop=1; //迭代记数变量
+ w) n0 q( b" t! f4 Iint 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; //变量个数
% d9 t2 ]9 l# p) \3 @; X) Sint num_st; //约束方程数
# P3 s% Z8 W  R' Fint num_ar=0; //人工变量个数
/ W, X4 k0 s9 i' k- R9 J: {, Jint arti[MAX]; //人工变量下标
int base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
<>void create(); //建立方程组
( f, V/ x# Z9 C. c% yvoid iterative(); //单纯型法迭代
! n/ B- X4 b# {2 B) R6 Tvoid output(); //输出结果
void banner(); //打印程序标题
void exchange(); //交换两阶段价值系数
' e: N% W) M$ v+ R+ E2 ovoid show(); //输出方程组</P>
<>void main() {
& P; p+ @$ P$ {( O int i,j,k;
banner();
create();
  L* k+ f! H3 y/ O+ W1 J' E0 M //保存原价值系数，转换为第一阶段价值系数
2 a' a3 q, U* T, ?+ Z: b- \ for(i=1;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
  if(k==1) temp_c=-1;
  else temp_c=0;
}
exchange(c,temp_c);</P>
$ [9 a* D& i* q% P4 i7 x" U' d<> printf("\n\n第一阶段问题为：\n\n");
3 d) z$ K- A% n, ~: ] show();
0 C/ N& o5 J' N( _/ a8 D* ^- Q  d! H3 \ step++;
( m/ H' j! D, W4 J printf("\n\n按回车开始第一阶段迭代");
# o! r6 L( S& d; M$ ~) x getchar();
getchar();
2 Q  L( u) _3 z4 u iterative();
$ {$ s# q! i' J/ Y  g4 [2 }! c: p if(status==-2) {
/ B! v( K' H" C9 C5 ]  puts("迭代超过限制次数强行终止！\n");
9 |3 }! k' [7 g7 u: L  A/ R  puts("\n按回车结束");
  getchar();
8 T$ i" G5 B: c' k- B  exit(0);
, a5 j2 c8 u. C }
( F+ Y+ e0 P: Y( |) F( Z output();</P>
, C# A; I" N5 M! A  @2 h3 i<> if(max!=0) {
$ Z) ^! ?& Q; N  puts("\n\n原问题无可行解。\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
}</P>
5 K3 o1 c3 @# r6 L<> //转换为第二阶段价值系数
exchange(c,temp_c);
6 R/ z* \3 ^4 W9 O- n: P //把人工变量列全设为0
for(i=1;i&lt;=num_ar;i++) {
: E# V# W: r! r! u, ?  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
}</P>
<> puts("\n\n第二阶段问题为：\n\n");
show();
puts("\n\n按回车开始第二阶段迭代");
" N0 i' M" C4 r getchar();
% \: K9 @  N- _3 m9 w) L4 [ iterative();
/ A3 L# w! s& t, n  \ switch(status) {
case 1:
/ u$ w7 [4 \/ {; R+ I* u- r5 d" Y  output();
  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
" u- c0 z2 k/ x  getchar();
  exit(0);
case 0:
) I- ?5 f" C* q& ?, @  puts("\n\n原问题为无界解。\n");
. d) ]& D: Z: R  puts("\n按回车结束");
  getchar();
7 c5 d4 o# ?* k  exit(0);
8 ?- ?+ y, D0 Q; o# @) D case -1:
$ a0 R" j8 ~" z* N) h9 G  output();
& V: @2 N4 p& U* \8 H  puts("\n\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
  v3 F' l  p5 v3 R8 n! G; ^  getchar();
5 G* p  H4 \; l" T2 R; U9 y1 H0 U( A  exit(0);
, m% M: S$ Y5 z8 I" j case -2:
  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
6 l3 N+ r& L: N& Y  getchar();
0 |* D# J5 h' j# X  exit(0);
8 z# C9 Z) a2 Z7 @6 R }//switch
}</P>
<>void banner() {
  I7 Y5 M" y$ s' r printf("\t\t****************************************\n");
printf("\t\t         单纯型法解线性规划问题\n");
printf("\t\t                         作者：Thunder\n");
; i! U, c8 F& e" k/ j$ {) E% h printf("\t\t****************************************\n");
printf("\n");
}</P>
<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
//为1的不显示系数1，-1系数只显示负号
& G) Q" B& n  y+ @2 m& G int i,j,k;
( ^- E2 H. F7 O+ ^0 q switch(step) {
case 1:
  printf("min z= ");
9 T; e+ k- a0 j/ K- V8 v4 i  printf("x[%d]",arti[1]);
- |7 D! ?/ J3 k" N- e2 a* T  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
  break;
& Q; u8 e8 M; P! ?% s case 2:
printf("max z= ");
' U: K' W( s8 N8 G+ `( z/ b printf("%lg x[%d]",c[1],1);
% L$ T4 C9 }* y( I, x for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
! k! v# [9 N; D  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
  else printf(" %lg x[%d]",c,i);
% g7 ]8 j/ b/ I+ ? }
break;
}</P>
<> printf("\nst:\n");
( \0 O+ c  ]4 R for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
7 p" s: U1 B( L5 _    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
  W1 T" s! L3 d( E, v( h    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
1 e) E5 B# l! P    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);
% M6 S9 a1 f3 ?4 J1 `/ o    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
( c- @- @6 J9 y! l    else printf(" %lg x[%d]",a[j],j);
1 r9 q( U2 g1 F9 N    k=1;
" H5 h1 k* N8 B   }
  h* p: s" J0 ]' y) b  v  }
# A4 C: s7 l7 X' |" e  printf(" == %lg\n",b);
}
printf(" x[1]~x[%d]&gt;=0",num_x);
) w" E! d" n! V) r}</P>
2 }5 J: ^4 @8 s7 C3 U( N<>void exchange() {
( o$ V9 C* s6 }/ E9 ` int i;
* S' ~1 o2 {/ \$ D1 G4 b double temp[MAX];
for(i=1;i&lt;=num_x;i++) {
; D+ o0 U* V3 D  temp=temp_c;
  temp_c=c;
6 t; {3 |4 t# ]3 h9 _0 q# D, n/ H  c=temp;
}
}</P>
) r- M; u/ B4 \; Q& m* l<>void create() {
9 A4 J, m0 o7 Q: h8 T, g8 r//输入方程组系数，每个方程输完后回显确认
2 q$ i9 m  b6 E' K7 ] int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
char confirm;

& n8 U& k' c2 r: B, N while(1) {
  printf("请选择：1、求最大值，2、求最小值：(1/2)");
! d) ?4 N- r6 S  scanf("%d",&amp;ma_mi);
  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
  else break;
}
# B" K) i/ K+ Q3 d) U" [* N
while(1) {
  printf("指定变量个数：");
; j6 ]+ K( a/ [! Q7 t! w. v, F  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);
1 x  O+ N1 F! r& Z9 x! y  }
  if(ma_mi==1) printf("max z= ");
+ j! S: Y$ E' }& _6 B3 e  else printf("min z= ");
  printf("%lg x[%d]",c[1],1);
; M! R3 ^# C( y  for(i=2;i&lt;=num_x;i++) {
   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
5 f# ?# ?: i7 l% }   else printf("  %lg x[%d]",c,i);
  }
printf("\n正确吗？：(y/n)");
# c. p* U1 [, A* U: h getchar();
confirm=getchar();
, m2 T$ z3 Z5 j* Y2 A  u  f if (confirm=='y') break;
else if(confirm=='n') continue;
}</P>
<> printf("输入约束方程组个数：");
- I: @/ J* s- B4 i6 ?+ ^* @ scanf("%d",&amp;num_st);
8 I9 W) i9 m. S4 j; j for(i=1;i&lt;=num_st;i++) {
  printf("st.%d:\n",i);
& U3 K& `; W+ \4 D1 P  while(1) {
   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
& a: h* a" x; y+ R" H+ Z   scanf("%d",&amp;re_st);
   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
   else break;
  }
6 \$ ]. v' X: d3 Q  Q  printf("输入技术系数：\n");
$ k" m2 s# V; o; g# ^  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
5 v4 y+ Q4 i, u7 g5 J$ s9 h! j   scanf("%lf",&amp;a[j]);
  }
) Y3 `& o  M* R, k1 z& }$ G, W  printf("输入资源拥有量：\nb%d=",i);
  scanf("%lf",&amp;b);
, F& z  m! `8 B9 [( }) A) o  
  printf("st.%i:\n",i);
6 L6 K6 [$ q- e9 [5 [  printf("%lg x[%d]",a[1],1);
  for(j=2;j&lt;=num_x;j++) {
& \8 h! |2 h+ N5 E   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
/ D/ m8 M" D+ S4 H( i* v   else printf(" %lg x[%d]",a[j],j);
, {( D1 \& h5 ]/ N  }
$ g' l3 ]+ R; e0 w# \  switch(re_st) {
   case 1: printf(" == %lg",b); break;
% }3 ]2 |8 m: ]! V   case 2: printf(" &gt;= %lg",b); break;
   case 3: printf(" &lt;= %lg",b); break;
  }</P>
( v9 E0 z; A' n. b9 z) h! d<>  while(1) {
   printf("\n正确吗？(y/n)");
   getchar();
  P& j) ~0 s+ `3 r" n, ~: F; ~   confirm=getchar();
9 a/ I' C: B3 i% x) c9 l   if (confirm=='y') break;
. ~/ E. i# P) L9 Z   else if(confirm=='n') {i-=1; break;}
  }
}</P>
7 y* [" Z$ H  X  w3 w+ |* C<>//显示输入的方程组
printf("\n原问题为：\n\n");
if(ma_mi==1) printf("max z= ");
3 Z' b  D# R3 ?! y+ I! `( h* p7 T else printf("min z= ");
$ W0 Y1 u: b. o+ L: {4 f printf("%lg x[%d]",c[1],1);
" [( z9 _7 A+ t+ J+ ~ for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
  r4 v( i* ^9 a& R/ z  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
" N% a" d/ D: _# l! B2 C  else printf(" %lg x[%d]",c,i);
}</P>
! i/ I/ w7 y4 \7 s% Y3 R<> printf("\nst:\n");
! Q3 F; a% Z0 @; T6 e8 j for(i=1;i&lt;=num_st;i++) {
- b, B7 V7 r. r1 W5 N4 P! Z  k=0;
  for(j=1;j&lt;=num_x;j++) {
* L# @" ?  E& g; \- a8 r   if(a[j]!=0) {
# q1 }7 e& w+ B1 _2 R# u    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);
; T" N. ]) C, ]; [" L    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);
! t3 k; G+ }# j" h) x( \    k=1;
& o: L0 C  D! z   }
  }
8 c; Z$ Y1 ]) f4 W( G+ J  switch(re_st) {
   case 1:
    printf(" == %lg\n",b);
    break;
   case 2:
, Z  b! E- q9 j8 u! z% R  h    printf(" &gt;= %lg\n",b);
8 h8 ?0 H3 ?8 q  @( j6 D% n5 z    break;
3 [' j3 Y$ s: g7 ]0 T  F! b   case 3:
    printf(" &lt;= %lg\n",b);
4 e; V9 o! q# e- Q1 ?& M    break;
  }
}
printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
<> tnum_x=num_x;
for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
  case 1:
1 X7 O7 m) Q6 J; f  case 3:
4 q! u/ V# y! ]" {   num_x+=1;
+ X* ?# z( L  D, O   break;
  case 2:
* A' B  R( `4 d$ @; M" Q   num_x+=2;
   break;
2 S4 i2 T6 R! J) f  }
7 j! B% V& j/ _% y: C1 X) e, i  e1 \ }</P>
<>//化为标准形式
if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
9 N% P- n9 ]" s4 `0 i" k% e# _ for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
. F' m, a! G8 T7 z1 ]7 ]9 E" i   case 1:
    num_addv++;
# e  y% }  v9 c; X1 _) m    num_ba++;
    num_ar++;
    c[tnum_x+num_addv]=0;
    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
- G/ V# b7 O( K' f    for(j=tnum_x+1;j&lt;=num_x;j++)
3 C+ |) C  T% G% M9 S$ _9 H4 u: a0 P. Y     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
' p1 U! B* H  [0 W    break;
8 E' E- `- W9 l  p9 [   case 2:
+ h6 E9 D4 b; B( S9 d5 F& o    num_addv++;
' d+ U. ]" G  \+ q  m    c[tnum_x+num_addv]=0;
7 C5 L# n6 i* J' n/ ^1 L- F. H    num_addv++;
    num_ba++;
    num_ar++;
    c[tnum_x+num_addv]=0;
    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
; @3 G; l, Y; S: p  H    for(j=tnum_x+1;j&lt;=num_x;j++)
5 n7 J8 E) \! y2 z3 p" N     if(j==tnum_x+num_addv-1) a[tnum_x+num_addv-1]=-1;
& o0 ]3 O5 {7 f. s- d% {8 d     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
    break;
   case 3:
    num_addv++;
" P7 I! O* a" ~8 S$ q0 B    num_ba++;
" w" H8 G* |* N! @8 u    c[tnum_x+num_addv]=0;
0 \& P: l; a6 N    base[num_ba]=tnum_x+num_addv;
5 D# x( t* O. S# A7 m) f    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# Q- y$ h; z% H8 {% n$ F    break;
5 y4 e" D- f4 [4 p& F7 G' q  }//switch
4 P: G6 J* {' R }//增加松弛变量、剩余变量、人工变量、确定基变量</P>
- n3 M& i* o- c0 ]# W! t  W1 Y) F<>//显示标准化后的方程组
printf("\n化为标准形式后：\n\n");
+ Y  M+ Y5 \' O$ D3 @- d3 c 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++) {
3 o' }1 C+ ?# {  C, Y, \  k=0;
* U" K& c) y, a- Q  for(j=1;j&lt;=num_ar;j++)
  if(i==arti[j]) k=1;
  if(k==1) printf(" -M x[%d]",i);
) l1 K3 I) Y8 B) n/ F8 Z& e7 n' N  else if(c==1) printf(" + x[%d]",i);
1 l) a4 T: p' f- {  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);
+ \  s/ v- _& I0 S# P9 R }</P>
+ {. r! D) a3 U  w8 |( Y<> printf("\nst:\n");
/ K+ H- K3 }3 v# B, F" l# i' K: l for(i=1;i&lt;=num_st;i++) {
  k=0;
3 m" D& U4 P( y1 q0 N* F1 Q  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
! I+ |: S$ y, H/ h8 g    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);
6 b6 z5 z4 }( d/ C    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
' Y, N! X+ j( q* `; K, ^7 \    else printf(" %lg x[%d]",a[j],j);
# ~% k2 f6 |& }3 @    k=1;
   }
% m2 T$ a, ~2 \+ Z! r  }
- u0 S/ v9 w3 Q7 l1 u0 j  printf(" == %lg\n",b);
* v& A9 L$ \$ @- n" [ }
printf(" x[1]~x[%d]&gt;=0",num_x);
. u( \/ z1 M  y( ^5 Q1 c& }! e}</P>
<>void iterative() {
6 D$ `5 f* G( G8 } int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
, U- @+ h/ k+ e' ^" L' Y int base_elem;
int base_out,base_in;
2 f- G4 v! T8 y8 i& d, ~ double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
<> printf("\n\n第%d次迭代：\n\n",stop);
. ]8 H; t$ G1 G- {. P$ G3 p- A for(i=1;i&lt;=num_st;i++) {
+ k) x! }* x* P/ |  printf("c%d=%lg\t",base,c[base]);
  printf("b%d=%lg\t",i,b);</P>
<>  switch(step) {
9 b; A3 L$ t" [  i& U3 D   case 1:
) I6 |( i) d+ d* u    for(j=1;j&lt;=num_x;j++)
% w2 b+ g& I& j    {
( f8 H9 P7 u  g( I* y     printf("a[%d][%d]=%lg\t",i,j,a[j]);
. h! e: N8 O+ N7 ~9 `* R3 N    }
    printf("\n");
- A( w1 N7 m. Z    break;
   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");
2 i7 e; Q: c, }$ i1 ]; ~    break;
9 {5 j! h9 g5 t# x4 f( M  }
}
2 F8 O: D4 l, D  L//求检验数sigma
for(i=1;i&lt;=num_x;i++) {
  sigma=c;
2 `* @4 j2 Z8 I" V5 D$ `% h% d  for(j=1;j&lt;=num_st;j++) sigma-=c[base[j]]*a[j];
$ k) ?5 W- j9 G% I' _' M/ c6 S  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
0 U5 M, P2 e  d- i1 ^  switch(step) {
   case 1:
    printf("sigma[%d]=%lg\t",i,sigma);
    break;
   case 2:
& S( |2 X5 _  {    k_a=0;
    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
: a% n0 V$ _, W! Q) K; b: f9 p    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
6 e) c6 H% @# c. Z$ J    break;
  }
}
putchar('\n');
//检验检验数sigma是否全小于等于0
5 `: f1 w( l2 M9 a$ K9 j( x' C8 W$ I& h k=0;
for(i=1;i&lt;=num_x;i++) {
  if(sigma&gt;0)
  k=1;
}
* D' S3 P. f0 P& u if(k==0) {
! W, e1 `* I" D7 g' f  //sigma是全小于等于0时，检查是否为无穷多最优解
2 n4 C0 m7 w9 Y& [  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 x% i2 o" C6 U: n( ~: ?0 Y3 F. J    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
   if(k_f==0) {status=-1; return;}
   }
4 C6 l3 A  `1 J( U3 e8 S, q  d( `  }
  status=1;
! k# d5 z( M5 d1 g  return;
) O1 W6 t4 E* m9 A! G7 Q+ z' } }
: x* Q  M* e( c$ V+ C- L: ]//检查是否为无界解
$ H3 y2 I6 e! ^5 E for(i=1;i&lt;=num_x;i++) {
+ D, T0 q% j& l  k_f=0;
( r, _6 r" s0 g9 A2 `; f, {  if(sigma&gt;0) {
  p# m! d: _$ k" R+ Q, E: g; P   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
   if(k_f!=1) {status=0; return;}
  }
}</P>
<>//确定换入变量
" w4 g: g/ O- f for(i=1;i&lt;=num_x;i++) {
  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
: ?2 n* \) e, L( G  if(k==0&amp;&amp;sigma&gt;0) temp=sigma-1;
( _: T" T- k- H9 t, y }//temp赋初值
+ e8 c; Z) ^5 s5 C8 H% W1 y: b 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)
   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
    base_in=i;
: k6 J% E: u0 \+ Q# M+ y    temp=sigma;
   }
& D/ v8 G8 f" _3 Z7 `* P  t; @ }</P>
. ]$ k: p4 }  t6 C<>//确定换出变量
; |* J' O8 z7 o. ~ for(i=1;i&lt;=num_st;i++)
  if(a[base_in]&gt;0) {
6 Z- ?! Q- @- d+ E   temp=b/a[base_in]+1;
   break;
  }//temp赋初值
; W9 i# K  o6 p for(i=1;i&lt;=num_st;i++) {
  P; i) g# ]2 _  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
2 I# ~+ i+ @* Y9 ]: Y   for(j=1;j&lt;=num_ar;j++)
    if(base==arti[j]) {
- o9 i, K$ u) x" [: l     base_out=base;
. w& ~% G; g/ u. f) G     base_elem=i;
! h( Z' i6 Q% ~! T* D. `     temp=b/a[base_in];
     break;
' n' m4 m) \# J8 U    }
  }//人工变量优先换出
  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];
2 d2 a! p3 G, X8 f  }
}</P>
<> printf(" 基变量：");
  c+ L" `8 p& R9 T. L for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
' b3 b/ H6 c3 e2 A. `3 l printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
//基变量变换，进行新方程初始化后迭代
" h( A, j. G* v/ o0 R! v4 H for(i=1;i&lt;=num_st;i++) {
  if(base==base_out) base=base_in;
2 R/ ]/ y9 e7 j' F }
+ ]% l' w1 L5 ^: b7 S1 e% _//初始化主元素行系数
value_be=a[base_elem][base_in];
( ?& a" D7 r. u* R& {5 j/ N b[base_elem]/=value_be;
1 ~1 Q5 [) A: d3 N& }; _2 B for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
6 S) d! P+ P& C& Q$ q- d" c6 g* F<> for(i=1;i&lt;=num_st;i++) {
8 C5 B0 L, G  l" W& z  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;
  b3 L+ k; _; u4 l  }
}
stop++;
  V" }" ~7 n7 F0 r9 ] if(stop&gt;STP) {status=-2; return;}
" K3 N& X3 d/ v+ G iterative();
}</P>
<>void output() {
int i,j;
) r& V! }, P& Y2 C0 e2 @: t1 g" P double X[MAX];
printf("\n结果如下：\n");
: k. v6 Z# A" M8 J printf("\nX=(");
for(i=1;i&lt;=num_x;i++) {
  for(j=1;j&lt;=num_st;j++)
# y. D7 s# Z/ N! G/ ~0 ~5 ~) N0 Y/ s: P5 X  if(i==base[j]) {X=b[j];break;}
0 G, x2 {0 {7 S( e% Z4 Q5 }# o  else X=0;
  printf("%lg ",X);
}
' C1 z6 \  x2 x2 _. v. y, ^ printf(")");
9 j6 Z' M" S4 W5 C for(i=1;i&lt;=num_x;i++) max+=c*X;
+ g/ V' k: W+ P, U: O( h if(ma_mi==1) printf("\nMax z= %lf\n",max);
3 {- }4 o: N( C* g9 T else printf("\nMin z= %lf\n",-max);
2 `( m. r0 }4 m; L7 a}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了