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 75 73 69 6e 67 20 6e 61 6d sert>..using nam
0140: 65 73 70 61 63 65 20 73 74 64 3b 0d 0a 74 79 70 espace std;..typ
0150: 65 64 65 66 20 6c 6f 6e 67 20 6c 6f 6e 67 20 4c edef long long L
0160: 4c 3b 0d 0a 74 79 70 65 64 65 66 20 63 6f 6d 70 L;..typedef comp
0170: 6c 65 78 3c 64 6f 75 62 6c 65 3e 20 43 4d 50 3b lex<double> CMP;
0180: 0d 0a 0d 0a 73 74 72 75 63 74 20 43 69 72 63 6c ....struct Circl
0190: 65 20 7b 0d 0a 09 43 4d 50 20 6f 3b 0d 0a 09 64 e {...CMP o;...d
01a0: 6f 75 62 6c 65 20 72 3b 0d 0a 7d 3b 0d 0a 0d 0a ouble r;..};....
01b0: 62 6f 6f 6c 20 6c 69 6e 65 5f 63 69 72 63 6c 65 bool line_circle
01c0: 28 43 4d 50 20 6c 61 2c 20 43 4d 50 20 6c 62 2c (CMP la, CMP lb,
01d0: 20 43 4d 50 20 63 2c 20 64 6f 75 62 6c 65 20 72 CMP c, double r
01e0: 2c 20 43 4d 50 2a 20 70 31 2c 20 43 4d 50 2a 20 , CMP* p1, CMP*
01f0: 70 32 29 0d 0a 7b 0d 0a 09 43 4d 50 20 76 20 3d p2)..{...CMP v =
0200: 20 28 6c 62 2d 6c 61 29 20 2f 20 61 62 73 28 6c (lb-la) / abs(l
0210: 62 2d 6c 61 29 3b 0d 0a 09 43 4d 50 20 6f 20 3d b-la);...CMP o =
0220: 20 28 63 2d 6c 61 29 20 2f 20 76 3b 0d 0a 09 69 (c-la) / v;...i
0230: 66 28 20 61 62 73 28 6f 2e 69 6d 61 67 28 29 29 f( abs(o.imag())
0240: 20 3e 20 72 20 29 0d 0a 09 09 72 65 74 75 72 6e > r )....return
0250: 20 66 61 6c 73 65 3b 0d 0a 09 64 6f 75 62 6c 65 false;...double
0260: 20 64 78 20 3d 20 73 71 72 74 28 72 2a 72 20 2d dx = sqrt(r*r -
0270: 20 6f 2e 69 6d 61 67 28 29 2a 6f 2e 69 6d 61 67 o.imag()*o.imag
0280: 28 29 29 3b 0d 0a 09 2a 70 31 20 3d 20 6c 61 20 ());...*p1 = la
0290: 2b 20 28 6f 2e 72 65 61 6c 28 29 2d 64 78 29 2a + (o.real()-dx)*
02a0: 76 3b 0d 0a 09 2a 70 32 20 3d 20 6c 61 20 2b 20 v;...*p2 = la +
02b0: 28 6f 2e 72 65 61 6c 28 29 2b 64 78 29 2a 76 3b (o.real()+dx)*v;
02c0: 0d 0a 09 72 65 74 75 72 6e 20 74 72 75 65 3b 0d ...return true;.
02d0: 0a 7d 0d 0a 0d 0a 63 6c 61 73 73 20 43 69 72 63 .}....class Circ
02e0: 75 73 54 65 6e 74 73 20 7b 20 70 75 62 6c 69 63 usTents { public
02f0: 3a 0d 0a 09 64 6f 75 62 6c 65 20 66 69 6e 64 4d :...double findM
0300: 61 78 69 6d 75 6d 44 69 73 74 61 6e 63 65 28 76 aximumDistance(v
0310: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 78 2c 20 76 ector <int> x, v
0320: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79 2c 20 76 ector <int> y, v
0330: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 72 29 0d 0a ector <int> r)..
0340: 09 7b 0d 0a 09 09 64 6f 75 62 6c 65 20 66 20 3d .{....double f =
0350: 20 72 5b 30 5d 3b 0d 0a 0d 0a 09 09 76 65 63 74 r[0];......vect
0360: 6f 72 3c 43 69 72 63 6c 65 3e 20 63 73 3b 0d 0a or<Circle> cs;..
0370: 09 09 66 6f 72 28 69 6e 74 20 69 3d 31 3b 20 69 ..for(int i=1; i
0380: 3c 78 2e 73 69 7a 65 28 29 3b 20 2b 2b 69 29 20 <x.size(); ++i)
0390: 7b 0d 0a 09 09 09 43 69 72 63 6c 65 20 79 6f 75 {.....Circle you
03a0: 20 3d 20 7b 43 4d 50 28 78 5b 31 5d 2d 78 5b 30 = {CMP(x[1]-x[0
03b0: 5d 2c 20 79 5b 31 5d 2d 79 5b 30 5d 29 2f 66 2c ], y[1]-y[0])/f,
03c0: 20 64 6f 75 62 6c 65 28 72 5b 31 5d 29 2f 66 7d double(r[1])/f}
03d0: 3b 0d 0a 09 09 09 63 73 2e 70 75 73 68 5f 62 61 ;.....cs.push_ba
03e0: 63 6b 28 79 6f 75 29 3b 0d 0a 09 09 7d 0d 0a 0d ck(you);....}...
03f0: 0a 09 09 72 65 74 75 72 6e 20 73 6f 6c 76 65 28 ...return solve(
0400: 63 73 29 20 2a 20 66 3b 0d 0a 09 7d 0d 0a 0d 0a cs) * f;...}....
0410: 09 64 6f 75 62 6c 65 20 73 6f 6c 76 65 28 63 6f .double solve(co
0420: 6e 73 74 20 76 65 63 74 6f 72 3c 43 69 72 63 6c nst vector<Circl
0430: 65 3e 26 20 63 73 29 0d 0a 09 7b 0d 0a 09 09 76 e>& cs)...{....v
0440: 65 63 74 6f 72 3c 64 6f 75 62 6c 65 3e 20 61 72 ector<double> ar
0450: 67 73 3b 0d 0a 0d 0a 09 09 66 6f 72 28 69 6e 74 gs;......for(int
0460: 20 69 3d 30 3b 20 69 3c 63 73 2e 73 69 7a 65 28 i=0; i<cs.size(
0470: 29 3b 20 2b 2b 69 29 0d 0a 09 09 66 6f 72 28 69 ); ++i)....for(i
0480: 6e 74 20 6b 3d 69 2b 31 3b 20 6b 3c 63 73 2e 73 nt k=i+1; k<cs.s
0490: 69 7a 65 28 29 3b 20 2b 2b 6b 29 20 7b 0d 0a 09 ize(); ++k) {...
04a0: 09 09 43 4d 50 20 70 20 3d 20 63 73 5b 69 5d 2e ..CMP p = cs[i].
04b0: 6f 3b 0d 0a 09 09 09 43 4d 50 20 71 20 3d 20 63 o;.....CMP q = c
04c0: 73 5b 6b 5d 2e 6f 3b 0d 0a 0d 0a 09 09 09 64 6f s[k].o;.......do
04d0: 75 62 6c 65 20 6c 65 6e 20 3d 20 61 62 73 28 71 uble len = abs(q
04e0: 2d 70 29 3b 0d 0a 09 09 09 64 6f 75 62 6c 65 20 -p);.....double
04f0: 6d 69 64 20 3d 20 28 6c 65 6e 20 2d 20 63 73 5b mid = (len - cs[
0500: 69 5d 2e 72 20 2d 20 63 73 5b 6b 5d 2e 72 29 2f i].r - cs[k].r)/
0510: 32 2b 63 73 5b 69 5d 2e 72 3b 0d 0a 09 09 09 43 2+cs[i].r;.....C
0520: 4d 50 20 72 20 3d 20 28 71 2d 70 29 2f 6c 65 6e MP r = (q-p)/len
0530: 2a 6d 69 64 3b 0d 0a 09 09 09 43 4d 50 20 64 20 *mid;.....CMP d
0540: 3d 20 28 71 2d 70 29 2f 6c 65 6e 2a 43 4d 50 28 = (q-p)/len*CMP(
0550: 30 2c 31 29 3b 0d 0a 09 09 09 2f 2f 2f 20 72 20 0,1);...../// r
0560: 2b 20 74 64 20 69 73 20 74 68 65 20 73 70 6c 69 + td is the spli
0570: 74 74 69 6e 67 20 6c 69 6e 65 2e 0d 0a 0d 0a 09 tting line......
0580: 09 09 43 4d 50 20 78 31 2c 20 78 32 3b 0d 0a 09 ..CMP x1, x2;...
0590: 09 09 69 66 28 20 6c 69 6e 65 5f 63 69 72 63 6c ..if( line_circl
05a0: 65 28 72 2c 20 72 2b 64 2c 20 43 4d 50 28 30 2c e(r, r+d, CMP(0,
05b0: 30 29 2c 20 31 2e 30 2c 20 26 78 31 2c 20 26 78 0), 1.0, &x1, &x
05c0: 32 29 20 29 20 7b 0d 0a 09 09 09 09 61 72 67 73 2) ) {......args
05d0: 2e 70 75 73 68 5f 62 61 63 6b 28 61 72 67 28 78 .push_back(arg(x
05e0: 31 29 29 3b 0d 0a 09 09 09 09 61 72 67 73 2e 70 1));......args.p
05f0: 75 73 68 5f 62 61 63 6b 28 61 72 67 28 78 32 29 ush_back(arg(x2)
0600: 29 3b 0d 0a 09 09 09 7d 0d 0a 09 09 7d 0d 0a 0d );.....}....}...
0610: 0a 09 09 69 66 28 61 72 67 73 2e 65 6d 70 74 79 ...if(args.empty
0620: 28 29 29 20 7b 0d 0a 09 09 09 61 72 67 73 2e 70 ()) {.....args.p
0630: 75 73 68 5f 62 61 63 6b 28 30 2e 30 29 3b 0d 0a ush_back(0.0);..
0640: 09 09 09 61 72 67 73 2e 70 75 73 68 5f 62 61 63 ...args.push_bac
0650: 6b 28 31 2e 30 29 3b 0d 0a 09 09 7d 0d 0a 09 09 k(1.0);....}....
0660: 73 6f 72 74 28 61 72 67 73 2e 62 65 67 69 6e 28 sort(args.begin(
0670: 29 2c 20 61 72 67 73 2e 65 6e 64 28 29 29 3b 0d ), args.end());.
0680: 0a 0d 0a 09 09 64 6f 75 62 6c 65 20 62 65 73 74 .....double best
0690: 20 3d 20 30 2e 30 3b 0d 0a 09 09 66 6f 72 28 69 = 0.0;....for(i
06a0: 6e 74 20 69 3d 30 3b 20 69 3c 61 72 67 73 2e 73 nt i=0; i<args.s
06b0: 69 7a 65 28 29 3b 20 2b 2b 69 29 20 7b 0d 0a 09 ize(); ++i) {...
06c0: 09 09 64 6f 75 62 6c 65 20 61 4c 20 3d 20 61 72 ..double aL = ar
06d0: 67 73 5b 69 5d 3b 0d 0a 09 09 09 64 6f 75 62 6c gs[i];.....doubl
06e0: 65 20 61 52 20 3d 20 28 69 2b 31 3d 3d 61 72 67 e aR = (i+1==arg
06f0: 73 2e 73 69 7a 65 28 29 20 3f 20 61 72 67 73 5b s.size() ? args[
0700: 30 5d 2b 4d 5f 50 49 2f 32 20 3a 20 61 72 67 73 0]+M_PI/2 : args
0710: 5b 69 2b 31 5d 29 3b 0d 0a 0d 0a 09 09 09 64 6f [i+1]);.......do
0720: 75 62 6c 65 20 61 43 20 3d 20 28 61 4c 2b 61 52 uble aC = (aL+aR
0730: 29 20 2f 20 32 3b 0d 0a 09 09 09 43 4d 50 20 70 ) / 2;.....CMP p
0740: 20 3d 20 70 6f 6c 61 72 28 31 2e 30 2c 20 61 43 = polar(1.0, aC
0750: 29 3b 0d 0a 0d 0a 09 09 09 69 6e 74 20 63 6c 6f );.......int clo
0760: 73 65 20 3d 20 2d 31 3b 20 64 6f 75 62 6c 65 20 se = -1; double
0770: 63 6c 6f 73 65 5f 64 20 3d 20 39 39 39 39 39 39 close_d = 999999
0780: 39 39 39 3b 0d 0a 09 09 09 66 6f 72 28 69 6e 74 999;.....for(int
0790: 20 6b 3d 30 3b 20 6b 3c 63 73 2e 73 69 7a 65 28 k=0; k<cs.size(
07a0: 29 3b 20 2b 2b 6b 29 20 7b 0d 0a 09 09 09 09 64 ); ++k) {......d
07b0: 6f 75 62 6c 65 20 64 20 3d 20 61 62 73 28 70 20 ouble d = abs(p
07c0: 2d 20 63 73 5b 6b 5d 2e 6f 29 20 2d 20 63 73 5b - cs[k].o) - cs[
07d0: 6b 5d 2e 72 3b 0d 0a 09 09 09 09 69 66 28 64 3c k].r;......if(d<
07e0: 63 6c 6f 73 65 5f 64 29 20 7b 63 6c 6f 73 65 3d close_d) {close=
07f0: 6b 2c 20 63 6c 6f 73 65 5f 64 3d 64 3b 7d 0d 0a k, close_d=d;}..
0800: 09 09 09 7d 0d 0a 0d 0a 09 09 09 43 4d 50 20 70 ...}.......CMP p
0810: 4c 20 3d 20 70 6f 6c 61 72 28 31 2e 30 2c 20 61 L = polar(1.0, a
0820: 4c 29 3b 0d 0a 09 09 09 43 4d 50 20 70 52 20 3d L);.....CMP pR =
0830: 20 70 6f 6c 61 72 28 31 2e 30 2c 20 61 52 29 3b polar(1.0, aR);
0840: 0d 0a 09 09 09 43 4d 50 20 76 20 3d 20 63 73 5b .....CMP v = cs[
0850: 63 6c 6f 73 65 5d 2e 6f 3b 0d 0a 09 09 09 70 4c close].o;.....pL
0860: 20 2f 3d 20 76 2f 61 62 73 28 76 29 3b 0d 0a 09 /= v/abs(v);...
0870: 09 09 70 52 20 2f 3d 20 76 2f 61 62 73 28 76 29 ..pR /= v/abs(v)
0880: 3b 0d 0a 09 09 09 64 6f 75 62 6c 65 20 61 61 4c ;.....double aaL
0890: 20 3d 20 61 72 67 28 70 4c 29 3b 0d 0a 09 09 09 = arg(pL);.....
08a0: 64 6f 75 62 6c 65 20 61 61 52 20 3d 20 61 72 67 double aaR = arg
08b0: 28 70 52 29 3b 0d 0a 09 09 09 64 6f 75 62 6c 65 (pR);.....double
08c0: 20 74 68 65 41 3b 0d 0a 09 09 09 69 66 28 61 61 theA;.....if(aa
08d0: 4c 20 3e 20 61 61 52 29 20 7b 0d 0a 09 09 09 09 L > aaR) {......
08e0: 74 68 65 41 20 3d 20 4d 5f 50 49 3b 0d 0a 09 09 theA = M_PI;....
08f0: 09 7d 20 65 6c 73 65 20 7b 0d 0a 09 09 09 09 74 .} else {......t
0900: 68 65 41 20 3d 20 28 61 62 73 28 61 61 4c 29 3e heA = (abs(aaL)>
0910: 61 62 73 28 61 61 52 29 20 3f 20 61 61 4c 20 3a abs(aaR) ? aaL :
0920: 20 61 61 52 29 3b 0d 0a 09 09 09 7d 0d 0a 0d 0a aaR);.....}....
0930: 09 09 09 2f 2f 20 64 69 73 74 61 6e 63 65 20 66 ...// distance f
0940: 72 6f 6d 20 70 6f 6c 61 72 28 31 2e 30 2c 20 74 rom polar(1.0, t
0950: 68 65 41 29 20 74 6f 20 43 69 72 63 28 28 61 62 heA) to Circ((ab
0960: 73 28 76 29 2c 30 29 2c 20 63 73 5b 63 6c 6f 73 s(v),0), cs[clos
0970: 65 5d 2e 72 29 0d 0a 09 09 09 2f 2f 20 2e 2e 2e e].r).....// ...
0980: 20 61 76 6f 69 64 69 6e 67 20 74 68 65 20 74 77 avoiding the tw
0990: 6f 20 63 69 72 63 6c 65 21 0d 0a 09 09 09 64 6f o circle!.....do
09a0: 75 62 6c 65 20 64 64 64 20 3d 20 30 3b 0d 0a 63 uble ddd = 0;..c
09b0: 65 72 72 3c 3c 70 6f 6c 61 72 28 31 2e 30 2c 20 err<<polar(1.0,
09c0: 74 68 65 41 29 3c 3c 65 6e 64 6c 3b 0d 0a 09 09 theA)<<endl;....
09d0: 09 62 65 73 74 20 3d 20 6d 69 6e 28 62 65 73 74 .best = min(best
09e0: 2c 20 64 64 64 29 3b 0d 0a 09 09 7d 0d 0a 09 09 , ddd);....}....
09f0: 72 65 74 75 72 6e 20 62 65 73 74 3b 0d 0a 09 7d return best;...}
0a00: 0d 0a 7d 3b 0d 0a 0d 0a 2f 2f 20 42 45 47 49 4e ..};....// BEGIN
0a10: 20 43 55 54 20 48 45 52 45 0d 0a 23 69 6e 63 6c CUT HERE..#incl
0a20: 75 64 65 20 3c 63 74 69 6d 65 3e 0d 0a 64 6f 75 ude <ctime>..dou
0a30: 62 6c 65 20 73 74 61 72 74 5f 74 69 6d 65 3b 20 ble start_time;
0a40: 73 74 72 69 6e 67 20 74 69 6d 65 72 28 29 0d 0a string timer()..
0a50: 20 7b 20 6f 73 74 72 69 6e 67 73 74 72 65 61 6d { ostringstream
0a60: 20 6f 73 3b 20 6f 73 20 3c 3c 20 22 20 28 22 20 os; os << " ("
0a70: 3c 3c 20 69 6e 74 28 28 63 6c 6f 63 6b 28 29 2d << int((clock()-
0a80: 73 74 61 72 74 5f 74 69 6d 65 29 2f 43 4c 4f 43 start_time)/CLOC
0a90: 4b 53 5f 50 45 52 5f 53 45 43 2a 31 30 30 30 29 KS_PER_SEC*1000)
0aa0: 20 3c 3c 20 22 20 6d 73 65 63 29 22 3b 20 72 65 << " msec)"; re
0ab0: 74 75 72 6e 20 6f 73 2e 73 74 72 28 29 3b 20 7d turn os.str(); }
0ac0: 0d 0a 74 65 6d 70 6c 61 74 65 3c 74 79 70 65 6e ..template<typen
0ad0: 61 6d 65 20 54 3e 20 6f 73 74 72 65 61 6d 26 20 ame T> ostream&
0ae0: 6f 70 65 72 61 74 6f 72 3c 3c 28 6f 73 74 72 65 operator<<(ostre
0af0: 61 6d 26 20 6f 73 2c 20 63 6f 6e 73 74 20 76 65 am& os, const ve
0b00: 63 74 6f 72 3c 54 3e 26 20 76 29 0d 0a 20 7b 20 ctor<T>& v).. {
0b10: 6f 73 20 3c 3c 20 22 7b 20 22 3b 0d 0a 20 20 20 os << "{ ";..
0b20: 66 6f 72 28 74 79 70 65 6e 61 6d 65 20 76 65 63 for(typename vec
0b30: 74 6f 72 3c 54 3e 3a 3a 63 6f 6e 73 74 5f 69 74 tor<T>::const_it
0b40: 65 72 61 74 6f 72 20 69 74 3d 76 2e 62 65 67 69 erator it=v.begi
0b50: 6e 28 29 3b 20 69 74 21 3d 76 2e 65 6e 64 28 29 n(); it!=v.end()
0b60: 3b 20 2b 2b 69 74 29 0d 0a 20 20 20 6f 73 20 3c ; ++it).. os <
0b70: 3c 20 27 5c 22 27 20 3c 3c 20 2a 69 74 20 3c 3c < '\"' << *it <<
0b80: 20 27 5c 22 27 20 3c 3c 20 28 69 74 2b 31 3d 3d '\"' << (it+1==
0b90: 76 2e 65 6e 64 28 29 20 3f 20 22 22 20 3a 20 22 v.end() ? "" : "
0ba0: 2c 20 22 29 3b 20 6f 73 20 3c 3c 20 22 20 7d 22 , "); os << " }"
0bb0: 3b 20 72 65 74 75 72 6e 20 6f 73 3b 20 7d 0d 0a ; return os; }..
0bc0: 76 6f 69 64 20 76 65 72 69 66 79 5f 63 61 73 65 void verify_case
0bd0: 28 63 6f 6e 73 74 20 64 6f 75 62 6c 65 26 20 45 (const double& E
0be0: 78 70 65 63 74 65 64 2c 20 63 6f 6e 73 74 20 64 xpected, const d
0bf0: 6f 75 62 6c 65 26 20 52 65 63 65 69 76 65 64 29 ouble& Received)
0c00: 20 7b 0d 0a 20 62 6f 6f 6c 20 6f 6b 20 3d 20 28 {.. bool ok = (
0c10: 61 62 73 28 45 78 70 65 63 74 65 64 20 2d 20 52 abs(Expected - R
0c20: 65 63 65 69 76 65 64 29 20 3c 20 31 65 2d 39 29 eceived) < 1e-9)
0c30: 3b 0d 0a 20 69 66 28 6f 6b 29 20 63 65 72 72 20 ;.. if(ok) cerr
0c40: 3c 3c 20 22 50 41 53 53 45 44 22 20 3c 3c 20 74 << "PASSED" << t
0c50: 69 6d 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b 20 imer() << endl;
0c60: 20 65 6c 73 65 20 7b 20 63 65 72 72 20 3c 3c 20 else { cerr <<
0c70: 22 46 41 49 4c 45 44 22 20 3c 3c 20 74 69 6d 65 "FAILED" << time
0c80: 72 28 29 20 3c 3c 20 65 6e 64 6c 3b 0d 0a 20 63 r() << endl;.. c
0c90: 65 72 72 20 3c 3c 20 22 5c 74 6f 3a 20 5c 22 22 err << "\to: \""
0ca0: 20 3c 3c 20 45 78 70 65 63 74 65 64 20 3c 3c 20 << Expected <<
0cb0: 27 5c 22 27 20 3c 3c 20 65 6e 64 6c 20 3c 3c 20 '\"' << endl <<
0cc0: 22 5c 74 78 3a 20 5c 22 22 20 3c 3c 20 52 65 63 "\tx: \"" << Rec
0cd0: 65 69 76 65 64 20 3c 3c 20 27 5c 22 27 20 3c 3c eived << '\"' <<
0ce0: 20 65 6e 64 6c 3b 20 7d 20 7d 0d 0a 23 64 65 66 endl; } }..#def
0cf0: 69 6e 65 20 43 41 53 45 28 4e 29 20 7b 63 65 72 ine CASE(N) {cer
0d00: 72 20 3c 3c 20 22 54 65 73 74 20 43 61 73 65 20 r << "Test Case
0d10: 23 22 20 3c 3c 20 4e 20 3c 3c 20 22 2e 2e 2e 22 #" << N << "..."
0d20: 20 3c 3c 20 66 6c 75 73 68 3b 20 73 74 61 72 74 << flush; start
0d30: 5f 74 69 6d 65 3d 63 6c 6f 63 6b 28 29 3b 0d 0a _time=clock();..
0d40: 23 64 65 66 69 6e 65 20 45 4e 44 09 20 76 65 72 #define END. ver
0d50: 69 66 79 5f 63 61 73 65 28 5f 2c 20 43 69 72 63 ify_case(_, Circ
0d60: 75 73 54 65 6e 74 73 28 29 2e 66 69 6e 64 4d 61 usTents().findMa
0d70: 78 69 6d 75 6d 44 69 73 74 61 6e 63 65 28 78 2c ximumDistance(x,
0d80: 20 79 2c 20 72 29 29 3b 7d 0d 0a 69 6e 74 20 6d y, r));}..int m
0d90: 61 69 6e 28 29 7b 0d 0a 0d 0a 43 41 53 45 28 30 ain(){....CASE(0
0da0: 29 0d 0a 09 69 6e 74 20 78 5f 5b 5d 20 3d 20 7b )...int x_[] = {
0db0: 30 2c 33 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 0,3};... vector
0dc0: 20 3c 69 6e 74 3e 20 78 28 78 5f 2c 20 78 5f 2b <int> x(x_, x_+
0dd0: 73 69 7a 65 6f 66 28 78 5f 29 2f 73 69 7a 65 6f sizeof(x_)/sizeo
0de0: 66 28 2a 78 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 f(*x_)); ...int
0df0: 79 5f 5b 5d 20 3d 20 7b 30 2c 30 7d 3b 0d 0a 09 y_[] = {0,0};...
0e00: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79 vector <int> y
0e10: 28 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f 66 28 79 (y_, y_+sizeof(y
0e20: 5f 29 2f 73 69 7a 65 6f 66 28 2a 79 5f 29 29 3b _)/sizeof(*y_));
0e30: 20 0d 0a 09 69 6e 74 20 72 5f 5b 5d 20 3d 20 7b ...int r_[] = {
0e40: 31 2c 31 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 1,1};... vector
0e50: 20 3c 69 6e 74 3e 20 72 28 72 5f 2c 20 72 5f 2b <int> r(r_, r_+
0e60: 73 69 7a 65 6f 66 28 72 5f 29 2f 73 69 7a 65 6f sizeof(r_)/sizeo
0e70: 66 28 2a 72 5f 29 29 3b 20 0d 0a 09 64 6f 75 62 f(*r_)); ...doub
0e80: 6c 65 20 5f 20 3d 20 33 2e 37 33 39 30 36 30 33 le _ = 3.7390603
0e90: 36 30 39 39 35 32 30 37 38 3b 20 0d 0a 45 4e 44 609952078; ..END
0ea0: 0d 0a 43 41 53 45 28 31 29 0d 0a 09 69 6e 74 20 ..CASE(1)...int
0eb0: 78 5f 5b 5d 20 3d 20 7b 30 2c 33 2c 2d 33 2c 33 x_[] = {0,3,-3,3
0ec0: 2c 2d 33 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 ,-3};... vector
0ed0: 20 3c 69 6e 74 3e 20 78 28 78 5f 2c 20 78 5f 2b <int> x(x_, x_+
0ee0: 73 69 7a 65 6f 66 28 78 5f 29 2f 73 69 7a 65 6f sizeof(x_)/sizeo
0ef0: 66 28 2a 78 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 f(*x_)); ...int
0f00: 79 5f 5b 5d 20 3d 20 7b 30 2c 33 2c 33 2c 2d 33 y_[] = {0,3,3,-3
0f10: 2c 2d 33 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 ,-3};... vector
0f20: 20 3c 69 6e 74 3e 20 79 28 79 5f 2c 20 79 5f 2b <int> y(y_, y_+
0f30: 73 69 7a 65 6f 66 28 79 5f 29 2f 73 69 7a 65 6f sizeof(y_)/sizeo
0f40: 66 28 2a 79 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 f(*y_)); ...int
0f50: 72 5f 5b 5d 20 3d 20 7b 31 2c 31 2c 31 2c 31 2c r_[] = {1,1,1,1,
0f60: 31 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 1};... vector <
0f70: 69 6e 74 3e 20 72 28 72 5f 2c 20 72 5f 2b 73 69 int> r(r_, r_+si
0f80: 7a 65 6f 66 28 72 5f 29 2f 73 69 7a 65 6f 66 28 zeof(r_)/sizeof(
0f90: 2a 72 5f 29 29 3b 20 0d 0a 09 64 6f 75 62 6c 65 *r_)); ...double
0fa0: 20 5f 20 3d 20 32 2e 36 30 35 35 35 31 32 37 35 _ = 2.605551275
0fb0: 34 36 33 39 38 38 37 3b 20 0d 0a 45 4e 44 0d 0a 4639887; ..END..
0fc0: 43 41 53 45 28 32 29 0d 0a 09 69 6e 74 20 78 5f CASE(2)...int x_
0fd0: 5b 5d 20 3d 20 7b 33 2c 37 2c 37 2c 37 2c 33 7d [] = {3,7,7,7,3}
0fe0: 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e ;... vector <in
0ff0: 74 3e 20 78 28 78 5f 2c 20 78 5f 2b 73 69 7a 65 t> x(x_, x_+size
1000: 6f 66 28 78 5f 29 2f 73 69 7a 65 6f 66 28 2a 78 of(x_)/sizeof(*x
1010: 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 79 5f 5b 5d _)); ...int y_[]
1020: 20 3d 20 7b 34 2c 36 2c 31 2c 2d 33 2c 30 7d 3b = {4,6,1,-3,0};
1030: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 ... vector <int
1040: 3e 20 79 28 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f > y(y_, y_+sizeo
1050: 66 28 79 5f 29 2f 73 69 7a 65 6f 66 28 2a 79 5f f(y_)/sizeof(*y_
1060: 29 29 3b 20 0d 0a 09 69 6e 74 20 72 5f 5b 5d 20 )); ...int r_[]
1070: 3d 20 7b 32 2c 32 2c 32 2c 31 2c 31 7d 3b 0d 0a = {2,2,2,1,1};..
1080: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 . vector <int>
1090: 72 28 72 5f 2c 20 72 5f 2b 73 69 7a 65 6f 66 28 r(r_, r_+sizeof(
10a0: 72 5f 29 2f 73 69 7a 65 6f 66 28 2a 72 5f 29 29 r_)/sizeof(*r_))
10b0: 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20 ; ...double _ =
10c0: 34 2e 33 32 36 34 34 35 39 30 39 39 36 32 30 37 4.32644590996207
10d0: 32 35 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 25; ..END..CASE(
10e0: 33 29 0d 0a 09 69 6e 74 20 78 5f 5b 5d 20 3d 20 3)...int x_[] =
10f0: 7b 31 30 2c 2d 31 7d 3b 0d 0a 09 20 20 76 65 63 {10,-1};... vec
1100: 74 6f 72 20 3c 69 6e 74 3e 20 78 28 78 5f 2c 20 tor <int> x(x_,
1110: 78 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 2f 73 69 x_+sizeof(x_)/si
1120: 7a 65 6f 66 28 2a 78 5f 29 29 3b 20 0d 0a 09 69 zeof(*x_)); ...i
1130: 6e 74 20 79 5f 5b 5d 20 3d 20 7b 30 2c 30 7d 3b nt y_[] = {0,0};
1140: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 ... vector <int
1150: 3e 20 79 28 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f > y(y_, y_+sizeo
1160: 66 28 79 5f 29 2f 73 69 7a 65 6f 66 28 2a 79 5f f(y_)/sizeof(*y_
1170: 29 29 3b 20 0d 0a 09 69 6e 74 20 72 5f 5b 5d 20 )); ...int r_[]
1180: 3d 20 7b 38 2c 32 7d 3b 0d 0a 09 20 20 76 65 63 = {8,2};... vec
1190: 74 6f 72 20 3c 69 6e 74 3e 20 72 28 72 5f 2c 20 tor <int> r(r_,
11a0: 72 5f 2b 73 69 7a 65 6f 66 28 72 5f 29 2f 73 69 r_+sizeof(r_)/si
11b0: 7a 65 6f 66 28 2a 72 5f 29 29 3b 20 0d 0a 09 64 zeof(*r_)); ...d
11c0: 6f 75 62 6c 65 20 5f 20 3d 20 32 34 2e 36 33 30 ouble _ = 24.630
11d0: 39 32 34 35 38 36 36 34 32 31 32 3b 20 0d 0a 45 92458664212; ..E
11e0: 4e 44 0d 0a 43 41 53 45 28 34 29 0d 0a 09 69 6e ND..CASE(4)...in
11f0: 74 20 78 5f 5b 5d 20 3d 20 7b 30 2c 34 2c 2d 34 t x_[] = {0,4,-4
1200: 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 };... vector <i
1210: 6e 74 3e 20 78 28 78 5f 2c 20 78 5f 2b 73 69 7a nt> x(x_, x_+siz
1220: 65 6f 66 28 78 5f 29 2f 73 69 7a 65 6f 66 28 2a eof(x_)/sizeof(*
1230: 78 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 79 5f 5b x_)); ...int y_[
1240: 5d 20 3d 20 7b 30 2c 34 2c 2d 34 7d 3b 0d 0a 09 ] = {0,4,-4};...
1250: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79 vector <int> y
1260: 28 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f 66 28 79 (y_, y_+sizeof(y
1270: 5f 29 2f 73 69 7a 65 6f 66 28 2a 79 5f 29 29 3b _)/sizeof(*y_));
1280: 20 0d 0a 09 69 6e 74 20 72 5f 5b 5d 20 3d 20 7b ...int r_[] = {
1290: 31 2c 31 2c 31 7d 3b 0d 0a 09 20 20 76 65 63 74 1,1,1};... vect
12a0: 6f 72 20 3c 69 6e 74 3e 20 72 28 72 5f 2c 20 72 or <int> r(r_, r
12b0: 5f 2b 73 69 7a 65 6f 66 28 72 5f 29 2f 73 69 7a _+sizeof(r_)/siz
12c0: 65 6f 66 28 2a 72 5f 29 29 3b 20 0d 0a 09 64 6f eof(*r_)); ...do
12d0: 75 62 6c 65 20 5f 20 3d 20 34 2e 37 34 35 34 37 uble _ = 4.74547
12e0: 34 39 36 33 36 37 35 31 33 33 3b 20 0d 0a 45 4e 4963675133; ..EN
12f0: 44 0d 0a 2f 2a 0d 0a 43 41 53 45 28 35 29 0d 0a D../*..CASE(5)..
1300: 09 69 6e 74 20 78 5f 5b 5d 20 3d 20 3b 0d 0a 09 .int x_[] = ;...
1310: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 78 vector <int> x
1320: 28 78 5f 2c 20 78 5f 2b 73 69 7a 65 6f 66 28 78 (x_, x_+sizeof(x
1330: 5f 29 2f 73 69 7a 65 6f 66 28 2a 78 5f 29 29 3b _)/sizeof(*x_));
1340: 20 0d 0a 09 69 6e 74 20 79 5f 5b 5d 20 3d 20 3b ...int y_[] = ;
1350: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 ... vector <int
1360: 3e 20 79 28 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f > y(y_, y_+sizeo
1370: 66 28 79 5f 29 2f 73 69 7a 65 6f 66 28 2a 79 5f f(y_)/sizeof(*y_
1380: 29 29 3b 20 0d 0a 09 69 6e 74 20 72 5f 5b 5d 20 )); ...int r_[]
1390: 3d 20 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c = ;... vector <
13a0: 69 6e 74 3e 20 72 28 72 5f 2c 20 72 5f 2b 73 69 int> r(r_, r_+si
13b0: 7a 65 6f 66 28 72 5f 29 2f 73 69 7a 65 6f 66 28 zeof(r_)/sizeof(
13c0: 2a 72 5f 29 29 3b 20 0d 0a 09 64 6f 75 62 6c 65 *r_)); ...double
13d0: 20 5f 20 3d 20 3b 20 0d 0a 45 4e 44 0d 0a 43 41 _ = ; ..END..CA
13e0: 53 45 28 36 29 0d 0a 09 69 6e 74 20 78 5f 5b 5d SE(6)...int x_[]
13f0: 20 3d 20 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 = ;... vector
1400: 3c 69 6e 74 3e 20 78 28 78 5f 2c 20 78 5f 2b 73 <int> x(x_, x_+s
1410: 69 7a 65 6f 66 28 78 5f 29 2f 73 69 7a 65 6f 66 izeof(x_)/sizeof
1420: 28 2a 78 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 79 (*x_)); ...int y
1430: 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 65 63 74 _[] = ;... vect
1440: 6f 72 20 3c 69 6e 74 3e 20 79 28 79 5f 2c 20 79 or <int> y(y_, y
1450: 5f 2b 73 69 7a 65 6f 66 28 79 5f 29 2f 73 69 7a _+sizeof(y_)/siz
1460: 65 6f 66 28 2a 79 5f 29 29 3b 20 0d 0a 09 69 6e eof(*y_)); ...in
1470: 74 20 72 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 t r_[] = ;... v
1480: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 72 28 72 5f ector <int> r(r_
1490: 2c 20 72 5f 2b 73 69 7a 65 6f 66 28 72 5f 29 2f , r_+sizeof(r_)/
14a0: 73 69 7a 65 6f 66 28 2a 72 5f 29 29 3b 20 0d 0a sizeof(*r_)); ..
14b0: 09 64 6f 75 62 6c 65 20 5f 20 3d 20 3b 20 0d 0a .double _ = ; ..
14c0: 45 4e 44 0d 0a 2a 2f 0d 0a 7d 0d 0a 2f 2f 20 45 END..*/..}..// E
14d0: 4e 44 20 43 55 54 20 48 45 52 45 0d 0a ND CUT HERE..