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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
) V4 [2 S6 z% Q4 u#include&lt;math.h&gt;
#define X 5
  v1 V. V" f7 {: ^! m#define Y 7
void xi_max(int *m2,int *mn1,float *c,int *is,int *ir,int *j0,float (*a)[X][Y])
{
* {* A" c- W: |  l/ b  int j;
4 [% X0 T5 @: r. I, j  *c=0;
- @# W4 t; j; F8 O3 {  for(j=1;j&lt;=*is;j++)
5 z% C+ d1 `: _( j4 ?    {
" x1 Q( Y( n5 b4 n7 ^% _      if((*a)[*ir][j]-*c&gt;0)
    {
4 ^, k7 @# ^/ s6 u      *c=(*a)[*ir][j];
% c* n8 z5 m2 z& i      *j0=j;
4 h0 {; e3 \2 d0 f( F8 I+ c    }
     }
}</P><>      /***************** 参数说明 **********************/
0 \7 P/ m( r7 g  ]$ e' Q      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
      /*   m2-m+2整个变量,(*a)[X][Y]存放初始数据       */
      /*   (*k)[]存放基变量脚标,(*x)存放基变量最优值   */
7 k6 S- C# A, D& [1 Q' u) O      /********* <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],
5 r3 {( L; S: p( H% p       int(*k)[],float(*x)[])
{
  int i,m1,mn,j0,i0,j;
6 F2 w: q' _" l( Y  float c,g;
  m1=m+1;
  mn=m+n;
  for(i=1;i&lt;=m;i++)
    (*k)=n+i;
+ q2 N% V9 u) q5 k% n3 m  jleap1:
  if(l1-1==0)
: X- ]0 @$ O; h+ M( r8 }' h6 a) d4 d     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  else
# [& B4 I+ E+ T6 @& N  ~! s/ H# H     {
5 [1 u# ~: h3 k1 S/ C# @$ `leap2:   if(l1-50==0)
" ]* q" X' ]$ Z# h  f8 c, k$ F      xi_max(&amp;m2,&amp;mn,&amp;c,&amp;n,&amp;m2,&amp;j0,a);
       else
/ u  ]# e# j; V! x      xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m2,&amp;j0,a);
' E! B1 g2 i' N% G& l     }
3 s/ G4 M; G4 _0 U4 J! f. P% T   c-=1e-8;
; W  Q" w: T% {8 \! y) ~( N! H   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
& J; s, |" r- V      {
    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
! C6 j: G  a5 c  k& m$ M' t    return(-1);
. D: h: V% o; a$ }1 O% Q$ [. I      }
   if((c&lt;0)&amp;&amp;(l1==1))
$ F+ \. W  P$ I       {
     l1=50;
     goto leap2;
! d" V: O# r% p       }
    if(c&lt;=0)
6 i7 M( f5 t3 h* e7 D       {
      for(i=1;i&lt;=m;i++)
4 j* _8 B. a3 k$ x; c7 O) D+ K% Q         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
      for(i=1;i&lt;=m;i++)
         printf("\ni=%d,x=%f\n",(*k),(*x));
3 c) A+ ]* s$ y4 [* z+ b  m) t5 O& s      printf("\nF=%f\n",(*a)[m2][mn1]);
      return 0;
5 L+ L; Q  V8 L( ?5 V    }
) {$ A( U$ R3 z* E5 i. x1 A$ M: A1 T: c     c=1e8;
     for(i=1;i&lt;=m;i++)
& _* Q$ Z; s5 l4 B9 h/ U7 J6 f       {
/ N$ ^. G% i1 K6 V     if((*a)[j0]&gt;1e-8)
, o8 d' `8 H& |: Y5 F0 s2 p+ \3 T        {
          g=(*a)[mn1]/(*a)[j0];
' y! s1 @+ Q5 b* \3 M; C          if(g-c&lt;0);
) i6 \" s8 {) s/ l         {
           c=g;
           i0=i;
2 F- [  F$ ]0 i          }
3 k! K% E! V% V6 q7 i& }1 D         }
    }
