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

lauqt1980612

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

lauqt1980612

最好是C语言的!

julysea

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

candid

<>不知以下这个程序满不满足你的要求</P><>#include&lt;stdio.h&gt;
#include&lt;math.h&gt;
#define X 5
4 y9 O* }2 K1 r#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)
0 `6 C, g3 Q/ M* e: g    {
' S/ R. w" i9 o      *c=(*a)[*ir][j];
& ^  x7 r4 f  m5 ]- n      *j0=j;
    }
     }
}</P><>      /***************** 参数说明 **********************/
2 P% \5 i- b, @" P3 q* S      /*   m_约束方程个数(基变量个数),n_非基变量个数   */
7 H3 K5 a! c1 r3 E      /*   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],
8 _/ D" A7 q: M! U8 I       int(*k)[],float(*x)[])
. K/ R( @' ~9 A$ }  g* S& Z/ z1 D: D{
  int i,m1,mn,j0,i0,j;
  float c,g;
- C* l  r# u, s1 [) m* ~) A* {  m1=m+1;
  mn=m+n;
* Y  D8 S' x, ^, ]  for(i=1;i&lt;=m;i++)
( }4 _/ S1 o+ k+ G    (*k)=n+i;
9 X# N/ d% I$ ?7 |. n. Y, Qleap1:
0 Z4 S5 N0 _/ W  if(l1-1==0)
     xi_max(&amp;m2,&amp;mn1,&amp;c,&amp;mn,&amp;m1,&amp;j0,a);
  else
     {
, K' n# E  P* L2 g0 t1 I9 u* dleap2:   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);
- E0 A/ A1 _5 o- m% g! |     }
" T3 n7 ~1 I3 S, C. J, R! F: \% [& X   c-=1e-8;
   if((c&lt;=0)&amp;&amp;(l1==1)&amp;&amp;((*a)[m1][mn1]-1e-8&gt;0))
- I" a' o2 v" I, d) D# @5 o      {
! I7 t  s( n2 i# k: c% W    printf("\n\t*********Not min&amp;&amp;No solution**********\n");
# q) F( i/ Y$ h9 ?    return(-1);
      }
9 t; o6 @1 g' @' S" @3 k7 g   if((c&lt;0)&amp;&amp;(l1==1))
       {
! l3 D# j6 H6 k6 c% K7 `     l1=50;
     goto leap2;
: @" F: {# t- i1 u       }
6 P' |: c6 N4 U    if(c&lt;=0)
       {
      for(i=1;i&lt;=m;i++)
         (*x)=(*a)[mn1];
         printf("\n\t********Optimal solution********\n");
      for(i=1;i&lt;=m;i++)
3 G6 r) X- r$ T$ j         printf("\ni=%d,x=%f\n",(*k),(*x));
6 ]) Q6 i. d5 U- {1 d# z, O      printf("\nF=%f\n",(*a)[m2][mn1]);
      return 0;
    }
     c=1e8;
