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

plgatc

<>/*************************************************************************
( g3 _- G1 f8 `                         表上作业法解运输问题        
% v9 n. a3 Y  K: S) U& T
7 |9 ^) @" p" P8 {3 n2 a" N, C8 P 编程环境：VC++6.0     
0 [- r- Z2 p, b$ ~+ T5 `  }  D 程序说明：
/ _& N/ X0 x7 c1 ?5 Y Vogel法求初始解，位势法检验是否最优，递归搜索闭回路改进解。
*************************************************************************/
& Y9 ^; ?; o2 O  M2 N#include &lt;stdio.h&gt;
2 Z- f9 }4 C$ H5 ?2 v, T6 C#include &lt;stdlib.h&gt;</P>
<>#define MAX 100
+ Y" `$ t/ v/ u8 {' d$ S: I#define M 999999</P>
) W* A3 a6 M! t, a<>int status; //1唯一最优解，-1无穷多最优解</P>
<>int num_a,num_b;
. \# }& q7 n4 n- F/ P0 ydouble a[MAX],b[MAX],temp_a[MAX],temp_b[MAX];
double u[MAX],v[MAX],pu[MAX],pv[MAX];
8 Z/ g+ k+ i# ndouble ab[MAX][MAX],temp_ab[MAX][MAX];
) m4 R6 K: p) c7 G0 I" Rdouble base[MAX][MAX];
: H- e4 B7 ^) d; D2 y) Gint pbase[MAX][MAX];
3 @2 i- h: I) {0 w4 w* U  K/ d3 `struct element {
1 V  f" `3 y2 Z) a( d int r;
# q; C- y9 h6 S% q# `. A1 V1 n int c;
+ I8 @& t$ n  }  y, q double value;
};
enum direction {mid,up,down,left,right};</P>
<>void create();
void banner();
double dabs(double d);
void findinit();
void computeuv();
int check();
int  circle(enum direction dir,struct element cir[],int *t);
void improve();
void output();</P>
<>void main() {
banner();
create();
printf("\n按回车确认后开始求解");
( e/ B! K* f. S: Q4 H: g) E; U4 v. b, N getchar();
8 f- P; v' y0 V8 R9 z3 W; U getchar();
2 _$ r- i, u' i6 j findinit();
improve();
switch(status) {
6 x4 H8 q" E9 M2 y0 A. ^9 k2 j9 t case 1:
  output();
! w3 I# l3 e/ |6 u- L4 A  puts("\n原问题有唯一最优解。\n");
  puts("\n按回车结束");
" J' M% h" ?& p9 j  getchar();
  exit(0);
* x- ~3 ^$ W' b& g1 F7 _ case -1:
  output();
  puts("\n原问题有无穷多最优解。\n");
' B3 Y& z, Q* F) l  puts("\n按回车结束");
  getchar();
  exit(0);
}
}</P>
<>void banner() {
, C; `/ q/ N3 E! U% H+ e, z7 K printf("\t\t****************************************\n");
printf("\t\t         表上作业法解运输问题\n");
% D) q, x/ w$ S* Q& I5 c printf("\t\t                              Thunder\n");
1 A4 M) F) }( O printf("\t\t****************************************\n");
$ s3 E# O' K4 ]5 v/ \8 v printf("\n");
7 o; T' W9 [$ ]! N% M- L}</P>
7 g# I0 T6 j$ n8 C8 y$ y/ ]<>void create() {
, U/ m, I+ y- b) N+ q3 i. ~. [ int i,j;
8 u1 F, O5 W$ A double sum_a,sum_b;
6 w3 O% v5 A6 v! J5 u( p char confirm;
" Q' i; X! u, o' J while(1) {
1 M. N8 Z/ Y  i* f& `  printf("请输入产地个数：");
5 S8 L- ?8 A. X2 o: Q9 K  scanf("%d",&amp;num_a);
0 b( v  L7 B( V$ m; F  for(i=1;i&lt;=num_a;i++) {
- {" E, |; ~7 m" H! T# Z& S   printf("A%d产量：",i);
* F; ?# m# ^; ~: b& Z7 p8 q/ Y   scanf("%lf",&amp;a);
  }
: N. v& ]9 Q; q8 E' o  l  for(i=1;i&lt;=num_a;i++) printf("A%d产量：%lg\t",i,a);
0 p# \5 i7 V6 P  printf("\n正确吗？：(y/n)");
  getchar();
; w: U7 }& z2 F. a) m+ J  confirm=getchar();
  if (confirm=='y') break;
& M) [6 S0 \5 O% D, M6 l  else if(confirm=='n') continue;
) g- k* x: z" C5 G% Z# G9 u; o }
" v, \" |9 H2 y0 | while(1) {
! f! x% `, y& P4 g  printf("\n请输入销地个数：");
  scanf("%d",&amp;num_b);
: v; n6 R6 |" M* }* z7 F. N' o. _  for(i=1;i&lt;=num_b;i++) {
   printf("B%d销量：",i);
) _7 c; y" [9 F+ g2 H7 P1 M* [   scanf("%lf",&amp;b);
0 |! t' E& B0 _' Q) g% w, o# o8 ~  }
" o9 Q1 O* b5 w$ o  for(i=1;i&lt;=num_b;i++) printf("B%d销量：%lg\t",i,b);
0 y9 w( X! X- @  printf("\n正确吗？：(y/n)");
  getchar();
. r  l4 H2 O" T# x8 M+ q  confirm=getchar();
: V1 d8 v2 {5 X& F! b. G  if (confirm=='y') break;
/ N) x3 b, w; L* O: `& m; l  else if(confirm=='n') continue;
}
putchar('\n');
for(i=1;i&lt;=num_a;i++) {
  for(j=1;j&lt;=num_b;j++) {
& [' R9 q; l& e( I4 @$ a   printf("A%d到B%d运价：",i,j);
/ x0 Y# _; Q3 L9 T$ u6 |8 a) w   scanf("%lf",&amp;ab[j]);
$ L3 P. O# X- r& n  }
# i; m* a& E/ o; Q! @# I  k  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--;
  putchar('\n');
# z- J/ g+ R8 r$ m/ F/ H7 p }
4 i; }9 ^; Z- W' q3 p //处理产销不平衡的情况
1 ]2 j# {& \( ~1 z sum_a=sum_b=0;
, M7 _# V7 U! V/ `, I for(i=1;i&lt;=num_a;i++) sum_a+=a;
for(i=1;i&lt;=num_b;i++) sum_b+=b;</P>
<> printf("总产量：%lg\t总销量：%lg",sum_a,sum_b);
if(sum_a==sum_b) printf("\t产销平衡。");
else if(sum_a&gt;sum_b) {
  printf("\t供大于求，增加假想销地B%d。\n",++num_b);
( a8 ?3 v! T# F  E  b[num_b]=sum_a-sum_b;
  for(i=1;i&lt;=num_a;i++) ab[num_b]=0;
" S: y* `% y/ d$ ?  e* z0 U0 \ }
5 f1 ^% V1 v: i5 A8 i else if(sum_a&lt;sum_b){
5 z6 e0 _  E" v9 i8 j7 I  printf("\t供不应求，增加假想产地A%d。\n",++num_a);
  a[num_a]=sum_b-sum_a;
% H4 ]# o  Q8 G/ {+ f+ j  for(i=1;i&lt;=num_b;i++) ab[num_a]=0;
2 {& ~/ W* B6 s" N }
//求解前的准备
0 t+ W. p. ^0 F7 H for(i=1;i&lt;=num_a;i++)
1 q2 o  _* [  K) Q2 i  for(j=1;j&lt;=num_b;j++)
0 ^% t* g- |- Z( j   base[j]=pbase[j]=0;
/ ^6 b7 U6 J, O0 v' ?! h3 {  ?( ~ for(i=1;i&lt;=num_a;i++) temp_a=a;
for(i=1;i&lt;=num_b;i++) temp_b=b;
7 J- j, u' u" E: @, C( Q9 h for(i=1;i&lt;=num_a;i++)
for(j=1;j&lt;=num_b;j++) temp_ab[j]=ab[j];
//回显问题
: v; G) t) {5 `1 w+ o- A printf("\n\n原问题为：\n\n");
   printf("产量：\n ");
for(i=1;i&lt;=num_a;i++) printf("A%d:%lg\t",i,a);
printf("\n运量：\n ");
for(i=1;i&lt;=num_b;i++) printf("B%d:%lg\t",i,b);
printf("\n运价为：\n");
for(i=1;i&lt;=num_a;i++) {
6 D$ G& V5 R( c4 H9 d! v# J9 Z  putchar(' ');
) Z3 {7 i% |- [5 e" v9 }  for(j=1;j&lt;=num_b;j++)
, }/ `. C' y& d8 {4 ^   printf("A%d-&gt;B%d:%lg\t",i,j,ab[j]);
  putchar('\n');
8 w3 z6 @# h' k8 B }
' X* Y% O2 [( ?+ ?$ a}</P>
<>double dabs(double d) {
if(d&gt;=0) return d;
else return -d;
# R; V- A, N- h$ P}</P>
<>void findinit() {
int i,j,k;
' U, }& d, \; ^- I double r_penum[MAX],c_penum[MAX];
struct largest {
  char rc;
0 p- S8 z7 t% |) n0 \+ F2 `  int num;
  double value;
};
struct largest lar1,lar2;
/ ]- g3 c6 {4 M9 |. W int r=0,c=0;
. l8 b: {- @0 u  O double a1,a2,b1,b2,temp;</P>
<> for(i=1;i&lt;=num_a;i++) {
  temp=temp_ab[1];
  for(j=1;j&lt;=num_b;j++)
9 Z+ F8 @: _. w$ ]" R4 J. R! X& C   if(temp_ab[j]&gt;temp) temp=temp_ab[j];
  a1=a2=temp;
8 b% H, |# R* T  for(j=1;j&lt;=num_b;j++)
$ J4 j* X# K; Y; M7 y" J   if(temp_ab[j]&lt;=a1) {
5 F  b: d; {, {. D1 b    a2=a1;
    a1=temp_ab[j];
: b' q4 I% G( i  Y  u6 p6 M   }
   else if(temp_ab[j]&lt;=a2) a2=temp_ab[j];
" j( r% A6 O- Y5 p/ [  r_penum=dabs(a1-a2);
, z' g8 d) R, S+ n }</P>
<> for(i=1;i&lt;=num_b;i++) {
' |4 f5 I' W# D2 l% A/ |  temp=temp_ab[1];
# T0 ^+ g7 y$ ?- Q  for(j=1;j&lt;=num_a;j++)
1 M! i; H+ U( F) g  @   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;
    b1=temp_ab[j];
! e$ X# }0 [3 _( A/ _   }
, g9 X) s- f; L( F$ G: K   else if(temp_ab[j]&lt;=b2) b2=temp_ab[j];
: R" b7 @5 _1 R: p  c_penum=dabs(b1-b2);
* E  F  N4 X+ Z! `9 F8 c }
/*
for(i=1;i&lt;=num_a;i++) printf("pa%d=%lg ",i,r_penum);
putchar('\n');
) [. s, m3 L( I( Z5 m; g for(i=1;i&lt;=num_b;i++) printf("pb%d=%lg ",i,c_penum);
! ~( H3 m! }3 ~9 `3 y7 c */
5 @4 E' [) `5 \3 K  L 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>
$ \  t! t9 u2 e4 n( b<> //考虑了有两个罚数相等的情况，大于两个相等只取其中两个进行比较
lar1.value=lar2.value=temp;
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';
   lar1.num=i;
6 m, X4 x% \( }, f  }# w   lar1.value=r_penum;
  }
  else if(r_penum&gt;=lar2.value) {
1 l- V4 p2 W# M2 _4 A6 d* g: w   lar2.rc='r';
   lar2.num=i;
   lar2.value=r_penum;
  x' u& r; h& r( r4 R  t/ ^, S% w  }</P>
. k# f& a. ^! p9 t/ q4 G: ]<> for(i=1;i&lt;=num_b;i++)
  if(c_penum&gt;=lar1.value) {
   lar2=lar1;
   lar1.rc='c';
   lar1.num=i;
  p% V9 V" G1 X% T! K   lar1.value=c_penum;
. F# O; Z9 c" |# N  }
  else if(c_penum&gt;=lar2.value) {
4 n7 P+ |7 l. R0 |; y9 F   lar2.rc='c';
   lar2.num=i;
   lar2.value=c_penum;
  }</P>
9 v7 B! z& n7 n; l<> if(lar1.value==lar2.value&amp;&amp;(lar1.num!=lar2.num||lar1.rc!=lar2.rc)) {
4 i9 f8 Y9 A: j7 T) ]  r6 U  if(lar1.rc=='r'&amp;&amp;lar2.rc=='r') {
   temp=temp_ab[lar1.num][1];
   r=lar1.num;
      c=1;
8 z6 H" {# k! b4 {$ g, n. T+ k   for(i=1;i&lt;=num_b;i++) {
    if(temp_ab[lar1.num]&lt;temp) {
     temp=temp_ab[lar1.num];
$ J! i8 A# m* w9 r, ]) ?& s$ N% l     r=lar1.num;
     c=i;
    }
; k1 O' _' B% l/ `    if(temp_ab[lar2.num]&lt;temp) {
     temp=temp_ab[lar2.num];
. a& C5 h. z/ @# p. s- E& r1 j     r=lar2.num;
     c=i;
5 ~! V3 h; `4 ]( `. e    }
   }
9 N) l/ j/ ?. C3 ^( R4 _  }
, ~# b7 }. H/ o, A8 _8 e1 Z  if(lar1.rc=='c'&amp;&amp;lar2.rc=='c') {
   temp=temp_ab[1][lar1.num];
& J7 b) J( _  E8 ?" K# J8 ]% o, W   r=1;
+ B; O% E# w+ K- W' D% D% S  k      c=lar1.num;
- T1 n$ W9 \  J  B7 {   for(i=1;i&lt;=num_a-1;i++) {
1 y5 e* Z2 x4 T" |    if(temp_ab[lar1.num]&lt;temp) {
: Q" |5 s# i* Y0 W+ D     temp=temp_ab[lar1.num];
& I4 C2 j& L( S4 b% t0 k     c=lar1.num;
     r=i;
    }
* n# D3 }0 Y: f    if(temp_ab[lar2.num]&lt;temp) {
     temp=temp_ab[lar2.num];
) J" G5 E( @- Q5 \( B! A     c=lar2.num;
     r=i;
    }
   }
  }
3 I3 s$ b, l( W; D6 J5 l. K  if(lar1.rc=='r'&amp;&amp;lar2.rc=='c') {
" U: n$ s: e- J3 D   temp=temp_ab[lar1.num][1];
" |6 z( E5 p+ ?  n$ p2 K4 O   r=lar1.num;
; x! l9 }% _* [- t! d      c=1;
5 M/ q! o  J, g  P4 v& g8 v   for(i=1;i&lt;=num_b;i++)
" [* c& e4 X9 E* ~    if(temp_ab[lar1.num]&lt;temp) {
     temp=temp_ab[lar1.num];
& i* B" `# s, k0 E: g     r=lar1.num;
; Q: w/ c9 N9 W& _5 e     c=i;
    }
   for(i=1;i&lt;=num_a;i++)
    if(temp_ab[lar2.num]&lt;temp) {
     temp=temp_ab[lar2.num];
     c=lar2.num;
1 `9 h' h- U+ u9 K" M  p/ X) |7 g     r=i;
    }
  }
  if(lar1.rc=='c'&amp;&amp;lar2.rc=='r') {
   temp=temp_ab[1][lar1.num];
& T  w9 F& X2 b   r=1;
      c=lar1.num;
# Y7 z2 A; M7 q   for(i=1;i&lt;=num_b;i++)
- e/ f8 z6 h7 E$ e$ k    if(temp_ab[lar2.num]&lt;temp) {
* d! Q- Z3 }4 |; }     temp=temp_ab[lar2.num];
6 s' ?$ e' [" }% [+ _     r=lar2.num;
/ P, H6 B" l) a/ X  O3 U; S/ s     c=i;
    }
   for(i=1;i&lt;=num_a;i++)
4 E9 T$ a; U) H3 p  `' i% o    if(temp_ab[lar1.num]&lt;temp) {
. a, a. d7 `# N- ~, a  U     temp=temp_ab[lar1.num];
     c=lar1.num;
& u2 ~8 |* U' I$ x0 r+ |/ i     r=i;
    }
  }
+ z+ X  V. H7 ~& [$ R  N8 [ }
# W* H1 ?* s" b  G2 m9 }- g( z else {
1 V9 g! z& I- v- j  if(lar1.rc=='r') {
   r=lar1.num;
   c=1;
   temp=temp_ab[lar1.num][1];
+ A" z: R6 X$ g, Z7 y8 u" p5 m   for(i=1;i&lt;=num_b;i++)
    if(temp_ab[lar1.num]&lt;temp) {
# f8 d5 @& r% Q3 p5 k* n" D0 G     c=i;
     temp=temp_ab[lar1.num];
6 o# j* d! d2 H    }
. Z  X) S. Q& |/ z  M+ ~* }! a  }
3 B, Q9 |& u/ a* }- c  else {
$ _1 ^+ {& g$ x- J# _1 R7 m) Q   r=1;
   c=lar1.num;
5 r5 M  W* G" z: ]2 g7 K! C   temp=temp_ab[1][lar1.num];
  @8 C- o# f4 t3 L, t! K   for(i=1;i&lt;=num_a;i++)
    if(temp_ab[lar1.num]&lt;temp) {
     r=i;
/ k" g0 d: U5 ^& L8 d! w) k     temp=temp_ab[lar1.num];
    }
: `3 m3 ]3 c' S4 E) G7 {- m" A* l  }
2 `- v% G! P" ^& a: J }
% X. b, _& `; v pbase[r][c]=1;
/ c' G0 a. U: _9 ]: T! Q if(temp_a[r]&gt;temp_b[c]) {
  base[r][c]=temp_b[c];
) s! b/ n( |( ~. W  temp_a[r]=temp_a[r]-temp_b[c];
  temp_b[c]=0;
  for(i=1;i&lt;=num_a;i++) temp_ab[c]=M;
+ E" w* j0 q* X9 u- X }
else if(temp_a[r]&lt;temp_b[c]) {
  base[r][c]=temp_a[r];
  temp_b[c]=temp_b[c]-temp_a[r];
  temp_a[r]=0;
  for(i=1;i&lt;=num_b;i++) temp_ab[r]=M;
}
3 E+ x' e, X: c else if(temp_a[r]==temp_b[c]) {
$ T- b+ h9 r8 z  E, @  base[r][c]=temp_a[r];
  temp_a[r]=temp_b[c]=0;
1 X; p7 d/ g6 t  for(i=1;i&lt;=num_a;i++) temp_ab[c]=M;
  for(i=1;i&lt;=num_b;i++) temp_ab[r]=M;
1 s/ R# s: {4 P1 l/ N  k=0;
; ~/ ~4 B) ^( i" q  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) {
4 I- f% m9 G5 R8 O+ V! i   k=0;
, ~2 ~5 L% X. k" f0 T; k   for(i=1;i&lt;=num_a;i++) if(pbase[c]!=1) {k=1;break;}
4 f3 X# d& G, ^1 P. p   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;
# s9 F& \% b; ?, `% E* Y  z  }
}
" ~- ?( N/ A' [8 u8 M& l/ m 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) return;
findinit();
}</P>
<>void computeuv() {
int i,j,k;
for(i=1;i&lt;=num_a;i++)
  for(j=1;j&lt;=num_b;j++)
3 M) u% f% p$ w- _+ m$ u( z   if(pbase[j]==1) {
! J  _& M: h, X5 ]# I$ C. ^1 |! D    if(pu==1) {v[j]=ab[j]-u; pv[j]=1;}
    if(pv[j]==1) {u=ab[j]-v[j]; pu=1;}
   }
% o( n1 T7 ?: x/ l& `1 ]; J( F4 X& L k=0;
1 d. B% A, W7 @4 a5 s. | for(i=1;i&lt;=num_a;i++) if(pu==0) k=1;
7 ~$ t: H$ ?6 [$ A# n: Z for(i=1;i&lt;=num_b;i++) if(pv==0) k=1;
if(k==1) computeuv();
}</P>
<>int check() {
  L: U" u& d5 @* B0 R. L int i,j,k,l;
k=l=0;
4 _/ I6 h# X+ B9 Y8 ?' ] for(i=1;i&lt;=num_a;i++)
  for(j=1;j&lt;=num_b;j++)
, d) `4 s, ~- _5 h7 t$ k   if(pbase[j]==0) {
) s3 J% o; ~/ X( S9 a    base[j]=ab[j]-u-v[j];
    if(base[j]&lt;0) k=1;
    else if(base[j]==0) l=1;
6 T) w1 N  p; X% I0 J9 f   }
% _+ s) E7 K3 ~$ U) g  p8 w1 W/* for(i=1;i&lt;=num_a;i++){
  for(j=1;j&lt;=num_b;j++) printf("base %lg\t",base[j]);
5 O! I0 ~) w+ I4 r' M$ Q0 ?  putchar('\n');}
. s7 L. S5 z" _, P* o+ | for(i=1;i&lt;=num_a;i++){
  for(j=1;j&lt;=num_b;j++) printf("pbase %d\t",pbase[j]);
  putchar('\n');}
, V5 n) Q$ U; _8 Y2 P*/
if(k==0&amp;&amp;l==0) {status=1; return 1;}
else if(k==0&amp;&amp;l!=0) {status=-1; return 1;}
, G1 A# O4 w! m else return 0;
2 @& y! V! s  C6 o! i  s9 c4 Q% d}</P>
5 b8 F6 L9 K+ C% T& C7 J# b<>int circle(enum direction dir,struct element cir[],int *t) {
int i,j,k;
/*
putchar('\n');
for(i=1;i&lt;=*t;i++) printf("%d: r%d c%d\t",i, cir.r,cir.c);
putchar('\n');
' P1 C  v. K' d# j  }( P$ s */
if((*t)!=1&amp;&amp;cir[(*t)].r==cir[1].r&amp;&amp;cir[*t].c==cir[1].c) {
/ ~1 N0 T" D  o0 t$ O  t--;
) k* a4 m7 I" I  return 1;
& m0 H! U2 h2 p( l: P0 ?8 D }
) _) @9 i; ?! g* ] if(dir!=down) {
* j9 n, S& n; V% U% B$ Y  k=0;
6 o& b: K  \6 L  Q: C. P4 {  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;}
; x( A6 N1 l! X1 s+ j  for(j=2;j&lt;=*t;j++) if(i==cir[j].r&amp;&amp;cir[*t].c==cir[j].c) k=0;
! ?( [$ Z( E- L  if(k==1) {
   if(dir==up) cir[*t].value=0;
   else cir[*t].value=1;
   (*t)++;
3 L) x4 T& e( V2 K1 N/ P   cir[*t].r=i;
# {, [6 i0 \2 I- Q   cir[*t].c=cir[(*t)-1].c;
   if(circle(up,cir,t)) return 1;
  }
}
% @% A, |  M2 [# q if(dir!=up) {
2 B  ^3 `3 I% R4 C( K/ i1 Z  k=0;
' E/ h9 h6 B0 t) n/ @$ a  for(i=cir[*t].r+1;i&lt;=num_a;i++)
7 v) x( B+ D2 b# G% G$ k! {0 ^   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) {
   if(dir==down) cir[*t].value=0;
   else cir[*t].value=1;
   (*t)++;
   cir[*t].r=i;
) H! I# I) L2 H   cir[*t].c=cir[(*t)-1].c;
   if(circle(down,cir,t)) return 1;
, b' F. k6 s- ?" d  }
}
if(dir!=right) {
- F9 M* y4 N; y. T4 n) e% S- l  k=0;
  for(i=cir[*t].c-1;i&gt;=1;i--)
& P1 m. H0 G" O& S* Q   if(pbase[cir[*t].r]==1||(cir[*t].r==cir[1].r&amp;&amp;i==cir[1].c)) {k=1; break;}
0 l/ o$ E* k$ t1 K  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) {
* Q5 e* P! h" n- J; w4 j! a   if(dir==left) cir[*t].value=0;
   else cir[*t].value=1;
7 A8 k3 I! k! e( t) J8 H   (*t)++;
) Y# S% B: `: _6 y! I6 N$ }& H   cir[*t].r=cir[(*t)-1].r;
4 y7 \& A* D* A& v$ y   cir[*t].c=i;
   if(circle(left,cir,t)) return 1;
/ S2 H$ d4 j4 b5 U% }. e* X  }
4 u$ B2 |+ `2 e, ^5 D; i }
if(dir!=left) {
  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;}
) r, b' ^$ X* n* g6 E$ q$ L# s  for(j=2;j&lt;=*t;j++) if(cir[*t].r==cir[j].r&amp;&amp;i==cir[j].c) k=0;
; E: w; {" m  A) R5 D+ J  if(k==1) {
0 i- p6 P/ T- o4 j   if(dir==right) cir[*t].value=0;
7 E+ y, S' k( {5 _   else cir[*t].value=1;
   (*t)++;
- q+ n- I$ c. n1 R   cir[*t].r=cir[(*t)-1].r;
+ S5 s. D4 [/ `# p   cir[*t].c=i;
' R7 u- _3 ~6 o/ k8 a   if(circle(right,cir,t)) return 1;
: B4 \8 C, L/ [6 L$ p! ^  }
  m% d- i9 B7 A }
( v% }! Z) a9 X3 r: M (*t)--;
return 0;
}</P>
<>void improve() {
int i,j,k,t;
8 Q" M' z. E$ R. i- `' p# Q. f struct element base_in,base_out,cir[2*MAX+1];</P>
  v' F5 P" U' Q- ]) r" V<> for(i=1;i&lt;=num_a;i++) pu=0;
  N# Q0 ]0 e9 q! d3 i8 _$ M3 b for(i=1;i&lt;=num_b;i++) pv=0;
