Hex Artifact Content
Not logged in

Artifact 7f6a0aa40df322bf38db67de0b3f3f7c6f23d5c0:


0000: 23 69 6e 63 6c 75 64 65 20 3c 69 6f 73 74 72 65  #include <iostre
0010: 61 6d 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 73  am>..#include <s
0020: 73 74 72 65 61 6d 3e 0d 0a 23 69 6e 63 6c 75 64  stream>..#includ
0030: 65 20 3c 69 6f 6d 61 6e 69 70 3e 0d 0a 23 69 6e  e <iomanip>..#in
0040: 63 6c 75 64 65 20 3c 76 65 63 74 6f 72 3e 0d 0a  clude <vector>..
0050: 23 69 6e 63 6c 75 64 65 20 3c 73 74 72 69 6e 67  #include <string
0060: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 6d 61 70  >..#include <map
0070: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 73 65 74  >..#include <set
0080: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 61 6c 67  >..#include <alg
0090: 6f 72 69 74 68 6d 3e 0d 0a 23 69 6e 63 6c 75 64  orithm>..#includ
00a0: 65 20 3c 6e 75 6d 65 72 69 63 3e 0d 0a 23 69 6e  e <numeric>..#in
00b0: 63 6c 75 64 65 20 3c 69 74 65 72 61 74 6f 72 3e  clude <iterator>
00c0: 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 66 75 6e 63  ..#include <func
00d0: 74 69 6f 6e 61 6c 3e 0d 0a 23 69 6e 63 6c 75 64  tional>..#includ
00e0: 65 20 3c 63 6f 6d 70 6c 65 78 3e 0d 0a 23 69 6e  e <complex>..#in
00f0: 63 6c 75 64 65 20 3c 71 75 65 75 65 3e 0d 0a 23  clude <queue>..#
0100: 69 6e 63 6c 75 64 65 20 3c 73 74 61 63 6b 3e 0d  include <stack>.
0110: 0a 23 69 6e 63 6c 75 64 65 20 3c 63 6d 61 74 68  .#include <cmath
0120: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 63 61 73  >..#include <cas
0130: 73 65 72 74 3e 0d 0a 23 69 6e 63 6c 75 64 65 20  sert>..#include 
0140: 3c 63 73 74 72 69 6e 67 3e 0d 0a 75 73 69 6e 67  <cstring>..using
0150: 20 6e 61 6d 65 73 70 61 63 65 20 73 74 64 3b 0d   namespace std;.
0160: 0a 74 79 70 65 64 65 66 20 6c 6f 6e 67 20 6c 6f  .typedef long lo
0170: 6e 67 20 4c 4c 3b 0d 0a 74 79 70 65 64 65 66 20  ng LL;..typedef 
0180: 63 6f 6d 70 6c 65 78 3c 64 6f 75 62 6c 65 3e 20  complex<double> 
0190: 43 4d 50 3b 0d 0a 0d 0a 63 6c 61 73 73 20 53 61  CMP;....class Sa
01a0: 76 65 54 68 65 54 72 65 65 73 20 7b 0d 0a 70 75  veTheTrees {..pu
01b0: 62 6c 69 63 3a 0d 0a 09 69 6e 74 20 6d 69 6e 69  blic:...int mini
01c0: 6d 75 6d 43 75 74 28 76 65 63 74 6f 72 20 3c 69  mumCut(vector <i
01d0: 6e 74 3e 20 78 2c 20 76 65 63 74 6f 72 20 3c 69  nt> x, vector <i
01e0: 6e 74 3e 20 79 2c 20 76 65 63 74 6f 72 20 3c 69  nt> y, vector <i
01f0: 6e 74 3e 20 68 29 20 0d 0a 09 7b 0d 0a 09 09 69  nt> h) ...{....i
0200: 6e 74 20 4e 20 3d 20 78 2e 73 69 7a 65 28 29 3b  nt N = x.size();
0210: 0d 0a 09 09 69 6e 74 20 61 6e 73 20 3d 20 4e 2d  ....int ans = N-
0220: 31 3b 0d 0a 09 09 66 6f 72 28 69 6e 74 20 69 31  1;....for(int i1
0230: 3d 30 3b 20 69 31 3c 4e 3b 20 2b 2b 69 31 29 0d  =0; i1<N; ++i1).
0240: 0a 09 09 66 6f 72 28 69 6e 74 20 69 32 3d 30 3b  ...for(int i2=0;
0250: 20 69 32 3c 4e 3b 20 2b 2b 69 32 29 0d 0a 09 09   i2<N; ++i2)....
0260: 66 6f 72 28 69 6e 74 20 6b 31 3d 30 3b 20 6b 31  for(int k1=0; k1
0270: 3c 4e 3b 20 2b 2b 6b 31 29 0d 0a 09 09 66 6f 72  <N; ++k1)....for
0280: 28 69 6e 74 20 6b 32 3d 30 3b 20 6b 32 3c 4e 3b  (int k2=0; k2<N;
0290: 20 2b 2b 6b 32 29 20 7b 0d 0a 09 09 09 69 6e 74   ++k2) {.....int
02a0: 20 6d 78 20 3d 20 6d 69 6e 28 78 5b 69 31 5d 2c   mx = min(x[i1],
02b0: 20 78 5b 69 32 5d 29 3b 0d 0a 09 09 09 69 6e 74   x[i2]);.....int
02c0: 20 4d 78 20 3d 20 6d 61 78 28 78 5b 69 31 5d 2c   Mx = max(x[i1],
02d0: 20 78 5b 69 32 5d 29 3b 0d 0a 09 09 09 69 6e 74   x[i2]);.....int
02e0: 20 6d 79 20 3d 20 6d 69 6e 28 79 5b 6b 31 5d 2c   my = min(y[k1],
02f0: 20 79 5b 6b 32 5d 29 3b 0d 0a 09 09 09 69 6e 74   y[k2]);.....int
0300: 20 4d 79 20 3d 20 6d 61 78 28 79 5b 6b 31 5d 2c   My = max(y[k1],
0310: 20 79 5b 6b 32 5d 29 3b 0d 0a 0d 0a 09 09 09 69   y[k2]);.......i
0320: 6e 74 20 68 68 5b 34 30 5d 2c 20 71 3d 30 3b 0d  nt hh[40], q=0;.
0330: 0a 09 09 09 69 6e 74 20 63 75 74 20 3d 20 30 2c  ....int cut = 0,
0340: 20 6e 65 65 64 20 3d 20 28 28 4d 78 2d 6d 78 29   need = ((Mx-mx)
0350: 20 2b 20 28 4d 79 2d 6d 79 29 29 2a 32 3b 0d 0a   + (My-my))*2;..
0360: 09 09 09 66 6f 72 28 69 6e 74 20 6a 3d 30 3b 20  ...for(int j=0; 
0370: 6a 3c 4e 3b 20 2b 2b 6a 29 0d 0a 09 09 09 09 69  j<N; ++j)......i
0380: 66 28 20 6d 78 3c 3d 78 5b 6a 5d 20 26 26 20 78  f( mx<=x[j] && x
0390: 5b 6a 5d 3c 3d 4d 78 20 26 26 20 6d 79 3c 3d 79  [j]<=Mx && my<=y
03a0: 5b 6a 5d 20 26 26 20 79 5b 6a 5d 3c 3d 4d 79 20  [j] && y[j]<=My 
03b0: 29 0d 0a 09 09 09 09 09 68 68 5b 71 2b 2b 5d 20  ).......hh[q++] 
03c0: 3d 20 68 5b 6a 5d 3b 0d 0a 09 09 09 09 65 6c 73  = h[j];......els
03d0: 65 20 7b 0d 0a 09 09 09 09 09 63 75 74 2b 2b 2c  e {.......cut++,
03e0: 20 6e 65 65 64 2d 3d 68 5b 6a 5d 3b 0d 0a 09 09   need-=h[j];....
03f0: 09 09 09 69 66 28 63 75 74 3e 3d 61 6e 73 29 20  ...if(cut>=ans) 
0400: 67 6f 74 6f 20 6e 65 78 74 3b 0d 0a 09 09 09 09  goto next;......
0410: 7d 0d 0a 09 09 09 69 66 28 20 6e 65 65 64 20 3e  }.....if( need >
0420: 20 30 20 29 0d 0a 09 09 09 7b 0d 0a 09 09 09 09   0 ).....{......
0430: 73 6f 72 74 28 20 26 68 68 5b 30 5d 2c 20 26 68  sort( &hh[0], &h
0440: 68 5b 71 5d 2c 20 67 72 65 61 74 65 72 3c 69 6e  h[q], greater<in
0450: 74 3e 28 29 20 29 3b 0d 0a 2f 2f 09 09 09 09 70  t>() );..//....p
0460: 61 72 74 69 61 6c 5f 73 75 6d 28 20 26 68 68 5b  artial_sum( &hh[
0470: 30 5d 2c 20 26 68 68 5b 71 5d 2c 20 26 68 68 5b  0], &hh[q], &hh[
0480: 30 5d 20 29 3b 0d 0a 2f 2f 09 09 09 09 63 75 74  0] );..//....cut
0490: 20 2b 3d 20 6c 6f 77 65 72 5f 62 6f 75 6e 64 28   += lower_bound(
04a0: 20 26 68 68 5b 30 5d 2c 20 26 68 68 5b 71 5d 2c   &hh[0], &hh[q],
04b0: 20 6e 65 65 64 20 29 20 2d 20 26 68 68 5b 30 5d   need ) - &hh[0]
04c0: 3b 0d 0a 09 09 09 09 66 6f 72 28 69 6e 74 20 70  ;......for(int p
04d0: 3d 30 3b 20 70 3c 71 3b 20 2b 2b 70 29 20 7b 0d  =0; p<q; ++p) {.
04e0: 0a 09 09 09 09 09 2b 2b 63 75 74 2c 20 6e 65 65  ......++cut, nee
04f0: 64 20 2d 3d 20 68 68 5b 70 5d 3b 0d 0a 09 09 09  d -= hh[p];.....
0500: 09 09 69 66 28 63 75 74 3e 3d 61 6e 73 29 20 67  ..if(cut>=ans) g
0510: 6f 74 6f 20 6e 65 78 74 3b 0d 0a 09 09 09 09 09  oto next;.......
0520: 69 66 28 6e 65 65 64 20 3c 3d 20 30 29 20 62 72  if(need <= 0) br
0530: 65 61 6b 3b 0d 0a 09 09 09 09 7d 0d 0a 09 09 09  eak;......}.....
0540: 7d 0d 0a 09 09 09 61 6e 73 20 3d 20 63 75 74 3b  }.....ans = cut;
0550: 0d 0a 09 09 6e 65 78 74 3a 3b 0d 0a 09 09 7d 0d  ....next:;....}.
0560: 0a 09 09 72 65 74 75 72 6e 20 61 6e 73 3b 0d 0a  ...return ans;..
0570: 09 7d 0d 0a 7d 3b 0d 0a 0d 0a 2f 2f 20 42 45 47  .}..};....// BEG
0580: 49 4e 20 43 55 54 20 48 45 52 45 0d 0a 23 69 6e  IN CUT HERE..#in
0590: 63 6c 75 64 65 20 3c 63 74 69 6d 65 3e 0d 0a 64  clude <ctime>..d
05a0: 6f 75 62 6c 65 20 73 74 61 72 74 5f 74 69 6d 65  ouble start_time
05b0: 3b 73 74 72 69 6e 67 20 74 69 6d 65 72 28 29 20  ;string timer() 
05c0: 7b 20 6f 73 74 72 69 6e 67 73 74 72 65 61 6d 20  { ostringstream 
05d0: 6f 73 3b 20 6f 73 20 3c 3c 20 22 20 28 22 20 3c  os; os << " (" <
05e0: 3c 20 69 6e 74 28 28 63 6c 6f 63 6b 28 29 2d 73  < int((clock()-s
05f0: 74 61 72 74 5f 74 69 6d 65 29 2f 43 4c 4f 43 4b  tart_time)/CLOCK
0600: 53 5f 50 45 52 5f 53 45 43 2a 31 30 30 30 29 20  S_PER_SEC*1000) 
0610: 3c 3c 20 22 20 6d 73 65 63 29 22 3b 20 72 65 74  << " msec)"; ret
0620: 75 72 6e 20 6f 73 2e 73 74 72 28 29 3b 20 7d 0d  urn os.str(); }.
0630: 0a 0d 0a 74 65 6d 70 6c 61 74 65 20 3c 74 79 70  ...template <typ
0640: 65 6e 61 6d 65 20 54 3e 20 73 74 72 69 6e 67 20  ename T> string 
0650: 70 72 69 6e 74 5f 61 72 72 61 79 28 63 6f 6e 73  print_array(cons
0660: 74 20 76 65 63 74 6f 72 3c 54 3e 20 26 56 29 20  t vector<T> &V) 
0670: 7b 20 6f 73 74 72 69 6e 67 73 74 72 65 61 6d 20  { ostringstream 
0680: 6f 73 3b 20 6f 73 20 3c 3c 20 22 7b 20 22 3b 20  os; os << "{ "; 
0690: 66 6f 72 20 28 74 79 70 65 6e 61 6d 65 20 76 65  for (typename ve
06a0: 63 74 6f 72 3c 54 3e 3a 3a 63 6f 6e 73 74 5f 69  ctor<T>::const_i
06b0: 74 65 72 61 74 6f 72 20 69 74 65 72 20 3d 20 56  terator iter = V
06c0: 2e 62 65 67 69 6e 28 29 3b 20 69 74 65 72 20 21  .begin(); iter !
06d0: 3d 20 56 2e 65 6e 64 28 29 3b 20 2b 2b 69 74 65  = V.end(); ++ite
06e0: 72 29 20 6f 73 20 3c 3c 20 27 5c 22 27 20 3c 3c  r) os << '\"' <<
06f0: 20 2a 69 74 65 72 20 3c 3c 20 22 5c 22 2c 22 3b   *iter << "\",";
0700: 20 6f 73 20 3c 3c 20 22 20 7d 22 3b 20 72 65 74   os << " }"; ret
0710: 75 72 6e 20 6f 73 2e 73 74 72 28 29 3b 20 7d 0d  urn os.str(); }.
0720: 0a 69 6e 74 20 76 65 72 69 66 79 5f 63 61 73 65  .int verify_case
0730: 28 63 6f 6e 73 74 20 69 6e 74 20 26 45 78 70 65  (const int &Expe
0740: 63 74 65 64 2c 20 63 6f 6e 73 74 20 69 6e 74 20  cted, const int 
0750: 26 52 65 63 65 69 76 65 64 29 20 7b 20 69 66 20  &Received) { if 
0760: 28 45 78 70 65 63 74 65 64 20 3d 3d 20 52 65 63  (Expected == Rec
0770: 65 69 76 65 64 29 20 63 65 72 72 20 3c 3c 20 22  eived) cerr << "
0780: 50 41 53 53 45 44 22 20 3c 3c 20 74 69 6d 65 72  PASSED" << timer
0790: 28 29 20 3c 3c 20 65 6e 64 6c 3b 20 65 6c 73 65  () << endl; else
07a0: 20 7b 20 63 65 72 72 20 3c 3c 20 22 46 41 49 4c   { cerr << "FAIL
07b0: 45 44 22 20 3c 3c 20 74 69 6d 65 72 28 29 20 3c  ED" << timer() <
07c0: 3c 20 65 6e 64 6c 3b 20 63 65 72 72 20 3c 3c 20  < endl; cerr << 
07d0: 22 5c 74 45 78 70 65 63 74 65 64 3a 20 5c 22 22  "\tExpected: \""
07e0: 20 3c 3c 20 45 78 70 65 63 74 65 64 20 3c 3c 20   << Expected << 
07f0: 27 5c 22 27 20 3c 3c 20 65 6e 64 6c 3b 20 63 65  '\"' << endl; ce
0800: 72 72 20 3c 3c 20 22 5c 74 52 65 63 65 69 76 65  rr << "\tReceive
0810: 64 3a 20 5c 22 22 20 3c 3c 20 52 65 63 65 69 76  d: \"" << Receiv
0820: 65 64 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 65 6e  ed << '\"' << en
0830: 64 6c 3b 20 7d 20 72 65 74 75 72 6e 20 30 3b 7d  dl; } return 0;}
0840: 0d 0a 0d 0a 74 65 6d 70 6c 61 74 65 3c 69 6e 74  ....template<int
0850: 20 4e 3e 20 73 74 72 75 63 74 20 43 61 73 65 5f   N> struct Case_
0860: 20 7b 20 43 61 73 65 5f 28 29 7b 73 74 61 72 74   { Case_(){start
0870: 5f 74 69 6d 65 3d 63 6c 6f 63 6b 28 29 3b 7d 20  _time=clock();} 
0880: 7d 3b 0d 0a 63 68 61 72 20 54 65 73 74 5f 28 2e  };..char Test_(.
0890: 2e 2e 29 3b 0d 0a 69 6e 74 20 54 65 73 74 5f 28  ..);..int Test_(
08a0: 43 61 73 65 5f 3c 30 3e 29 20 7b 0d 0a 09 69 6e  Case_<0>) {...in
08b0: 74 20 78 5f 5b 5d 20 3d 20 7b 31 2c 20 32 2c 20  t x_[] = {1, 2, 
08c0: 38 2c 20 39 2c 20 35 7d 3b 0d 0a 09 20 20 76 65  8, 9, 5};...  ve
08d0: 63 74 6f 72 20 3c 69 6e 74 3e 20 78 28 78 5f 2c  ctor <int> x(x_,
08e0: 20 78 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 2f 73   x_+sizeof(x_)/s
08f0: 69 7a 65 6f 66 28 2a 78 5f 29 29 3b 20 0d 0a 09  izeof(*x_)); ...
0900: 69 6e 74 20 79 5f 5b 5d 20 3d 20 7b 31 2c 20 38  int y_[] = {1, 8
0910: 2c 20 32 2c 20 39 2c 20 35 7d 3b 0d 0a 09 20 20  , 2, 9, 5};...  
0920: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79 28 79  vector <int> y(y
0930: 5f 2c 20 79 5f 2b 73 69 7a 65 6f 66 28 79 5f 29  _, y_+sizeof(y_)
0940: 2f 73 69 7a 65 6f 66 28 2a 79 5f 29 29 3b 20 0d  /sizeof(*y_)); .
0950: 0a 09 69 6e 74 20 68 5f 5b 5d 20 3d 20 7b 31 2c  ..int h_[] = {1,
0960: 20 31 2c 20 31 2c 20 31 2c 20 33 32 7d 3b 0d 0a   1, 1, 1, 32};..
0970: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20  .  vector <int> 
0980: 68 28 68 5f 2c 20 68 5f 2b 73 69 7a 65 6f 66 28  h(h_, h_+sizeof(
0990: 68 5f 29 2f 73 69 7a 65 6f 66 28 2a 68 5f 29 29  h_)/sizeof(*h_))
09a0: 3b 20 0d 0a 09 69 6e 74 20 52 65 74 56 61 6c 20  ; ...int RetVal 
09b0: 3d 20 31 3b 20 0d 0a 09 72 65 74 75 72 6e 20 76  = 1; ...return v
09c0: 65 72 69 66 79 5f 63 61 73 65 28 52 65 74 56 61  erify_case(RetVa
09d0: 6c 2c 20 53 61 76 65 54 68 65 54 72 65 65 73 28  l, SaveTheTrees(
09e0: 29 2e 6d 69 6e 69 6d 75 6d 43 75 74 28 78 2c 20  ).minimumCut(x, 
09f0: 79 2c 20 68 29 29 3b 20 7d 0d 0a 69 6e 74 20 54  y, h)); }..int T
0a00: 65 73 74 5f 28 43 61 73 65 5f 3c 31 3e 29 20 7b  est_(Case_<1>) {
0a10: 0d 0a 09 69 6e 74 20 78 5f 5b 5d 20 3d 20 7b 31  ...int x_[] = {1
0a20: 2c 20 32 2c 20 38 2c 20 39 2c 20 35 7d 3b 0d 0a  , 2, 8, 9, 5};..
0a30: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20  .  vector <int> 
0a40: 78 28 78 5f 2c 20 78 5f 2b 73 69 7a 65 6f 66 28  x(x_, x_+sizeof(
0a50: 78 5f 29 2f 73 69 7a 65 6f 66 28 2a 78 5f 29 29  x_)/sizeof(*x_))
0a60: 3b 20 0d 0a 09 69 6e 74 20 79 5f 5b 5d 20 3d 20  ; ...int y_[] = 
0a70: 7b 31 2c 20 38 2c 20 32 2c 20 39 2c 20 35 7d 3b  {1, 8, 2, 9, 5};
0a80: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74  ...  vector <int
0a90: 3e 20 79 28 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f  > y(y_, y_+sizeo
0aa0: 66 28 79 5f 29 2f 73 69 7a 65 6f 66 28 2a 79 5f  f(y_)/sizeof(*y_
0ab0: 29 29 3b 20 0d 0a 09 69 6e 74 20 68 5f 5b 5d 20  )); ...int h_[] 
0ac0: 3d 20 7b 31 2c 20 31 2c 20 31 2c 20 31 2c 20 32  = {1, 1, 1, 1, 2
0ad0: 37 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c  7};...  vector <
0ae0: 69 6e 74 3e 20 68 28 68 5f 2c 20 68 5f 2b 73 69  int> h(h_, h_+si
0af0: 7a 65 6f 66 28 68 5f 29 2f 73 69 7a 65 6f 66 28  zeof(h_)/sizeof(
0b00: 2a 68 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 52 65  *h_)); ...int Re
0b10: 74 56 61 6c 20 3d 20 32 3b 20 0d 0a 09 72 65 74  tVal = 2; ...ret
0b20: 75 72 6e 20 76 65 72 69 66 79 5f 63 61 73 65 28  urn verify_case(
0b30: 52 65 74 56 61 6c 2c 20 53 61 76 65 54 68 65 54  RetVal, SaveTheT
0b40: 72 65 65 73 28 29 2e 6d 69 6e 69 6d 75 6d 43 75  rees().minimumCu
0b50: 74 28 78 2c 20 79 2c 20 68 29 29 3b 20 7d 0d 0a  t(x, y, h)); }..
0b60: 69 6e 74 20 54 65 73 74 5f 28 43 61 73 65 5f 3c  int Test_(Case_<
0b70: 32 3e 29 20 7b 0d 0a 09 69 6e 74 20 78 5f 5b 5d  2>) {...int x_[]
0b80: 20 3d 20 7b 31 2c 20 32 2c 20 38 2c 20 39 2c 20   = {1, 2, 8, 9, 
0b90: 35 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c  5};...  vector <
0ba0: 69 6e 74 3e 20 78 28 78 5f 2c 20 78 5f 2b 73 69  int> x(x_, x_+si
0bb0: 7a 65 6f 66 28 78 5f 29 2f 73 69 7a 65 6f 66 28  zeof(x_)/sizeof(
0bc0: 2a 78 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 79 5f  *x_)); ...int y_
0bd0: 5b 5d 20 3d 20 7b 31 2c 20 38 2c 20 32 2c 20 39  [] = {1, 8, 2, 9
0be0: 2c 20 35 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72  , 5};...  vector
0bf0: 20 3c 69 6e 74 3e 20 79 28 79 5f 2c 20 79 5f 2b   <int> y(y_, y_+
0c00: 73 69 7a 65 6f 66 28 79 5f 29 2f 73 69 7a 65 6f  sizeof(y_)/sizeo
0c10: 66 28 2a 79 5f 29 29 3b 20 0d 0a 09 69 6e 74 20  f(*y_)); ...int 
0c20: 68 5f 5b 5d 20 3d 20 7b 34 2c 20 34 2c 20 34 2c  h_[] = {4, 4, 4,
0c30: 20 34 2c 20 34 7d 3b 0d 0a 09 20 20 76 65 63 74   4, 4};...  vect
0c40: 6f 72 20 3c 69 6e 74 3e 20 68 28 68 5f 2c 20 68  or <int> h(h_, h
0c50: 5f 2b 73 69 7a 65 6f 66 28 68 5f 29 2f 73 69 7a  _+sizeof(h_)/siz
0c60: 65 6f 66 28 2a 68 5f 29 29 3b 20 0d 0a 09 69 6e  eof(*h_)); ...in
0c70: 74 20 52 65 74 56 61 6c 20 3d 20 33 3b 20 0d 0a  t RetVal = 3; ..
0c80: 09 72 65 74 75 72 6e 20 76 65 72 69 66 79 5f 63  .return verify_c
0c90: 61 73 65 28 52 65 74 56 61 6c 2c 20 53 61 76 65  ase(RetVal, Save
0ca0: 54 68 65 54 72 65 65 73 28 29 2e 6d 69 6e 69 6d  TheTrees().minim
0cb0: 75 6d 43 75 74 28 78 2c 20 79 2c 20 68 29 29 3b  umCut(x, y, h));
0cc0: 20 7d 0d 0a 69 6e 74 20 54 65 73 74 5f 28 43 61   }..int Test_(Ca
0cd0: 73 65 5f 3c 33 3e 29 20 7b 0d 0a 09 69 6e 74 20  se_<3>) {...int 
0ce0: 78 5f 5b 5d 20 3d 20 7b 31 2c 20 32 2c 20 38 2c  x_[] = {1, 2, 8,
0cf0: 20 39 2c 20 35 7d 3b 0d 0a 09 20 20 76 65 63 74   9, 5};...  vect
0d00: 6f 72 20 3c 69 6e 74 3e 20 78 28 78 5f 2c 20 78  or <int> x(x_, x
0d10: 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 2f 73 69 7a  _+sizeof(x_)/siz
0d20: 65 6f 66 28 2a 78 5f 29 29 3b 20 0d 0a 09 69 6e  eof(*x_)); ...in
0d30: 74 20 79 5f 5b 5d 20 3d 20 7b 31 2c 20 38 2c 20  t y_[] = {1, 8, 
0d40: 32 2c 20 39 2c 20 35 7d 3b 0d 0a 09 20 20 76 65  2, 9, 5};...  ve
0d50: 63 74 6f 72 20 3c 69 6e 74 3e 20 79 28 79 5f 2c  ctor <int> y(y_,
0d60: 20 79 5f 2b 73 69 7a 65 6f 66 28 79 5f 29 2f 73   y_+sizeof(y_)/s
0d70: 69 7a 65 6f 66 28 2a 79 5f 29 29 3b 20 0d 0a 09  izeof(*y_)); ...
0d80: 69 6e 74 20 68 5f 5b 5d 20 3d 20 7b 33 2c 20 33  int h_[] = {3, 3
0d90: 2c 20 33 2c 20 33 2c 20 33 7d 3b 0d 0a 09 20 20  , 3, 3, 3};...  
0da0: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 68 28 68  vector <int> h(h
0db0: 5f 2c 20 68 5f 2b 73 69 7a 65 6f 66 28 68 5f 29  _, h_+sizeof(h_)
0dc0: 2f 73 69 7a 65 6f 66 28 2a 68 5f 29 29 3b 20 0d  /sizeof(*h_)); .
0dd0: 0a 09 69 6e 74 20 52 65 74 56 61 6c 20 3d 20 34  ..int RetVal = 4
0de0: 3b 20 0d 0a 09 72 65 74 75 72 6e 20 76 65 72 69  ; ...return veri
0df0: 66 79 5f 63 61 73 65 28 52 65 74 56 61 6c 2c 20  fy_case(RetVal, 
0e00: 53 61 76 65 54 68 65 54 72 65 65 73 28 29 2e 6d  SaveTheTrees().m
0e10: 69 6e 69 6d 75 6d 43 75 74 28 78 2c 20 79 2c 20  inimumCut(x, y, 
0e20: 68 29 29 3b 20 7d 0d 0a 69 6e 74 20 54 65 73 74  h)); }..int Test
0e30: 5f 28 43 61 73 65 5f 3c 34 3e 29 20 7b 0d 0a 09  _(Case_<4>) {...
0e40: 69 6e 74 20 78 5f 5b 5d 20 3d 20 7b 33 2c 20 31  int x_[] = {3, 1
0e50: 30 2c 20 36 2c 20 32 2c 20 38 2c 20 37 2c 20 35  0, 6, 2, 8, 7, 5
0e60: 2c 20 34 2c 20 39 2c 20 31 7d 3b 0d 0a 09 20 20  , 4, 9, 1};...  
0e70: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 78 28 78  vector <int> x(x
0e80: 5f 2c 20 78 5f 2b 73 69 7a 65 6f 66 28 78 5f 29  _, x_+sizeof(x_)
0e90: 2f 73 69 7a 65 6f 66 28 2a 78 5f 29 29 3b 20 0d  /sizeof(*x_)); .
0ea0: 0a 09 69 6e 74 20 79 5f 5b 5d 20 3d 20 7b 34 2c  ..int y_[] = {4,
0eb0: 20 31 30 2c 20 36 2c 20 39 2c 20 38 2c 20 37 2c   10, 6, 9, 8, 7,
0ec0: 20 33 2c 20 32 2c 20 35 2c 20 31 7d 3b 0d 0a 09   3, 2, 5, 1};...
0ed0: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79    vector <int> y
0ee0: 28 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f 66 28 79  (y_, y_+sizeof(y
0ef0: 5f 29 2f 73 69 7a 65 6f 66 28 2a 79 5f 29 29 3b  _)/sizeof(*y_));
0f00: 20 0d 0a 09 69 6e 74 20 68 5f 5b 5d 20 3d 20 7b   ...int h_[] = {
0f10: 34 2c 20 32 2c 20 32 2c 20 34 2c 20 34 2c 20 35  4, 2, 2, 4, 4, 5
0f20: 2c 20 34 2c 20 34 2c 20 32 2c 20 33 7d 3b 0d 0a  , 4, 4, 2, 3};..
0f30: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20  .  vector <int> 
0f40: 68 28 68 5f 2c 20 68 5f 2b 73 69 7a 65 6f 66 28  h(h_, h_+sizeof(
0f50: 68 5f 29 2f 73 69 7a 65 6f 66 28 2a 68 5f 29 29  h_)/sizeof(*h_))
0f60: 3b 20 0d 0a 09 69 6e 74 20 52 65 74 56 61 6c 20  ; ...int RetVal 
0f70: 3d 20 35 3b 20 0d 0a 09 72 65 74 75 72 6e 20 76  = 5; ...return v
0f80: 65 72 69 66 79 5f 63 61 73 65 28 52 65 74 56 61  erify_case(RetVa
0f90: 6c 2c 20 53 61 76 65 54 68 65 54 72 65 65 73 28  l, SaveTheTrees(
0fa0: 29 2e 6d 69 6e 69 6d 75 6d 43 75 74 28 78 2c 20  ).minimumCut(x, 
0fb0: 79 2c 20 68 29 29 3b 20 7d 0d 0a 69 6e 74 20 54  y, h)); }..int T
0fc0: 65 73 74 5f 28 43 61 73 65 5f 3c 35 3e 29 20 7b  est_(Case_<5>) {
0fd0: 0d 0a 09 69 6e 74 20 78 5f 5b 5d 20 3d 20 7b 35  ...int x_[] = {5
0fe0: 35 2c 20 36 37 2c 20 31 30 30 2c 20 33 38 2c 20  5, 67, 100, 38, 
0ff0: 38 30 2c 20 39 38 2c 20 34 37 2c 20 35 38 2c 20  80, 98, 47, 58, 
1000: 36 31 2c 20 33 33 7d 3b 0d 0a 09 20 20 76 65 63  61, 33};...  vec
1010: 74 6f 72 20 3c 69 6e 74 3e 20 78 28 78 5f 2c 20  tor <int> x(x_, 
1020: 78 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 2f 73 69  x_+sizeof(x_)/si
1030: 7a 65 6f 66 28 2a 78 5f 29 29 3b 20 0d 0a 09 69  zeof(*x_)); ...i
1040: 6e 74 20 79 5f 5b 5d 20 3d 20 7b 38 37 2c 20 31  nt y_[] = {87, 1
1050: 37 2c 20 32 30 2c 20 37 2c 20 35 37 2c 20 31 39  7, 20, 7, 57, 19
1060: 2c 20 32 33 2c 20 36 38 2c 20 32 37 2c 20 33 39  , 23, 68, 27, 39
1070: 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69  };...  vector <i
1080: 6e 74 3e 20 79 28 79 5f 2c 20 79 5f 2b 73 69 7a  nt> y(y_, y_+siz
1090: 65 6f 66 28 79 5f 29 2f 73 69 7a 65 6f 66 28 2a  eof(y_)/sizeof(*
10a0: 79 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 68 5f 5b  y_)); ...int h_[
10b0: 5d 20 3d 20 7b 36 2c 20 32 2c 20 38 2c 20 38 2c  ] = {6, 2, 8, 8,
10c0: 20 33 31 2c 20 32 35 2c 20 32 33 2c 20 31 39 2c   31, 25, 23, 19,
10d0: 20 34 35 2c 20 34 7d 3b 0d 0a 09 20 20 76 65 63   45, 4};...  vec
10e0: 74 6f 72 20 3c 69 6e 74 3e 20 68 28 68 5f 2c 20  tor <int> h(h_, 
10f0: 68 5f 2b 73 69 7a 65 6f 66 28 68 5f 29 2f 73 69  h_+sizeof(h_)/si
1100: 7a 65 6f 66 28 2a 68 5f 29 29 3b 20 0d 0a 09 69  zeof(*h_)); ...i
1110: 6e 74 20 52 65 74 56 61 6c 20 3d 20 36 3b 20 0d  nt RetVal = 6; .
1120: 0a 09 72 65 74 75 72 6e 20 76 65 72 69 66 79 5f  ..return verify_
1130: 63 61 73 65 28 52 65 74 56 61 6c 2c 20 53 61 76  case(RetVal, Sav
1140: 65 54 68 65 54 72 65 65 73 28 29 2e 6d 69 6e 69  eTheTrees().mini
1150: 6d 75 6d 43 75 74 28 78 2c 20 79 2c 20 68 29 29  mumCut(x, y, h))
1160: 3b 20 7d 0d 0a 69 6e 74 20 54 65 73 74 5f 28 43  ; }..int Test_(C
1170: 61 73 65 5f 3c 36 3e 29 20 7b 0d 0a 09 69 6e 74  ase_<6>) {...int
1180: 20 78 5f 5b 5d 20 3d 20 7b 31 2c 32 2c 33 2c 34   x_[] = {1,2,3,4
1190: 2c 35 2c 36 2c 37 2c 38 2c 39 2c 31 30 2c 31 31  ,5,6,7,8,9,10,11
11a0: 2c 31 32 2c 31 33 2c 31 34 2c 31 35 2c 31 36 2c  ,12,13,14,15,16,
11b0: 31 37 2c 31 38 2c 31 39 2c 32 30 2c 32 31 2c 32  17,18,19,20,21,2
11c0: 32 2c 32 33 2c 32 34 2c 32 35 2c 32 36 2c 32 37  2,23,24,25,26,27
11d0: 2c 32 38 2c 32 39 2c 33 30 2c 33 31 2c 33 32 2c  ,28,29,30,31,32,
11e0: 33 33 2c 33 34 2c 33 35 2c 33 36 2c 33 37 2c 33  33,34,35,36,37,3
11f0: 38 2c 33 39 2c 34 30 7d 3b 0d 0a 09 20 20 76 65  8,39,40};...  ve
1200: 63 74 6f 72 20 3c 69 6e 74 3e 20 78 28 78 5f 2c  ctor <int> x(x_,
1210: 20 78 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 2f 73   x_+sizeof(x_)/s
1220: 69 7a 65 6f 66 28 2a 78 5f 29 29 3b 20 0d 0a 09  izeof(*x_)); ...
1230: 69 6e 74 20 79 5f 5b 5d 20 3d 20 7b 31 2c 32 2c  int y_[] = {1,2,
1240: 33 2c 34 2c 35 2c 36 2c 37 2c 38 2c 39 2c 31 30  3,4,5,6,7,8,9,10
1250: 2c 31 31 2c 31 32 2c 31 33 2c 31 34 2c 31 35 2c  ,11,12,13,14,15,
1260: 31 36 2c 31 37 2c 31 38 2c 31 39 2c 32 30 2c 32  16,17,18,19,20,2
1270: 31 2c 32 32 2c 32 33 2c 32 34 2c 32 35 2c 32 36  1,22,23,24,25,26
1280: 2c 32 37 2c 32 38 2c 32 39 2c 33 30 2c 33 31 2c  ,27,28,29,30,31,
1290: 33 32 2c 33 33 2c 33 34 2c 33 35 2c 33 36 2c 33  32,33,34,35,36,3
12a0: 37 2c 33 38 2c 33 39 2c 34 30 7d 3b 0d 0a 09 20  7,38,39,40};... 
12b0: 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79 28   vector <int> y(
12c0: 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f 66 28 79 5f  y_, y_+sizeof(y_
12d0: 29 2f 73 69 7a 65 6f 66 28 2a 79 5f 29 29 3b 20  )/sizeof(*y_)); 
12e0: 0d 0a 09 69 6e 74 20 68 5f 5b 5d 20 3d 20 7b 31  ...int h_[] = {1
12f0: 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31  ,1,1,1,1,1,1,1,1
1300: 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31  ,1,1,1,1,1,1,1,1
1310: 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31  ,1,1,1,1,1,1,1,1
1320: 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31  ,1,1,1,1,1,1,1,1
1330: 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31 7d 3b  ,1,1,1,1,1,1,1};
1340: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74  ...  vector <int
1350: 3e 20 68 28 68 5f 2c 20 68 5f 2b 73 69 7a 65 6f  > h(h_, h_+sizeo
1360: 66 28 68 5f 29 2f 73 69 7a 65 6f 66 28 2a 68 5f  f(h_)/sizeof(*h_
1370: 29 29 3b 20 0d 0a 09 69 6e 74 20 52 65 74 56 61  )); ...int RetVa
1380: 6c 20 3d 20 36 3b 20 0d 0a 09 72 65 74 75 72 6e  l = 6; ...return
1390: 20 76 65 72 69 66 79 5f 63 61 73 65 28 52 65 74   verify_case(Ret
13a0: 56 61 6c 2c 20 53 61 76 65 54 68 65 54 72 65 65  Val, SaveTheTree
13b0: 73 28 29 2e 6d 69 6e 69 6d 75 6d 43 75 74 28 78  s().minimumCut(x
13c0: 2c 20 79 2c 20 68 29 29 3b 20 7d 0d 0a 0d 0a 74  , y, h)); }....t
13d0: 65 6d 70 6c 61 74 65 3c 69 6e 74 20 4e 3e 20 76  emplate<int N> v
13e0: 6f 69 64 20 52 75 6e 5f 28 29 20 7b 20 63 65 72  oid Run_() { cer
13f0: 72 20 3c 3c 20 22 54 65 73 74 20 43 61 73 65 20  r << "Test Case 
1400: 23 22 20 3c 3c 20 4e 20 3c 3c 20 22 2e 2e 2e 22  #" << N << "..."
1410: 20 3c 3c 20 66 6c 75 73 68 3b 20 54 65 73 74 5f   << flush; Test_
1420: 28 43 61 73 65 5f 3c 4e 3e 28 29 29 3b 20 52 75  (Case_<N>()); Ru
1430: 6e 5f 3c 73 69 7a 65 6f 66 28 54 65 73 74 5f 28  n_<sizeof(Test_(
1440: 43 61 73 65 5f 3c 4e 2b 31 3e 28 29 29 29 3d 3d  Case_<N+1>()))==
1450: 31 20 3f 20 2d 31 20 3a 20 4e 2b 31 3e 28 29 3b  1 ? -1 : N+1>();
1460: 20 7d 0d 0a 74 65 6d 70 6c 61 74 65 3c 3e 20 20   }..template<>  
1470: 20 20 20 20 76 6f 69 64 20 52 75 6e 5f 3c 2d 31      void Run_<-1
1480: 3e 28 29 20 7b 7d 0d 0a 69 6e 74 20 6d 61 69 6e  >() {}..int main
1490: 28 29 20 7b 20 52 75 6e 5f 3c 30 3e 28 29 3b 20  () { Run_<0>(); 
14a0: 7d 0d 0a 2f 2f 20 45 4e 44 20 43 55 54 20 48 45  }..// END CUT HE
14b0: 52 45 0d 0a 0d 0a                                RE....