个人写的运输问题算法，还望指教

plgatc

<>/*************************************************************************
8 t4 P! u  ?: n9 o8 V                         表上作业法解运输问题        

编程环境：VC++6.0     
程序说明：
" e$ K/ G4 n9 ^6 E% c Vogel法求初始解，位势法检验是否最优，递归搜索闭回路改进解。
  E& H3 v" y, v0 `9 D *************************************************************************/
, ^1 u+ X3 A& E( m#include &lt;stdio.h&gt;
#include &lt;stdlib.h&gt;</P>
<>#define MAX 100
#define M 999999</P>
) [/ ?: `3 M9 e! A) d8 }4 e! S9 K<>int status; //1唯一最优解，-1无穷多最优解</P>
<>int num_a,num_b;
" ]. a+ p6 [- m: {* Sdouble a[MAX],b[MAX],temp_a[MAX],temp_b[MAX];
double u[MAX],v[MAX],pu[MAX],pv[MAX];
double ab[MAX][MAX],temp_ab[MAX][MAX];
double base[MAX][MAX];
int pbase[MAX][MAX];
3 k0 X/ `3 j* x- z5 estruct element {
int r;
int c;
/ w$ a6 b# k% j& t# ~ double value;
+ F! S) d5 u# b  O, `};
: ~2 X  e5 j+ B- `: D  benum direction {mid,up,down,left,right};</P>
+ p* i. a3 S. `  I<>void create();
void banner();
double dabs(double d);
3 K/ I. o  n& w% Fvoid findinit();
void computeuv();
int check();
int  circle(enum direction dir,struct element cir[],int *t);
) H: w* A: o" _8 D9 e5 `% Wvoid improve();
9 B: [# t+ @$ I% xvoid output();</P>
<>void main() {
9 w! w$ y4 L5 n, N3 _7 T8 G6 T banner();
create();
8 j; b' p" G0 U7 ` printf("\n按回车确认后开始求解");
4 x( _: S. z: z) F/ y' q- k getchar();
getchar();
7 a: ]% |! W+ h' s  P findinit();
improve();
switch(status) {
case 1:
0 m& ^' j1 Q. @  x  output();
  puts("\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
! u2 K0 b' y* I- }' Z0 F) {  getchar();
  exit(0);
case -1:
  output();
1 [3 h- s% |( r  puts("\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
/ ~( [: M) U8 [. y. X  getchar();
  exit(0);
}
7 Y& @! k% t6 Q) L  C$ I+ G8 B! p}</P>
<>void banner() {
5 e6 a9 }7 d$ d. Y printf("\t\t****************************************\n");
" @# e$ S6 D/ e& d3 i printf("\t\t         表上作业法解运输问题\n");
: v) t. U+ |9 R) Z- U: M printf("\t\t                              Thunder\n");
printf("\t\t****************************************\n");
printf("\n");
}</P>
<>void create() {
int i,j;
3 W$ t7 D( n5 G/ z4 e double sum_a,sum_b;
$ z3 n% h4 p, n4 m6 S char confirm;
while(1) {
6 {; v- @! X9 i1 e" V% m; K, d  printf("请输入产地个数：");
" p: \1 N( p% ~  scanf("%d",&amp;num_a);
  for(i=1;i&lt;=num_a;i++) {
   printf("A%d产量：",i);
   scanf("%lf",&amp;a);
  }
  for(i=1;i&lt;=num_a;i++) printf("A%d产量：%lg\t",i,a);
  printf("\n正确吗？：(y/n)");
5 ?5 Z5 f; r4 p* }! d  getchar();
; A# n# L4 s# X7 ^& g  confirm=getchar();
  if (confirm=='y') break;
  else if(confirm=='n') continue;
}
while(1) {
  printf("\n请输入销地个数：");
  scanf("%d",&amp;num_b);
  for(i=1;i&lt;=num_b;i++) {
7 u! o# T6 Z8 q. E, F   printf("B%d销量：",i);
   scanf("%lf",&amp;b);
9 H2 i+ s& N+ p6 w6 P, j  }
  for(i=1;i&lt;=num_b;i++) printf("B%d销量：%lg\t",i,b);
" B, E6 F3 [+ T$ u5 O  printf("\n正确吗？：(y/n)");
  getchar();
  confirm=getchar();
. s% x) Y$ W! k( T+ s/ G  if (confirm=='y') break;
  else if(confirm=='n') continue;
6 R: ~; @! @% I  D' k% l1 X }
putchar('\n');
for(i=1;i&lt;=num_a;i++) {
  ~% r' Y/ L" w  for(j=1;j&lt;=num_b;j++) {
! t+ ~1 @$ l( G' o! N4 q- \1 Q   printf("A%d到B%d运价：",i,j);
   scanf("%lf",&amp;ab[j]);
  }
$ h4 ^& f0 }: W  for(j=1;j&lt;=num_b;j++) printf("A%d到B%d：%lg\t",i,j,ab[j]);
  printf("\n正确吗？：(y/n)");
  getchar();
  confirm=getchar();
  if(confirm=='n') i--;
. x- e2 k, W5 o3 x9 M  R4 \4 Z7 d  putchar('\n');
}
$ Z% f. u4 E8 O! J, C) F' P" l //处理产销不平衡的情况
5 s8 t6 b3 N& u, L  m+ D4 r( g sum_a=sum_b=0;
1 w" A9 j2 g1 A" @2 s7 R for(i=1;i&lt;=num_a;i++) sum_a+=a;
for(i=1;i&lt;=num_b;i++) sum_b+=b;</P>
) h. G5 C* w: E4 G5 t( J8 q! U6 p<> printf("总产量：%lg\t总销量：%lg",sum_a,sum_b);
" _# f  W+ R8 x; D' }; T if(sum_a==sum_b) printf("\t产销平衡。");
else if(sum_a&gt;sum_b) {
  printf("\t供大于求，增加假想销地B%d。\n",++num_b);
7 T; e# A6 R& H, U0 ?" i" p  b[num_b]=sum_a-sum_b;
% U; @" Q  X+ a2 [* n  for(i=1;i&lt;=num_a;i++) ab[num_b]=0;
: K' F% P8 B# R! t6 M8 w! x }
/ O5 u7 _: x0 o% I else if(sum_a&lt;sum_b){
  printf("\t供不应求，增加假想产地A%d。\n",++num_a);
6 d. k: a' ?/ ^  [+ H  a[num_a]=sum_b-sum_a;
  for(i=1;i&lt;=num_b;i++) ab[num_a]=0;
}
//求解前的准备
* J; W+ S# s1 N* A, h6 b& J- b% J for(i=1;i&lt;=num_a;i++)
  for(j=1;j&lt;=num_b;j++)
7 n( ~' v  Z# Y! F   base[j]=pbase[j]=0;
for(i=1;i&lt;=num_a;i++) temp_a=a;
5 c* i; D. d) [/ e9 D$ o0 q' J for(i=1;i&lt;=num_b;i++) temp_b=b;
  F2 X/ g* ]/ `6 j- L" x) P for(i=1;i&lt;=num_a;i++)
for(j=1;j&lt;=num_b;j++) temp_ab[j]=ab[j];
( ^5 C' }1 |* X //回显问题
printf("\n\n原问题为：\n\n");
; C9 ^9 n: M+ O; N2 X$ o   printf("产量：\n ");
for(i=1;i&lt;=num_a;i++) printf("A%d:%lg\t",i,a);
% P! r1 E1 i# j2 l4 r printf("\n运量：\n ");
" I6 O" A9 `, N4 Z' H+ L for(i=1;i&lt;=num_b;i++) printf("B%d:%lg\t",i,b);
printf("\n运价为：\n");
- J- a9 _) |  Z6 j% M7 t7 j6 U for(i=1;i&lt;=num_a;i++) {
1 k6 Y0 F' m) [& O# S( o  putchar(' ');
  for(j=1;j&lt;=num_b;j++)
   printf("A%d-&gt;B%d:%lg\t",i,j,ab[j]);
  putchar('\n');
}
}</P>
# d8 W9 o  {" Z, Y' P  k% |7 Z9 D<>double dabs(double d) {
if(d&gt;=0) return d;
else return -d;
}</P>
6 q  }9 W; I$ a0 A# J4 F# E/ n<>void findinit() {
. F$ z5 |+ X$ f int i,j,k;
) w" _& k2 U! @. f* l9 l% Z double r_penum[MAX],c_penum[MAX];
# X7 Z2 V  L3 p1 T- {7 y3 V. [8 M2 E struct largest {
; \& {# c1 P! f8 o' C/ r  char rc;
9 t* K1 D: S+ q, V8 q7 W  int num;
  double value;
' }; s7 B+ X* o' [4 d; r };
% Z3 ^' I3 d* {  O/ B% h struct largest lar1,lar2;
5 B' n0 `' p+ e  g+ w; z int r=0,c=0;
0 g2 Z) i" g' B& o double a1,a2,b1,b2,temp;</P>
8 O( u1 s; b: Z9 V# r<> for(i=1;i&lt;=num_a;i++) {
$ x, p! |6 T( F. a% t  temp=temp_ab[1];
  for(j=1;j&lt;=num_b;j++)
   if(temp_ab[j]&gt;temp) temp=temp_ab[j];
  a1=a2=temp;
, ]' Q' ~6 N, H# K& A  for(j=1;j&lt;=num_b;j++)
6 ]5 w7 X  _& q& {8 J5 A# e   if(temp_ab[j]&lt;=a1) {
    a2=a1;
2 W4 S: H8 M7 K! J( T5 [, g    a1=temp_ab[j];
5 T8 e5 h6 I7 s7 R3 r4 `4 c# H/ C   }
' D: {  ]; W- v   else if(temp_ab[j]&lt;=a2) a2=temp_ab[j];
  r_penum=dabs(a1-a2);
8 z  F& c- q  z' h2 y3 z) ] }</P>
8 T4 v2 D& Y9 S; a<> for(i=1;i&lt;=num_b;i++) {
  temp=temp_ab[1];
- N  X0 D; e" S  for(j=1;j&lt;=num_a;j++)
   if(temp_ab[j]&gt;temp) temp=temp_ab[j];
  b1=b2=temp;
  for(j=1;j&lt;=num_a;j++)  
   if(temp_ab[j]&lt;=b1) {
    b2=b1;
2 ]3 c. L4 o/ }6 M  |! E: i7 d( d2 V    b1=temp_ab[j];
. ~7 C+ U+ B9 ^- w. p& q   }
   else if(temp_ab[j]&lt;=b2) b2=temp_ab[j];
! K" T7 [+ F$ j& j+ J  c_penum=dabs(b1-b2);
}
; x% S& d/ N7 g0 Y5 ~- r) v" S /*
; p8 k6 _7 M$ w+ { for(i=1;i&lt;=num_a;i++) printf("pa%d=%lg ",i,r_penum);
5 b8 f: V+ g& `& s putchar('\n');
! T# s6 o  E  C" Z/ K3 q0 c; m for(i=1;i&lt;=num_b;i++) printf("pb%d=%lg ",i,c_penum);
*/
temp=r_penum[1];
for(i=1;i&lt;=num_a;i++) if(r_penum&lt;temp) temp=r_penum;
for(i=1;i&lt;=num_b;i++) if(c_penum&lt;temp) temp=c_penum;</P>
9 Q. e$ w% c" h<> //考虑了有两个罚数相等的情况，大于两个相等只取其中两个进行比较
lar1.value=lar2.value=temp;
  X# @4 @; Q0 ]" n& k- h- ~ lar1.num=lar2.num=1;
lar1.rc=lar2.rc='r';
for(i=1;i&lt;=num_a;i++)
' G$ [9 n0 Z' c  V2 M! r  if(r_penum&gt;=lar1.value) {
   lar2=lar1;
   lar1.rc='r';
8 x3 A& Y) G# c8 E   lar1.num=i;
   lar1.value=r_penum;
  }
! d' g/ T% S3 {8 W2 v  else if(r_penum&gt;=lar2.value) {
   lar2.rc='r';
% E7 ], E# }+ L" a8 `8 h( q: O( S   lar2.num=i;
1 b% A0 _" u3 K# f  W$ ^  L- t   lar2.value=r_penum;
8 R2 k+ _' f+ f3 v$ w5 D! G' U  }</P>
& E; }3 l" i$ m8 b/ E2 V9 b8 k<> for(i=1;i&lt;=num_b;i++)
: j/ A8 E8 f  y2 z" U2 y  if(c_penum&gt;=lar1.value) {
* ]. p! D1 H, w  U   lar2=lar1;
   lar1.rc='c';