1 F, f- l" v2 L0 N( l+ `      if(c==1e+8)
9 _3 q5 {/ x5 s- l$ p; I3 h$ W    {
: I# v7 D* I9 y# d( \      printf("\n\t*******Lp no solution********\n");
9 U% [; s' F8 j9 N1 H8 ]& A      return -3;
  O3 \. ]9 G  [; J# p+ m     }
/ n, K% M' H8 h6 B2 P6 t4 S9 s       (*k)[i0]=j0;
: _6 G! _, F% n- D8 \# ~       for(j=1;j&lt;=mn1;j++)
# J( B2 T6 `9 n( o  ]1 Q) q# H3 L      {
+ V. b2 f1 r* s6 }- G6 h" a        if((j==j0)||(l1==50)&amp;&amp;(n&lt;j)&amp;&amp;(j&lt;mn1))
5 _4 j- l) p  `( o' ~! E9 D           continue;
         g=(*a)[i0][j]/(*a)[i0][j0];
4 Q. ^# K( ~. t4 r3 a6 ?) [; _           (*a)[i0][j]=g;
' t3 G# n7 s. @: r. U; G9 P, m* u         for(i=1;i&lt;=m2;i++)
           {
. A+ W$ P/ B4 b! G% n+ @  T: e6 u         if((i==i0)||(!(l1==1)&amp;(i==m1)))
- c+ \8 i% d" `4 n2 y            continue;
         (*a)[j]=(*a)[j]-(*a)[j0]*g;
           }
       }
       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
       (*a)[i0][j0]=1;
       goto leap1;
}  
main()
, e! C8 Z1 v+ `+ Y{
  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},
7 l. G$ o- j4 j( _                 {0,0,0, 0,  0,0,0},
' _# N# u# `- U                 {0,69,0,144,0,0,300}};
/ }) q1 M& i3 K2 y$ A1 p1 D" \2 |  int k[3];
# y* K5 b. h4 ~  float x[3];
  clrscr();
  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
}
1 [" b1 h) |) r# }% n5 b</P>

plgatc

<>/*************************************************************************
! b/ h9 H+ h9 o+ I3 I' _                    单纯型法解线性规划问题（两阶段法）        
& }- c' r2 N9 X; z
编程环境：VC++6.0     
方程组输入说明：
' }* f% F, G" V) L 变量非负，按提示输入相关参数。
+ {' L. e* J, R$ ]# e: c. M *************************************************************************/
#include &lt;stdio.h&gt;
#include &lt;stdlib.h&gt;
#define MAX 100
& k% D) p8 p5 d#define STP 100</P>
<>int stop=1; //迭代记数变量
% j  R) z5 `: n1 V3 A) wint status; //iterative迭代返回值：1唯一最优，0无界解，-1无穷多最优解 -2迭代超过限制次数
0 x+ W4 s$ w$ H7 U- Y0 j9 \1 sint 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]; //人工变量下标
3 @$ X1 W  N! y- h4 O7 u/ vint base[MAX]; //基变量下标
' _1 ^% D- ?9 i# B9 |" Kint ma_mi; //1为求最大值，2为求最小值</P>
<>void create(); //建立方程组
void iterative(); //单纯型法迭代
5 r: S. V3 [* U' n& o4 j, uvoid output(); //输出结果
void banner(); //打印程序标题
void exchange(); //交换两阶段价值系数
: M, Q5 S& \- K* gvoid show(); //输出方程组</P>
3 l! z; B' [5 ]- _' J<>void main() {
, k$ _$ I0 e1 h int i,j,k;
banner();
create();
//保存原价值系数，转换为第一阶段价值系数
6 ^, |2 \4 C1 g5 r# S2 l for(i=1;i&lt;=num_x;i++) {
, |, b3 K/ r& k- x  j% h5 s  k=0;
# }$ j7 @) ~1 U- r2 L6 x6 P  for(j=1;j&lt;=num_ar;j++) if(i==arti[j]) k=1;
  if(k==1) temp_c=-1;
  else temp_c=0;
- h! p" \6 Z/ n4 X3 ^1 Q }
: P1 |3 e0 r* q exchange(c,temp_c);</P>
7 I9 G8 o" d# P6 f0 W3 d<> printf("\n\n第一阶段问题为：\n\n");
show();
step++;
printf("\n\n按回车开始第一阶段迭代");
4 v. m5 {$ w3 n! M( |# W. Q/ T5 j* C! u getchar();
# w% b. d9 E. g) W( @' L# v getchar();
iterative();
2 \9 @0 c1 ]7 `% W9 s! C if(status==-2) {
- k( x6 N& b2 g5 c. K7 }  puts("迭代超过限制次数强行终止！\n");
. h* p0 p' d: l  puts("\n按回车结束");
, ]% `4 c- z) A  getchar();
# O' z( r" i. {( A  exit(0);
* f" P! L, N# |# c% m$ X }
output();</P>
# G* k# g( _) P- c) @& u<> if(max!=0) {
  puts("\n\n原问题无可行解。\n");
  puts("\n按回车结束");
  getchar();
; b0 N/ ~" K7 V5 k- i! _# G  exit(0);
- s1 x5 g* q, B; B }</P>
<> //转换为第二阶段价值系数
exchange(c,temp_c);
/ h5 Q% L7 \1 F3 l% |6 ~4 { //把人工变量列全设为0
for(i=1;i&lt;=num_ar;i++) {
  c[arti]=0;
, L' C; P  X4 E  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
}</P>
<> puts("\n\n第二阶段问题为：\n\n");
. u7 c4 O, q- N0 ]$ a show();
" L5 H/ ^: u! N) N* G! f puts("\n\n按回车开始第二阶段迭代");
getchar();
8 q# a+ W4 W# K/ ?" w" H iterative();
switch(status) {
: Y) q3 O6 [$ V" ~6 t case 1:
  output();
  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
9 X! P5 Z; @, u5 V  getchar();
7 D& Z! f" Y" ~& s4 I  exit(0);
' b) R1 r5 P* |4 O+ A, C: b6 b1 v case 0:
$ x6 @! ?; O+ o  puts("\n\n原问题为无界解。\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
case -1:
  output();
  puts("\n\n原问题有无穷多最优解。\n");
8 Q4 R* m* x, r- f. d7 {  puts("\n按回车结束");
7 V* j9 X) J* Q$ E) f  getchar();
  exit(0);
case -2:
3 A# V: o( i' e2 l6 \. S, ]" |. G  puts("迭代超过限制次数强行终止！\n");
  puts("\n按回车结束");
  getchar();
  exit(0);
5 Y" {8 q/ Y2 D- B9 e9 m! B }//switch
+ V4 C$ f6 s% T+ G. b' g}</P>
4 Y  n9 T6 H/ c% ]+ l- ^2 O% x<>void banner() {
printf("\t\t****************************************\n");
% F7 [/ V+ \) \, o9 e: q printf("\t\t         单纯型法解线性规划问题\n");
7 `0 V8 z# {8 y% p% `6 S printf("\t\t                         作者：Thunder\n");
" y) Y, h/ r6 v: C4 x: v printf("\t\t****************************************\n");
printf("\n");
* b/ _4 U/ ?% T- X}</P>
) w1 S0 A# b& p) x, K) F. |* h<>void show() {
//对方程组以自然的格式输出，系数为零的x不显示
9 L6 }+ d+ B2 F$ o//为1的不显示系数1，-1系数只显示负号
int i,j,k;
switch(step) {
. F( \% T1 w3 I. b case 1:
4 m6 P; T+ m+ H: |: V- ?  printf("min z= ");
  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++) {
  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);
}
break;
}</P>
<> printf("\nst:\n");
; _5 h' M( W; A for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
) k+ q9 Q1 q1 ^( S7 M    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);
" S% E, a) }! z  B    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);
3 G0 T: c! K2 X  d) b+ F* }6 |4 Y    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;
5 I9 B, z% P6 z% J$ u   }
9 p+ Q- c  e- G  }
  printf(" == %lg\n",b);