* W( K7 R3 c! r) i# {     for(i=1;i&lt;=m;i++)
       {
     if((*a)[j0]&gt;1e-8)
        {
          g=(*a)[mn1]/(*a)[j0];
          if(g-c&lt;0);
# m+ K$ d8 l  Y6 Z' _         {
           c=g;
           i0=i;
5 `0 c# \$ N( z+ ~          }
: \3 p) {( H- s- J- k- R4 Q6 _6 w) p         }
    }
      if(c==1e+8)
    {
      printf("\n\t*******Lp no solution********\n");
      return -3;
( W, l8 f7 l& a8 {3 |     }
       (*k)[i0]=j0;
# q& }$ E) A" W9 J( f       for(j=1;j&lt;=mn1;j++)
8 \& n- M( I# I% v: C/ y      {
        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)))
1 x$ d1 D7 |% x- N6 U/ p6 s9 Z8 j) U            continue;
         (*a)[j]=(*a)[j]-(*a)[j0]*g;
/ L# A8 L3 a- l' S$ @  J6 B           }
       }
       for(i=1;i&lt;=m2;i++)
      (*a)[j0]=0;
       (*a)[i0][j0]=1;
       goto leap1;
}  
main()
{
: H/ P9 a0 _2 j. L  float a[5][7]={{0,0,0, 0,  0,0,0},
- Q5 C$ s2 b7 _  J                 {0,3,-4,3,  1,0,12},
" K' d! O) p4 Q( J& I1 S                 {0,3,0, 6,  0,1,12},
! t( |8 s, }  V8 A                 {0,0,0, 0,  0,0,0},
" a) {; Q0 h, p                 {0,69,0,144,0,0,300}};
* N" K! r% g3 m0 i* f* N+ n  int k[3];
  float x[3];
  clrscr();
  xi_sm(2,3,4,6,0,&amp;a,&amp;k,&amp;x);
! Q  f( ^3 N, C! E) Y6 v' {( ^# t}
$ n: n8 F# t7 B: \$ C0 ^* }* S</P>

plgatc

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

编程环境：VC++6.0     
! V% y5 c* i+ T. E 方程组输入说明：
变量非负，按提示输入相关参数。
/ M6 g5 G9 x, `( R5 W& v *************************************************************************/
#include &lt;stdio.h&gt;
$ j& @# X; q% F#include &lt;stdlib.h&gt;
#define MAX 100
0 \7 [2 i$ n, \& F#define STP 100</P>
<>int stop=1; //迭代记数变量
( D/ K/ x# P' c% l# |int 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; //变量个数
: \6 V. b4 d' d+ Z- ^7 r! Oint num_st; //约束方程数
int num_ar=0; //人工变量个数
int arti[MAX]; //人工变量下标
1 d3 V- a1 b0 I: w' @" o; Aint base[MAX]; //基变量下标
int ma_mi; //1为求最大值，2为求最小值</P>
, x3 d5 O3 G, R<>void create(); //建立方程组
void iterative(); //单纯型法迭代
- H' X; v% }/ k1 k2 c2 U2 bvoid output(); //输出结果
void banner(); //打印程序标题
void exchange(); //交换两阶段价值系数
8 Z$ |" Y: x5 {. y: N: z' ~void show(); //输出方程组</P>
<>void main() {
int i,j,k;
banner();
- H( U  x) d* u9 |$ s. Z- @ create();
//保存原价值系数，转换为第一阶段价值系数
for(i=1;i&lt;=num_x;i++) {
6 ^7 z* O  t: T0 T+ I! D# \  k=0;
& ]# G7 e  i$ m  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>
<> printf("\n\n第一阶段问题为：\n\n");
: W% w6 x9 _: z' k show();
step++;
printf("\n\n按回车开始第一阶段迭代");
' ^6 ~* k2 M  C8 }$ _$ x, f getchar();
3 O; |/ q9 ^2 z. g$ s3 j% y getchar();
iterative();
) i, }6 I# U+ h5 z% E  c if(status==-2) {
  puts("迭代超过限制次数强行终止！\n");
. g0 ?$ X5 Q3 u% ^7 [  puts("\n按回车结束");
' R: |+ ?3 @3 f# {+ f+ z, ?  `  getchar();
" M# W- z7 |, P$ J$ D+ s  exit(0);
/ r9 I+ q2 x- k9 m7 m9 |! G- N! D. V6 |% B }
1 |+ o- @$ H, m# \, } output();</P>
<> if(max!=0) {
' N; f$ X# o, v7 N% X# c  puts("\n\n原问题无可行解。\n");
$ t. D* R3 [- `) R  puts("\n按回车结束");
  getchar();
( k5 }: O; h  ]; |  exit(0);
# V3 r' s5 `0 X- `1 Y }</P>
% }. i$ {* V) l) y% p- {, |, R<> //转换为第二阶段价值系数
* t6 Q: O5 V+ D$ H% S exchange(c,temp_c);
8 z" O/ u4 |& }$ {& Q( C% `. j# d$ f //把人工变量列全设为0
for(i=1;i&lt;=num_ar;i++) {
: h* \( O2 S, X( y  c[arti]=0;
  for(j=1;j&lt;=num_st;j++) a[j][arti]=0;
' U7 ~9 M' l- e' \0 C! W# V# S }</P>
% w1 X) n9 \% ~" C1 z! h/ K<> puts("\n\n第二阶段问题为：\n\n");
show();
% F& a4 g' B% Y0 v2 R0 i puts("\n\n按回车开始第二阶段迭代");
9 b+ V! e4 b' }. ]) o1 J6 I: p getchar();
iterative();
switch(status) {
case 1:
: O3 y# f/ S9 E# t$ i/ H; F. c  output();
! O" r( F* d2 a2 E  puts("\n\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
" b+ q& ^& |: `! j" J8 N$ h& m1 w. L  getchar();
, v: B$ f& f) j  B  exit(0);
case 0:
  puts("\n\n原问题为无界解。\n");
  puts("\n按回车结束");
0 I" ?8 Q" f& \  A7 o8 z  getchar();
& d1 b# t- `: {  exit(0);
case -1:
  output();
  puts("\n\n原问题有无穷多最优解。\n");
, j% @/ d! V: t, a7 S8 b  puts("\n按回车结束");
  getchar();
' ?& D- @% m# j6 F- X  exit(0);
case -2:
  puts("迭代超过限制次数强行终止！\n");
& o; V( r6 D$ {, C/ o( h  puts("\n按回车结束");
( V# P0 S4 g- e4 A* H  getchar();
) i* P% T! b% {  exit(0);
}//switch
$ ?0 j: ?6 ^  x- _}</P>
<>void banner() {
1 Z; U6 \; o" _. m3 c printf("\t\t****************************************\n");
8 o' G# s, A5 a9 \* h printf("\t\t         单纯型法解线性规划问题\n");
printf("\t\t                         作者：Thunder\n");
, S- p# J( i2 }0 j" ?, ]& ^ printf("\t\t****************************************\n");
printf("\n");
2 f- F' O- k7 W0 F" L}</P>
$ H( M6 x" e$ t4 I4 n0 _<>void show() {
3 f' i* ?- U- B' W1 B8 F# w& b" }3 u//对方程组以自然的格式输出，系数为零的x不显示
& A  e8 y( S, |5 W//为1的不显示系数1，-1系数只显示负号
' |. n0 o. j4 U% K; B6 f int i,j,k;
switch(step) {
case 1:
  printf("min z= ");
  printf("x[%d]",arti[1]);
5 ?& H2 B& M1 [! O* M  for(i=2;i&lt;=num_ar;i++) printf(" + x[%d]",arti);
6 e2 C, `2 c; I0 j# K. |  break;
case 2:
printf("max z= ");
printf("%lg x[%d]",c[1],1);
0 P9 q# _, k+ p( R for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
  else if(c==-1) printf(" - x[%d]",i);
' e0 I4 F. k- V  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
; c0 u4 |  M9 A0 d  else printf(" %lg x[%d]",c,i);
' s+ F) `0 h! J! M1 ? }
break;
}</P>
<> printf("\nst:\n");
for(i=1;i&lt;=num_st;i++) {
" Z! V$ O* g; C  k=0;
  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
; n: I$ x$ i' s) w0 i! \) z# y2 S    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);
' i, ^8 Z  b( k    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);
. b) V; e3 h- X' r" \/ [    else printf(" %lg x[%d]",a[j],j);
  x7 Q. _$ Z; G4 @6 O. }6 k    k=1;
   }
! p; ~, H4 v/ L  }
, x' M* o# }$ C' E  g1 i  printf(" == %lg\n",b);
9 O% |* R$ q4 v" {. p }
printf(" x[1]~x[%d]&gt;=0",num_x);
' M8 Y4 _* ^. n0 \; s! {- E}</P>
# k& m8 Z# o) ~. P0 j<>void exchange() {
' J+ e! t9 {9 h) t: Z" S. D7 s8 m7 ~ int i;
  m; Y1 k: w# d0 L* X& x" M double temp[MAX];
. o! x0 `4 a! g. } for(i=1;i&lt;=num_x;i++) {
# w  J4 z9 W" s) a# D6 q  temp=temp_c;
& d; _& z# j1 M+ C4 h2 G  temp_c=c;
- O. f4 w+ t+ Q2 a2 _# r  c=temp;
- Y& b9 n4 z5 K4 K7 E* x }
# T' O. T2 E  d9 V0 O6 h& q}</P>
<>void create() {
  Z9 K5 {6 ?/ Z$ `* E//输入方程组系数，每个方程输完后回显确认
int i,j,k,re_st[MAX],tnum_x,num_addv=0,num_ba=0;
9 G+ F6 P% }2 T char confirm;

: @) k; R5 ^- m9 O- W% \- r9 h while(1) {
  printf("请选择：1、求最大值，2、求最小值：(1/2)");
  scanf("%d",&amp;ma_mi);
  if(ma_mi!=1&amp;&amp;ma_mi!=2) printf("输入错误，重新选择。");
- y5 X* C3 A# o# Y! `4 U6 Z, `  else break;
}
5 p0 w3 _  D* a8 A
while(1) {
: k' ?! C1 v- Z; b+ f# [  printf("指定变量个数：");
) B4 {$ A# R4 X1 U. m  scanf("%d",&amp;num_x);
; @* T+ g3 L& f$ |3 e5 K) X  printf("输入价值系数c1-c%d:\n",num_x);
  for(i=1;i&lt;=num_x;i++) {   
   printf("c%d=",i);
9 @5 J& b0 g$ S3 q' w   scanf("%lf",&amp;c);
  }
  if(ma_mi==1) printf("max z= ");
  else printf("min z= ");
; d- b( F$ C9 n, f  printf("%lg x[%d]",c[1],1);
1 o, i/ a1 y4 h/ G5 d: T/ s8 }  for(i=2;i&lt;=num_x;i++) {
, r9 J. M3 F  T. t% i" F6 R2 X   if(c&gt;=0) printf(" +%lg x[%d]",c,i);
4 f3 G, j- V" {: B- n. ~$ Z   else printf("  %lg x[%d]",c,i);
  }
& f" n6 }2 z" v( ~9 \  p' l printf("\n正确吗？：(y/n)");
getchar();
confirm=getchar();
* F$ d7 }# ~' a. [% K if (confirm=='y') break;
9 [$ L. |/ [. |& f3 I( F) d* r1 t else if(confirm=='n') continue;
. ~1 q2 ?6 {6 O* k& G }</P>
<> printf("输入约束方程组个数：");
; X4 }* ]: W/ S scanf("%d",&amp;num_st);
for(i=1;i&lt;=num_st;i++) {
$ Z8 ~' Z3 k  e9 c: U  printf("st.%d:\n",i);
# V$ \: F$ [& @" K- x: e  while(1) {
. K. c  o7 R6 P( @( V* [   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("输入错误，请重新选择。");
8 e/ A8 _4 F4 v2 k- @) B& }   else break;
  }
  printf("输入技术系数：\n");
  for(j=1;j&lt;=num_x;j++) {  
   printf("a%d=",j);
0 m# p* B! J0 w; [- U& O   scanf("%lf",&amp;a[j]);
  }
3 t! B+ T- K' S  H0 r: H! T  printf("输入资源拥有量：\nb%d=",i);
  scanf("%lf",&amp;b);
% _& Z( Y3 S! {& v' m. M5 T  
; i+ P5 f$ x1 W1 _  L  printf("st.%i:\n",i);
0 B) Z' ^: [0 v  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);
! p9 T6 q; \4 b: m' i6 [. s5 T   else printf(" %lg x[%d]",a[j],j);
- |6 O( i6 X5 p4 o& F5 z  }
( h" l! c( g, _) C  switch(re_st) {
4 t* A4 A5 ^) L( f/ S* P, l   case 1: printf(" == %lg",b); break;
# F( R. k/ M6 w. j8 Q- S   case 2: printf(" &gt;= %lg",b); break;
   case 3: printf(" &lt;= %lg",b); break;
  }</P>
<>  while(1) {
7 G9 ~. @- l3 \' c8 H3 m3 k   printf("\n正确吗？(y/n)");
   getchar();
   confirm=getchar();
   if (confirm=='y') break;
# P: Q5 G) s5 h; x$ V# f   else if(confirm=='n') {i-=1; break;}
: N. [: J4 z5 e. p' @  }
# ?, A- U  K; M& O/ n7 s }</P>
<>//显示输入的方程组
* K* U% s0 e4 g  z% ? printf("\n原问题为：\n\n");
if(ma_mi==1) printf("max z= ");
else printf("min z= ");
- i( {( B# x: { printf("%lg x[%d]",c[1],1);
for(i=2;i&lt;=num_x;i++) {
  if(c==1) printf(" + x[%d]",i);
2 @, N( w" r. z/ w# G  else if(c==-1) printf(" - x[%d]",i);
& t7 X, {( X: a! D0 j$ T8 c, }  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
; w8 Q9 x3 a8 {3 p! }) Y  else printf(" %lg x[%d]",c,i);
9 R" {  |. v7 W/ n" m1 f }</P>
' W$ Y+ O. M2 {<> printf("\nst:\n");
5 o: k" L, R- |- y2 X0 A0 @ for(i=1;i&lt;=num_st;i++) {
4 \5 F0 p, l4 ?4 _  k=0;
' @" V) ?- S( l: ]+ e7 W2 j* Q  for(j=1;j&lt;=num_x;j++) {
   if(a[j]!=0) {
. X; q- i$ Z: X  @  H0 j3 |    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);
& {( w" K% q  G% L( M5 l9 r    k=1;
   }
* X* V7 g5 \7 i  @+ `  }
1 Y6 M0 u& o5 F) d) ~  switch(re_st) {
" ^* l* e* c' V/ x0 o  I6 P& Y   case 1:
    printf(" == %lg\n",b);
8 h7 m% m+ U$ H6 d1 i    break;
   case 2:
    printf(" &gt;= %lg\n",b);
5 C3 S+ y5 E" _# P3 s. |    break;
   case 3:
, x1 C% Y& b# {9 J" O3 i& U    printf(" &lt;= %lg\n",b);
    break;
6 b4 t+ g+ y" |4 a  }
. M, S, O; O4 B; R }
, b6 [$ `% {8 } printf(" x[1]~x[%d]&gt;=0\n",num_x);</P>
. Q/ @- e9 K. G0 o2 @<> tnum_x=num_x;
. t. h0 j+ H! y# n" a for(i=1;i&lt;=num_st;i++) {
  switch(re_st) {
3 y' t( n) v# d! w: \4 q  case 1:
  case 3:
. M9 V$ C: @* W# s1 ^" z8 o   num_x+=1;
   break;
  case 2:
   num_x+=2;
' Z& _  I/ m! W   break;
  }
}</P>
; N5 V" |. e# h6 q/ r7 }- k<>//化为标准形式
! e0 R7 u) a4 H( n* b& h) Y9 t if(ma_mi==2) for(i=1;i&lt;=tnum_x;i++) c*=-1; //求最小值时，系数变相反数
" D5 E7 t( p  N: C" q for(i=1;i&lt;=num_st;i++) {
. p0 \' f& u' y; d! ]  switch(re_st) {
   case 1:
    num_addv++;
    num_ba++;
5 D: Q8 H8 B3 ~. A6 z    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++)
8 U, n0 s3 r) T* t     if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
     else a[j]=0;
    break;
   case 2:
    num_addv++;
    c[tnum_x+num_addv]=0;
    num_addv++;
    num_ba++;
    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;
     else if(j==tnum_x+num_addv) a[tnum_x+num_addv]=1;
2 }" g& H2 ?& F: q     else a[j]=0;
    break;
$ W1 d" z2 E' u# L: k4 E  L& C* ~   case 3:
9 u* U( s- s! @3 K5 y  F    num_addv++;
1 m6 ?( Q! D# p& {2 i. ~* B    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;
" S0 ]. M# ~6 g& S     else a[j]=0;
4 p/ Q4 i3 K" S' W) C& E+ c& J    break;
& J' f$ O  u& W. a% b9 D7 ]3 k* N% u4 t  }//switch
9 B& s% U, q; q. B5 D }//增加松弛变量、剩余变量、人工变量、确定基变量</P>
<>//显示标准化后的方程组
/ }7 Y9 ]4 Q, b; y- P+ d- \: z printf("\n化为标准形式后：\n\n");
if(ma_mi==1) printf("max z= ");
else printf("max z'= ");
printf("%lg x[%d]",c[1],1);
6 r& X, m+ l0 c* r, o  m, N7 g for(i=2;i&lt;=num_x;i++) {
0 l* w# y$ H- ^, h  k=0;
  for(j=1;j&lt;=num_ar;j++)
  if(i==arti[j]) k=1;
* A9 n! u# c  N9 i) N5 Z' P5 q& `  if(k==1) printf(" -M x[%d]",i);
( r) B8 x, M& T. U  else if(c==1) printf(" + x[%d]",i);
0 V8 `6 B3 y' ?  else if(c==-1) printf(" - x[%d]",i);
: ~3 @4 Z, b* Q1 ^. u  else if(c&gt;=0) printf(" +%lg x[%d]",c,i);
4 }: d+ N6 X! ~6 Y  else printf(" %lg x[%d]",c,i);
}</P>
<> printf("\nst:\n");
8 Z0 L: G, k6 o# } for(i=1;i&lt;=num_st;i++) {
5 U9 H0 y( H" ?4 u8 r' I* c6 t  k=0;
' r% x0 h' \$ Q: F8 `# B, l  for(j=1;j&lt;=num_x;j++) {
  ]1 j  A8 U  g   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 R, n; I! Q" q( l( x! o    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);
0 A6 U; l6 i/ g$ V  F+ g    else if(a[j]&gt;=0&amp;&amp;k==0) printf(" %lg x[%d]",a[j],j);
; t2 y* d- I4 R    else printf(" %lg x[%d]",a[j],j);
    k=1;
   }
  }
  printf(" == %lg\n",b);
) T% @9 }5 W  s, p# H }
printf(" x[1]~x[%d]&gt;=0",num_x);
, ]" C* K6 w' D}</P>
<>void iterative() {
int i,j,k,k_a,k_f,l; //k_a,k_f值为0或1，记录当前下标在arti[]或base[]里的搜索结果
* y8 {, m1 M' i, D7 p4 w+ a int base_elem;
int base_out,base_in;
double sigma[MAX],temp;
double value_be; //高斯消元里保存主元素值</P>
<> printf("\n\n第%d次迭代：\n\n",stop);
6 o2 d1 p- {+ O$ Z for(i=1;i&lt;=num_st;i++) {
" {0 `" E! y$ u, k3 }  printf("c%d=%lg\t",base,c[base]);
  printf("b%d=%lg\t",i,b);</P>
<>  switch(step) {
: ^% s% ~' p! |4 l. H% [   case 1:
    for(j=1;j&lt;=num_x;j++)
3 z/ h* S& q5 N0 Y/ t/ d; J$ d5 E    {
     printf("a[%d][%d]=%lg\t",i,j,a[j]);
    }
1 q/ J6 o4 y2 Z4 C6 \    printf("\n");
2 }5 h7 ?9 C  K0 h4 @7 n    break;
   case 2:
3 K% E7 l2 d% M5 b  j% }    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]);
    }
1 S6 g0 O4 e4 I4 Q& y1 b    printf("\n");
    break;
  }
}
//求检验数sigma
for(i=1;i&lt;=num_x;i++) {
  sigma=c;
0 g; F1 V  B: i) S5 X( U) |0 v  for(j=1;j&lt;=num_st;j++) sigma-=c[base[j]]*a[j];
* p, s+ o2 R8 n) e' o! r' E  for(j=1;j&lt;=num_st;j++) if(i==base[j]) sigma=0;
. B  m  x8 g1 f. E0 R  Q! y7 E- {  switch(step) {
7 H9 w1 T; ^( g1 t& f8 T* e( {9 f   case 1:
$ C6 P* W1 h" ^* {. `, w    printf("sigma[%d]=%lg\t",i,sigma);
" o4 U! Z" u+ r9 W    break;
& B4 G0 p6 n' R5 R: v   case 2:
# ^) s8 \4 P2 R' g' B4 Y. ^0 z    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 C7 n0 ~% V  c* `' |5 |3 T! I, _    break;
9 m) e. p6 Q: H$ B: k- ?  }
/ x! l3 s5 e) R  x8 k, S }
putchar('\n');
0 q) L# H( D% Q$ x- h//检验检验数sigma是否全小于等于0
k=0;
for(i=1;i&lt;=num_x;i++) {
3 z& I  G# |: Z0 c4 G  if(sigma&gt;0)
2 E2 U0 Q- R" o6 t  k=1;
( ]5 ~+ l2 `4 ]) l# _8 L  V! B( q }
if(k==0) {
  //sigma是全小于等于0时，检查是否为无穷多最优解
, ~- f" D( z) T! o  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 M/ b: G8 G" B  ~. f$ ]    for(j=1;j&lt;=num_st;j++) if(i==base[j]) k_f=1;
2 a) N4 F! J" d; X$ O   if(k_f==0) {status=-1; return;}
& s  Q, n3 x; S9 i) a% @   }
! K+ g$ C  @* U/ l  }
1 Q3 Y4 _, n# ~) v. X  status=1;
- k- I" X( m7 Y: o! ]; i+ U  return;
' G( U; I1 @" G  W5 T }
//检查是否为无界解
; i- I9 v7 g6 Q% O+ T( E% ?  P for(i=1;i&lt;=num_x;i++) {
  k_f=0;
  if(sigma&gt;0) {
   for(j=1;j&lt;=num_st;j++) if(a[j]&gt;0) k_f=1;
   if(k_f!=1) {status=0; return;}
* p  a1 |6 C. m  }
3 U& `9 ?! b8 W+ U- b4 v' p }</P>
<>//确定换入变量
# ^) `8 d( X/ z3 W; h  I for(i=1;i&lt;=num_x;i++) {
% M# @  l" p6 q6 c* q8 i  k=0;
9 C" C0 {0 V' P8 T7 t3 W, ^$ j  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;
, n! E) i! X& n- `6 u }//temp赋初值
5 ]3 [/ F" W  _: J0 O, e7 ^ for(i=1;i&lt;=num_x;i++) {
  k=0;
* p5 q9 n$ x6 _( ^. R) z6 q& |  for(j=1;j&lt;=num_st;j++) if(i==base[j]) k=1;
  if(k==0)
* n1 h; s! q" P0 L  M   if(sigma&gt;temp&amp;&amp;sigma&gt;0) {
' b9 U  f8 c, b/ ]7 x9 H7 H    base_in=i;
    temp=sigma;
   }
3 A/ p2 Q: j4 s2 d+ P. j9 ]) I }</P>
<>//确定换出变量
for(i=1;i&lt;=num_st;i++)
. x3 S' M4 b3 P5 U+ ]$ b  if(a[base_in]&gt;0) {
1 I' u0 O- Q5 \& k$ ^8 T   temp=b/a[base_in]+1;
   break;
$ |1 G  D7 W" u' P  }//temp赋初值
2 k3 b# K! \; ?* | 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++)
" a. W, N$ r; h    if(base==arti[j]) {
7 o: _  m9 M, l2 C2 k     base_out=base;
     base_elem=i;
1 ^% n6 j% K& Q5 G     temp=b/a[base_in];
     break;
    }
7 ?( t5 t6 J, T$ m) y8 M  h9 e( h6 |  }//人工变量优先换出
  if(b/a[base_in]&lt;temp&amp;&amp;a[base_in]&gt;0) {
   base_out=base;
# D, f7 v3 r) U' C" P   base_elem=i;
   temp=b/a[base_in];
, S- |+ _3 l5 ?  M' Y  }
}</P>
" M1 M! }4 D3 y0 Z4 w  `! U4 T<> printf(" 基变量：");
for(i=1;i&lt;=num_st;i++) printf("x[%d] ",base);
printf("换入变量：x[%d] 换出变量：x[%d]",base_in,base_out);
# N  h+ I$ Q5 n//基变量变换，进行新方程初始化后迭代
for(i=1;i&lt;=num_st;i++) {
: D5 a% C7 e% v, E! Q/ ?! t! D: n9 v: ?/ e  if(base==base_out) base=base_in;
}
6 O& `& j- c# @//初始化主元素行系数
value_be=a[base_elem][base_in];
b[base_elem]/=value_be;
for(i=1;i&lt;=num_x;i++) a[base_elem]/=value_be;</P>
7 ~0 x5 S( ?+ Y, L4 n<> for(i=1;i&lt;=num_st;i++) {
  if(i!=base_elem) {
, w+ o" }: s& m; t/ u   b-=b[base_elem]*a[base_in];
   value_be=a[base_in];
6 D3 A/ l" l1 o   for(j=1;j&lt;=num_x;j++) a[j]-=a[base_elem][j]*value_be;
' X+ X. D: |9 b' t. I1 o9 G) V  }
}
stop++;
& m" X; O( C/ l% w' h# s if(stop&gt;STP) {status=-2; return;}
iterative();
2 w( d$ X" ^8 ~}</P>
<>void output() {
% X4 i8 H/ w6 c$ I: j+ { int i,j;
  G3 A8 T7 q% x double X[MAX];
; [( ^0 {7 g4 E- O+ ^) e printf("\n结果如下：\n");
printf("\nX=(");
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;
( J2 ~" o7 [; P9 Q  printf("%lg ",X);
( h1 j0 j1 o: o0 g/ Y }
3 N& H( U. a7 j! [" S! X printf(")");
6 q  f0 C% [: E3 s% [1 O 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？帅呆了