7 |: T& Z% g: W( x6 p   lar1.num=i;
& u4 M0 m% F, U  j   lar1.value=c_penum;
  }
! P2 |1 V( L8 I  else if(c_penum&gt;=lar2.value) {
   lar2.rc='c';
   lar2.num=i;
$ @  ~6 G2 U2 W  E   lar2.value=c_penum;
  }</P>
: d' r- W" D0 l4 L$ Z3 I6 x) e3 o<> if(lar1.value==lar2.value&amp;&amp;(lar1.num!=lar2.num||lar1.rc!=lar2.rc)) {
  if(lar1.rc=='r'&amp;&amp;lar2.rc=='r') {
   temp=temp_ab[lar1.num][1];
: O. p  p, m9 q1 x5 ]* v   r=lar1.num;
      c=1;
, h+ E& e0 y: \  h7 n  d   for(i=1;i&lt;=num_b;i++) {
    if(temp_ab[lar1.num]&lt;temp) {
     temp=temp_ab[lar1.num];
7 ~8 ]& H# p. [6 H( E& k( d     r=lar1.num;
+ N7 a2 z% t( e3 H% A- W/ j  Z+ j" g     c=i;
. ?! e& z$ ~8 J5 g) f8 s) ^    }
    if(temp_ab[lar2.num]&lt;temp) {
     temp=temp_ab[lar2.num];
     r=lar2.num;
# E& z+ _, p0 m$ {: c& w) R6 [     c=i;
    }
   }