1 J- R9 v5 y7 i% j" @  k u[1]=0; pu[1]=1;
computeuv();
if(check()) return;
for(i=1;i&lt;=num_a;i++)
  for(j=1;j&lt;=num_b;j++)
   if(pbase[j]==0) break;
8 S8 U4 E( A6 [4 j+ |$ O base_in.r=i;
base_in.c=j;
base_in.value=base[j];</P>
/ D! j% e+ G$ |) ?- A' _4 p$ b! k# H<> for(i=1;i&lt;=num_a;i++)
2 ~2 X' `. ~1 E, i- W  for(j=1;j&lt;=num_b;j++)
- z5 q& I9 F" M- ~9 \   if(pbase[j]==0&amp;&amp;base[j]&lt;base_in.value) {
" U* l, e) d: I; Y' h    base_in.r=i;
    base_in.c=j;
    base_in.value=base[j];
( _. X, u) ^# o- K# y   }</P>
<> t=1;
cir[t].r=base_in.r;
cir[t].c=base_in.c;
cir[t].value=1;
9 j, x0 t' }; m# R' _, D if(!circle(mid,cir,&amp;t)) {
3 i4 n) @4 V" l. n" A  printf("程序出错！按回车结束");
  getchar();
! M& w; G& L% a& l% A4 X8 c  exit(0);
2 ^, v( S( z1 C }
& I- w9 a& D4 e+ c6 y8 _. [: g 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>
<> for(i=2;i&lt;=t;i++) if(cir.value==1) break;
$ L2 H: j- U' G6 u' A/ A$ F1 m6 v base_out=cir;
base_out.value=base[cir.r][cir.c];
k=1;
, _5 B0 U; r- [7 `- L3 G for(i=1;i&lt;=t;i++)
( L, X/ V$ v- O7 c% A* R4 U& t  if(k%2==0&amp;&amp;cir.value==1&amp;&amp;base[cir.r][cir.c]&lt;base_out.value) {
   base_out=cir;
   base_out.value=base[cir.r][cir.c];
, p2 e" y9 x: _) G   k++;
# n2 B* T7 L9 d# c. @  }
4 E' }" k* _& J  else if(cir.value==1) k++;
5 }2 M3 Q3 r# \ base[cir[1].r][cir[1].c]=0;
1 A+ Z1 ~  w" } k=1;
4 O6 A; k' a8 ?( {& v0 _ for(i=1;i&lt;=t;i++) {
8 ?  F2 B$ O7 U1 f7 G) _, O3 n5 Q  if(k%2==1&amp;&amp;cir.value==1) {
   base[cir.r][cir.c]+=base_out.value;
0 c8 M4 y3 R# p0 q  Q9 Q) J4 ]   k++;
3 l5 @6 @& T/ u  M5 L& p+ j  }
  ~, K9 f6 Z. }$ R+ C2 v1 W" ^  v) ~  else if(k%2==0&amp;&amp;cir.value==1) {
   base[cir.r][cir.c]-=base_out.value;
7 c) O  H  _. h' r/ V9 i$ @7 R( X! X   k++;
, n2 R- ~) l: \+ P: _  }
+ r$ ]4 Q2 _: ^ }
pbase[base_out.r][base_out.c]=0;
  L+ u( ^$ d4 S% M2 z pbase[base_in.r][base_in.c]=1;</P>
<> improve();
}</P>
<>void output() {
" ]! a/ J: ?+ [1 \" J% r% q, f+ D int i,j;
7 k) Q+ m( i, v! q; K/ [ double sum=0;
printf("\n运量为：\n");
& W, z; J7 J! }; I+ V: d2 k0 @0 \ for(i=1;i&lt;=num_a;i++) {
% c* I; A9 l* r1 A# J/ o  putchar(' ');
  for(j=1;j&lt;=num_b;j++)
  J# C( f6 e/ A4 O   if(pbase[j]==1) {
# t; {* m2 {* K% U/ L, {) U    printf("A%d-&gt;B%d:%lg\t",i,j,base[j]);
    sum+=base[j]*ab[j];
: \- o/ U# I+ E: Q& Y   }
   else printf("A%d-&gt;B%d:0\t",i,j);
/ ]4 w0 a* q. P; V7 I$ q  putchar('\n');
! I2 |0 o0 C% U }
9 ~6 T" T* k# X printf("\n最低总运费为：%lg\n",sum);
}
</P>

zhongguohelin

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

enrique08

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

hzlhm

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

liunengwu

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