3 n+ b) @. M4 U- z7 g }
printf(" x[1]~x[%d]&gt;=0",num_x);
' `, L6 r9 H# g- a) H- |) K}</P>
<>void exchange() {
int i;
) t, p: k" ~5 G' E  ~ double temp[MAX];
for(i=1;i&lt;=num_x;i++) {
  temp=temp_c;
  temp_c=c;
  c=temp;
}
& ?+ T  H; n6 x" Y) v}</P>
% B0 C1 ?+ i( ?  t8 p& }<>void create() {
//输入方程组系数，每个方程输完后回显确认
int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
char confirm;
7 M5 J8 u, o9 y# V! ~
0 ^7 T+ F6 ]3 o, z" a# ~ while(1) {
5 X& e4 _# i7 X) P1 I$ h  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  |9 N; R( k% g1 w3 m0 q  scanf("%d",&amp;ma_mi);
  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
; w- X# X& a3 T" h& A  else break;
4 x2 E9 `  Z* c1 k$ ] }
+ _% M8 X6 b3 _
while(1) {
( b: {& J$ i# K" H3 V6 q  printf("指定变量个数：");
  scanf("%d",&amp;num_x);
  printf("输入价值系数c1-c%d:\n",num_x);
7 X+ S' ~- N/ M, A  for(i=1;i&lt;=num_x;i++) {   
' c$ c! `6 x& B" [' g   printf("c%d=",i);
   scanf("%lf",&amp;c);
  }
  if(ma_mi==1) printf("max z= ");
  else printf("min z= ");