1 f9 q+ q$ E. d" k+ u3 W+ {  }
  if(lar1.rc=='c'&amp;&amp;lar2.rc=='c') {
   temp=temp_ab[1][lar1.num];
   r=1;
8 C- I4 R5 u3 _3 e2 v! z0 _      c=lar1.num;
   for(i=1;i&lt;=num_a-1;i++) {
    if(temp_ab[lar1.num]&lt;temp) {
& T. L/ x/ M1 @- i     temp=temp_ab[lar1.num];
     c=lar1.num;
     r=i;
    }
    if(temp_ab[lar2.num]&lt;temp) {
     temp=temp_ab[lar2.num];
9 k, {6 x. S& k1 @/ m) `1 ^7 U, c     c=lar2.num;
     r=i;
" {6 K- K% m' z- R# b    }
   }
4 a$ U. G3 b) N  }
  if(lar1.rc=='r'&amp;&amp;lar2.rc=='c') {
2 l4 d* y1 F8 V7 V   temp=temp_ab[lar1.num][1];
0 B1 @4 [( K: E2 y' P$ S   r=lar1.num;
4 Z$ T+ u$ R, S: Q      c=1;
   for(i=1;i&lt;=num_b;i++)
8 t0 \4 l" c0 A& ], [4 I0 g  B4 g    if(temp_ab[lar1.num]&lt;temp) {
     temp=temp_ab[lar1.num];
     r=lar1.num;
% g, {: G, D5 i     c=i;
! I1 Q2 u( w# d5 E! z9 e+ [    }
/ r% m2 F; H) e: `  O8 L* d   for(i=1;i&lt;=num_a;i++)
    if(temp_ab[lar2.num]&lt;temp) {
, r. s/ ?$ Q# h  m6 F( b. V+ ^+ n     temp=temp_ab[lar2.num];
' d$ V) g' `# e4 s5 [7 \     c=lar2.num;
" q9 d( p/ Q9 R* e     r=i;
9 U& P' o: T. W0 q9 T# ]5 n    }
2 d. m# a9 |; }/ @5 i  }
  if(lar1.rc=='c'&amp;&amp;lar2.rc=='r') {
   temp=temp_ab[1][lar1.num];
+ G% h! m2 e# D( @1 R/ {, D   r=1;
4 I: ]& [& Q5 B4 j      c=lar1.num;
9 \: K( K$ K0 _   for(i=1;i&lt;=num_b;i++)
    if(temp_ab[lar2.num]&lt;temp) {
3 J" @! S; q4 J0 B. H& ]. i% ]     temp=temp_ab[lar2.num];
. W9 v. X6 `) c4 k; \- ^6 O; S     r=lar2.num;
: k4 n. ?) t. \5 F- {     c=i;
    }
   for(i=1;i&lt;=num_a;i++)
    if(temp_ab[lar1.num]&lt;temp) {
# U- {9 D% H" e" G4 C     temp=temp_ab[lar1.num];
     c=lar1.num;
     r=i;
$ A+ z1 @# F8 O) l* b8 S+ O    }
; O: f, B' \( y% I# r  }
}
else {
  if(lar1.rc=='r') {
   r=lar1.num;
   c=1;
   temp=temp_ab[lar1.num][1];
   for(i=1;i&lt;=num_b;i++)
* _- d0 M: m/ c) I    if(temp_ab[lar1.num]&lt;temp) {
% `3 Q" c% }% `. w3 d  V4 X# H     c=i;
     temp=temp_ab[lar1.num];
    }
  Z! e% `1 E' U1 O- h  }
  else {
   r=1;
) ~1 V' x7 h" ?   c=lar1.num;
$ X4 V- G4 E3 X; B* f- l* j5 X   temp=temp_ab[1][lar1.num];
4 {' P% {: m) i: T9 ]   for(i=1;i&lt;=num_a;i++)
    if(temp_ab[lar1.num]&lt;temp) {
     r=i;
6 X2 Q# @4 O) b! J7 n7 U     temp=temp_ab[lar1.num];
    }
9 c: J) b6 c, g; G9 ~4 X! X/ w  }
- y8 V% y/ a8 R/ \# b! K; U }
pbase[r][c]=1;
if(temp_a[r]&gt;temp_b[c]) {
: X( x" Q- R/ N( F, X8 s  base[r][c]=temp_b[c];
  temp_a[r]=temp_a[r]-temp_b[c];
  temp_b[c]=0;
% X) w1 J/ H( t  D  for(i=1;i&lt;=num_a;i++) temp_ab[c]=M;
6 [1 k- \  [  W: _2 l }
+ |9 e- d) K6 G else if(temp_a[r]&lt;temp_b[c]) {
  base[r][c]=temp_a[r];
. M& u9 y* I4 w+ _  temp_b[c]=temp_b[c]-temp_a[r];
0 n' x( I$ f  ^% R; Y4 H" T  temp_a[r]=0;
  for(i=1;i&lt;=num_b;i++) temp_ab[r]=M;
}
else if(temp_a[r]==temp_b[c]) {
  base[r][c]=temp_a[r];
9 b, x, j1 k& |# i0 a7 O  temp_a[r]=temp_b[c]=0;
3 p. v9 ~! s* f  h1 \  for(i=1;i&lt;=num_a;i++) temp_ab[c]=M;
3 o  C) \/ k0 ^) D/ Y) B0 f  for(i=1;i&lt;=num_b;i++) temp_ab[r]=M;
  J! H2 T( O" x  ^$ d0 v' c2 `2 u  k=0;
: L9 `0 s* ~: Z  for(i=1;i&lt;=num_a;i++) if(temp_a!=0) k=1;
  for(i=1;i&lt;=num_b;i++) if(temp_b!=0) k=1;
) X9 X1 g# |( _* j" @  if(k!=0) {
   k=0;
8 ^: L/ W1 @% Y' K! ^& a& `8 a   for(i=1;i&lt;=num_a;i++) if(pbase[c]!=1) {k=1;break;}
& c7 z2 Y+ f( n! v; v( H   for(j=1;j&lt;=num_b;j++) if(pbase[r][j]!=1) {k=2;break;}
   if(k==1) pbase[c]=1;
   else if(k==2) pbase[r][j]=1;
) w4 m5 q( B1 g: x! F  }
}
k=0;
for(i=1;i&lt;=num_a;i++) if(temp_a!=0) k=1;
0 _: e7 x- Y  ~1 M for(i=1;i&lt;=num_b;i++) if(temp_b!=0) k=1;
) K/ J" Q# \" ]; r: B if(k==0) return;
! Z! n8 `/ m& X1 i$ r findinit();
}</P>
9 u; @5 G3 }) D% D. h<>void computeuv() {
* d6 O( R* H! X, E* `5 B int i,j,k;
% k7 n$ V. {% y& n  k; V0 l3 X) P for(i=1;i&lt;=num_a;i++)
6 h) F; B  ^3 M* _/ t, B  for(j=1;j&lt;=num_b;j++)
   if(pbase[j]==1) {
, n' ]# j7 w" ]& s2 c. a  P% F; {1 m    if(pu==1) {v[j]=ab[j]-u; pv[j]=1;}
) {4 M: l  O; v- H! L; O  n    if(pv[j]==1) {u=ab[j]-v[j]; pu=1;}
   }
k=0;
0 z- a* M: L7 b- L! y for(i=1;i&lt;=num_a;i++) if(pu==0) k=1;
* W( m6 ~0 p0 J$ X& o5 f- m% Y for(i=1;i&lt;=num_b;i++) if(pv==0) k=1;
) J5 \( S8 c! y9 x" ]+ j if(k==1) computeuv();
}</P>
) o' @- g: F1 o<>int check() {
int i,j,k,l;
2 m- o1 N% Q/ H( Z k=l=0;
for(i=1;i&lt;=num_a;i++)
  for(j=1;j&lt;=num_b;j++)
   if(pbase[j]==0) {
    base[j]=ab[j]-u-v[j];
    if(base[j]&lt;0) k=1;
- @, d) J' T8 X0 q6 ]    else if(base[j]==0) l=1;
   }
