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

plgatc

<>/*************************************************************************
                         表上作业法解运输问题        

; ?. m+ v2 j' P" {& v) ~( W6 H9 f 编程环境：VC++6.0     
5 K. W7 k' I! x5 t5 O/ f; b: A7 |* t: P 程序说明：
Vogel法求初始解，位势法检验是否最优，递归搜索闭回路改进解。
*************************************************************************/
# c4 k) ~, M2 z9 |# r#include &lt;stdio.h&gt;
#include &lt;stdlib.h&gt;</P>
<>#define MAX 100
  r3 r9 p* B5 r+ N7 t8 t#define M 999999</P>
- @0 E( K% j3 y: X. J# D<>int status; //1唯一最优解，-1无穷多最优解</P>
$ i  U/ [( B* g3 A6 ^* l<>int num_a,num_b;
7 T& f8 y9 o8 @# Y4 J$ x6 \# Edouble a[MAX],b[MAX],temp_a[MAX],temp_b[MAX];
5 J) g; Q; c& e0 F, A* Ndouble u[MAX],v[MAX],pu[MAX],pv[MAX];
, }" U- O# |, F3 f3 r# }: Kdouble ab[MAX][MAX],temp_ab[MAX][MAX];
double base[MAX][MAX];
: [/ g4 r% T  Eint pbase[MAX][MAX];
  f# l  J  N3 K9 A1 M0 c& Bstruct element {
  r& o0 |' i4 e3 d( x. a2 v' F int r;
int c;
0 V: d% s* B2 I5 l/ r double value;
' z! D* Z1 ^8 \% F5 I. ^};
enum direction {mid,up,down,left,right};</P>
<>void create();
void banner();
double dabs(double d);
+ G' j3 f" V% X" h# m$ o% t, ~void findinit();
8 `9 M) {" w) V4 f; R3 m# G8 p3 M7 l# ovoid computeuv();
* y* U, ?; k. [  V/ mint check();
' [: o1 i) R5 D/ d3 U- g6 d8 Hint  circle(enum direction dir,struct element cir[],int *t);
void improve();
% [. s$ L7 Z# r$ ^7 w3 e2 @" b) Yvoid output();</P>
<>void main() {
banner();
6 `: g) ]5 ?! d9 z' C+ D" B3 R create();
printf("\n按回车确认后开始求解");
getchar();
getchar();
findinit();
% m# C! U: l/ \  w: t' v+ g+ X improve();
- j1 |2 T8 C, B3 V! H, E- ]+ U switch(status) {
case 1:
! o: ^) k% U' t  p, A  output();
  puts("\n原问题有唯一最优解。\n");
  ~3 g2 U9 N+ M/ `  puts("\n按回车结束");
  getchar();
  exit(0);
! }) N) x4 b, W" M* V% x case -1:
  output();
  puts("\n原问题有无穷多最优解。\n");
  puts("\n按回车结束");
$ o) A5 Y$ d! ~7 {0 d- u  getchar();
4 |; l; k" |% n) i: B- ?  exit(0);
: N- R5 x: v1 ` }
}</P>
<>void banner() {
printf("\t\t****************************************\n");
2 j% z0 O3 q/ B7 f  i' {0 ?8 | printf("\t\t         表上作业法解运输问题\n");
5 w# }1 M9 y6 U6 b printf("\t\t                              Thunder\n");
printf("\t\t****************************************\n");
& P$ w. ^3 T+ H4 |& v* K% E4 ?* i7 B printf("\n");
}</P>
<>void create() {
% W' d  x5 ?6 f/ r5 n1 v, n& L9 D int i,j;
3 A2 ]6 N$ J: M4 t double sum_a,sum_b;
char confirm;
while(1) {
: L$ I) b# R# s& X; [# g( N  printf("请输入产地个数：");
5 O  ?/ b' ]: x' x/ W5 K  scanf("%d",&amp;num_a);
2 P, {# S2 o2 C3 e* @  for(i=1;i&lt;=num_a;i++) {
% H( n8 N9 g& N" \8 J   printf("A%d产量：",i);
# f  X+ F" c# v$ [2 V- g$ p   scanf("%lf",&amp;a);
$ k* z: [/ ^+ ^0 O7 p  }
( c. f0 A, L: f5 N  for(i=1;i&lt;=num_a;i++) printf("A%d产量：%lg\t",i,a);
" @+ P2 x5 i- J% c! i3 t/ v  y* W! M  printf("\n正确吗？：(y/n)");
  getchar();
  confirm=getchar();
3 e' B" r; c, c- t, N) U# N  if (confirm=='y') break;
  else if(confirm=='n') continue;
* q) y5 D& o1 q7 z8 T& f( s }
( T2 C) I, L0 R, C$ N  X) t8 v while(1) {
  printf("\n请输入销地个数：");
. P" g. b6 H, c* a8 _. e, I+ Q: U  scanf("%d",&amp;num_b);
4 h1 A2 q+ ^4 {) n5 U7 u- P  for(i=1;i&lt;=num_b;i++) {
   printf("B%d销量：",i);
3 i  E" \7 _" G   scanf("%lf",&amp;b);
1 `8 {9 L" G/ i  m6 V  }
( H, s! o" j. [5 u: w4 {0 r& V4 G/ B  for(i=1;i&lt;=num_b;i++) printf("B%d销量：%lg\t",i,b);
  printf("\n正确吗？：(y/n)");
9 M/ b4 t6 e) A1 t& F  getchar();
  confirm=getchar();
  if (confirm=='y') break;
  else if(confirm=='n') continue;
3 T2 `' o9 L5 @- [* N, L: Y/ x }
  W* F9 X$ u2 v# d* z: v. | putchar('\n');
for(i=1;i&lt;=num_a;i++) {
  for(j=1;j&lt;=num_b;j++) {
- Q7 m# p4 N; \1 j5 v$ Q1 n   printf("A%d到B%d运价：",i,j);
   scanf("%lf",&amp;ab[j]);
  }
  for(j=1;j&lt;=num_b;j++) printf("A%d到B%d：%lg\t",i,j,ab[j]);
  printf("\n正确吗？：(y/n)");
  getchar();
$ t4 w1 S, r8 T3 o& g/ x  confirm=getchar();
; ^: B0 W# Z. f2 L7 ]/ b  if(confirm=='n') i--;
& r( ~) W- R' U. m2 A; k- _  putchar('\n');
}
5 ~7 B, I  u% v1 s( Z; x- I //处理产销不平衡的情况
sum_a=sum_b=0;
for(i=1;i&lt;=num_a;i++) sum_a+=a;
for(i=1;i&lt;=num_b;i++) sum_b+=b;</P>
2 [# `. C8 r7 y- |) h9 O$ d7 P0 h<> printf("总产量：%lg\t总销量：%lg",sum_a,sum_b);
if(sum_a==sum_b) printf("\t产销平衡。");
3 s  g9 r4 I# m; ?) @" R3 x* Y else if(sum_a&gt;sum_b) {
5 Z$ Z' y6 _- q  printf("\t供大于求，增加假想销地B%d。\n",++num_b);
  b[num_b]=sum_a-sum_b;
/ c' V! N5 S7 C+ a3 Y  S3 m! V) T  for(i=1;i&lt;=num_a;i++) ab[num_b]=0;
- }8 |  ?3 K+ W9 u' U6 U }
3 x( `1 g6 t( A( n4 G; [ else if(sum_a&lt;sum_b){
3 t- s4 ]) X5 E" w* Z2 e- A2 @  printf("\t供不应求，增加假想产地A%d。\n",++num_a);
  a[num_a]=sum_b-sum_a;
: r: L3 [2 V- y& C/ d$ k  for(i=1;i&lt;=num_b;i++) ab[num_a]=0;
}
//求解前的准备
for(i=1;i&lt;=num_a;i++)
9 J- C& n6 y  W9 j  for(j=1;j&lt;=num_b;j++)
   base[j]=pbase[j]=0;
for(i=1;i&lt;=num_a;i++) temp_a=a;
for(i=1;i&lt;=num_b;i++) temp_b=b;
for(i=1;i&lt;=num_a;i++)
! l- {0 v* G+ k. b/ N7 U8 t9 L for(j=1;j&lt;=num_b;j++) temp_ab[j]=ab[j];
//回显问题
' k5 }3 I+ G( ]" }/ Z; N printf("\n\n原问题为：\n\n");
# x- l- ^$ _6 Y, S8 A   printf("产量：\n ");
for(i=1;i&lt;=num_a;i++) printf("A%d:%lg\t",i,a);
6 ~9 f, T$ U! G  Q6 z+ o( x, l& p printf("\n运量：\n ");
for(i=1;i&lt;=num_b;i++) printf("B%d:%lg\t",i,b);
printf("\n运价为：\n");
3 p- d! U% k$ P2 s- S& v& U3 i" ` for(i=1;i&lt;=num_a;i++) {
  putchar(' ');
6 }1 H" H7 x" }- p4 q  for(j=1;j&lt;=num_b;j++)
9 h) P0 Q2 M4 g9 |* B- y   printf("A%d-&gt;B%d:%lg\t",i,j,ab[j]);
  putchar('\n');
- S- Q5 E- E) U+ X- `7 i; ` }
}</P>
; L5 P- W( @- i0 Z" t, }; o<>double dabs(double d) {
if(d&gt;=0) return d;
else return -d;
}</P>
<>void findinit() {
0 b& {: T1 v* j* ?! q0 D* T2 n8 d int i,j,k;
; v6 H( c& v- ~( K double r_penum[MAX],c_penum[MAX];
: G- D' n4 J* `% f+ z struct largest {
- z8 W8 D5 y7 [2 y9 `) U* s4 j+ [' ~2 D  char rc;
$ N& u! [! _, a# P& I  int num;
  double value;
};
3 w# x0 F" A3 [1 c7 H% X struct largest lar1,lar2;
int r=0,c=0;
double a1,a2,b1,b2,temp;</P>
# {7 z0 X) N$ C7 `' @<> for(i=1;i&lt;=num_a;i++) {
  temp=temp_ab[1];
% a8 \7 D9 a6 h/ t  for(j=1;j&lt;=num_b;j++)
$ f% _$ j) I* P+ W  }3 M   if(temp_ab[j]&gt;temp) temp=temp_ab[j];
" K; u9 u7 I$ y. p4 g  a1=a2=temp;
  for(j=1;j&lt;=num_b;j++)
, s7 a- J9 g6 T) W5 B( ]1 q   if(temp_ab[j]&lt;=a1) {
$ W6 `9 w: H2 |9 ^    a2=a1;
4 s/ y4 }. Q1 @8 o" A  r    a1=temp_ab[j];
5 y3 o* o5 Y( o2 u( b   }
   else if(temp_ab[j]&lt;=a2) a2=temp_ab[j];
  r_penum=dabs(a1-a2);
" _# h) m5 n. m7 x" Q9 H }</P>
<> for(i=1;i&lt;=num_b;i++) {
% M9 W& S. K! A: [- D5 O' _  temp=temp_ab[1];
! X; o3 X' [9 K& k$ j3 E' H  for(j=1;j&lt;=num_a;j++)
   if(temp_ab[j]&gt;temp) temp=temp_ab[j];
8 ?* G3 f/ J6 T  b1=b2=temp;
3 ]* @" p. v; e2 h3 Z  for(j=1;j&lt;=num_a;j++)  
- h$ L  |7 o' ]2 ^3 x+ O; |1 _   if(temp_ab[j]&lt;=b1) {
$ z- j; X0 |3 F0 N6 B( d/ i: |    b2=b1;
& F4 l/ l( O& D$ Y    b1=temp_ab[j];
   }
   else if(temp_ab[j]&lt;=b2) b2=temp_ab[j];
  c_penum=dabs(b1-b2);
7 p- u, h" R6 ~5 x. K4 B* p }
4 n4 c% `' q% F* h /*
for(i=1;i&lt;=num_a;i++) printf("pa%d=%lg ",i,r_penum);
putchar('\n');
for(i=1;i&lt;=num_b;i++) printf("pb%d=%lg ",i,c_penum);
0 a: d% L9 O. Z& t/ b# M2 [9 _ */
temp=r_penum[1];
for(i=1;i&lt;=num_a;i++) if(r_penum&lt;temp) temp=r_penum;
: v! i3 H0 B2 |% `( ?5 I! c! e for(i=1;i&lt;=num_b;i++) if(c_penum&lt;temp) temp=c_penum;</P>
, a) D  l8 y% Z4 e- B<> //考虑了有两个罚数相等的情况，大于两个相等只取其中两个进行比较
lar1.value=lar2.value=temp;
+ y  g1 V5 @+ M+ ^" G( z lar1.num=lar2.num=1;
lar1.rc=lar2.rc='r';
for(i=1;i&lt;=num_a;i++)
  if(r_penum&gt;=lar1.value) {
   lar2=lar1;
   lar1.rc='r';
% x+ m% Q) {; [- ?9 f: X   lar1.num=i;
   lar1.value=r_penum;
  }
  else if(r_penum&gt;=lar2.value) {
$ L; f# u7 i/ q2 o2 L* q6 g$ d   lar2.rc='r';
4 y' k* n, u. n' B7 i( o" T8 d   lar2.num=i;
4 r+ k2 u/ X- i8 |" e6 O   lar2.value=r_penum;
  }</P>
<> for(i=1;i&lt;=num_b;i++)
  if(c_penum&gt;=lar1.value) {
   lar2=lar1;
' D1 Y' Q, }* z# |( F& S+ f9 P   lar1.rc='c';
; S" @1 b5 k; J% V& h  }7 Q; @   lar1.num=i;
. {1 r1 q$ Y4 f3 {   lar1.value=c_penum;
* x* |! |% n3 S' _! F  }
+ }# H1 x: J# Z7 X  else if(c_penum&gt;=lar2.value) {
& w' `7 Q  S/ x1 d( m* K  O5 K   lar2.rc='c';
0 b7 E! _1 r% c/ E3 K6 e8 G   lar2.num=i;
   lar2.value=c_penum;
  }</P>
% c( \) Q5 O0 V5 a  [5 l<> if(lar1.value==lar2.value&amp;&amp;(lar1.num!=lar2.num||lar1.rc!=lar2.rc)) {
  if(lar1.rc=='r'&amp;&amp;lar2.rc=='r') {
- _: ^* O# X  ?* C+ {- \+ q   temp=temp_ab[lar1.num][1];
   r=lar1.num;
      c=1;
) G. ]2 |" T) _% i  o2 T   for(i=1;i&lt;=num_b;i++) {
    if(temp_ab[lar1.num]&lt;temp) {
     temp=temp_ab[lar1.num];
     r=lar1.num;
     c=i;
    }
    if(temp_ab[lar2.num]&lt;temp) {
8 q# ]2 K+ s9 z. f- ~% O8 X     temp=temp_ab[lar2.num];
% }. g0 o# ^' o& v. \* k     r=lar2.num;
     c=i;
8 x/ K3 W  m. t: N. {" |* s    }
   }
  }
  if(lar1.rc=='c'&amp;&amp;lar2.rc=='c') {
; d4 o3 }) D' f6 \7 C   temp=temp_ab[1][lar1.num];
   r=1;
      c=lar1.num;
   for(i=1;i&lt;=num_a-1;i++) {
    if(temp_ab[lar1.num]&lt;temp) {
: h6 K% Q  q0 G1 W* g# t  E     temp=temp_ab[lar1.num];
     c=lar1.num;
     r=i;
    }
    if(temp_ab[lar2.num]&lt;temp) {
0 a6 U3 z* S" W7 U0 j     temp=temp_ab[lar2.num];
     c=lar2.num;
     r=i;
    }
2 @* W* x. L0 W& Z6 O  ^. D   }
$ n7 v2 n. i5 i- k! P  }
. s( I. n  p$ z& B( [# w  if(lar1.rc=='r'&amp;&amp;lar2.rc=='c') {
: |, {0 L3 z- B  f8 E8 j   temp=temp_ab[lar1.num][1];
) W' A& Z$ J* g5 n1 G" e5 {2 m$ [   r=lar1.num;
      c=1;
9 Q8 z5 L/ t7 r: x8 _   for(i=1;i&lt;=num_b;i++)
    if(temp_ab[lar1.num]&lt;temp) {
     temp=temp_ab[lar1.num];
2 h2 s/ ^- }6 l9 O1 b     r=lar1.num;
     c=i;
    }
# G1 H9 O# |5 g; H   for(i=1;i&lt;=num_a;i++)
    if(temp_ab[lar2.num]&lt;temp) {
     temp=temp_ab[lar2.num];
& E6 p7 i8 ^$ P: C$ V8 M     c=lar2.num;
; `  h8 s" M7 w, a% x     r=i;
) Z  f/ u% ^' X    }
  }
8 ?/ E% h+ ]5 c" h% _  if(lar1.rc=='c'&amp;&amp;lar2.rc=='r') {
9 k. Z' M" F: G. J- |0 `   temp=temp_ab[1][lar1.num];
   r=1;
      c=lar1.num;
   for(i=1;i&lt;=num_b;i++)
5 t6 G: u, F0 C) a    if(temp_ab[lar2.num]&lt;temp) {
     temp=temp_ab[lar2.num];
4 c0 ?- j& Y2 t4 E# y" i     r=lar2.num;
( }1 ?4 [; K/ Z  S0 w     c=i;
" P, i; ~1 y% \- K    }
% g, d2 q: ?+ @* F# D- Z   for(i=1;i&lt;=num_a;i++)
    if(temp_ab[lar1.num]&lt;temp) {
" D. o' C* q1 P# L  Z     temp=temp_ab[lar1.num];
     c=lar1.num;
' p: M$ x/ Q9 S( U% R3 r     r=i;
5 N9 ~" E' W# I0 n) {; I6 y" ^    }
  }
4 |7 [+ i4 |, P$ ?8 ^" P: Y }
else {
  if(lar1.rc=='r') {
% z2 C, t# Q2 D1 U   r=lar1.num;
, }2 m3 q& H' s- P7 _   c=1;
' {! }0 t8 c8 z( g+ Z   temp=temp_ab[lar1.num][1];
   for(i=1;i&lt;=num_b;i++)
    if(temp_ab[lar1.num]&lt;temp) {
% O4 |9 M- h9 L; U7 V+ G8 J/ [     c=i;
2 ?' ?% L$ M3 M) H2 K$ G; I! E- t     temp=temp_ab[lar1.num];
    }
  }
8 i' T6 N8 k7 R: L. f  else {
. o7 x& B3 t0 _9 `! v! x7 k   r=1;
; a5 k& M- l$ y( c% S1 A' B0 u) C   c=lar1.num;
- @$ }$ m) }! W5 x+ i   temp=temp_ab[1][lar1.num];
   for(i=1;i&lt;=num_a;i++)
    if(temp_ab[lar1.num]&lt;temp) {
     r=i;
( x! B) y" v7 H* E& j     temp=temp_ab[lar1.num];
    }
  }
}
: r& J" z) t, Q4 b" M, Z. G pbase[r][c]=1;
if(temp_a[r]&gt;temp_b[c]) {
$ i# d2 J% i4 ?+ {4 A8 s4 i4 C  base[r][c]=temp_b[c];
  temp_a[r]=temp_a[r]-temp_b[c];
% d/ Q# J& w% r  temp_b[c]=0;
( }# T9 q3 {9 V  o  for(i=1;i&lt;=num_a;i++) temp_ab[c]=M;
}
else if(temp_a[r]&lt;temp_b[c]) {
7 Y4 J4 K0 k; P7 Z" r  base[r][c]=temp_a[r];
  temp_b[c]=temp_b[c]-temp_a[r];
" {0 D9 P  {/ O6 ~  temp_a[r]=0;
! K$ ]" p3 @3 z7 T; Z' g$ B  A  for(i=1;i&lt;=num_b;i++) temp_ab[r]=M;
}
else if(temp_a[r]==temp_b[c]) {
3 a9 W/ o$ W  h* d! Q# E  base[r][c]=temp_a[r];
7 [( m) o- G7 ]* g! [7 Q, c  temp_a[r]=temp_b[c]=0;
  v' _9 P" {5 q+ {8 X2 p4 Y0 h  for(i=1;i&lt;=num_a;i++) temp_ab[c]=M;
  for(i=1;i&lt;=num_b;i++) temp_ab[r]=M;
. f! v! C$ N' u4 p  k=0;
  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;
  if(k!=0) {
   k=0;
! E! Z: h1 ~: A0 B0 ?- I, R   for(i=1;i&lt;=num_a;i++) if(pbase[c]!=1) {k=1;break;}
   for(j=1;j&lt;=num_b;j++) if(pbase[r][j]!=1) {k=2;break;}
. m0 y3 M$ B, {) u4 n8 @7 A   if(k==1) pbase[c]=1;
   else if(k==2) pbase[r][j]=1;
  }
}
: Q9 q1 C. [# q/ x, F  ^ k=0;
8 g2 f8 Z+ c& S8 j6 w9 R3 } for(i=1;i&lt;=num_a;i++) if(temp_a!=0) k=1;
% |  \* `  @' E* p0 C for(i=1;i&lt;=num_b;i++) if(temp_b!=0) k=1;
if(k==0) return;
findinit();
1 ?0 e5 h2 ^/ r* @' `# v}</P>
1 m" q) s- N1 Y0 T6 }<>void computeuv() {
int i,j,k;
for(i=1;i&lt;=num_a;i++)
( W7 Z* o* ^, D+ ?9 g$ o  for(j=1;j&lt;=num_b;j++)
   if(pbase[j]==1) {
6 u4 r) {5 m* D& _; u( j    if(pu==1) {v[j]=ab[j]-u; pv[j]=1;}
- B$ l" S5 G. _+ Y. |. Z1 p9 ?    if(pv[j]==1) {u=ab[j]-v[j]; pu=1;}
   }
k=0;
for(i=1;i&lt;=num_a;i++) if(pu==0) k=1;
for(i=1;i&lt;=num_b;i++) if(pv==0) k=1;
4 Y$ ]0 `. I8 Q if(k==1) computeuv();
}</P>
* e1 A2 \' H8 K: n<>int check() {
" f0 L' D$ Q/ p! J  d5 q int i,j,k,l;
5 ]) B2 A, \. ^: m4 H  k 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;
    else if(base[j]==0) l=1;
   }
/* for(i=1;i&lt;=num_a;i++){
" ^0 O+ x  l" L) @. ]: ?. n  for(j=1;j&lt;=num_b;j++) printf("base %lg\t",base[j]);
  putchar('\n');}
for(i=1;i&lt;=num_a;i++){
  for(j=1;j&lt;=num_b;j++) printf("pbase %d\t",pbase[j]);
2 H  [, a9 W* K: ~, k& t4 @  putchar('\n');}
7 J  ~" |! G. U, Y/ D" u. Y8 v*/
if(k==0&amp;&amp;l==0) {status=1; return 1;}
else if(k==0&amp;&amp;l!=0) {status=-1; return 1;}
8 w2 I  H% w  @* Z) P: ^- ]) g else return 0;
}</P>
<>int circle(enum direction dir,struct element cir[],int *t) {
int i,j,k;
$ W3 i1 }% j  P! `+ @, @ /*
9 O& E) M5 s' y4 r, J4 n putchar('\n');
8 q$ l* [( v# j- S for(i=1;i&lt;=*t;i++) printf("%d: r%d c%d\t",i, cir.r,cir.c);
5 K3 b/ s: k* i8 Q9 A+ [$ i putchar('\n');
*/
% L9 U7 Z. t/ t" P+ B* Q if((*t)!=1&amp;&amp;cir[(*t)].r==cir[1].r&amp;&amp;cir[*t].c==cir[1].c) {
) E. Y1 B* k- m, a! o+ i  t--;
: u! O9 L0 b6 o. g: `2 F  return 1;
1 V. f$ o& Q! B, I7 o$ k8 ^6 Y. { }
if(dir!=down) {
  k=0;
  for(i=cir[*t].r-1;i&gt;=1;i--)
, w+ n! C, U1 H7 }) v   if(pbase[cir[*t].c]==1||(i==cir[1].r&amp;&amp;cir[*t].c==cir[1].c)) {k=1; break;}
. K9 V6 V$ j9 d6 k" O& {  ^2 {& v  for(j=2;j&lt;=*t;j++) if(i==cir[j].r&amp;&amp;cir[*t].c==cir[j].c) k=0;
3 P0 S4 G8 B! T1 D3 P2 J  if(k==1) {
$ j) M/ u, q% z' ~$ w5 F& I/ @   if(dir==up) cir[*t].value=0;
   else cir[*t].value=1;
$ S. }& Y2 b) T* x; F0 x   (*t)++;
   cir[*t].r=i;
   cir[*t].c=cir[(*t)-1].c;
7 B8 G& u0 y  F: p   if(circle(up,cir,t)) return 1;
  }
5 e3 c# q( Q6 o# N: H9 w }
if(dir!=up) {
  k=0;
2 f2 v' s  G$ R* [  for(i=cir[*t].r+1;i&lt;=num_a;i++)
   if(pbase[cir[*t].c]==1||(i==cir[1].r&amp;&amp;cir[*t].c==cir[1].c)) {k=1; break;}
* a, b1 j# q: i6 P: F7 V& K  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) {
   if(dir==down) cir[*t].value=0;
4 k6 V: m4 R2 R/ j6 w" X   else cir[*t].value=1;
   (*t)++;
   cir[*t].r=i;
2 {" F! e4 o! y9 x   cir[*t].c=cir[(*t)-1].c;
+ \# j% j) i8 H/ B; v# A   if(circle(down,cir,t)) return 1;
  }
}
if(dir!=right) {
  k=0;
  for(i=cir[*t].c-1;i&gt;=1;i--)
- r# u0 D7 R: q* z2 e   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;
  if(k==1) {
1 w/ ]/ l" y4 R5 k( m7 G$ e  {   if(dir==left) cir[*t].value=0;
  N$ P! v: l& `   else cir[*t].value=1;
  q" i# X( c7 z   (*t)++;
   cir[*t].r=cir[(*t)-1].r;
   cir[*t].c=i;
   if(circle(left,cir,t)) return 1;
* k0 Q1 l9 t6 H6 ^  }
}
if(dir!=left) {
" U$ x7 Z* u1 b7 `. j  k=0;
  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;}
  for(j=2;j&lt;=*t;j++) if(cir[*t].r==cir[j].r&amp;&amp;i==cir[j].c) k=0;
6 i& m- a3 r  J- V% H9 }5 f! ^  if(k==1) {
! B. b. j4 V5 o5 i5 J   if(dir==right) cir[*t].value=0;
6 {( E) J: B7 a: _  y! I8 d/ |   else cir[*t].value=1;
   (*t)++;
   cir[*t].r=cir[(*t)-1].r;
" {5 c+ |; ^& @( j* H   cir[*t].c=i;
   if(circle(right,cir,t)) return 1;
! l+ I0 D! Y$ h- b6 C3 f3 A2 r  }
9 I0 b+ d8 g/ y7 z! B }
(*t)--;
& U7 N' @7 g* m  e- x2 o/ _ return 0;
}</P>
- ]5 h( b  ~% G# H  y<>void improve() {
int i,j,k,t;
* Q' X# X+ {) z' m struct element base_in,base_out,cir[2*MAX+1];</P>
3 D* `' x2 Q. I+ T<> for(i=1;i&lt;=num_a;i++) pu=0;
for(i=1;i&lt;=num_b;i++) pv=0;
8 N- ^: n( ]4 R7 ] u[1]=0; pu[1]=1;
' k( a4 Q7 K( U$ H! c& d computeuv();
if(check()) return;
for(i=1;i&lt;=num_a;i++)
/ |( p. w" m5 T4 T' {  for(j=1;j&lt;=num_b;j++)
   if(pbase[j]==0) break;
base_in.r=i;
. i  [) @2 }  i8 u9 B8 g/ l4 j# j+ E base_in.c=j;
, W& w$ t! e0 R( { base_in.value=base[j];</P>
<> for(i=1;i&lt;=num_a;i++)
  for(j=1;j&lt;=num_b;j++)
   if(pbase[j]==0&amp;&amp;base[j]&lt;base_in.value) {
    base_in.r=i;
    base_in.c=j;
    base_in.value=base[j];
$ p( O' i' G) W& _1 L   }</P>
<> t=1;
; b9 o( ?8 n( ~- o cir[t].r=base_in.r;
cir[t].c=base_in.c;
cir[t].value=1;
1 v  A' P% G3 K( Z  R" [" J if(!circle(mid,cir,&amp;t)) {
6 g( M  K: w& P8 g8 o6 O  printf("程序出错！按回车结束");
  getchar();
* F- I$ m- A6 M3 n/ k# C  exit(0);
# o+ f* o0 z3 S* n( T }
t--;
// for(i=1;i&lt;=t;i++) printf("%d:r%d c%d v%lg\t",i,cir.r,cir.c,cir.value);
) E2 E' D8 [1 X// putchar('\n');</P>
<> for(i=2;i&lt;=t;i++) if(cir.value==1) break;
/ H0 Y" `2 U+ Y3 f( H7 N base_out=cir;
base_out.value=base[cir.r][cir.c];
k=1;
for(i=1;i&lt;=t;i++)
2 }, x, Q& ]6 u2 z4 l  if(k%2==0&amp;&amp;cir.value==1&amp;&amp;base[cir.r][cir.c]&lt;base_out.value) {
: ~9 ?% {1 v; v( a6 T/ [5 i5 x9 X   base_out=cir;
3 W3 ?; U8 E( e9 b4 p7 M   base_out.value=base[cir.r][cir.c];
   k++;
  }
6 z! `* }2 R8 B* q$ p4 |6 w  else if(cir.value==1) k++;
base[cir[1].r][cir[1].c]=0;
* R$ _& z/ o( ] k=1;
6 @+ {  |: @9 ^& u5 J7 M6 `; \ for(i=1;i&lt;=t;i++) {
  t( P( L' s) @" j5 T6 M8 k# L9 `  if(k%2==1&amp;&amp;cir.value==1) {
$ T# _. F( T: R( }+ ?4 f* @   base[cir.r][cir.c]+=base_out.value;
   k++;
  }
9 Z) h! X6 c8 |2 }5 A  else if(k%2==0&amp;&amp;cir.value==1) {
   base[cir.r][cir.c]-=base_out.value;
   k++;
# d/ J  Y: ?: c& T/ C. z  }
8 n3 v$ k. v& _: b; s, _- a }
( I% z; _! q2 C  @2 O" Q pbase[base_out.r][base_out.c]=0;
pbase[base_in.r][base_in.c]=1;</P>
<> improve();
/ r; D& j( h6 |! C}</P>
<>void output() {
$ J$ s" A. A. C int i,j;
double sum=0;
- w- T# |* B, f' C' ~3 { printf("\n运量为：\n");
for(i=1;i&lt;=num_a;i++) {
6 B% i, z3 D* E6 }  putchar(' ');
# L" s; @- q* j+ A, O& P  for(j=1;j&lt;=num_b;j++)
   if(pbase[j]==1) {
    printf("A%d-&gt;B%d:%lg\t",i,j,base[j]);
0 ]- F& y+ I) l    sum+=base[j]*ab[j];
   }
4 Z. M1 j+ d, U4 X   else printf("A%d-&gt;B%d:0\t",i,j);
  putchar('\n');
}
printf("\n最低总运费为：%lg\n",sum);
- }/ ^( G, c( j' b/ j}
</P>

zhongguohelin

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

enrique08

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

hzlhm

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

liunengwu

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