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 54 68 CMP;....class Th
01a0: 65 4e 65 77 48 6f 75 73 65 44 69 76 4f 6e 65 20 eNewHouseDivOne
01b0: 7b 20 70 75 62 6c 69 63 3a 0d 0a 09 64 6f 75 62 { public:...doub
01c0: 6c 65 20 64 69 73 74 61 6e 63 65 28 76 65 63 74 le distance(vect
01d0: 6f 72 20 3c 69 6e 74 3e 20 78 2c 20 76 65 63 74 or <int> x, vect
01e0: 6f 72 20 3c 69 6e 74 3e 20 79 2c 20 69 6e 74 20 or <int> y, int
01f0: 6b 29 20 0d 0a 09 7b 0d 0a 09 09 64 6f 75 62 6c k) ...{....doubl
0200: 65 20 6d 69 6e 44 20 3d 20 31 65 2b 39 39 3b 0d e minD = 1e+99;.
0210: 0a 09 09 66 6f 72 28 64 6f 75 62 6c 65 20 58 3d ...for(double X=
0220: 2d 35 30 3b 20 58 3c 3d 2b 35 30 3b 20 58 2b 3d -50; X<=+50; X+=
0230: 30 2e 35 29 0d 0a 09 09 09 66 6f 72 28 64 6f 75 0.5).....for(dou
0240: 62 6c 65 20 59 3d 2d 35 30 3b 20 59 3c 3d 2b 35 ble Y=-50; Y<=+5
0250: 30 3b 20 59 2b 3d 30 2e 35 29 0d 0a 09 09 09 7b 0; Y+=0.5).....{
0260: 0d 0a 09 09 09 09 76 65 63 74 6f 72 3c 64 6f 75 ......vector<dou
0270: 62 6c 65 3e 20 64 3b 0d 0a 09 09 09 09 66 6f 72 ble> d;......for
0280: 28 69 6e 74 20 69 3d 30 3b 20 69 3c 78 2e 73 69 (int i=0; i<x.si
0290: 7a 65 28 29 3b 20 2b 2b 69 29 0d 0a 09 09 09 09 ze(); ++i)......
02a0: 09 64 2e 70 75 73 68 5f 62 61 63 6b 28 20 61 62 .d.push_back( ab
02b0: 73 28 78 5b 69 5d 2d 58 29 20 2b 20 61 62 73 28 s(x[i]-X) + abs(
02c0: 79 5b 69 5d 2d 59 29 20 29 3b 0d 0a 09 09 09 09 y[i]-Y) );......
02d0: 6e 74 68 5f 65 6c 65 6d 65 6e 74 28 64 2e 62 65 nth_element(d.be
02e0: 67 69 6e 28 29 2c 20 64 2e 62 65 67 69 6e 28 29 gin(), d.begin()
02f0: 2b 6b 2d 31 2c 20 64 2e 65 6e 64 28 29 29 3b 0d +k-1, d.end());.
0300: 0a 09 09 09 09 6d 69 6e 44 20 3d 20 6d 69 6e 28 .....minD = min(
0310: 6d 69 6e 44 2c 20 64 5b 6b 2d 31 5d 29 3b 0d 0a minD, d[k-1]);..
0320: 09 09 09 7d 0d 0a 09 09 72 65 74 75 72 6e 20 6d ...}....return m
0330: 69 6e 44 3b 0d 0a 09 7d 0d 0a 7d 3b 0d 0a 0d 0a inD;...}..};....
0340: 2f 2f 20 42 45 47 49 4e 20 43 55 54 20 48 45 52 // BEGIN CUT HER
0350: 45 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 63 74 69 E..#include <cti
0360: 6d 65 3e 0d 0a 64 6f 75 62 6c 65 20 73 74 61 72 me>..double star
0370: 74 5f 74 69 6d 65 3b 20 73 74 72 69 6e 67 20 74 t_time; string t
0380: 69 6d 65 72 28 29 0d 0a 20 7b 20 6f 73 74 72 69 imer().. { ostri
0390: 6e 67 73 74 72 65 61 6d 20 6f 73 3b 20 6f 73 20 ngstream os; os
03a0: 3c 3c 20 22 20 28 22 20 3c 3c 20 69 6e 74 28 28 << " (" << int((
03b0: 63 6c 6f 63 6b 28 29 2d 73 74 61 72 74 5f 74 69 clock()-start_ti
03c0: 6d 65 29 2f 43 4c 4f 43 4b 53 5f 50 45 52 5f 53 me)/CLOCKS_PER_S
03d0: 45 43 2a 31 30 30 30 29 20 3c 3c 20 22 20 6d 73 EC*1000) << " ms
03e0: 65 63 29 22 3b 20 72 65 74 75 72 6e 20 6f 73 2e ec)"; return os.
03f0: 73 74 72 28 29 3b 20 7d 0d 0a 74 65 6d 70 6c 61 str(); }..templa
0400: 74 65 3c 74 79 70 65 6e 61 6d 65 20 54 3e 20 6f te<typename T> o
0410: 73 74 72 65 61 6d 26 20 6f 70 65 72 61 74 6f 72 stream& operator
0420: 3c 3c 28 6f 73 74 72 65 61 6d 26 20 6f 73 2c 20 <<(ostream& os,
0430: 63 6f 6e 73 74 20 76 65 63 74 6f 72 3c 54 3e 26 const vector<T>&
0440: 20 76 29 0d 0a 20 7b 20 6f 73 20 3c 3c 20 22 7b v).. { os << "{
0450: 20 22 3b 0d 0a 20 20 20 66 6f 72 28 74 79 70 65 ";.. for(type
0460: 6e 61 6d 65 20 76 65 63 74 6f 72 3c 54 3e 3a 3a name vector<T>::
0470: 63 6f 6e 73 74 5f 69 74 65 72 61 74 6f 72 20 69 const_iterator i
0480: 74 3d 76 2e 62 65 67 69 6e 28 29 3b 20 69 74 21 t=v.begin(); it!
0490: 3d 76 2e 65 6e 64 28 29 3b 20 2b 2b 69 74 29 0d =v.end(); ++it).
04a0: 0a 20 20 20 6f 73 20 3c 3c 20 27 5c 22 27 20 3c . os << '\"' <
04b0: 3c 20 2a 69 74 20 3c 3c 20 27 5c 22 27 20 3c 3c < *it << '\"' <<
04c0: 20 28 69 74 2b 31 3d 3d 76 2e 65 6e 64 28 29 20 (it+1==v.end()
04d0: 3f 20 22 22 20 3a 20 22 2c 20 22 29 3b 20 6f 73 ? "" : ", "); os
04e0: 20 3c 3c 20 22 20 7d 22 3b 20 72 65 74 75 72 6e << " }"; return
04f0: 20 6f 73 3b 20 7d 0d 0a 76 6f 69 64 20 76 65 72 os; }..void ver
0500: 69 66 79 5f 63 61 73 65 28 63 6f 6e 73 74 20 64 ify_case(const d
0510: 6f 75 62 6c 65 26 20 45 78 70 65 63 74 65 64 2c ouble& Expected,
0520: 20 63 6f 6e 73 74 20 64 6f 75 62 6c 65 26 20 52 const double& R
0530: 65 63 65 69 76 65 64 29 20 7b 0d 0a 20 64 6f 75 eceived) {.. dou
0540: 62 6c 65 20 64 69 66 66 20 3d 20 45 78 70 65 63 ble diff = Expec
0550: 74 65 64 20 2d 20 52 65 63 65 69 76 65 64 3b 20 ted - Received;
0560: 69 66 20 28 64 69 66 66 20 3c 20 30 29 20 64 69 if (diff < 0) di
0570: 66 66 20 3d 20 2d 64 69 66 66 3b 20 62 6f 6f 6c ff = -diff; bool
0580: 20 6f 6b 20 3d 20 28 64 69 66 66 20 3c 20 31 65 ok = (diff < 1e
0590: 2d 39 29 3b 0d 0a 20 69 66 28 6f 6b 29 20 63 65 -9);.. if(ok) ce
05a0: 72 72 20 3c 3c 20 22 50 41 53 53 45 44 22 20 3c rr << "PASSED" <
05b0: 3c 20 74 69 6d 65 72 28 29 20 3c 3c 20 65 6e 64 < timer() << end
05c0: 6c 3b 20 20 65 6c 73 65 20 7b 20 63 65 72 72 20 l; else { cerr
05d0: 3c 3c 20 22 46 41 49 4c 45 44 22 20 3c 3c 20 74 << "FAILED" << t
05e0: 69 6d 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b 0d imer() << endl;.
05f0: 0a 20 63 65 72 72 20 3c 3c 20 22 5c 74 6f 3a 20 . cerr << "\to:
0600: 5c 22 22 20 3c 3c 20 45 78 70 65 63 74 65 64 20 \"" << Expected
0610: 3c 3c 20 27 5c 22 27 20 3c 3c 20 65 6e 64 6c 20 << '\"' << endl
0620: 3c 3c 20 22 5c 74 78 3a 20 5c 22 22 20 3c 3c 20 << "\tx: \"" <<
0630: 52 65 63 65 69 76 65 64 20 3c 3c 20 27 5c 22 27 Received << '\"'
0640: 20 3c 3c 20 65 6e 64 6c 3b 20 7d 20 7d 0d 0a 23 << endl; } }..#
0650: 64 65 66 69 6e 65 20 43 41 53 45 28 4e 29 20 7b define CASE(N) {
0660: 63 65 72 72 20 3c 3c 20 22 54 65 73 74 20 43 61 cerr << "Test Ca
0670: 73 65 20 23 22 20 3c 3c 20 4e 20 3c 3c 20 22 2e se #" << N << ".
0680: 2e 2e 22 20 3c 3c 20 66 6c 75 73 68 3b 20 73 74 .." << flush; st
0690: 61 72 74 5f 74 69 6d 65 3d 63 6c 6f 63 6b 28 29 art_time=clock()
06a0: 3b 0d 0a 23 64 65 66 69 6e 65 20 45 4e 44 09 20 ;..#define END.
06b0: 76 65 72 69 66 79 5f 63 61 73 65 28 5f 2c 20 54 verify_case(_, T
06c0: 68 65 4e 65 77 48 6f 75 73 65 44 69 76 4f 6e 65 heNewHouseDivOne
06d0: 28 29 2e 64 69 73 74 61 6e 63 65 28 78 2c 20 79 ().distance(x, y
06e0: 2c 20 6b 29 29 3b 7d 0d 0a 69 6e 74 20 6d 61 69 , k));}..int mai
06f0: 6e 28 29 7b 0d 0a 0d 0a 43 41 53 45 28 30 29 0d n(){....CASE(0).
0700: 0a 09 69 6e 74 20 78 5f 5b 5d 20 3d 20 7b 2d 31 ..int x_[] = {-1
0710: 2c 20 2d 31 2c 20 31 2c 20 31 7d 3b 0d 0a 09 20 , -1, 1, 1};...
0720: 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 78 28 vector <int> x(
0730: 78 5f 2c 20 78 5f 2b 73 69 7a 65 6f 66 28 78 5f x_, x_+sizeof(x_
0740: 29 2f 73 69 7a 65 6f 66 28 2a 78 5f 29 29 3b 20 )/sizeof(*x_));
0750: 0d 0a 09 69 6e 74 20 79 5f 5b 5d 20 3d 20 7b 2d ...int y_[] = {-
0760: 31 2c 20 31 2c 20 2d 31 2c 20 31 7d 3b 0d 0a 09 1, 1, -1, 1};...
0770: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79 vector <int> y
0780: 28 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f 66 28 79 (y_, y_+sizeof(y
0790: 5f 29 2f 73 69 7a 65 6f 66 28 2a 79 5f 29 29 3b _)/sizeof(*y_));
07a0: 20 0d 0a 09 69 6e 74 20 6b 20 3d 20 33 3b 20 0d ...int k = 3; .
07b0: 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20 32 2e 30 ..double _ = 2.0
07c0: 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 31 29 ; ..END..CASE(1)
07d0: 0d 0a 09 69 6e 74 20 78 5f 5b 5d 20 3d 20 7b 2d ...int x_[] = {-
07e0: 31 2c 20 2d 31 2c 20 31 2c 20 31 7d 3b 0d 0a 09 1, -1, 1, 1};...
07f0: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 78 vector <int> x
0800: 28 78 5f 2c 20 78 5f 2b 73 69 7a 65 6f 66 28 78 (x_, x_+sizeof(x
0810: 5f 29 2f 73 69 7a 65 6f 66 28 2a 78 5f 29 29 3b _)/sizeof(*x_));
0820: 20 0d 0a 09 69 6e 74 20 79 5f 5b 5d 20 3d 20 7b ...int y_[] = {
0830: 2d 31 2c 20 31 2c 20 2d 31 2c 20 31 7d 3b 0d 0a -1, 1, -1, 1};..
0840: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 . vector <int>
0850: 79 28 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f 66 28 y(y_, y_+sizeof(
0860: 79 5f 29 2f 73 69 7a 65 6f 66 28 2a 79 5f 29 29 y_)/sizeof(*y_))
0870: 3b 20 0d 0a 09 69 6e 74 20 6b 20 3d 20 32 3b 20 ; ...int k = 2;
0880: 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20 31 2e ...double _ = 1.
0890: 30 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 32 0; ..END..CASE(2
08a0: 29 0d 0a 09 69 6e 74 20 78 5f 5b 5d 20 3d 20 7b )...int x_[] = {
08b0: 34 2c 20 34 2c 20 34 2c 20 34 2c 20 34 2c 20 33 4, 4, 4, 4, 4, 3
08c0: 2c 20 33 2c 20 35 2c 20 35 7d 3b 0d 0a 09 20 20 , 3, 5, 5};...
08d0: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 78 28 78 vector <int> x(x
08e0: 5f 2c 20 78 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 _, x_+sizeof(x_)
08f0: 2f 73 69 7a 65 6f 66 28 2a 78 5f 29 29 3b 20 0d /sizeof(*x_)); .
0900: 0a 09 69 6e 74 20 79 5f 5b 5d 20 3d 20 7b 37 2c ..int y_[] = {7,
0910: 20 37 2c 20 37 2c 20 34 2c 20 34 2c 20 35 2c 20 7, 7, 4, 4, 5,
0920: 36 2c 20 35 2c 20 36 7d 3b 0d 0a 09 20 20 76 65 6, 5, 6};... ve
0930: 63 74 6f 72 20 3c 69 6e 74 3e 20 79 28 79 5f 2c ctor <int> y(y_,
0940: 20 79 5f 2b 73 69 7a 65 6f 66 28 79 5f 29 2f 73 y_+sizeof(y_)/s
0950: 69 7a 65 6f 66 28 2a 79 5f 29 29 3b 20 0d 0a 09 izeof(*y_)); ...
0960: 69 6e 74 20 6b 20 3d 20 39 3b 20 0d 0a 09 64 6f int k = 9; ...do
0970: 75 62 6c 65 20 5f 20 3d 20 31 2e 35 3b 20 0d 0a uble _ = 1.5; ..
0980: 45 4e 44 0d 0a 43 41 53 45 28 33 29 0d 0a 09 69 END..CASE(3)...i
0990: 6e 74 20 78 5f 5b 5d 20 3d 20 7b 33 30 2c 20 2d nt x_[] = {30, -
09a0: 31 35 2c 20 32 34 2c 20 2d 32 33 2c 20 34 33 2c 15, 24, -23, 43,
09b0: 20 2d 38 2c 20 2d 36 2c 20 2d 34 37 2c 20 32 33 -8, -6, -47, 23
09c0: 2c 20 2d 31 31 2c 20 34 33 2c 20 36 2c 20 2d 31 , -11, 43, 6, -1
09d0: 38 2c 20 34 34 2c 20 2d 34 36 2c 20 31 36 2c 20 8, 44, -46, 16,
09e0: 33 32 2c 20 33 31 2c 20 31 33 2c 20 39 2c 20 32 32, 31, 13, 9, 2
09f0: 32 2c 20 32 35 2c 20 34 2c 20 2d 34 34 2c 20 33 2, 25, 4, -44, 3
0a00: 38 2c 20 2d 34 31 2c 20 2d 32 30 2c 20 34 31 2c 8, -41, -20, 41,
0a10: 20 2d 33 35 2c 20 37 2c 20 32 35 2c 20 33 39 2c -35, 7, 25, 39,
0a20: 20 2d 33 36 2c 20 2d 32 30 2c 20 2d 35 2c 20 2d -36, -20, -5, -
0a30: 33 38 2c 20 2d 31 35 2c 20 2d 32 32 2c 20 30 7d 38, -15, -22, 0}
0a40: 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e ;... vector <in
0a50: 74 3e 20 78 28 78 5f 2c 20 78 5f 2b 73 69 7a 65 t> x(x_, x_+size
0a60: 6f 66 28 78 5f 29 2f 73 69 7a 65 6f 66 28 2a 78 of(x_)/sizeof(*x
0a70: 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 79 5f 5b 5d _)); ...int y_[]
0a80: 20 3d 20 7b 2d 34 35 2c 20 2d 37 2c 20 2d 33 33 = {-45, -7, -33
0a90: 2c 20 33 31 2c 20 2d 38 2c 20 2d 33 33 2c 20 2d , 31, -8, -33, -
0aa0: 32 30 2c 20 2d 31 34 2c 20 2d 35 30 2c 20 2d 34 20, -14, -50, -4
0ab0: 38 2c 20 2d 33 31 2c 20 33 35 2c 20 2d 32 34 2c 8, -31, 35, -24,
0ac0: 20 2d 33 31 2c 20 2d 31 31 2c 20 34 31 2c 20 2d -31, -11, 41, -
0ad0: 34 31 2c 20 2d 31 31 2c 20 34 36 2c 20 2d 31 2c 41, -11, 46, -1,
0ae0: 20 2d 35 2c 20 35 2c 20 2d 33 39 2c 20 2d 32 36 -5, 5, -39, -26
0af0: 2c 20 2d 34 30 2c 20 2d 39 2c 20 31 36 2c 20 33 , -40, -9, 16, 3
0b00: 38 2c 20 2d 33 30 2c 20 33 34 2c 20 34 36 2c 20 8, -30, 34, 46,
0b10: 2d 31 37 2c 20 2d 32 37 2c 20 33 33 2c 20 2d 33 -17, -27, 33, -3
0b20: 38 2c 20 32 38 2c 20 34 36 2c 20 2d 31 36 2c 20 8, 28, 46, -16,
0b30: 2d 34 36 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 -46};... vector
0b40: 20 3c 69 6e 74 3e 20 79 28 79 5f 2c 20 79 5f 2b <int> y(y_, y_+
0b50: 73 69 7a 65 6f 66 28 79 5f 29 2f 73 69 7a 65 6f sizeof(y_)/sizeo
0b60: 66 28 2a 79 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 f(*y_)); ...int
0b70: 6b 20 3d 20 31 33 3b 20 0d 0a 09 64 6f 75 62 6c k = 13; ...doubl
0b80: 65 20 5f 20 3d 20 33 32 2e 30 3b 20 0d 0a 45 4e e _ = 32.0; ..EN
0b90: 44 0d 0a 0d 0a 7d 0d 0a 2f 2f 20 45 4e 44 20 43 D....}..// END C
0ba0: 55 54 20 48 45 52 45 0d 0a UT HERE..