/* for(i=1;i&lt;=num_a;i++){
  for(j=1;j&lt;=num_b;j++) printf("base %lg\t",base[j]);
  putchar('\n');}
, }4 a5 l+ E2 A, N; ], x for(i=1;i&lt;=num_a;i++){
9 u2 A% g* e* |  n  for(j=1;j&lt;=num_b;j++) printf("pbase %d\t",pbase[j]);
6 h8 Q% Y. N1 F$ q  putchar('\n');}
1 X( ^6 F1 W$ ?- M*/
if(k==0&amp;&amp;l==0) {status=1; return 1;}
else if(k==0&amp;&amp;l!=0) {status=-1; return 1;}
else return 0;
* O; e0 O0 b3 U. z}</P>
<>int circle(enum direction dir,struct element cir[],int *t) {
- Q# V6 h9 p- [ int i,j,k;
/*
putchar('\n');
& Y9 }" c& \: w& s* U0 R5 R; X4 J for(i=1;i&lt;=*t;i++) printf("%d: r%d c%d\t",i, cir.r,cir.c);
% u; V9 z5 I: D- Q putchar('\n');
, t) s2 N7 N" h& I0 ] */
( |, A, G$ K* Y3 h* V- X8 ]+ O if((*t)!=1&amp;&amp;cir[(*t)].r==cir[1].r&amp;&amp;cir[*t].c==cir[1].c) {
: w6 ?4 }! c" r4 ]; ^  t--;
  return 1;
! n! P) L5 C% H9 x: I1 N }
if(dir!=down) {
! s' [9 E  A' G  k=0;
  for(i=cir[*t].r-1;i&gt;=1;i--)
   if(pbase[cir[*t].c]==1||(i==cir[1].r&amp;&amp;cir[*t].c==cir[1].c)) {k=1; break;}
, _$ z3 D- m" l& n' {  for(j=2;j&lt;=*t;j++) if(i==cir[j].r&amp;&amp;cir[*t].c==cir[j].c) k=0;
  if(k==1) {
, f/ O: j% j5 Y& _9 ^+ R   if(dir==up) cir[*t].value=0;
   else cir[*t].value=1;
- k9 w8 m* e4 G! I& f$ [   (*t)++;
   cir[*t].r=i;
   cir[*t].c=cir[(*t)-1].c;
   if(circle(up,cir,t)) return 1;
/ T1 w( ^7 O% U1 O" n- f  }
}
- t1 J9 f: F% w3 `' S if(dir!=up) {
  k=0;
  for(i=cir[*t].r+1;i&lt;=num_a;i++)
  {/ r. k! ]1 |  c   if(pbase[cir[*t].c]==1||(i==cir[1].r&amp;&amp;cir[*t].c==cir[1].c)) {k=1; break;}
  for(j=2;j&lt;=*t;j++) if(i==cir[j].r&amp;&amp;cir[*t].c==cir[j].c) k=0;
  if(k==1) {
* h+ g9 I. ]8 w0 n7 Z   if(dir==down) cir[*t].value=0;
1 ^! {! Y# v% [8 \* Z0 w$ v   else cir[*t].value=1;
   (*t)++;
   cir[*t].r=i;
/ {" L. p8 f* Y9 ~. S   cir[*t].c=cir[(*t)-1].c;
   if(circle(down,cir,t)) return 1;
3 f* D) V% I" C% i5 H/ K. f  }
}
; l0 A( M! F! C5 O( K if(dir!=right) {
  k=0;
* n9 F" w0 L7 c6 x/ G! h  for(i=cir[*t].c-1;i&gt;=1;i--)
2 L) `) E/ l0 _- [: o* V2 f   if(pbase[cir[*t].r]==1||(cir[*t].r==cir[1].r&amp;&amp;i==cir[1].c)) {k=1; break;}
  for(j=2;j&lt;=*t;j++) if(cir[*t].r==cir[j].r&amp;&amp;i==cir[j].c) k=0;
9 ~3 h, a7 A9 ?0 `6 O8 E  if(k==1) {
   if(dir==left) cir[*t].value=0;
   else cir[*t].value=1;
4 N2 n$ t* X$ C( i$ q: N5 g   (*t)++;
) {& C/ B( z# h, ?   cir[*t].r=cir[(*t)-1].r;
   cir[*t].c=i;
# W1 \: G' p: @/ U   if(circle(left,cir,t)) return 1;
  }
}
; \& ]1 i# ^6 G8 ` if(dir!=left) {
  k=0;
3 g0 X( [7 F. c8 K  for(i=cir[*t].c+1;i&lt;=num_b;i++)
   if(pbase[cir[*t].r]==1||(cir[*t].r==cir[1].r&amp;&amp;i==cir[1].c)) {k=1; break;}
* I7 L% R; E& N% [  for(j=2;j&lt;=*t;j++) if(cir[*t].r==cir[j].r&amp;&amp;i==cir[j].c) k=0;
0 m% T9 m$ l: c: s# B- T% `  if(k==1) {
% Z" N8 G+ b: P   if(dir==right) cir[*t].value=0;
# H. B7 S, O3 N  C4 [* [   else cir[*t].value=1;
0 i6 \" \* @* i, ]* j+ F: j   (*t)++;
   cir[*t].r=cir[(*t)-1].r;
   cir[*t].c=i;
   if(circle(right,cir,t)) return 1;
  }
% x0 k& t  J; u7 ?/ R9 ~ }
. q# m" o7 S. D (*t)--;
return 0;
. @& ]1 d2 F- m+ W/ G4 c" }& p# `}</P>
<>void improve() {
int i,j,k,t;
# P; }. e+ U' n8 u/ h3 p6 c struct element base_in,base_out,cir[2*MAX+1];</P>
<> for(i=1;i&lt;=num_a;i++) pu=0;
4 J& u0 _: I7 U8 M for(i=1;i&lt;=num_b;i++) pv=0;
u[1]=0; pu[1]=1;
computeuv();
0 y+ u" @9 r4 @! x# I& ^& k9 h if(check()) return;
8 U9 E" b* h/ ?* s' K$ O- {( E! K; c for(i=1;i&lt;=num_a;i++)
  for(j=1;j&lt;=num_b;j++)
   if(pbase[j]==0) break;
base_in.r=i;
2 D, d: H/ B0 ] base_in.c=j;
base_in.value=base[j];</P>
/ f+ z+ N; Z; |7 C& c2 `3 ?% t6 i<> for(i=1;i&lt;=num_a;i++)
  for(j=1;j&lt;=num_b;j++)
" U% e+ k5 H- h4 L& d9 ?7 K   if(pbase[j]==0&amp;&amp;base[j]&lt;base_in.value) {
4 ~( ~. @  n$ f: _# J' p* q    base_in.r=i;
3 u1 K: L; h+ ^/ P. s2 e) c" I- L    base_in.c=j;
    base_in.value=base[j];
   }</P>
$ N$ s- S& [  a9 `) |2 s) \$ @- S<> t=1;
cir[t].r=base_in.r;
cir[t].c=base_in.c;
cir[t].value=1;
if(!circle(mid,cir,&amp;t)) {
  printf("程序出错！按回车结束");
  getchar();
  exit(0);
2 u, O5 i& q2 ^$ u/ o }
t--;
// for(i=1;i&lt;=t;i++) printf("%d:r%d c%d v%lg\t",i,cir.r,cir.c,cir.value);
// putchar('\n');</P>
6 M# ^. x8 r! N" b. F$ X9 J<> for(i=2;i&lt;=t;i++) if(cir.value==1) break;
base_out=cir;
base_out.value=base[cir.r][cir.c];
k=1;
! g4 H$ I9 x$ B6 |* ~$ T for(i=1;i&lt;=t;i++)
  if(k%2==0&amp;&amp;cir.value==1&amp;&amp;base[cir.r][cir.c]&lt;base_out.value) {
   base_out=cir;
6 q5 P8 F4 T9 X4 ^$ R   base_out.value=base[cir.r][cir.c];
   k++;
! W: k6 H& u( l2 X4 V1 C/ p  }
$ Z  f+ U- V& q9 \; ~  else if(cir.value==1) k++;
& a! f& d3 i- N  I base[cir[1].r][cir[1].c]=0;
k=1;
" A/ d, a2 l, _3 V+ B& N% ? for(i=1;i&lt;=t;i++) {
% {7 }2 d5 f& c7 w  E  if(k%2==1&amp;&amp;cir.value==1) {
9 C) `  }7 n) r" H, m   base[cir.r][cir.c]+=base_out.value;
   k++;
( g9 `: W7 {' c8 Q2 x  }
  else if(k%2==0&amp;&amp;cir.value==1) {
   base[cir.r][cir.c]-=base_out.value;
   k++;
) n7 {) }' T  p. v. v6 X7 X. D5 R; Q  }
}
& M& a( w3 W$ t8 ~9 x pbase[base_out.r][base_out.c]=0;
0 [3 P6 f0 u" j- C/ b pbase[base_in.r][base_in.c]=1;</P>
<> improve();
}</P>
# ]; o7 _, I6 o" D# O5 _<>void output() {
int i,j;
double sum=0;
* w. ~( b( e. P# `* r printf("\n运量为：\n");
for(i=1;i&lt;=num_a;i++) {
  putchar(' ');
  for(j=1;j&lt;=num_b;j++)
- T9 m* K0 r$ ]: @. L& c   if(pbase[j]==1) {
    printf("A%d-&gt;B%d:%lg\t",i,j,base[j]);
    sum+=base[j]*ab[j];
   }
% \, ~/ G* k$ ~6 @   else printf("A%d-&gt;B%d:0\t",i,j);
  putchar('\n');
}
& E  h5 r9 O  E. N printf("\n最低总运费为：%lg\n",sum);
  q% t+ m* x' I5 x1 ]+ f}
</P>

liunengwu

好！顶顶！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！

hzlhm

乱了让人不想看了。。。。。。。。

enrique08

太乱了吧！！！！！！！！！！

zhongguohelin

good!   就是希望下次把格式稍稍注意下，这样可读性太差了呀~