; [$ R3 y" _+ G- Q8 m  printf("%lg x[%d]",c[1],1);
) _7 c! j' a8 ]% A  Y7 ^4 k+ V  for(i=2;i&lt;=num_x;i++) {
   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
. t+ c$ W1 H" y& _2 ?- u  m1 q   else printf("  %lg x[%d]",c,i);
3 D) B) x8 p" r( c, m, L2 A. S  }
6 s& c9 N; y4 l' l3 F8 o printf("\n正确吗？：(y/n)");
( t6 I) N5 V/ i/ T( Q' {8 j( k getchar();
confirm=getchar();
/ P8 G% L& `. u# W! ]$ Y if (confirm=='y') break;
  S* k/ Y& `& V; m6 R9 h- b else if(confirm=='n') continue;
}</P>
<> printf("输入约束方程组个数：");
# u# O4 z  i! }) z! p0 u scanf("%d",&amp;num_st);
1 `7 k& w2 ?5 e: O' |( m: s& ^ for(i=1;i&lt;=num_st;i++) {
  printf("st.%d:\n",i);
1 j7 B* {/ d5 n9 d( D# o. B  while(1) {
   printf("请选择：1、==，2、&gt;=，3、&lt;= ：(1/2/3)");
   scanf("%d",&amp;re_st);
& q0 c* t0 \- x) l   if(re_st!=1&amp;&amp;re_st!=2&amp;&amp;re_st!=3) printf("输入错误，请重新选择。");
8 N8 ~7 v5 ~7 h: M' b( g5 E   else break;
. f/ d% N7 c  A, w. f  }
  printf("输入技术系数：\n");
  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
   scanf("%lf",&amp;a[j]);
  }
  O2 J& L% L+ z5 P+ C$ ^" y1 N  printf("输入资源拥有量：\nb%d=",i);
+ a1 K, H" t4 O0 }( \8 [  scanf("%lf",&amp;b);
7 P% w. _8 F/ B7 u1 B8 }! V; i  
1 e7 P3 y( V  J: O0 x) b  printf("st.%i:\n",i);
0 v0 _9 ~/ |) `- a) ^: y/ r" M  printf("%lg x[%d]",a[1],1);
  for(j=2;j&lt;=num_x;j++) {
5 |9 B' O+ y& R! k# i) b/ L/ t" z   if(a[j]&gt;=0) printf(" +%lg x[%d]",a[j],j);
   else printf(" %lg x[%d]",a[j],j);
# T: N2 B% ~+ ^' s, O' Q# L# x  }
9 I2 T. Z$ |" x3 s) a, l/ l, {  switch(re_st) {
* t( M# N4 i2 T- I' P2 c5 S   case 1: printf(" == %lg",b); break;
   case 2: printf(" &gt;= %lg",b); break;
! z8 ^& u( x8 i+ \   case 3: printf(" &lt;= %lg",b); break;
# b+ ^  |6 Y9 R$ \+ t  g5 ^  }</P>
4 p0 T% M7 P8 v<>  while(1) {
   printf("\n正确吗？(y/n)");
4 ~1 W& {) u& w9 a: Y" Z* Q   getchar();
   confirm=getchar();
   if (confirm=='y') break;
   else if(confirm=='n') {i-=1; break;}
! G% z- t5 A0 M( y  Z, T  }
5 Z  i) b; K" v" v8 n4 p- V }</P>
<>//显示输入的方程组
( X+ S  }* c2 N* r( v) v printf("\n原问题为：\n\n");
if(ma_mi==1) printf("max z= ");
, m* h  E( J* W8 P( [ else printf("min z= ");
+ m+ T7 F* J5 b! C: }8 K6 L; ~ 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);
+ ^) G5 C1 W( r4 A" M& N0 T3 e5 S  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
  else printf(" %lg x[%d]",c,i);
. B* G' o; ]) ]5 k }</P>
<> printf("\nst:\n");
5 @( U7 \. F7 D6 K for(i=1;i&lt;=num_st;i++) {
  k=0;
  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
; s7 Y+ P  n4 [    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
  P  S5 ?2 g( }+ r% }. d* w; {: n8 L    else if(a[j]==-1) printf(" - x[%d]",j);
" I# ?' n8 r6 O+ D- 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);
6 @' h0 S# ^# p    k=1;
   }
  }
" ]5 y2 u0 l/ G# ?6 k  switch(re_st) {
' L8 S; ?" Q5 }9 U; y   case 1:
    printf(" == %lg\n",b);
5 b& ?3 o0 }1 e  U) Y- u' q    break;
   case 2:
    printf(" &gt;= %lg\n",b);
    break;
   case 3:
    printf(" &lt;= %lg\n",b);
" w, V4 |9 i2 k7 b$ O( x    break;
  }
}
5 n( p3 z% w' T+ S2 o8 N' Q& B) C, C4 J printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
6 l* d, R9 c9 }0 \( O<> tnum_x=num_x;
( R, a# C$ D! I' W' x2 i for(i=1;i&lt;=num_st;i++) {
3 ~9 @% ~6 p2 \! B  switch(re_st) {
) E" s9 }+ }# {! O; Z  case 1:
  case 3:
   num_x+=1;
   break;
  case 2:
   num_x+=2;
   break;
0 H1 S- w6 K7 d1 p: Z  }
}</P>
7 O2 z9 H5 H# A7 `0 ~' N<>//化为标准形式
% Q, t, ^5 U* q% A if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
, V; J" [. ?- s4 ]+ u for(i=1;i&lt;=num_st;i++) {
& Z' E3 k; o* j* K! q  switch(re_st) {
1 }. x+ `3 \7 r; m   case 1:
    num_addv++;
9 c' f4 L1 v6 h( V3 b( ?" Y    num_ba++;
    num_ar++;
; R2 L/ ~" q& R( Z0 d    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) a[tnum_x+num_addv]=1;
     else a[j]=0;
$ R, E+ K0 b" T. Z6 O+ X    break;
- H& v4 \" _. J. [4 `1 i% i2 ~. R   case 2:
    num_addv++;
