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 42 61 CMP;....class Ba
01a0: 72 62 61 72 69 61 6e 49 6e 76 61 73 69 6f 6e 32 rbarianInvasion2
01b0: 20 7b 20 70 75 62 6c 69 63 3a 0d 0a 09 64 6f 75 { public:...dou
01c0: 62 6c 65 20 6d 69 6e 69 6d 75 6d 54 69 6d 65 28 ble minimumTime(
01d0: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 62 6f 75 vector <int> bou
01e0: 6e 64 61 72 79 58 2c 20 76 65 63 74 6f 72 20 3c ndaryX, vector <
01f0: 69 6e 74 3e 20 62 6f 75 6e 64 61 72 79 59 2c 20 int> boundaryY,
0200: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 63 69 74 vector <int> cit
0210: 79 58 2c 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e yX, vector <int>
0220: 20 63 69 74 79 59 29 20 0d 0a 09 7b 0d 0a 09 09 cityY) ...{....
0230: 76 65 63 74 6f 72 3c 43 4d 50 3e 20 62 2c 20 63 vector<CMP> b, c
0240: 3b 0d 0a 09 09 66 6f 72 28 69 6e 74 20 69 3d 30 ;....for(int i=0
0250: 3b 20 69 3c 62 6f 75 6e 64 61 72 79 58 2e 73 69 ; i<boundaryX.si
0260: 7a 65 28 29 3b 20 2b 2b 69 29 0d 0a 09 09 09 62 ze(); ++i).....b
0270: 2e 70 75 73 68 5f 62 61 63 6b 28 43 4d 50 28 62 .push_back(CMP(b
0280: 6f 75 6e 64 61 72 79 58 5b 69 5d 2c 20 62 6f 75 oundaryX[i], bou
0290: 6e 64 61 72 79 59 5b 69 5d 29 29 3b 0d 0a 09 09 ndaryY[i]));....
02a0: 66 6f 72 28 69 6e 74 20 69 3d 30 3b 20 69 3c 63 for(int i=0; i<c
02b0: 69 74 79 58 2e 73 69 7a 65 28 29 3b 20 2b 2b 69 ityX.size(); ++i
02c0: 29 0d 0a 09 09 09 63 2e 70 75 73 68 5f 62 61 63 ).....c.push_bac
02d0: 6b 28 43 4d 50 28 63 69 74 79 58 5b 69 5d 2c 20 k(CMP(cityX[i],
02e0: 63 69 74 79 59 5b 69 5d 29 29 3b 0d 0a 0d 0a 09 cityY[i]));.....
02f0: 09 64 6f 75 62 6c 65 20 61 6e 73 77 65 72 20 3d .double answer =
0300: 20 30 3b 0d 0a 09 09 66 6f 72 28 69 6e 74 20 69 0;....for(int i
0310: 3d 30 3b 20 69 3c 62 2e 73 69 7a 65 28 29 3b 20 =0; i<b.size();
0320: 2b 2b 69 29 0d 0a 09 09 09 61 6e 73 77 65 72 20 ++i).....answer
0330: 3d 20 6d 61 78 28 61 6e 73 77 65 72 2c 20 6d 69 = max(answer, mi
0340: 6e 5f 6d 61 78 28 62 5b 69 5d 2c 20 62 5b 28 69 n_max(b[i], b[(i
0350: 2b 31 29 25 62 2e 73 69 7a 65 28 29 5d 2c 20 63 +1)%b.size()], c
0360: 29 29 3b 0d 0a 09 09 72 65 74 75 72 6e 20 61 6e ));....return an
0370: 73 77 65 72 3b 0d 0a 09 7d 0d 0a 0d 0a 09 64 6f swer;...}.....do
0380: 75 62 6c 65 20 6d 69 6e 5f 6d 61 78 28 43 4d 50 uble min_max(CMP
0390: 20 61 2c 20 43 4d 50 20 62 2c 20 76 65 63 74 6f a, CMP b, vecto
03a0: 72 3c 43 4d 50 3e 20 63 73 29 0d 0a 09 7b 0d 0a r<CMP> cs)...{..
03b0: 09 09 2f 2f 20 72 65 74 75 72 6e 20 6d 61 78 5b ..// return max[
03c0: 70 20 6f 6e 20 61 2d 62 5d 20 6d 69 6e 5b 63 20 p on a-b] min[c
03d0: 69 6e 20 63 73 5d 20 7c 63 2d 70 7c 0d 0a 0d 0a in cs] |c-p|....
03e0: 09 09 64 6f 75 62 6c 65 20 66 61 63 74 6f 72 20 ..double factor
03f0: 3d 20 61 62 73 28 62 2d 61 29 3b 0d 0a 09 09 7b = abs(b-a);....{
0400: 0d 0a 09 09 09 62 20 2d 3d 20 61 3b 0d 0a 09 09 .....b -= a;....
0410: 09 66 6f 72 28 69 6e 74 20 69 3d 30 3b 20 69 3c .for(int i=0; i<
0420: 63 73 2e 73 69 7a 65 28 29 3b 20 2b 2b 69 29 0d cs.size(); ++i).
0430: 0a 09 09 09 09 63 73 5b 69 5d 20 2d 3d 20 61 3b .....cs[i] -= a;
0440: 0d 0a 09 09 09 66 6f 72 28 69 6e 74 20 69 3d 30 .....for(int i=0
0450: 3b 20 69 3c 63 73 2e 73 69 7a 65 28 29 3b 20 2b ; i<cs.size(); +
0460: 2b 69 29 0d 0a 09 09 09 09 63 73 5b 69 5d 20 2f +i)......cs[i] /
0470: 3d 20 62 3b 0d 0a 09 09 7d 0d 0a 0d 0a 09 09 2f = b;....}....../
0480: 2f 20 72 65 74 75 72 6e 20 66 61 63 74 6f 72 20 / return factor
0490: 2a 20 6d 61 78 5b 70 20 6f 6e 20 30 2d 31 5d 20 * max[p on 0-1]
04a0: 6d 69 6e 5b 63 20 69 6e 20 63 73 5d 20 7c 63 2d min[c in cs] |c-
04b0: 70 7c 0d 0a 09 09 76 65 63 74 6f 72 3c 64 6f 75 p|....vector<dou
04c0: 62 6c 65 3e 20 73 73 3b 0d 0a 09 09 66 6f 72 28 ble> ss;....for(
04d0: 69 6e 74 20 69 3d 30 3b 20 69 3c 63 73 2e 73 69 int i=0; i<cs.si
04e0: 7a 65 28 29 3b 20 2b 2b 69 29 0d 0a 09 09 09 66 ze(); ++i).....f
04f0: 6f 72 28 69 6e 74 20 6a 3d 69 2b 31 3b 20 6a 3c or(int j=i+1; j<
0500: 63 73 2e 73 69 7a 65 28 29 3b 20 2b 2b 6a 29 0d cs.size(); ++j).
0510: 0a 09 09 09 09 69 66 28 20 63 73 5b 69 5d 2e 72 .....if( cs[i].r
0520: 65 61 6c 28 29 20 21 3d 20 63 73 5b 6a 5d 2e 72 eal() != cs[j].r
0530: 65 61 6c 28 29 20 29 0d 0a 09 09 09 09 09 73 73 eal() ).......ss
0540: 2e 70 75 73 68 5f 62 61 63 6b 28 20 73 65 63 74 .push_back( sect
0550: 70 6f 69 6e 74 28 63 73 5b 69 5d 2c 20 63 73 5b point(cs[i], cs[
0560: 6a 5d 29 20 29 3b 0d 0a 09 09 73 6f 72 74 28 73 j]) );....sort(s
0570: 73 2e 62 65 67 69 6e 28 29 2c 20 73 73 2e 65 6e s.begin(), ss.en
0580: 64 28 29 29 3b 0d 0a 09 09 64 6f 75 62 6c 65 20 d());....double
0590: 66 6d 69 6e 5f 6d 61 78 20 3d 20 30 2e 30 3b 0d fmin_max = 0.0;.
05a0: 0a 0d 0a 09 09 64 6f 75 62 6c 65 20 7a 20 3d 20 .....double z =
05b0: 30 2e 30 3b 0d 0a 09 09 66 6f 72 28 69 6e 74 20 0.0;....for(int
05c0: 69 3d 30 3b 20 69 3c 73 73 2e 73 69 7a 65 28 29 i=0; i<ss.size()
05d0: 3b 20 2b 2b 69 29 20 7b 0d 0a 09 09 09 64 6f 75 ; ++i) {.....dou
05e0: 62 6c 65 20 74 20 3d 20 6d 69 6e 28 31 2e 30 2c ble t = min(1.0,
05f0: 20 73 73 5b 69 5d 29 3b 0d 0a 09 09 09 69 66 28 ss[i]);.....if(
0600: 20 74 20 3e 20 7a 20 29 20 7b 0d 0a 09 09 09 09 t > z ) {......
0610: 2f 2f 20 5b 7a 2c 20 74 29 0d 0a 09 09 09 09 7b // [z, t)......{
0620: 0d 0a 09 09 09 09 09 64 6f 75 62 6c 65 20 70 74 .......double pt
0630: 20 3d 20 28 7a 2b 74 29 2f 32 3b 0d 0a 09 09 09 = (z+t)/2;.....
0640: 09 09 69 6e 74 20 6d 69 6e 4b 20 3d 20 30 3b 0d ..int minK = 0;.
0650: 0a 09 09 09 09 09 66 6f 72 28 69 6e 74 20 6b 3d ......for(int k=
0660: 31 3b 20 6b 3c 63 73 2e 73 69 7a 65 28 29 3b 20 1; k<cs.size();
0670: 2b 2b 6b 29 0d 0a 09 09 09 09 09 09 69 66 28 20 ++k)........if(
0680: 61 62 73 28 63 73 5b 6d 69 6e 4b 5d 2d 70 74 29 abs(cs[minK]-pt)
0690: 20 3e 20 61 62 73 28 63 73 5b 6b 5d 2d 70 74 29 > abs(cs[k]-pt)
06a0: 20 29 0d 0a 09 09 09 09 09 09 09 6d 69 6e 4b 20 ).........minK
06b0: 3d 20 6b 3b 0d 0a 09 09 09 09 09 64 6f 75 62 6c = k;.......doubl
06c0: 65 20 66 6d 69 6e 20 3d 20 66 75 72 74 68 65 73 e fmin = furthes
06d0: 74 28 7a 2c 20 74 2c 20 63 73 5b 6d 69 6e 4b 5d t(z, t, cs[minK]
06e0: 29 3b 0d 0a 09 09 09 09 09 66 6d 69 6e 5f 6d 61 );.......fmin_ma
06f0: 78 20 3d 20 6d 61 78 28 66 6d 69 6e 5f 6d 61 78 x = max(fmin_max
0700: 2c 20 66 6d 69 6e 29 3b 0d 0a 09 09 09 09 7d 0d , fmin);......}.
0710: 0a 09 09 09 09 7a 20 3d 20 74 3b 0d 0a 09 09 09 .....z = t;.....
0720: 7d 0d 0a 09 09 09 69 66 28 20 7a 3e 3d 31 2e 30 }.....if( z>=1.0
0730: 20 29 0d 0a 09 09 09 09 62 72 65 61 6b 3b 0d 0a )......break;..
0740: 09 09 7d 0d 0a 09 09 69 66 28 20 7a 20 3c 20 31 ..}....if( z < 1
0750: 2e 30 20 29 20 7b 0d 0a 09 09 09 64 6f 75 62 6c .0 ) {.....doubl
0760: 65 20 20 74 20 3d 20 31 2e 30 3b 0d 0a 09 09 09 e t = 1.0;.....
0770: 2f 2f 20 5b 7a 2c 20 31 2e 30 29 0d 0a 09 09 09 // [z, 1.0).....
0780: 7b 0d 0a 09 09 09 09 64 6f 75 62 6c 65 20 70 74 {......double pt
0790: 20 3d 20 28 7a 2b 74 29 2f 32 3b 0d 0a 09 09 09 = (z+t)/2;.....
07a0: 09 69 6e 74 20 6d 69 6e 4b 20 3d 20 30 3b 0d 0a .int minK = 0;..
07b0: 09 09 09 09 66 6f 72 28 69 6e 74 20 6b 3d 31 3b ....for(int k=1;
07c0: 20 6b 3c 63 73 2e 73 69 7a 65 28 29 3b 20 2b 2b k<cs.size(); ++
07d0: 6b 29 0d 0a 09 09 09 09 09 69 66 28 20 61 62 73 k).......if( abs
07e0: 28 63 73 5b 6d 69 6e 4b 5d 2d 70 74 29 20 3e 20 (cs[minK]-pt) >
07f0: 61 62 73 28 63 73 5b 6b 5d 2d 70 74 29 20 29 0d abs(cs[k]-pt) ).
0800: 0a 09 09 09 09 09 09 6d 69 6e 4b 20 3d 20 6b 3b .......minK = k;
0810: 0d 0a 09 09 09 09 64 6f 75 62 6c 65 20 66 6d 69 ......double fmi
0820: 6e 20 3d 20 66 75 72 74 68 65 73 74 28 7a 2c 20 n = furthest(z,
0830: 74 2c 20 63 73 5b 6d 69 6e 4b 5d 29 3b 0d 0a 09 t, cs[minK]);...
0840: 09 09 09 66 6d 69 6e 5f 6d 61 78 20 3d 20 6d 61 ...fmin_max = ma
0850: 78 28 66 6d 69 6e 5f 6d 61 78 2c 20 66 6d 69 6e x(fmin_max, fmin
0860: 29 3b 0d 0a 09 09 09 7d 0d 0a 09 09 7d 0d 0a 09 );.....}....}...
0870: 09 63 65 72 72 20 3c 3c 20 66 61 63 74 6f 72 2a .cerr << factor*
0880: 66 6d 69 6e 5f 6d 61 78 20 3c 3c 20 65 6e 64 6c fmin_max << endl
0890: 3b 0d 0a 09 09 72 65 74 75 72 6e 20 66 61 63 74 ;....return fact
08a0: 6f 72 20 2a 20 66 6d 69 6e 5f 6d 61 78 3b 0d 0a or * fmin_max;..
08b0: 09 7d 0d 0a 0d 0a 09 64 6f 75 62 6c 65 20 73 65 .}.....double se
08c0: 63 74 70 6f 69 6e 74 28 43 4d 50 20 61 2c 20 43 ctpoint(CMP a, C
08d0: 4d 50 20 62 29 0d 0a 09 7b 0d 0a 09 09 43 4d 50 MP b)...{....CMP
08e0: 20 63 20 3d 20 28 61 2b 62 29 2f 32 2e 30 3b 0d c = (a+b)/2.0;.
08f0: 0a 09 09 43 4d 50 20 64 69 72 20 3d 20 28 61 2d ...CMP dir = (a-
0900: 62 29 2a 43 4d 50 28 30 2c 31 29 3b 0d 0a 09 09 b)*CMP(0,1);....
0910: 64 69 72 20 2f 3d 20 61 62 73 28 64 69 72 29 3b dir /= abs(dir);
0920: 0d 0a 09 09 72 65 74 75 72 6e 20 28 63 20 2b 20 ....return (c +
0930: 63 2e 69 6d 61 67 28 29 20 2f 20 64 69 72 2e 69 c.imag() / dir.i
0940: 6d 61 67 28 29 20 2a 20 64 69 72 29 2e 72 65 61 mag() * dir).rea
0950: 6c 28 29 3b 0d 0a 09 7d 0d 0a 0d 0a 09 64 6f 75 l();...}.....dou
0960: 62 6c 65 20 66 75 72 74 68 65 73 74 28 64 6f 75 ble furthest(dou
0970: 62 6c 65 20 61 2c 20 64 6f 75 62 6c 65 20 62 2c ble a, double b,
0980: 20 43 4d 50 20 63 29 0d 0a 09 7b 0d 0a 09 09 69 CMP c)...{....i
0990: 66 28 20 63 2e 72 65 61 6c 28 29 20 3e 20 28 61 f( c.real() > (a
09a0: 2b 62 29 2f 32 20 29 0d 0a 09 09 09 72 65 74 75 +b)/2 ).....retu
09b0: 72 6e 20 61 62 73 28 63 20 2d 20 61 29 3b 0d 0a rn abs(c - a);..
09c0: 09 09 65 6c 73 65 0d 0a 09 09 09 72 65 74 75 72 ..else.....retur
09d0: 6e 20 61 62 73 28 63 20 2d 20 62 29 3b 0d 0a 09 n abs(c - b);...
09e0: 7d 0d 0a 7d 3b 0d 0a 0d 0a 2f 2f 20 42 45 47 49 }..};....// BEGI
09f0: 4e 20 43 55 54 20 48 45 52 45 0d 0a 23 69 6e 63 N CUT HERE..#inc
0a00: 6c 75 64 65 20 3c 63 74 69 6d 65 3e 0d 0a 64 6f lude <ctime>..do
0a10: 75 62 6c 65 20 73 74 61 72 74 5f 74 69 6d 65 3b uble start_time;
0a20: 20 73 74 72 69 6e 67 20 74 69 6d 65 72 28 29 0d string timer().
0a30: 0a 20 7b 20 6f 73 74 72 69 6e 67 73 74 72 65 61 . { ostringstrea
0a40: 6d 20 6f 73 3b 20 6f 73 20 3c 3c 20 22 20 28 22 m os; os << " ("
0a50: 20 3c 3c 20 69 6e 74 28 28 63 6c 6f 63 6b 28 29 << int((clock()
0a60: 2d 73 74 61 72 74 5f 74 69 6d 65 29 2f 43 4c 4f -start_time)/CLO
0a70: 43 4b 53 5f 50 45 52 5f 53 45 43 2a 31 30 30 30 CKS_PER_SEC*1000
0a80: 29 20 3c 3c 20 22 20 6d 73 65 63 29 22 3b 20 72 ) << " msec)"; r
0a90: 65 74 75 72 6e 20 6f 73 2e 73 74 72 28 29 3b 20 eturn os.str();
0aa0: 7d 0d 0a 74 65 6d 70 6c 61 74 65 3c 74 79 70 65 }..template<type
0ab0: 6e 61 6d 65 20 54 3e 20 6f 73 74 72 65 61 6d 26 name T> ostream&
0ac0: 20 6f 70 65 72 61 74 6f 72 3c 3c 28 6f 73 74 72 operator<<(ostr
0ad0: 65 61 6d 26 20 6f 73 2c 20 63 6f 6e 73 74 20 76 eam& os, const v
0ae0: 65 63 74 6f 72 3c 54 3e 26 20 76 29 0d 0a 20 7b ector<T>& v).. {
0af0: 20 6f 73 20 3c 3c 20 22 7b 20 22 3b 0d 0a 20 20 os << "{ ";..
0b00: 20 66 6f 72 28 74 79 70 65 6e 61 6d 65 20 76 65 for(typename ve
0b10: 63 74 6f 72 3c 54 3e 3a 3a 63 6f 6e 73 74 5f 69 ctor<T>::const_i
0b20: 74 65 72 61 74 6f 72 20 69 74 3d 76 2e 62 65 67 terator it=v.beg
0b30: 69 6e 28 29 3b 20 69 74 21 3d 76 2e 65 6e 64 28 in(); it!=v.end(
0b40: 29 3b 20 2b 2b 69 74 29 0d 0a 20 20 20 6f 73 20 ); ++it).. os
0b50: 3c 3c 20 27 5c 22 27 20 3c 3c 20 2a 69 74 20 3c << '\"' << *it <
0b60: 3c 20 27 5c 22 27 20 3c 3c 20 28 69 74 2b 31 3d < '\"' << (it+1=
0b70: 3d 76 2e 65 6e 64 28 29 20 3f 20 22 22 20 3a 20 =v.end() ? "" :
0b80: 22 2c 20 22 29 3b 20 6f 73 20 3c 3c 20 22 20 7d ", "); os << " }
0b90: 22 3b 20 72 65 74 75 72 6e 20 6f 73 3b 20 7d 0d "; return os; }.
0ba0: 0a 76 6f 69 64 20 76 65 72 69 66 79 5f 63 61 73 .void verify_cas
0bb0: 65 28 63 6f 6e 73 74 20 64 6f 75 62 6c 65 26 20 e(const double&
0bc0: 45 78 70 65 63 74 65 64 2c 20 63 6f 6e 73 74 20 Expected, const
0bd0: 64 6f 75 62 6c 65 26 20 52 65 63 65 69 76 65 64 double& Received
0be0: 29 20 7b 0d 0a 20 62 6f 6f 6c 20 6f 6b 20 3d 20 ) {.. bool ok =
0bf0: 28 61 62 73 28 45 78 70 65 63 74 65 64 20 2d 20 (abs(Expected -
0c00: 52 65 63 65 69 76 65 64 29 20 3c 20 31 65 2d 39 Received) < 1e-9
0c10: 29 3b 0d 0a 20 69 66 28 6f 6b 29 20 63 65 72 72 );.. if(ok) cerr
0c20: 20 3c 3c 20 22 50 41 53 53 45 44 22 20 3c 3c 20 << "PASSED" <<
0c30: 74 69 6d 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b timer() << endl;
0c40: 20 20 65 6c 73 65 20 7b 20 63 65 72 72 20 3c 3c else { cerr <<
0c50: 20 22 46 41 49 4c 45 44 22 20 3c 3c 20 74 69 6d "FAILED" << tim
0c60: 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b 0d 0a 20 er() << endl;..
0c70: 63 65 72 72 20 3c 3c 20 22 5c 74 6f 3a 20 5c 22 cerr << "\to: \"
0c80: 22 20 3c 3c 20 45 78 70 65 63 74 65 64 20 3c 3c " << Expected <<
0c90: 20 27 5c 22 27 20 3c 3c 20 65 6e 64 6c 20 3c 3c '\"' << endl <<
0ca0: 20 22 5c 74 78 3a 20 5c 22 22 20 3c 3c 20 52 65 "\tx: \"" << Re
0cb0: 63 65 69 76 65 64 20 3c 3c 20 27 5c 22 27 20 3c ceived << '\"' <
0cc0: 3c 20 65 6e 64 6c 3b 20 7d 20 7d 0d 0a 23 64 65 < endl; } }..#de
0cd0: 66 69 6e 65 20 43 41 53 45 28 4e 29 20 7b 63 65 fine CASE(N) {ce
0ce0: 72 72 20 3c 3c 20 22 54 65 73 74 20 43 61 73 65 rr << "Test Case
0cf0: 20 23 22 20 3c 3c 20 4e 20 3c 3c 20 22 2e 2e 2e #" << N << "...
0d00: 22 20 3c 3c 20 66 6c 75 73 68 3b 20 73 74 61 72 " << flush; star
0d10: 74 5f 74 69 6d 65 3d 63 6c 6f 63 6b 28 29 3b 0d t_time=clock();.
0d20: 0a 23 64 65 66 69 6e 65 20 45 4e 44 09 20 76 65 .#define END. ve
0d30: 72 69 66 79 5f 63 61 73 65 28 5f 2c 20 42 61 72 rify_case(_, Bar
0d40: 62 61 72 69 61 6e 49 6e 76 61 73 69 6f 6e 32 28 barianInvasion2(
0d50: 29 2e 6d 69 6e 69 6d 75 6d 54 69 6d 65 28 62 6f ).minimumTime(bo
0d60: 75 6e 64 61 72 79 58 2c 20 62 6f 75 6e 64 61 72 undaryX, boundar
0d70: 79 59 2c 20 63 69 74 79 58 2c 20 63 69 74 79 59 yY, cityX, cityY
0d80: 29 29 3b 7d 0d 0a 69 6e 74 20 6d 61 69 6e 28 29 ));}..int main()
0d90: 7b 0d 0a 0d 0a 43 41 53 45 28 30 29 0d 0a 09 69 {....CASE(0)...i
0da0: 6e 74 20 62 6f 75 6e 64 61 72 79 58 5f 5b 5d 20 nt boundaryX_[]
0db0: 3d 20 7b 30 2c 32 2c 32 2c 30 7d 3b 0d 0a 09 20 = {0,2,2,0};...
0dc0: 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 62 6f vector <int> bo
0dd0: 75 6e 64 61 72 79 58 28 62 6f 75 6e 64 61 72 79 undaryX(boundary
0de0: 58 5f 2c 20 62 6f 75 6e 64 61 72 79 58 5f 2b 73 X_, boundaryX_+s
0df0: 69 7a 65 6f 66 28 62 6f 75 6e 64 61 72 79 58 5f izeof(boundaryX_
0e00: 29 2f 73 69 7a 65 6f 66 28 2a 62 6f 75 6e 64 61 )/sizeof(*bounda
0e10: 72 79 58 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 62 ryX_)); ...int b
0e20: 6f 75 6e 64 61 72 79 59 5f 5b 5d 20 3d 20 7b 30 oundaryY_[] = {0
0e30: 2c 30 2c 32 2c 32 7d 3b 0d 0a 09 20 20 76 65 63 ,0,2,2};... vec
0e40: 74 6f 72 20 3c 69 6e 74 3e 20 62 6f 75 6e 64 61 tor <int> bounda
0e50: 72 79 59 28 62 6f 75 6e 64 61 72 79 59 5f 2c 20 ryY(boundaryY_,
0e60: 62 6f 75 6e 64 61 72 79 59 5f 2b 73 69 7a 65 6f boundaryY_+sizeo
0e70: 66 28 62 6f 75 6e 64 61 72 79 59 5f 29 2f 73 69 f(boundaryY_)/si
0e80: 7a 65 6f 66 28 2a 62 6f 75 6e 64 61 72 79 59 5f zeof(*boundaryY_
0e90: 29 29 3b 20 0d 0a 09 69 6e 74 20 63 69 74 79 58 )); ...int cityX
0ea0: 5f 5b 5d 20 3d 20 7b 31 7d 3b 0d 0a 09 20 20 76 _[] = {1};... v
0eb0: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 63 69 74 79 ector <int> city
0ec0: 58 28 63 69 74 79 58 5f 2c 20 63 69 74 79 58 5f X(cityX_, cityX_
0ed0: 2b 73 69 7a 65 6f 66 28 63 69 74 79 58 5f 29 2f +sizeof(cityX_)/
0ee0: 73 69 7a 65 6f 66 28 2a 63 69 74 79 58 5f 29 29 sizeof(*cityX_))
0ef0: 3b 20 0d 0a 09 69 6e 74 20 63 69 74 79 59 5f 5b ; ...int cityY_[
0f00: 5d 20 3d 20 7b 31 7d 3b 0d 0a 09 20 20 76 65 63 ] = {1};... vec
0f10: 74 6f 72 20 3c 69 6e 74 3e 20 63 69 74 79 59 28 tor <int> cityY(
0f20: 63 69 74 79 59 5f 2c 20 63 69 74 79 59 5f 2b 73 cityY_, cityY_+s
0f30: 69 7a 65 6f 66 28 63 69 74 79 59 5f 29 2f 73 69 izeof(cityY_)/si
0f40: 7a 65 6f 66 28 2a 63 69 74 79 59 5f 29 29 3b 20 zeof(*cityY_));
0f50: 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20 31 2e ...double _ = 1.
0f60: 34 31 34 32 31 33 35 36 32 33 37 33 30 38 38 3b 414213562373088;
0f70: 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 31 29 0d ..END..CASE(1).
0f80: 0a 09 69 6e 74 20 62 6f 75 6e 64 61 72 79 58 5f ..int boundaryX_
0f90: 5b 5d 20 3d 20 7b 30 2c 33 2c 33 2c 30 7d 3b 0d [] = {0,3,3,0};.
0fa0: 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e .. vector <int>
0fb0: 20 62 6f 75 6e 64 61 72 79 58 28 62 6f 75 6e 64 boundaryX(bound
0fc0: 61 72 79 58 5f 2c 20 62 6f 75 6e 64 61 72 79 58 aryX_, boundaryX
0fd0: 5f 2b 73 69 7a 65 6f 66 28 62 6f 75 6e 64 61 72 _+sizeof(boundar
0fe0: 79 58 5f 29 2f 73 69 7a 65 6f 66 28 2a 62 6f 75 yX_)/sizeof(*bou
0ff0: 6e 64 61 72 79 58 5f 29 29 3b 20 0d 0a 09 69 6e ndaryX_)); ...in
1000: 74 20 62 6f 75 6e 64 61 72 79 59 5f 5b 5d 20 3d t boundaryY_[] =
1010: 20 7b 30 2c 30 2c 33 2c 33 7d 3b 0d 0a 09 20 20 {0,0,3,3};...
1020: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 62 6f 75 vector <int> bou
1030: 6e 64 61 72 79 59 28 62 6f 75 6e 64 61 72 79 59 ndaryY(boundaryY
1040: 5f 2c 20 62 6f 75 6e 64 61 72 79 59 5f 2b 73 69 _, boundaryY_+si
1050: 7a 65 6f 66 28 62 6f 75 6e 64 61 72 79 59 5f 29 zeof(boundaryY_)
1060: 2f 73 69 7a 65 6f 66 28 2a 62 6f 75 6e 64 61 72 /sizeof(*boundar
1070: 79 59 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 63 69 yY_)); ...int ci
1080: 74 79 58 5f 5b 5d 20 3d 20 7b 31 7d 3b 0d 0a 09 tyX_[] = {1};...
1090: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 63 vector <int> c
10a0: 69 74 79 58 28 63 69 74 79 58 5f 2c 20 63 69 74 ityX(cityX_, cit
10b0: 79 58 5f 2b 73 69 7a 65 6f 66 28 63 69 74 79 58 yX_+sizeof(cityX
10c0: 5f 29 2f 73 69 7a 65 6f 66 28 2a 63 69 74 79 58 _)/sizeof(*cityX
10d0: 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 63 69 74 79 _)); ...int city
10e0: 59 5f 5b 5d 20 3d 20 7b 31 7d 3b 0d 0a 09 20 20 Y_[] = {1};...
10f0: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 63 69 74 vector <int> cit
1100: 79 59 28 63 69 74 79 59 5f 2c 20 63 69 74 79 59 yY(cityY_, cityY
1110: 5f 2b 73 69 7a 65 6f 66 28 63 69 74 79 59 5f 29 _+sizeof(cityY_)
1120: 2f 73 69 7a 65 6f 66 28 2a 63 69 74 79 59 5f 29 /sizeof(*cityY_)
1130: 29 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 3d ); ...double _ =
1140: 20 32 2e 38 32 38 34 32 37 31 32 34 37 34 36 31 2.8284271247461
1150: 34 38 35 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 485; ..END..CASE
1160: 28 32 29 0d 0a 09 69 6e 74 20 62 6f 75 6e 64 61 (2)...int bounda
1170: 72 79 58 5f 5b 5d 20 3d 20 7b 30 2c 33 2c 33 2c ryX_[] = {0,3,3,
1180: 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 0};... vector <
1190: 69 6e 74 3e 20 62 6f 75 6e 64 61 72 79 58 28 62 int> boundaryX(b
11a0: 6f 75 6e 64 61 72 79 58 5f 2c 20 62 6f 75 6e 64 oundaryX_, bound
11b0: 61 72 79 58 5f 2b 73 69 7a 65 6f 66 28 62 6f 75 aryX_+sizeof(bou
11c0: 6e 64 61 72 79 58 5f 29 2f 73 69 7a 65 6f 66 28 ndaryX_)/sizeof(
11d0: 2a 62 6f 75 6e 64 61 72 79 58 5f 29 29 3b 20 0d *boundaryX_)); .
11e0: 0a 09 69 6e 74 20 62 6f 75 6e 64 61 72 79 59 5f ..int boundaryY_
11f0: 5b 5d 20 3d 20 7b 30 2c 30 2c 33 2c 33 7d 3b 0d [] = {0,0,3,3};.
1200: 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e .. vector <int>
1210: 20 62 6f 75 6e 64 61 72 79 59 28 62 6f 75 6e 64 boundaryY(bound
1220: 61 72 79 59 5f 2c 20 62 6f 75 6e 64 61 72 79 59 aryY_, boundaryY
1230: 5f 2b 73 69 7a 65 6f 66 28 62 6f 75 6e 64 61 72 _+sizeof(boundar
1240: 79 59 5f 29 2f 73 69 7a 65 6f 66 28 2a 62 6f 75 yY_)/sizeof(*bou
1250: 6e 64 61 72 79 59 5f 29 29 3b 20 0d 0a 09 69 6e ndaryY_)); ...in
1260: 74 20 63 69 74 79 58 5f 5b 5d 20 3d 20 7b 31 2c t cityX_[] = {1,
1270: 32 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 2};... vector <
1280: 69 6e 74 3e 20 63 69 74 79 58 28 63 69 74 79 58 int> cityX(cityX
1290: 5f 2c 20 63 69 74 79 58 5f 2b 73 69 7a 65 6f 66 _, cityX_+sizeof
12a0: 28 63 69 74 79 58 5f 29 2f 73 69 7a 65 6f 66 28 (cityX_)/sizeof(
12b0: 2a 63 69 74 79 58 5f 29 29 3b 20 0d 0a 09 69 6e *cityX_)); ...in
12c0: 74 20 63 69 74 79 59 5f 5b 5d 20 3d 20 7b 32 2c t cityY_[] = {2,
12d0: 31 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 1};... vector <
12e0: 69 6e 74 3e 20 63 69 74 79 59 28 63 69 74 79 59 int> cityY(cityY
12f0: 5f 2c 20 63 69 74 79 59 5f 2b 73 69 7a 65 6f 66 _, cityY_+sizeof
1300: 28 63 69 74 79 59 5f 29 2f 73 69 7a 65 6f 66 28 (cityY_)/sizeof(
1310: 2a 63 69 74 79 59 5f 29 29 3b 20 0d 0a 09 64 6f *cityY_)); ...do
1320: 75 62 6c 65 20 5f 20 3d 20 32 2e 32 33 36 30 36 uble _ = 2.23606
1330: 37 39 37 37 34 39 39 37 37 32 3b 20 0d 0a 45 4e 7977499772; ..EN
1340: 44 0d 0a 43 41 53 45 28 33 29 0d 0a 09 69 6e 74 D..CASE(3)...int
1350: 20 62 6f 75 6e 64 61 72 79 58 5f 5b 5d 20 3d 20 boundaryX_[] =
1360: 7b 30 2c 34 30 2c 34 30 2c 30 7d 3b 0d 0a 09 20 {0,40,40,0};...
1370: 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 62 6f vector <int> bo
1380: 75 6e 64 61 72 79 58 28 62 6f 75 6e 64 61 72 79 undaryX(boundary
1390: 58 5f 2c 20 62 6f 75 6e 64 61 72 79 58 5f 2b 73 X_, boundaryX_+s
13a0: 69 7a 65 6f 66 28 62 6f 75 6e 64 61 72 79 58 5f izeof(boundaryX_
13b0: 29 2f 73 69 7a 65 6f 66 28 2a 62 6f 75 6e 64 61 )/sizeof(*bounda
13c0: 72 79 58 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 62 ryX_)); ...int b
13d0: 6f 75 6e 64 61 72 79 59 5f 5b 5d 20 3d 20 7b 30 oundaryY_[] = {0
13e0: 2c 30 2c 34 30 2c 34 30 7d 3b 0d 0a 09 20 20 76 ,0,40,40};... v
13f0: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 62 6f 75 6e ector <int> boun
1400: 64 61 72 79 59 28 62 6f 75 6e 64 61 72 79 59 5f daryY(boundaryY_
1410: 2c 20 62 6f 75 6e 64 61 72 79 59 5f 2b 73 69 7a , boundaryY_+siz
1420: 65 6f 66 28 62 6f 75 6e 64 61 72 79 59 5f 29 2f eof(boundaryY_)/
1430: 73 69 7a 65 6f 66 28 2a 62 6f 75 6e 64 61 72 79 sizeof(*boundary
1440: 59 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 63 69 74 Y_)); ...int cit
1450: 79 58 5f 5b 5d 20 3d 20 7b 31 2c 32 2c 33 31 2c yX_[] = {1,2,31,
1460: 32 2c 31 35 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 2,15};... vecto
1470: 72 20 3c 69 6e 74 3e 20 63 69 74 79 58 28 63 69 r <int> cityX(ci
1480: 74 79 58 5f 2c 20 63 69 74 79 58 5f 2b 73 69 7a tyX_, cityX_+siz
1490: 65 6f 66 28 63 69 74 79 58 5f 29 2f 73 69 7a 65 eof(cityX_)/size
14a0: 6f 66 28 2a 63 69 74 79 58 5f 29 29 3b 20 0d 0a of(*cityX_)); ..
14b0: 09 69 6e 74 20 63 69 74 79 59 5f 5b 5d 20 3d 20 .int cityY_[] =
14c0: 7b 31 2c 32 2c 33 2c 33 2c 32 34 7d 3b 0d 0a 09 {1,2,3,3,24};...
14d0: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 63 vector <int> c
14e0: 69 74 79 59 28 63 69 74 79 59 5f 2c 20 63 69 74 ityY(cityY_, cit
14f0: 79 59 5f 2b 73 69 7a 65 6f 66 28 63 69 74 79 59 yY_+sizeof(cityY
1500: 5f 29 2f 73 69 7a 65 6f 66 28 2a 63 69 74 79 59 _)/sizeof(*cityY
1510: 5f 29 29 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 5f _)); ...double _
1520: 20 3d 20 33 38 2e 30 35 37 34 38 31 35 33 35 35 = 38.0574815355
1530: 31 39 39 34 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 1994; ..END..CAS
1540: 45 28 34 29 0d 0a 09 69 6e 74 20 62 6f 75 6e 64 E(4)...int bound
1550: 61 72 79 58 5f 5b 5d 20 3d 20 7b 30 2c 31 32 34 aryX_[] = {0,124
1560: 2c 2d 36 2c 2d 31 32 30 2c 2d 33 30 30 7d 3b 0d ,-6,-120,-300};.
1570: 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e .. vector <int>
1580: 20 62 6f 75 6e 64 61 72 79 58 28 62 6f 75 6e 64 boundaryX(bound
1590: 61 72 79 58 5f 2c 20 62 6f 75 6e 64 61 72 79 58 aryX_, boundaryX
15a0: 5f 2b 73 69 7a 65 6f 66 28 62 6f 75 6e 64 61 72 _+sizeof(boundar
15b0: 79 58 5f 29 2f 73 69 7a 65 6f 66 28 2a 62 6f 75 yX_)/sizeof(*bou
15c0: 6e 64 61 72 79 58 5f 29 29 3b 20 0d 0a 09 69 6e ndaryX_)); ...in
15d0: 74 20 62 6f 75 6e 64 61 72 79 59 5f 5b 5d 20 3d t boundaryY_[] =
15e0: 20 7b 30 2c 31 32 35 2c 31 34 30 2c 31 33 37 2c {0,125,140,137,
15f0: 2d 31 30 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f -100};... vecto
1600: 72 20 3c 69 6e 74 3e 20 62 6f 75 6e 64 61 72 79 r <int> boundary
1610: 59 28 62 6f 75 6e 64 61 72 79 59 5f 2c 20 62 6f Y(boundaryY_, bo
1620: 75 6e 64 61 72 79 59 5f 2b 73 69 7a 65 6f 66 28 undaryY_+sizeof(
1630: 62 6f 75 6e 64 61 72 79 59 5f 29 2f 73 69 7a 65 boundaryY_)/size
1640: 6f 66 28 2a 62 6f 75 6e 64 61 72 79 59 5f 29 29 of(*boundaryY_))
1650: 3b 20 0d 0a 09 69 6e 74 20 63 69 74 79 58 5f 5b ; ...int cityX_[
1660: 5d 20 3d 20 7b 31 30 2c 31 30 2c 31 30 2c 31 30 ] = {10,10,10,10
1670: 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 };... vector <i
1680: 6e 74 3e 20 63 69 74 79 58 28 63 69 74 79 58 5f nt> cityX(cityX_
1690: 2c 20 63 69 74 79 58 5f 2b 73 69 7a 65 6f 66 28 , cityX_+sizeof(
16a0: 63 69 74 79 58 5f 29 2f 73 69 7a 65 6f 66 28 2a cityX_)/sizeof(*
16b0: 63 69 74 79 58 5f 29 29 3b 20 0d 0a 09 69 6e 74 cityX_)); ...int
16c0: 20 63 69 74 79 59 5f 5b 5d 20 3d 20 7b 35 30 2c cityY_[] = {50,
16d0: 35 31 2c 35 32 2c 32 31 7d 3b 0d 0a 09 20 20 76 51,52,21};... v
16e0: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 63 69 74 79 ector <int> city
16f0: 59 28 63 69 74 79 59 5f 2c 20 63 69 74 79 59 5f Y(cityY_, cityY_
1700: 2b 73 69 7a 65 6f 66 28 63 69 74 79 59 5f 29 2f +sizeof(cityY_)/
1710: 73 69 7a 65 6f 66 28 2a 63 69 74 79 59 5f 29 29 sizeof(*cityY_))
1720: 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20 ; ...double _ =
1730: 33 33 32 2e 37 37 37 37 30 33 35 38 30 30 32 37 332.777703580027
1740: 38 33 3b 20 0d 0a 45 4e 44 0d 0a 2f 2a 0d 0a 43 83; ..END../*..C
1750: 41 53 45 28 35 29 0d 0a 09 69 6e 74 20 62 6f 75 ASE(5)...int bou
1760: 6e 64 61 72 79 58 5f 5b 5d 20 3d 20 3b 0d 0a 09 ndaryX_[] = ;...
1770: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 62 vector <int> b
1780: 6f 75 6e 64 61 72 79 58 28 62 6f 75 6e 64 61 72 oundaryX(boundar
1790: 79 58 5f 2c 20 62 6f 75 6e 64 61 72 79 58 5f 2b yX_, boundaryX_+
17a0: 73 69 7a 65 6f 66 28 62 6f 75 6e 64 61 72 79 58 sizeof(boundaryX
17b0: 5f 29 2f 73 69 7a 65 6f 66 28 2a 62 6f 75 6e 64 _)/sizeof(*bound
17c0: 61 72 79 58 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 aryX_)); ...int
17d0: 62 6f 75 6e 64 61 72 79 59 5f 5b 5d 20 3d 20 3b boundaryY_[] = ;
17e0: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 ... vector <int
17f0: 3e 20 62 6f 75 6e 64 61 72 79 59 28 62 6f 75 6e > boundaryY(boun
1800: 64 61 72 79 59 5f 2c 20 62 6f 75 6e 64 61 72 79 daryY_, boundary
1810: 59 5f 2b 73 69 7a 65 6f 66 28 62 6f 75 6e 64 61 Y_+sizeof(bounda
1820: 72 79 59 5f 29 2f 73 69 7a 65 6f 66 28 2a 62 6f ryY_)/sizeof(*bo
1830: 75 6e 64 61 72 79 59 5f 29 29 3b 20 0d 0a 09 69 undaryY_)); ...i
1840: 6e 74 20 63 69 74 79 58 5f 5b 5d 20 3d 20 3b 0d nt cityX_[] = ;.
1850: 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e .. vector <int>
1860: 20 63 69 74 79 58 28 63 69 74 79 58 5f 2c 20 63 cityX(cityX_, c
1870: 69 74 79 58 5f 2b 73 69 7a 65 6f 66 28 63 69 74 ityX_+sizeof(cit
1880: 79 58 5f 29 2f 73 69 7a 65 6f 66 28 2a 63 69 74 yX_)/sizeof(*cit
1890: 79 58 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 63 69 yX_)); ...int ci
18a0: 74 79 59 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 tyY_[] = ;... v
18b0: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 63 69 74 79 ector <int> city
18c0: 59 28 63 69 74 79 59 5f 2c 20 63 69 74 79 59 5f Y(cityY_, cityY_
18d0: 2b 73 69 7a 65 6f 66 28 63 69 74 79 59 5f 29 2f +sizeof(cityY_)/
18e0: 73 69 7a 65 6f 66 28 2a 63 69 74 79 59 5f 29 29 sizeof(*cityY_))
18f0: 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20 ; ...double _ =
1900: 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 36 29 ; ..END..CASE(6)
1910: 0d 0a 09 69 6e 74 20 62 6f 75 6e 64 61 72 79 58 ...int boundaryX
1920: 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 65 63 74 _[] = ;... vect
1930: 6f 72 20 3c 69 6e 74 3e 20 62 6f 75 6e 64 61 72 or <int> boundar
1940: 79 58 28 62 6f 75 6e 64 61 72 79 58 5f 2c 20 62 yX(boundaryX_, b
1950: 6f 75 6e 64 61 72 79 58 5f 2b 73 69 7a 65 6f 66 oundaryX_+sizeof
1960: 28 62 6f 75 6e 64 61 72 79 58 5f 29 2f 73 69 7a (boundaryX_)/siz
1970: 65 6f 66 28 2a 62 6f 75 6e 64 61 72 79 58 5f 29 eof(*boundaryX_)
1980: 29 3b 20 0d 0a 09 69 6e 74 20 62 6f 75 6e 64 61 ); ...int bounda
1990: 72 79 59 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 ryY_[] = ;... v
19a0: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 62 6f 75 6e ector <int> boun
19b0: 64 61 72 79 59 28 62 6f 75 6e 64 61 72 79 59 5f daryY(boundaryY_
19c0: 2c 20 62 6f 75 6e 64 61 72 79 59 5f 2b 73 69 7a , boundaryY_+siz
19d0: 65 6f 66 28 62 6f 75 6e 64 61 72 79 59 5f 29 2f eof(boundaryY_)/
19e0: 73 69 7a 65 6f 66 28 2a 62 6f 75 6e 64 61 72 79 sizeof(*boundary
19f0: 59 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 63 69 74 Y_)); ...int cit
1a00: 79 58 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 65 yX_[] = ;... ve
1a10: 63 74 6f 72 20 3c 69 6e 74 3e 20 63 69 74 79 58 ctor <int> cityX
1a20: 28 63 69 74 79 58 5f 2c 20 63 69 74 79 58 5f 2b (cityX_, cityX_+
1a30: 73 69 7a 65 6f 66 28 63 69 74 79 58 5f 29 2f 73 sizeof(cityX_)/s
1a40: 69 7a 65 6f 66 28 2a 63 69 74 79 58 5f 29 29 3b izeof(*cityX_));
1a50: 20 0d 0a 09 69 6e 74 20 63 69 74 79 59 5f 5b 5d ...int cityY_[]
1a60: 20 3d 20 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 = ;... vector
1a70: 3c 69 6e 74 3e 20 63 69 74 79 59 28 63 69 74 79 <int> cityY(city
1a80: 59 5f 2c 20 63 69 74 79 59 5f 2b 73 69 7a 65 6f Y_, cityY_+sizeo
1a90: 66 28 63 69 74 79 59 5f 29 2f 73 69 7a 65 6f 66 f(cityY_)/sizeof
1aa0: 28 2a 63 69 74 79 59 5f 29 29 3b 20 0d 0a 09 64 (*cityY_)); ...d
1ab0: 6f 75 62 6c 65 20 5f 20 3d 20 3b 20 0d 0a 45 4e ouble _ = ; ..EN
1ac0: 44 0d 0a 2a 2f 0d 0a 7d 0d 0a 2f 2f 20 45 4e 44 D..*/..}..// END
1ad0: 20 43 55 54 20 48 45 52 45 0d 0a CUT HERE..