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 74 75 70 6c 65 3e 0d 0a 75 73 69 6e 67 20 6e <tuple>..using n
0150: 61 6d 65 73 70 61 63 65 20 73 74 64 3b 0d 0a 74 amespace std;..t
0160: 79 70 65 64 65 66 20 6c 6f 6e 67 20 6c 6f 6e 67 ypedef long long
0170: 20 4c 4c 3b 0d 0a 74 79 70 65 64 65 66 20 63 6f LL;..typedef co
0180: 6d 70 6c 65 78 3c 64 6f 75 62 6c 65 3e 20 43 4d mplex<double> CM
0190: 50 3b 0d 0a 0d 0a 63 6c 61 73 73 20 54 68 65 47 P;....class TheG
01a0: 72 69 64 44 69 76 4f 6e 65 20 7b 20 70 75 62 6c ridDivOne { publ
01b0: 69 63 3a 0d 0a 09 69 6e 74 20 66 69 6e 64 28 76 ic:...int find(v
01c0: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 78 2c 20 76 ector <int> x, v
01d0: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79 2c 20 69 ector <int> y, i
01e0: 6e 74 20 6b 29 0d 0a 09 7b 0d 0a 09 09 65 6e 75 nt k)...{....enu
01f0: 6d 20 7b 58 2c 20 59 7d 3b 0d 0a 09 09 74 79 70 m {X, Y};....typ
0200: 65 64 65 66 20 74 75 70 6c 65 3c 69 6e 74 2c 69 edef tuple<int,i
0210: 6e 74 3e 20 50 6f 69 6e 74 3b 0d 0a 0d 0a 09 09 nt> Point;......
0220: 73 65 74 3c 50 6f 69 6e 74 3e 20 62 6c 6f 63 6b set<Point> block
0230: 65 64 3b 0d 0a 09 09 6d 61 70 3c 69 6e 74 2c 20 ed;....map<int,
0240: 73 65 74 3c 69 6e 74 3e 3e 20 78 32 79 3b 0d 0a set<int>> x2y;..
0250: 09 09 78 32 79 5b 30 5d 2e 65 6d 70 6c 61 63 65 ..x2y[0].emplace
0260: 28 30 29 3b 0d 0a 09 09 66 6f 72 28 69 6e 74 20 (0);....for(int
0270: 69 3d 30 3b 20 69 3c 78 2e 73 69 7a 65 28 29 3b i=0; i<x.size();
0280: 20 2b 2b 69 29 20 7b 0d 0a 09 09 09 62 6c 6f 63 ++i) {.....bloc
0290: 6b 65 64 2e 65 6d 70 6c 61 63 65 28 78 5b 69 5d ked.emplace(x[i]
02a0: 2c 20 79 5b 69 5d 29 3b 0d 0a 09 09 09 78 32 79 , y[i]);.....x2y
02b0: 5b 78 5b 69 5d 2d 31 5d 2e 65 6d 70 6c 61 63 65 [x[i]-1].emplace
02c0: 28 79 5b 69 5d 2d 31 29 3b 0d 0a 09 09 09 78 32 (y[i]-1);.....x2
02d0: 79 5b 78 5b 69 5d 2d 31 5d 2e 65 6d 70 6c 61 63 y[x[i]-1].emplac
02e0: 65 28 79 5b 69 5d 2b 30 29 3b 0d 0a 09 09 09 78 e(y[i]+0);.....x
02f0: 32 79 5b 78 5b 69 5d 2d 31 5d 2e 65 6d 70 6c 61 2y[x[i]-1].empla
0300: 63 65 28 79 5b 69 5d 2b 31 29 3b 0d 0a 09 09 09 ce(y[i]+1);.....
0310: 78 32 79 5b 78 5b 69 5d 2b 30 5d 2e 65 6d 70 6c x2y[x[i]+0].empl
0320: 61 63 65 28 79 5b 69 5d 2d 31 29 3b 0d 0a 09 09 ace(y[i]-1);....
0330: 09 78 32 79 5b 78 5b 69 5d 2b 30 5d 2e 65 6d 70 .x2y[x[i]+0].emp
0340: 6c 61 63 65 28 79 5b 69 5d 2b 31 29 3b 0d 0a 09 lace(y[i]+1);...
0350: 09 09 78 32 79 5b 78 5b 69 5d 2b 31 5d 2e 65 6d ..x2y[x[i]+1].em
0360: 70 6c 61 63 65 28 79 5b 69 5d 2d 31 29 3b 0d 0a place(y[i]-1);..
0370: 09 09 09 78 32 79 5b 78 5b 69 5d 2b 31 5d 2e 65 ...x2y[x[i]+1].e
0380: 6d 70 6c 61 63 65 28 79 5b 69 5d 2b 30 29 3b 0d mplace(y[i]+0);.
0390: 0a 09 09 09 78 32 79 5b 78 5b 69 5d 2b 31 5d 2e ....x2y[x[i]+1].
03a0: 65 6d 70 6c 61 63 65 28 79 5b 69 5d 2b 31 29 3b emplace(y[i]+1);
03b0: 0d 0a 09 09 7d 0d 0a 0d 0a 09 09 6d 61 70 3c 69 ....}......map<i
03c0: 6e 74 2c 20 73 65 74 3c 69 6e 74 3e 3e 20 78 32 nt, set<int>> x2
03d0: 79 5f 61 75 67 20 3d 20 78 32 79 3b 0d 0a 09 09 y_aug = x2y;....
03e0: 66 6f 72 28 61 75 74 6f 20 69 74 3d 78 32 79 2e for(auto it=x2y.
03f0: 62 65 67 69 6e 28 29 3b 20 69 74 21 3d 78 32 79 begin(); it!=x2y
0400: 2e 65 6e 64 28 29 3b 20 2b 2b 69 74 29 20 7b 0d .end(); ++it) {.
0410: 0a 09 09 09 61 75 74 6f 20 6b 74 20 3d 20 69 74 ....auto kt = it
0420: 3b 20 2b 2b 6b 74 3b 0d 0a 09 09 09 69 66 28 6b ; ++kt;.....if(k
0430: 74 20 21 3d 20 78 32 79 2e 65 6e 64 28 29 29 20 t != x2y.end())
0440: 7b 0d 0a 09 09 09 09 66 6f 72 28 69 6e 74 20 79 {......for(int y
0450: 3a 20 69 74 2d 3e 73 65 63 6f 6e 64 29 0d 0a 09 : it->second)...
0460: 09 09 09 09 78 32 79 5f 61 75 67 5b 6b 74 2d 3e ....x2y_aug[kt->
0470: 66 69 72 73 74 5d 2e 65 6d 70 6c 61 63 65 28 79 first].emplace(y
0480: 29 3b 0d 0a 09 09 09 7d 0d 0a 09 09 09 69 66 28 );.....}.....if(
0490: 69 74 20 21 3d 20 78 32 79 2e 62 65 67 69 6e 28 it != x2y.begin(
04a0: 29 29 20 7b 0d 0a 09 09 09 09 61 75 74 6f 20 6b )) {......auto k
04b0: 74 20 3d 20 69 74 3b 20 2d 2d 6b 74 3b 0d 0a 09 t = it; --kt;...
04c0: 09 09 09 66 6f 72 28 69 6e 74 20 79 3a 20 69 74 ...for(int y: it
04d0: 2d 3e 73 65 63 6f 6e 64 29 0d 0a 09 09 09 09 09 ->second).......
04e0: 78 32 79 5f 61 75 67 5b 6b 74 2d 3e 66 69 72 73 x2y_aug[kt->firs
04f0: 74 5d 2e 65 6d 70 6c 61 63 65 28 79 29 3b 0d 0a t].emplace(y);..
0500: 09 09 09 7d 0d 0a 09 09 7d 0d 0a 0d 0a 09 09 76 ...}....}......v
0510: 65 63 74 6f 72 3c 50 6f 69 6e 74 3e 20 70 74 73 ector<Point> pts
0520: 3b 0d 0a 09 09 66 6f 72 28 61 75 74 6f 26 20 78 ;....for(auto& x
0530: 79 73 3a 20 78 32 79 5f 61 75 67 29 20 7b 0d 0a ys: x2y_aug) {..
0540: 09 09 09 69 6e 74 20 78 20 3d 20 78 79 73 2e 66 ...int x = xys.f
0550: 69 72 73 74 3b 0d 0a 09 09 09 66 6f 72 28 69 6e irst;.....for(in
0560: 74 20 79 3a 20 78 79 73 2e 73 65 63 6f 6e 64 29 t y: xys.second)
0570: 20 7b 0d 0a 09 09 09 09 50 6f 69 6e 74 20 71 28 {......Point q(
0580: 78 2c 79 29 3b 0d 0a 09 09 09 09 69 66 28 21 62 x,y);......if(!b
0590: 6c 6f 63 6b 65 64 2e 63 6f 75 6e 74 28 71 29 29 locked.count(q))
05a0: 0d 0a 09 09 09 09 09 70 74 73 2e 65 6d 70 6c 61 .......pts.empla
05b0: 63 65 5f 62 61 63 6b 28 71 29 3b 0d 0a 09 09 09 ce_back(q);.....
05c0: 7d 0d 0a 09 09 7d 0d 0a 0d 0a 09 09 74 79 70 65 }....}......type
05d0: 64 65 66 20 69 6e 74 20 56 65 72 74 3b 0d 0a 09 def int Vert;...
05e0: 09 63 6f 6e 73 74 20 69 6e 74 20 4e 20 3d 20 70 .const int N = p
05f0: 74 73 2e 73 69 7a 65 28 29 3b 0d 0a 09 09 63 6f ts.size();....co
0600: 6e 73 74 20 69 6e 74 20 53 20 3d 20 73 74 64 3a nst int S = std:
0610: 3a 66 69 6e 64 28 70 74 73 2e 62 65 67 69 6e 28 :find(pts.begin(
0620: 29 2c 20 70 74 73 2e 65 6e 64 28 29 2c 20 50 6f ), pts.end(), Po
0630: 69 6e 74 28 30 2c 30 29 29 20 2d 20 70 74 73 2e int(0,0)) - pts.
0640: 62 65 67 69 6e 28 29 3b 0d 0a 09 09 74 79 70 65 begin();....type
0650: 64 65 66 20 69 6e 74 20 43 6f 73 74 3b 0d 0a 09 def int Cost;...
0660: 09 74 79 70 65 64 65 66 20 70 61 69 72 3c 43 6f .typedef pair<Co
0670: 73 74 2c 56 65 72 74 3e 20 45 64 67 65 3b 0d 0a st,Vert> Edge;..
0680: 09 09 74 79 70 65 64 65 66 20 76 65 63 74 6f 72 ..typedef vector
0690: 3c 45 64 67 65 3e 20 45 64 67 65 73 3b 0d 0a 09 <Edge> Edges;...
06a0: 09 74 79 70 65 64 65 66 20 76 65 63 74 6f 72 3c .typedef vector<
06b0: 45 64 67 65 73 3e 20 47 72 61 70 68 3b 0d 0a 0d Edges> Graph;...
06c0: 0a 09 09 47 72 61 70 68 20 47 28 4e 29 3b 0d 0a ...Graph G(N);..
06d0: 09 09 66 6f 72 28 69 6e 74 20 69 3d 30 3b 20 69 ..for(int i=0; i
06e0: 3c 4e 3b 20 2b 2b 69 29 0d 0a 09 09 66 6f 72 28 <N; ++i)....for(
06f0: 69 6e 74 20 6b 3d 69 2b 31 3b 20 6b 3c 4e 3b 20 int k=i+1; k<N;
0700: 2b 2b 6b 29 20 7b 0d 0a 09 09 09 62 6f 6f 6c 20 ++k) {.....bool
0710: 63 6f 6e 6e 65 63 74 65 64 20 3d 20 74 72 75 65 connected = true
0720: 3b 0d 0a 09 09 09 69 6e 74 20 78 31 20 3d 20 67 ;.....int x1 = g
0730: 65 74 3c 58 3e 28 70 74 73 5b 69 5d 29 3b 0d 0a et<X>(pts[i]);..
0740: 09 09 09 69 6e 74 20 78 32 20 3d 20 67 65 74 3c ...int x2 = get<
0750: 58 3e 28 70 74 73 5b 6b 5d 29 3b 0d 0a 09 09 09 X>(pts[k]);.....
0760: 69 6e 74 20 79 31 20 3d 20 67 65 74 3c 59 3e 28 int y1 = get<Y>(
0770: 70 74 73 5b 69 5d 29 3b 0d 0a 09 09 09 69 6e 74 pts[i]);.....int
0780: 20 79 32 20 3d 20 67 65 74 3c 59 3e 28 70 74 73 y2 = get<Y>(pts
0790: 5b 6b 5d 29 3b 0d 0a 09 09 09 66 6f 72 28 50 6f [k]);.....for(Po
07a0: 69 6e 74 20 62 3a 20 62 6c 6f 63 6b 65 64 29 20 int b: blocked)
07b0: 7b 0d 0a 09 09 09 09 69 66 28 6d 69 6e 28 79 31 {......if(min(y1
07c0: 2c 79 32 29 3c 3d 67 65 74 3c 59 3e 28 62 29 20 ,y2)<=get<Y>(b)
07d0: 26 26 20 67 65 74 3c 59 3e 28 62 29 3c 3d 6d 61 && get<Y>(b)<=ma
07e0: 78 28 79 31 2c 79 32 29 20 26 26 0d 0a 09 09 09 x(y1,y2) &&.....
07f0: 09 20 20 20 6d 69 6e 28 78 31 2c 78 32 29 3c 3d . min(x1,x2)<=
0800: 67 65 74 3c 58 3e 28 62 29 20 26 26 20 67 65 74 get<X>(b) && get
0810: 3c 58 3e 28 62 29 3c 3d 6d 61 78 28 78 31 2c 78 <X>(b)<=max(x1,x
0820: 32 29 29 20 7b 0d 0a 09 09 09 09 09 63 6f 6e 6e 2)) {.......conn
0830: 65 63 74 65 64 20 3d 20 66 61 6c 73 65 3b 0d 0a ected = false;..
0840: 09 09 09 09 09 62 72 65 61 6b 3b 0d 0a 09 09 09 .....break;.....
0850: 09 7d 0d 0a 09 09 09 7d 0d 0a 09 09 09 2f 2f 20 .}.....}.....//
0860: 49 66 20 6e 6f 74 20 62 6c 6f 63 6b 65 64 20 74 If not blocked t
0870: 68 65 6e 20 61 64 64 20 69 3c 2d 2d 3e 6b 20 77 hen add i<-->k w
0880: 69 74 68 20 4d 61 6e 68 61 74 74 61 6e 20 63 6f ith Manhattan co
0890: 73 74 0d 0a 09 09 09 69 66 28 63 6f 6e 6e 65 63 st.....if(connec
08a0: 74 65 64 29 20 7b 0d 0a 09 09 09 09 47 5b 69 5d ted) {......G[i]
08b0: 2e 65 6d 70 6c 61 63 65 5f 62 61 63 6b 28 61 62 .emplace_back(ab
08c0: 73 28 78 31 2d 78 32 29 2b 61 62 73 28 79 31 2d s(x1-x2)+abs(y1-
08d0: 79 32 29 2c 20 6b 29 3b 0d 0a 09 09 09 09 47 5b y2), k);......G[
08e0: 6b 5d 2e 65 6d 70 6c 61 63 65 5f 62 61 63 6b 28 k].emplace_back(
08f0: 61 62 73 28 78 31 2d 78 32 29 2b 61 62 73 28 79 abs(x1-x2)+abs(y
0900: 31 2d 79 32 29 2c 20 69 29 3b 0d 0a 09 09 09 7d 1-y2), i);.....}
0910: 0d 0a 09 09 7d 0d 0a 0d 0a 09 09 2f 2f 20 44 69 ....}......// Di
0920: 6a 6b 73 74 72 61 0d 0a 09 09 70 72 69 6f 72 69 jkstra....priori
0930: 74 79 5f 71 75 65 75 65 3c 45 64 67 65 2c 20 76 ty_queue<Edge, v
0940: 65 63 74 6f 72 3c 45 64 67 65 3e 2c 20 67 72 65 ector<Edge>, gre
0950: 61 74 65 72 3c 45 64 67 65 3e 3e 20 51 3b 0d 0a ater<Edge>> Q;..
0960: 09 09 51 2e 65 6d 70 6c 61 63 65 28 30 2c 20 53 ..Q.emplace(0, S
0970: 29 3b 0d 0a 09 09 76 65 63 74 6f 72 3c 43 6f 73 );....vector<Cos
0980: 74 3e 20 44 28 4e 2c 20 2d 31 29 3b 0d 0a 09 09 t> D(N, -1);....
0990: 77 68 69 6c 65 28 21 51 2e 65 6d 70 74 79 28 29 while(!Q.empty()
09a0: 29 20 7b 0d 0a 09 09 09 43 6f 73 74 20 63 20 3d ) {.....Cost c =
09b0: 20 51 2e 74 6f 70 28 29 2e 66 69 72 73 74 3b 0d Q.top().first;.
09c0: 0a 09 09 09 56 65 72 74 20 76 20 3d 20 51 2e 74 ....Vert v = Q.t
09d0: 6f 70 28 29 2e 73 65 63 6f 6e 64 3b 0d 0a 09 09 op().second;....
09e0: 09 51 2e 70 6f 70 28 29 3b 0d 0a 09 09 09 69 66 .Q.pop();.....if
09f0: 28 44 5b 76 5d 21 3d 2d 31 29 0d 0a 09 09 09 09 (D[v]!=-1)......
0a00: 63 6f 6e 74 69 6e 75 65 3b 0d 0a 09 09 09 44 5b continue;.....D[
0a10: 76 5d 20 3d 20 63 3b 0d 0a 0d 0a 09 09 09 66 6f v] = c;.......fo
0a20: 72 28 61 75 74 6f 20 65 3a 20 47 5b 76 5d 29 20 r(auto e: G[v])
0a30: 7b 0d 0a 09 09 09 09 43 6f 73 74 20 63 32 20 3d {......Cost c2 =
0a40: 20 63 20 2b 20 65 2e 66 69 72 73 74 3b 0d 0a 09 c + e.first;...
0a50: 09 09 09 56 65 72 74 20 71 20 3d 20 65 2e 73 65 ...Vert q = e.se
0a60: 63 6f 6e 64 3b 0d 0a 09 09 09 09 69 66 28 44 5b cond;......if(D[
0a70: 71 5d 3d 3d 2d 31 29 0d 0a 09 09 09 09 09 51 2e q]==-1).......Q.
0a80: 65 6d 70 6c 61 63 65 28 63 32 2c 20 71 29 3b 0d emplace(c2, q);.
0a90: 0a 09 09 09 7d 0d 0a 09 09 7d 0d 0a 0d 0a 09 09 ....}....}......
0aa0: 69 6e 74 20 62 65 73 74 20 3d 20 30 3b 0d 0a 09 int best = 0;...
0ab0: 09 66 6f 72 28 56 65 72 74 20 76 3d 30 3b 20 76 .for(Vert v=0; v
0ac0: 3c 4e 3b 20 2b 2b 76 29 20 7b 0d 0a 09 09 09 69 <N; ++v) {.....i
0ad0: 66 28 44 5b 76 5d 20 21 3d 20 2d 31 20 26 26 20 f(D[v] != -1 &&
0ae0: 44 5b 76 5d 3c 3d 6b 29 20 7b 0d 0a 09 09 09 09 D[v]<=k) {......
0af0: 2f 2f 20 74 68 65 20 70 6f 69 6e 74 0d 0a 09 09 // the point....
0b00: 09 09 62 65 73 74 20 3d 20 6d 61 78 28 62 65 73 ..best = max(bes
0b10: 74 2c 20 67 65 74 3c 58 3e 28 70 74 73 5b 76 5d t, get<X>(pts[v]
0b20: 29 29 3b 0d 0a 0d 0a 0d 0a 09 09 09 09 2f 2f 20 ));..........//
0b30: 69 66 20 6e 6f 74 20 62 6c 6f 63 6b 65 64 2e 2e if not blocked..
0b40: 2e 0d 0a 09 09 09 09 50 6f 69 6e 74 20 72 69 67 .......Point rig
0b50: 68 74 28 67 65 74 3c 58 3e 28 70 74 73 5b 76 5d ht(get<X>(pts[v]
0b60: 29 2b 31 2c 20 67 65 74 3c 59 3e 28 70 74 73 5b )+1, get<Y>(pts[
0b70: 76 5d 29 29 3b 0d 0a 09 09 09 09 69 66 28 62 6c v]));......if(bl
0b80: 6f 63 6b 65 64 2e 63 6f 75 6e 74 28 72 69 67 68 ocked.count(righ
0b90: 74 29 29 0d 0a 09 09 09 09 09 63 6f 6e 74 69 6e t)).......contin
0ba0: 75 65 3b 0d 0a 0d 0a 09 09 09 09 2f 2f 20 67 6f ue;........// go
0bb0: 20 6d 6f 72 65 20 65 61 73 74 0d 0a 09 09 09 09 more east......
0bc0: 62 6f 6f 6c 20 73 6b 69 70 20 3d 20 66 61 6c 73 bool skip = fals
0bd0: 65 3b 0d 0a 09 09 09 09 66 6f 72 28 45 64 67 65 e;......for(Edge
0be0: 20 65 3a 20 47 5b 76 5d 29 20 7b 0d 0a 09 09 09 e: G[v]) {.....
0bf0: 09 09 56 65 72 74 20 75 20 3d 20 65 2e 73 65 63 ..Vert u = e.sec
0c00: 6f 6e 64 3b 0d 0a 09 09 09 09 09 69 66 28 67 65 ond;.......if(ge
0c10: 74 3c 59 3e 28 70 74 73 5b 75 5d 29 20 3d 3d 20 t<Y>(pts[u]) ==
0c20: 67 65 74 3c 59 3e 28 70 74 73 5b 76 5d 29 20 26 get<Y>(pts[v]) &
0c30: 26 20 67 65 74 3c 58 3e 28 70 74 73 5b 75 5d 29 & get<X>(pts[u])
0c40: 20 3e 20 67 65 74 3c 58 3e 28 70 74 73 5b 76 5d > get<X>(pts[v]
0c50: 29 0d 0a 09 09 09 09 09 09 26 26 20 44 5b 75 5d )........&& D[u]
0c60: 21 3d 2d 31 20 26 26 20 44 5b 75 5d 3c 3d 6b 29 !=-1 && D[u]<=k)
0c70: 20 7b 0d 0a 09 09 09 09 09 09 73 6b 69 70 20 3d {........skip =
0c80: 20 74 72 75 65 3b 0d 0a 09 09 09 09 09 7d 0d 0a true;.......}..
0c90: 09 09 09 09 7d 0d 0a 09 09 09 09 69 66 28 21 73 ....}......if(!s
0ca0: 6b 69 70 29 20 7b 0d 0a 09 09 09 09 09 69 6e 74 kip) {.......int
0cb0: 20 73 65 63 20 3d 20 6b 2d 44 5b 76 5d 3b 0d 0a sec = k-D[v];..
0cc0: 09 09 09 09 09 62 65 73 74 20 3d 20 6d 61 78 28 .....best = max(
0cd0: 62 65 73 74 2c 20 67 65 74 3c 58 3e 28 70 74 73 best, get<X>(pts
0ce0: 5b 76 5d 29 20 2b 20 73 65 63 29 3b 0d 0a 09 09 [v]) + sec);....
0cf0: 09 09 7d 0d 0a 09 09 09 7d 0d 0a 09 09 7d 0d 0a ..}.....}....}..
0d00: 0d 0a 09 09 72 65 74 75 72 6e 20 62 65 73 74 3b ....return best;
0d10: 0d 0a 09 7d 0d 0a 7d 3b 0d 0a 0d 0a 2f 2f 20 42 ...}..};....// B
0d20: 45 47 49 4e 20 43 55 54 20 48 45 52 45 0d 0a 23 EGIN CUT HERE..#
0d30: 69 6e 63 6c 75 64 65 20 3c 63 74 69 6d 65 3e 0d include <ctime>.
0d40: 0a 64 6f 75 62 6c 65 20 73 74 61 72 74 5f 74 69 .double start_ti
0d50: 6d 65 3b 20 73 74 72 69 6e 67 20 74 69 6d 65 72 me; string timer
0d60: 28 29 0d 0a 20 7b 20 6f 73 74 72 69 6e 67 73 74 ().. { ostringst
0d70: 72 65 61 6d 20 6f 73 3b 20 6f 73 20 3c 3c 20 22 ream os; os << "
0d80: 20 28 22 20 3c 3c 20 69 6e 74 28 28 63 6c 6f 63 (" << int((cloc
0d90: 6b 28 29 2d 73 74 61 72 74 5f 74 69 6d 65 29 2f k()-start_time)/
0da0: 43 4c 4f 43 4b 53 5f 50 45 52 5f 53 45 43 2a 31 CLOCKS_PER_SEC*1
0db0: 30 30 30 29 20 3c 3c 20 22 20 6d 73 65 63 29 22 000) << " msec)"
0dc0: 3b 20 72 65 74 75 72 6e 20 6f 73 2e 73 74 72 28 ; return os.str(
0dd0: 29 3b 20 7d 0d 0a 74 65 6d 70 6c 61 74 65 3c 74 ); }..template<t
0de0: 79 70 65 6e 61 6d 65 20 54 3e 20 6f 73 74 72 65 ypename T> ostre
0df0: 61 6d 26 20 6f 70 65 72 61 74 6f 72 3c 3c 28 6f am& operator<<(o
0e00: 73 74 72 65 61 6d 26 20 6f 73 2c 20 63 6f 6e 73 stream& os, cons
0e10: 74 20 76 65 63 74 6f 72 3c 54 3e 26 20 76 29 0d t vector<T>& v).
0e20: 0a 20 7b 20 6f 73 20 3c 3c 20 22 7b 20 22 3b 0d . { os << "{ ";.
0e30: 0a 20 20 20 66 6f 72 28 74 79 70 65 6e 61 6d 65 . for(typename
0e40: 20 76 65 63 74 6f 72 3c 54 3e 3a 3a 63 6f 6e 73 vector<T>::cons
0e50: 74 5f 69 74 65 72 61 74 6f 72 20 69 74 3d 76 2e t_iterator it=v.
0e60: 62 65 67 69 6e 28 29 3b 20 69 74 21 3d 76 2e 65 begin(); it!=v.e
0e70: 6e 64 28 29 3b 20 2b 2b 69 74 29 0d 0a 20 20 20 nd(); ++it)..
0e80: 6f 73 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 2a 69 os << '\"' << *i
0e90: 74 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 28 69 74 t << '\"' << (it
0ea0: 2b 31 3d 3d 76 2e 65 6e 64 28 29 20 3f 20 22 22 +1==v.end() ? ""
0eb0: 20 3a 20 22 2c 20 22 29 3b 20 6f 73 20 3c 3c 20 : ", "); os <<
0ec0: 22 20 7d 22 3b 20 72 65 74 75 72 6e 20 6f 73 3b " }"; return os;
0ed0: 20 7d 0d 0a 76 6f 69 64 20 76 65 72 69 66 79 5f }..void verify_
0ee0: 63 61 73 65 28 63 6f 6e 73 74 20 69 6e 74 26 20 case(const int&
0ef0: 45 78 70 65 63 74 65 64 2c 20 63 6f 6e 73 74 20 Expected, const
0f00: 69 6e 74 26 20 52 65 63 65 69 76 65 64 29 20 7b int& Received) {
0f10: 0d 0a 20 62 6f 6f 6c 20 6f 6b 20 3d 20 28 45 78 .. bool ok = (Ex
0f20: 70 65 63 74 65 64 20 3d 3d 20 52 65 63 65 69 76 pected == Receiv
0f30: 65 64 29 3b 0d 0a 20 69 66 28 6f 6b 29 20 63 65 ed);.. if(ok) ce
0f40: 72 72 20 3c 3c 20 22 50 41 53 53 45 44 22 20 3c rr << "PASSED" <
0f50: 3c 20 74 69 6d 65 72 28 29 20 3c 3c 20 65 6e 64 < timer() << end
0f60: 6c 3b 20 20 65 6c 73 65 20 7b 20 63 65 72 72 20 l; else { cerr
0f70: 3c 3c 20 22 46 41 49 4c 45 44 22 20 3c 3c 20 74 << "FAILED" << t
0f80: 69 6d 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b 0d imer() << endl;.
0f90: 0a 20 63 65 72 72 20 3c 3c 20 22 5c 74 6f 3a 20 . cerr << "\to:
0fa0: 5c 22 22 20 3c 3c 20 45 78 70 65 63 74 65 64 20 \"" << Expected
0fb0: 3c 3c 20 27 5c 22 27 20 3c 3c 20 65 6e 64 6c 20 << '\"' << endl
0fc0: 3c 3c 20 22 5c 74 78 3a 20 5c 22 22 20 3c 3c 20 << "\tx: \"" <<
0fd0: 52 65 63 65 69 76 65 64 20 3c 3c 20 27 5c 22 27 Received << '\"'
0fe0: 20 3c 3c 20 65 6e 64 6c 3b 20 7d 20 7d 0d 0a 23 << endl; } }..#
0ff0: 64 65 66 69 6e 65 20 43 41 53 45 28 4e 29 20 7b define CASE(N) {
1000: 63 65 72 72 20 3c 3c 20 22 54 65 73 74 20 43 61 cerr << "Test Ca
1010: 73 65 20 23 22 20 3c 3c 20 4e 20 3c 3c 20 22 2e se #" << N << ".
1020: 2e 2e 22 20 3c 3c 20 66 6c 75 73 68 3b 20 73 74 .." << flush; st
1030: 61 72 74 5f 74 69 6d 65 3d 63 6c 6f 63 6b 28 29 art_time=clock()
1040: 3b 0d 0a 23 64 65 66 69 6e 65 20 45 4e 44 09 20 ;..#define END.
1050: 76 65 72 69 66 79 5f 63 61 73 65 28 5f 2c 20 54 verify_case(_, T
1060: 68 65 47 72 69 64 44 69 76 4f 6e 65 28 29 2e 66 heGridDivOne().f
1070: 69 6e 64 28 78 2c 20 79 2c 20 6b 29 29 3b 7d 0d ind(x, y, k));}.
1080: 0a 69 6e 74 20 6d 61 69 6e 28 29 7b 0d 0a 0d 0a .int main(){....
1090: 43 41 53 45 28 30 29 0d 0a 09 69 6e 74 20 78 5f CASE(0)...int x_
10a0: 5b 5d 20 3d 20 7b 31 2c 31 2c 31 2c 31 7d 3b 0d [] = {1,1,1,1};.
10b0: 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e .. vector <int>
10c0: 20 78 28 78 5f 2c 20 78 5f 2b 73 69 7a 65 6f 66 x(x_, x_+sizeof
10d0: 28 78 5f 29 2f 73 69 7a 65 6f 66 28 2a 78 5f 29 (x_)/sizeof(*x_)
10e0: 29 3b 20 0d 0a 09 69 6e 74 20 79 5f 5b 5d 20 3d ); ...int y_[] =
10f0: 20 7b 2d 32 2c 2d 31 2c 30 2c 31 7d 3b 0d 0a 09 {-2,-1,0,1};...
1100: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79 vector <int> y
1110: 28 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f 66 28 79 (y_, y_+sizeof(y
1120: 5f 29 2f 73 69 7a 65 6f 66 28 2a 79 5f 29 29 3b _)/sizeof(*y_));
1130: 20 0d 0a 09 69 6e 74 20 6b 20 3d 20 34 3b 20 0d ...int k = 4; .
1140: 0a 09 69 6e 74 20 5f 20 3d 20 32 3b 20 0d 0a 45 ..int _ = 2; ..E
1150: 4e 44 0d 0a 43 41 53 45 28 31 29 0d 0a 09 69 6e ND..CASE(1)...in
1160: 74 20 78 5f 5b 5d 20 3d 20 7b 2d 31 2c 20 30 2c t x_[] = {-1, 0,
1170: 20 30 2c 20 31 7d 3b 0d 0a 09 20 20 76 65 63 74 0, 1};... vect
1180: 6f 72 20 3c 69 6e 74 3e 20 78 28 78 5f 2c 20 78 or <int> x(x_, x
1190: 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 2f 73 69 7a _+sizeof(x_)/siz
11a0: 65 6f 66 28 2a 78 5f 29 29 3b 20 0d 0a 09 69 6e eof(*x_)); ...in
11b0: 74 20 79 5f 5b 5d 20 3d 20 7b 30 2c 20 2d 31 2c t y_[] = {0, -1,
11c0: 20 31 2c 20 30 7d 3b 0d 0a 09 20 20 76 65 63 74 1, 0};... vect
11d0: 6f 72 20 3c 69 6e 74 3e 20 79 28 79 5f 2c 20 79 or <int> y(y_, y
11e0: 5f 2b 73 69 7a 65 6f 66 28 79 5f 29 2f 73 69 7a _+sizeof(y_)/siz
11f0: 65 6f 66 28 2a 79 5f 29 29 3b 20 0d 0a 09 69 6e eof(*y_)); ...in
1200: 74 20 6b 20 3d 20 39 3b 20 0d 0a 09 69 6e 74 20 t k = 9; ...int
1210: 5f 20 3d 20 30 3b 20 0d 0a 45 4e 44 0d 0a 43 41 _ = 0; ..END..CA
1220: 53 45 28 32 29 0d 0a 09 76 65 63 74 6f 72 20 3c SE(2)...vector <
1230: 69 6e 74 3e 20 78 3b 20 0d 0a 09 76 65 63 74 6f int> x; ...vecto
1240: 72 20 3c 69 6e 74 3e 20 79 3b 20 0d 0a 09 69 6e r <int> y; ...in
1250: 74 20 6b 20 3d 20 31 30 30 30 3b 20 0d 0a 09 69 t k = 1000; ...i
1260: 6e 74 20 5f 20 3d 20 31 30 30 30 3b 20 0d 0a 45 nt _ = 1000; ..E
1270: 4e 44 0d 0a 43 41 53 45 28 33 29 0d 0a 09 69 6e ND..CASE(3)...in
1280: 74 20 78 5f 5b 5d 20 3d 20 7b 31 2c 30 2c 30 2c t x_[] = {1,0,0,
1290: 2d 31 2c 2d 31 2c 2d 32 2c 2d 32 2c 2d 33 2c 2d -1,-1,-2,-2,-3,-
12a0: 33 2c 2d 34 2c 2d 34 7d 3b 0d 0a 09 20 20 76 65 3,-4,-4};... ve
12b0: 63 74 6f 72 20 3c 69 6e 74 3e 20 78 28 78 5f 2c ctor <int> x(x_,
12c0: 20 78 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 2f 73 x_+sizeof(x_)/s
12d0: 69 7a 65 6f 66 28 2a 78 5f 29 29 3b 20 0d 0a 09 izeof(*x_)); ...
12e0: 69 6e 74 20 79 5f 5b 5d 20 3d 20 7b 30 2c 2d 31 int y_[] = {0,-1
12f0: 2c 31 2c 2d 32 2c 32 2c 2d 33 2c 33 2c 2d 34 2c ,1,-2,2,-3,3,-4,
1300: 34 2c 2d 35 2c 35 7d 3b 0d 0a 09 20 20 76 65 63 4,-5,5};... vec
1310: 74 6f 72 20 3c 69 6e 74 3e 20 79 28 79 5f 2c 20 tor <int> y(y_,
1320: 79 5f 2b 73 69 7a 65 6f 66 28 79 5f 29 2f 73 69 y_+sizeof(y_)/si
1330: 7a 65 6f 66 28 2a 79 5f 29 29 3b 20 0d 0a 09 69 zeof(*y_)); ...i
1340: 6e 74 20 6b 20 3d 20 34 37 3b 20 0d 0a 09 69 6e nt k = 47; ...in
1350: 74 20 5f 20 3d 20 33 31 3b 20 0d 0a 45 4e 44 0d t _ = 31; ..END.
1360: 0a 2f 2a 0d 0a 43 41 53 45 28 34 29 0d 0a 09 69 ./*..CASE(4)...i
1370: 6e 74 20 78 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 nt x_[] = ;...
1380: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 78 28 78 vector <int> x(x
1390: 5f 2c 20 78 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 _, x_+sizeof(x_)
13a0: 2f 73 69 7a 65 6f 66 28 2a 78 5f 29 29 3b 20 0d /sizeof(*x_)); .
13b0: 0a 09 69 6e 74 20 79 5f 5b 5d 20 3d 20 3b 0d 0a ..int y_[] = ;..
13c0: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 . vector <int>
13d0: 79 28 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f 66 28 y(y_, y_+sizeof(
13e0: 79 5f 29 2f 73 69 7a 65 6f 66 28 2a 79 5f 29 29 y_)/sizeof(*y_))
13f0: 3b 20 0d 0a 09 69 6e 74 20 6b 20 3d 20 3b 20 0d ; ...int k = ; .
1400: 0a 09 69 6e 74 20 5f 20 3d 20 3b 20 0d 0a 45 4e ..int _ = ; ..EN
1410: 44 0d 0a 43 41 53 45 28 35 29 0d 0a 09 69 6e 74 D..CASE(5)...int
1420: 20 78 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 65 x_[] = ;... ve
1430: 63 74 6f 72 20 3c 69 6e 74 3e 20 78 28 78 5f 2c ctor <int> x(x_,
1440: 20 78 5f 2b 73 69 7a 65 6f 66 28 78 5f 29 2f 73 x_+sizeof(x_)/s
1450: 69 7a 65 6f 66 28 2a 78 5f 29 29 3b 20 0d 0a 09 izeof(*x_)); ...
1460: 69 6e 74 20 79 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 int y_[] = ;...
1470: 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 79 28 vector <int> y(
1480: 79 5f 2c 20 79 5f 2b 73 69 7a 65 6f 66 28 79 5f y_, y_+sizeof(y_
1490: 29 2f 73 69 7a 65 6f 66 28 2a 79 5f 29 29 3b 20 )/sizeof(*y_));
14a0: 0d 0a 09 69 6e 74 20 6b 20 3d 20 3b 20 0d 0a 09 ...int k = ; ...
14b0: 69 6e 74 20 5f 20 3d 20 3b 20 0d 0a 45 4e 44 0d int _ = ; ..END.
14c0: 0a 2a 2f 0d 0a 7d 0d 0a 2f 2f 20 45 4e 44 20 43 .*/..}..// END C
14d0: 55 54 20 48 45 52 45 0d 0a UT HERE..