$ ?# |' |- p) P1 }1 l  z2 |" i    c[tnum_x+num_addv]=0;
    num_addv++;
    num_ba++;
    num_ar++;
4 y! @- ^' x$ }( Q    c[tnum_x+num_addv]=0;
    base[num_ba]=arti[num_ar]=tnum_x+num_addv;
# @+ l" T3 U3 c' o" y6 d# i    for(j=tnum_x+1;j&lt;=num_x;j++)
0 {* u$ Q2 W9 l2 a& |% C     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;
; z6 u$ D5 F9 d1 I  [6 X* T     else a[j]=0;
    break;
2 k: f" o6 `9 t0 G  q8 E& _( w   case 3:
$ S( D. R4 T1 P# K, j2 D& ]5 M! E( @    num_addv++;
3 I4 P, Y- |+ d5 P    num_ba++;
    c[tnum_x+num_addv]=0;
. t9 P3 d' I% L0 a; S    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;
3 g. |. e8 i- l* e; q) \     else a[j]=0;
    break;
  }//switch
' y( h$ F  d4 e5 p9 u* M  c7 e. u }//增加松弛变量、剩余变量、人工变量、确定基变量</P>
# U4 F- D' ^6 ?0 N+ \3 E6 d<>//显示标准化后的方程组
printf("\n化为标准形式后：\n\n");
if(ma_mi==1) printf("max z= ");
0 X, e6 ~, e# w( D  V: S) d else printf("max z'= ");
printf("%lg x[%d]",c[1],1);
2 e4 Q& U; l, P; x7 E$ H" Y! A for(i=2;i&lt;=num_x;i++) {
7 m4 w; M& F: L  O$ g: q  k=0;
  for(j=1;j&lt;=num_ar;j++)
: k) J& w( K4 H( y! l  if(i==arti[j]) k=1;
  if(k==1) printf(" -M x[%d]",i);
  else if(c==1) printf(" + x[%d]",i);
' s  y! ~3 n5 \7 r) r( _, U: C+ I  else if(c==-1) printf(" - x[%d]",i);
+ F' t, e# E3 z, w( o$ z! D' B  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
  else printf(" %lg x[%d]",c,i);
$ j, p5 Q. C6 F3 w2 N4 P3 U* L }</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
# K! L& a' W9 Z) K! b8 p  k=0;
7 K7 i& |+ J3 h% ~/ o! F  for(j=1;j&lt;=num_x;j++) {
1 K1 n) p) b2 l& \* j& S   if(a[j]!=0) {
    if(a[j]==1&amp;&amp;k!=0) printf(" + x[%d]",j);
* j# W/ a3 N8 R* j7 G% O    else if(a[j]==1&amp;&amp;k==0) printf("  x[%d]",j);
# F! a4 }9 ^- ^- h1 w    else if(a[j]==-1) printf(" - x[%d]",j);
# d0 c) e* ]# i- j6 T    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;
- H5 H+ T8 t, i8 o$ X; Z* T   }
& s8 b  {3 u0 }$ l7 |$ w  }
  printf(" == %lg\n",b);
9 T! z' m+ x. c3 X( D: P( A }
  @: c, `- E" B+ M printf(" x[1]~x[%d]&gt;=0",num_x);
}</P>
0 E1 [" O9 r3 Y0 n3 j<>void iterative() {
int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
5 Z$ A( ]5 h" \* Z int base_elem;
int base_out,base_in;
double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
<> 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>
, C% p( W7 L! H0 C: M$ z: |6 `# i4 Q<>  switch(step) {
4 [4 @! y/ l3 a7 \) Y& C   case 1:
2 |- j  D' X/ X% ?# G. u    for(j=1;j&lt;=num_x;j++)
    {
     printf("a[%d][%d]=%lg\t",i,j,a[j]);
) ^1 ]: A$ X7 b; n8 ]    }
; c3 R; E/ P# S! x2 p    printf("\n");
3 m! Q( V. ?/ d- ]    break;
   case 2:
    for(j=1;j&lt;=num_x;j++) {
2 m* |3 v2 ~% c: m: Y     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;
  }
}
* c0 [' h: v2 e8 K//求检验数sigma
for(i=1;i&lt;=num_x;i++) {
  G6 I- t4 [8 ~2 T% Q$ Q# m  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;
  switch(step) {
+ H6 Y+ o! A9 f) S   case 1:
    printf("sigma[%d]=%lg\t",i,sigma);
    break;
' q/ W6 ?5 _- \' U   case 2:
$ g, U/ m" K; H2 B  {0 D    k_a=0;
    for(l=1;l&lt;=num_ar;l++) if(i==arti[l]) k_a=1;
7 |9 m* e. Q% i9 H  A- Y    if(k_a!=1) printf("sigma[%d]=%lg\t",i,sigma);
8 H3 B& r: L$ A# P9 f. M( u    break;
: R' J2 ^+ p* S5 A9 K" |  }
}
putchar('\n');
//检验检验数sigma是否全小于等于0
$ S) d/ ^' K! o k=0;
for(i=1;i&lt;=num_x;i++) {
/ j. K: c/ A5 g6 h/ e1 t  if(sigma&gt;0)
  k=1;
% s+ q. s2 l( H# f, P1 ?" ?/ p }
9 O& S5 Y! {) V4 T8 m; E0 ?; f( g8 V if(k==0) {
1 E- Y( |6 W% ^+ k( j3 |6 B+ [( ^3 K; o  //sigma是全小于等于0时，检查是否为无穷多最优解
+ L  m' @5 c5 p- O' }: X  for(i=1;i&lt;=num_x;i++) {
/ f: t% i% B# \2 O+ D- @9 e( F   k_f=k_a=0;
   for(j=1;j&lt;=num_ar;j++)
   if(i==arti[j]) k_a=1;
( M9 D9 z' M* i% `4 W9 c   if(sigma==0&amp;&amp;k_a!=1) {
: }; d# F6 \  m0 s* R    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
   if(k_f==0) {status=-1; return;}
   }
  }
  status=1;
  return;
3 u5 f" w  ~7 \. _# u  h0 `0 X }
7 ~& ?. P* _  q, n8 D! m1 [//检查是否为无界解
) |/ E+ p/ n" @$ F$ n2 L  |+ Q$ I for(i=1;i&lt;=num_x;i++) {
  k_f=0;
: [. D! g8 y7 d' L% B  if(sigma&gt;0) {
' [6 v. m* M- f8 g+ _6 q   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
# l0 C$ `" W# ]   if(k_f!=1) {status=0; return;}
2 e/ U! D3 d) O+ J6 W  }
}</P>
<>//确定换入变量
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&amp;&amp;sigma&gt;0) temp=sigma-1;
3 t+ ]) g4 W: F/ S+ u* W0 k }//temp赋初值
for(i=1;i&lt;=num_x;i++) {
; Z" d" V) [) t2 ?! _, M7 J  k=0;
  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  if(k==0)
3 e3 A& \0 J) X% C( J2 m( h  B: H   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
' l' m# e9 g. C1 E* e# f! D" X    base_in=i;
    temp=sigma;
   }
}</P>
. b: w# [# U. F# t* m<>//确定换出变量
; H0 |7 K/ x: U! v% a+ D+ q for(i=1;i&lt;=num_st;i++)
3 Z- @5 k3 P8 _0 h  if(a[base_in]&gt;0) {
   temp=b/a[base_in]+1;
   break;
( ^5 {' p" u. L1 x2 S  }//temp赋初值
for(i=1;i&lt;=num_st;i++) {
3 K# l% S9 Z) W8 ~5 V0 {/ |3 x  if(b/a[base_in]&lt;=temp&amp;&amp;a[base_in]&gt;0) {
$ r  |: M* |: A. W% D% m   for(j=1;j&lt;=num_ar;j++)
    if(base==arti[j]) {
     base_out=base;
0 F- Z8 e: P  X     base_elem=i;
! l9 Q3 q" i+ v' ]5 k     temp=b/a[base_in];
' Q- Y( |) U+ z4 K7 ~. h     break;
: Y; x! n( f- K  I    }
  }//人工变量优先换出
: ~7 d1 D2 K9 Q  E7 @0 [  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
   base_out=base;
   base_elem=i;
; F. w' O# G  @7 k$ \: G   temp=b/a[base_in];
( t3 r1 p* H# J) \" h; g9 x" B  }
4 B. e" G" H9 N" \ }</P>
+ x" i7 h2 A. D$ e: q<> printf(" 基变量：");
0 n1 e1 I! E# }; H# J' {! |: M1 K for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
' V3 |  v2 `* Y5 P9 | printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
  \) D' H0 T( n' z7 m6 d//基变量变换，进行新方程初始化后迭代
for(i=1;i&lt;=num_st;i++) {
  if(base==base_out) base=base_in;
$ U$ ?. |0 U# ~# O1 i, {4 N0 b! m! ] }
//初始化主元素行系数
2 N/ _, t% y2 n: f5 g4 z value_be=a[base_elem][base_in];
' W$ r) V: ~5 N4 `4 b" e7 s4 t9 u: W b[base_elem]/=value_be;
8 F+ H% @8 a3 Y8 {7 O. j' a8 p 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) {
   b-=b[base_elem]*a[base_in];
- C$ t" ]. z+ Q- e; x5 V   value_be=a[base_in];
' q1 R+ \5 ]5 o   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
  }
$ ]4 `0 O0 r. t4 z' E: v }
stop++;
4 L* f& |6 f& }& l3 D if(stop&gt;STP) {status=-2; return;}
iterative();
}</P>
/ \2 z5 K- }1 ]1 i$ J3 H9 j<>void output() {
int i,j;
% f6 L% y$ A" M2 Y3 B double X[MAX];
. [( U& y8 ]- |" A printf("\n结果如下：\n");
' X4 u5 L5 y. _! \8 S! d# ] printf("\nX=(");
* o; [% p7 l* T5 V1 A2 { for(i=1;i&lt;=num_x;i++) {
% J" ?$ P# H7 c0 E  for(j=1;j&lt;=num_st;j++)
' J  J: \% F: E) G. d3 G3 h# B5 ~7 e  if(i==base[j]) {X=b[j];break;}
# k  E8 I# I2 k, T  else X=0;
  printf("%lg ",X);
}
' n0 H) b! \$ _7 G* k printf(")");
- M0 g$ f. C1 F% u 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);
! U8 X' ~- H* _3 q( [}</P>

sg47

好详细！

yqm10507

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

tw1982

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

tw1982

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

winzipftp

是你吗？thunder？帅呆了