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 74 65 6d 70 6c 61 74 65 3c 74 P;....template<t
01a0: 79 70 65 6e 61 6d 65 20 54 3e 0d 0a 63 6c 61 73 ypename T>..clas
01b0: 73 20 49 64 47 65 6e 0d 0a 7b 0d 0a 09 6d 61 70 s IdGen..{...map
01c0: 3c 54 2c 20 69 6e 74 3e 20 76 32 69 64 5f 3b 0d <T, int> v2id_;.
01d0: 0a 09 76 65 63 74 6f 72 3c 54 3e 20 20 20 69 64 ..vector<T> id
01e0: 32 76 5f 3b 0d 0a 70 75 62 6c 69 63 3a 0d 0a 09 2v_;..public:...
01f0: 69 6e 74 20 76 32 69 64 28 63 6f 6e 73 74 20 54 int v2id(const T
0200: 26 20 76 29 20 7b 0d 0a 09 09 69 66 28 20 21 76 & v) {....if( !v
0210: 32 69 64 5f 2e 63 6f 75 6e 74 28 76 29 20 29 20 2id_.count(v) )
0220: 7b 20 76 32 69 64 5f 5b 76 5d 20 3d 20 73 69 7a { v2id_[v] = siz
0230: 65 28 29 3b 20 69 64 32 76 5f 2e 70 75 73 68 5f e(); id2v_.push_
0240: 62 61 63 6b 28 76 29 3b 20 7d 0d 0a 09 09 72 65 back(v); }....re
0250: 74 75 72 6e 20 76 32 69 64 5f 5b 76 5d 3b 0d 0a turn v2id_[v];..
0260: 09 7d 0d 0a 09 63 6f 6e 73 74 20 54 26 20 69 64 .}...const T& id
0270: 32 76 28 69 6e 74 20 69 29 20 63 6f 6e 73 74 20 2v(int i) const
0280: 7b 20 72 65 74 75 72 6e 20 69 64 32 76 5f 5b 69 { return id2v_[i
0290: 5d 3b 20 7d 0d 0a 09 69 6e 74 20 73 69 7a 65 28 ]; }...int size(
02a0: 29 20 63 6f 6e 73 74 20 7b 20 72 65 74 75 72 6e ) const { return
02b0: 20 69 64 32 76 5f 2e 73 69 7a 65 28 29 3b 20 7d id2v_.size(); }
02c0: 0d 0a 7d 3b 0d 0a 0d 0a 74 65 6d 70 6c 61 74 65 ..};....template
02d0: 3c 74 79 70 65 6e 61 6d 65 20 56 65 72 74 3d 70 <typename Vert=p
02e0: 61 69 72 3c 69 6e 74 2c 69 6e 74 3e 2c 20 74 79 air<int,int>, ty
02f0: 70 65 6e 61 6d 65 20 46 6c 6f 77 3d 4c 4c 2c 20 pename Flow=LL,
0300: 69 6e 74 20 4e 56 3d 32 30 34 38 3e 0d 0a 63 6c int NV=2048>..cl
0310: 61 73 73 20 4d 61 78 46 6c 6f 77 0d 0a 7b 0d 0a ass MaxFlow..{..
0320: 09 49 64 47 65 6e 3c 56 65 72 74 3e 20 69 64 67 .IdGen<Vert> idg
0330: 65 6e 3b 0d 0a 09 76 65 63 74 6f 72 3c 69 6e 74 en;...vector<int
0340: 3e 20 47 5b 4e 56 5d 3b 0d 0a 09 46 6c 6f 77 20 > G[NV];...Flow
0350: 46 5b 4e 56 5d 5b 4e 56 5d 3b 0d 0a 0d 0a 70 75 F[NV][NV];....pu
0360: 62 6c 69 63 3a 0d 0a 09 76 6f 69 64 20 61 64 64 blic:...void add
0370: 45 64 67 65 28 20 56 65 72 74 20 73 5f 2c 20 56 Edge( Vert s_, V
0380: 65 72 74 20 74 5f 2c 20 46 6c 6f 77 20 66 20 29 ert t_, Flow f )
0390: 0d 0a 09 7b 0d 0a 09 09 63 6f 6e 73 74 20 69 6e ...{....const in
03a0: 74 20 73 20 3d 20 69 64 67 65 6e 2e 76 32 69 64 t s = idgen.v2id
03b0: 28 73 5f 29 2c 20 74 20 3d 20 69 64 67 65 6e 2e (s_), t = idgen.
03c0: 76 32 69 64 28 74 5f 29 3b 0d 0a 09 09 47 5b 73 v2id(t_);....G[s
03d0: 5d 2e 70 75 73 68 5f 62 61 63 6b 28 74 29 3b 0d ].push_back(t);.
03e0: 0a 09 09 47 5b 74 5d 2e 70 75 73 68 5f 62 61 63 ...G[t].push_bac
03f0: 6b 28 73 29 3b 0d 0a 09 09 46 5b 73 5d 5b 74 5d k(s);....F[s][t]
0400: 20 3d 20 66 3b 0d 0a 09 09 46 5b 74 5d 5b 73 5d = f;....F[t][s]
0410: 20 3d 20 30 3b 0d 0a 09 7d 0d 0a 0d 0a 09 46 6c = 0;...}.....Fl
0420: 6f 77 20 63 61 6c 63 28 20 56 65 72 74 20 73 5f ow calc( Vert s_
0430: 2c 20 56 65 72 74 20 74 5f 20 29 0d 0a 09 7b 0d , Vert t_ )...{.
0440: 0a 09 09 63 6f 6e 73 74 20 69 6e 74 20 53 20 3d ...const int S =
0450: 20 69 64 67 65 6e 2e 76 32 69 64 28 73 5f 29 2c idgen.v2id(s_),
0460: 20 44 20 3d 20 69 64 67 65 6e 2e 76 32 69 64 28 D = idgen.v2id(
0470: 74 5f 29 3b 0d 0a 09 09 66 6f 72 28 20 46 6c 6f t_);....for( Flo
0480: 77 20 74 6f 74 61 6c 3d 30 20 3b 3b 20 29 20 7b w total=0 ;; ) {
0490: 0d 0a 09 09 09 2f 2f 20 44 6f 20 42 46 53 20 61 .....// Do BFS a
04a0: 6e 64 20 63 6f 6d 70 75 74 65 20 74 68 65 20 6c nd compute the l
04b0: 65 76 65 6c 20 66 6f 72 20 65 61 63 68 20 6e 6f evel for each no
04c0: 64 65 2e 0d 0a 09 09 09 69 6e 74 20 4c 56 5b 4e de......int LV[N
04d0: 56 5d 20 3d 20 7b 30 7d 3b 0d 0a 09 09 09 76 65 V] = {0};.....ve
04e0: 63 74 6f 72 3c 69 6e 74 3e 20 51 28 31 2c 20 53 ctor<int> Q(1, S
04f0: 29 3b 0d 0a 09 09 09 66 6f 72 28 69 6e 74 20 6c );.....for(int l
0500: 76 3d 31 3b 20 21 51 2e 65 6d 70 74 79 28 29 3b v=1; !Q.empty();
0510: 20 2b 2b 6c 76 29 20 7b 0d 0a 09 09 09 09 76 65 ++lv) {......ve
0520: 63 74 6f 72 3c 69 6e 74 3e 20 51 32 3b 0d 0a 09 ctor<int> Q2;...
0530: 09 09 09 66 6f 72 28 73 69 7a 65 5f 74 20 69 3d ...for(size_t i=
0540: 30 3b 20 69 21 3d 51 2e 73 69 7a 65 28 29 3b 20 0; i!=Q.size();
0550: 2b 2b 69 29 20 7b 0d 0a 09 09 09 09 09 63 6f 6e ++i) {.......con
0560: 73 74 20 76 65 63 74 6f 72 3c 69 6e 74 3e 26 20 st vector<int>&
0570: 6e 65 20 3d 20 47 5b 51 5b 69 5d 5d 3b 0d 0a 09 ne = G[Q[i]];...
0580: 09 09 09 09 66 6f 72 28 73 69 7a 65 5f 74 20 6a ....for(size_t j
0590: 3d 30 3b 20 6a 21 3d 6e 65 2e 73 69 7a 65 28 29 =0; j!=ne.size()
05a0: 3b 20 2b 2b 6a 29 0d 0a 09 09 09 09 09 09 69 66 ; ++j)........if
05b0: 28 20 46 5b 51 5b 69 5d 5d 5b 6e 65 5b 6a 5d 5d ( F[Q[i]][ne[j]]
05c0: 20 26 26 20 21 4c 56 5b 6e 65 5b 6a 5d 5d 20 26 && !LV[ne[j]] &
05d0: 26 20 6e 65 5b 6a 5d 21 3d 53 20 29 0d 0a 09 09 & ne[j]!=S )....
05e0: 09 09 09 09 09 4c 56 5b 6e 65 5b 6a 5d 5d 3d 6c .....LV[ne[j]]=l
05f0: 76 2c 20 51 32 2e 70 75 73 68 5f 62 61 63 6b 28 v, Q2.push_back(
0600: 6e 65 5b 6a 5d 29 3b 0d 0a 09 09 09 09 7d 0d 0a ne[j]);......}..
0610: 09 09 09 09 51 2e 73 77 61 70 28 51 32 29 3b 0d ....Q.swap(Q2);.
0620: 0a 09 09 09 7d 0d 0a 0d 0a 09 09 09 2f 2f 20 44 ....}.......// D
0630: 65 73 74 69 6e 61 74 69 6f 6e 20 69 73 20 6e 6f estination is no
0640: 77 20 75 6e 72 65 61 63 68 61 62 6c 65 2e 20 44 w unreachable. D
0650: 6f 6e 65 2e 0d 0a 09 09 09 69 66 28 20 21 4c 56 one......if( !LV
0660: 5b 44 5d 20 29 0d 0a 09 09 09 09 72 65 74 75 72 [D] )......retur
0670: 6e 20 74 6f 74 61 6c 3b 0d 0a 0d 0a 09 09 09 2f n total;......./
0680: 2f 20 49 74 65 72 61 74 69 6e 67 20 44 46 53 2e / Iterating DFS.
0690: 0d 0a 09 09 09 62 6f 6f 6c 20 62 6c 6f 63 6b 65 .....bool blocke
06a0: 64 5b 4e 56 5d 20 3d 20 7b 7d 3b 0d 0a 09 09 09 d[NV] = {};.....
06b0: 74 6f 74 61 6c 20 2b 3d 20 64 69 6e 69 63 5f 64 total += dinic_d
06c0: 66 73 28 20 53 2c 20 44 2c 20 4c 56 2c 20 30 78 fs( S, D, LV, 0x
06d0: 37 66 66 66 66 66 66 66 2c 20 62 6c 6f 63 6b 65 7fffffff, blocke
06e0: 64 20 29 3b 0d 0a 09 09 7d 0d 0a 09 7d 0d 0a 0d d );....}...}...
06f0: 0a 70 72 69 76 61 74 65 3a 0d 0a 09 46 6c 6f 77 .private:...Flow
0700: 20 64 69 6e 69 63 5f 64 66 73 28 20 69 6e 74 20 dinic_dfs( int
0710: 76 2c 20 69 6e 74 20 44 2c 20 69 6e 74 20 4c 56 v, int D, int LV
0720: 5b 5d 2c 20 46 6c 6f 77 20 66 6c 6f 77 5f 69 6e [], Flow flow_in
0730: 2c 20 62 6f 6f 6c 20 62 6c 6f 63 6b 65 64 5b 5d , bool blocked[]
0740: 20 29 0d 0a 09 7b 0d 0a 09 09 46 6c 6f 77 20 66 )...{....Flow f
0750: 6c 6f 77 5f 6f 75 74 20 3d 20 30 3b 0d 0a 09 09 low_out = 0;....
0760: 66 6f 72 28 73 69 7a 65 5f 74 20 69 3d 30 3b 20 for(size_t i=0;
0770: 69 21 3d 47 5b 76 5d 2e 73 69 7a 65 28 29 3b 20 i!=G[v].size();
0780: 2b 2b 69 29 20 7b 0d 0a 09 09 09 69 6e 74 20 75 ++i) {.....int u
0790: 20 3d 20 47 5b 76 5d 5b 69 5d 3b 0d 0a 09 09 09 = G[v][i];.....
07a0: 69 66 28 20 4c 56 5b 76 5d 2b 31 3d 3d 4c 56 5b if( LV[v]+1==LV[
07b0: 75 5d 20 26 26 20 46 5b 76 5d 5b 75 5d 20 29 20 u] && F[v][u] )
07c0: 7b 0d 0a 09 09 09 09 46 6c 6f 77 20 66 20 3d 20 {......Flow f =
07d0: 6d 69 6e 28 66 6c 6f 77 5f 69 6e 2d 66 6c 6f 77 min(flow_in-flow
07e0: 5f 6f 75 74 2c 20 46 5b 76 5d 5b 75 5d 29 3b 0d _out, F[v][u]);.
07f0: 0a 09 09 09 09 69 66 28 20 75 3d 3d 44 20 7c 7c .....if( u==D ||
0800: 20 21 62 6c 6f 63 6b 65 64 5b 75 5d 20 26 26 20 !blocked[u] &&
0810: 28 66 3d 64 69 6e 69 63 5f 64 66 73 28 75 2c 44 (f=dinic_dfs(u,D
0820: 2c 4c 56 2c 66 2c 62 6c 6f 63 6b 65 64 29 29 3e ,LV,f,blocked))>
0830: 30 20 29 20 7b 0d 0a 09 09 09 09 09 46 5b 76 5d 0 ) {.......F[v]
0840: 5b 75 5d 20 20 2d 3d 20 66 3b 0d 0a 09 09 09 09 [u] -= f;......
0850: 09 46 5b 75 5d 5b 76 5d 20 20 2b 3d 20 66 3b 0d .F[u][v] += f;.
0860: 0a 09 09 09 09 09 66 6c 6f 77 5f 6f 75 74 20 2b ......flow_out +
0870: 3d 20 66 3b 0d 0a 09 09 09 09 09 69 66 28 20 66 = f;.......if( f
0880: 6c 6f 77 5f 69 6e 20 3d 3d 20 66 6c 6f 77 5f 6f low_in == flow_o
0890: 75 74 20 29 20 72 65 74 75 72 6e 20 66 6c 6f 77 ut ) return flow
08a0: 5f 6f 75 74 3b 0d 0a 09 09 09 09 7d 0d 0a 09 09 _out;......}....
08b0: 09 7d 0d 0a 09 09 7d 0d 0a 09 09 62 6c 6f 63 6b .}....}....block
08c0: 65 64 5b 76 5d 20 3d 20 28 66 6c 6f 77 5f 6f 75 ed[v] = (flow_ou
08d0: 74 3d 3d 30 29 3b 0d 0a 09 09 72 65 74 75 72 6e t==0);....return
08e0: 20 66 6c 6f 77 5f 6f 75 74 3b 0d 0a 09 7d 0d 0a flow_out;...}..
08f0: 7d 3b 0d 0a 0d 0a 63 6c 61 73 73 20 53 69 6e 67 };....class Sing
0900: 69 6e 67 20 7b 20 70 75 62 6c 69 63 3a 0d 0a 09 ing { public:...
0910: 69 6e 74 20 73 6f 6c 76 65 28 69 6e 74 20 4e 2c int solve(int N,
0920: 20 69 6e 74 20 6c 6f 77 2c 20 69 6e 74 20 68 69 int low, int hi
0930: 67 68 2c 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e gh, vector <int>
0940: 20 70 69 74 63 68 29 0d 0a 09 7b 0d 0a 09 09 6d pitch)...{....m
0950: 61 70 3c 70 61 69 72 3c 69 6e 74 2c 69 6e 74 3e ap<pair<int,int>
0960: 2c 20 69 6e 74 3e 20 65 64 67 65 3b 0d 0a 09 09 , int> edge;....
0970: 66 6f 72 28 69 6e 74 20 69 3d 30 3b 20 69 2b 31 for(int i=0; i+1
0980: 3c 70 69 74 63 68 2e 73 69 7a 65 28 29 3b 20 2b <pitch.size(); +
0990: 2b 69 29 20 7b 0d 0a 09 09 09 69 6e 74 20 70 30 +i) {.....int p0
09a0: 20 3d 20 6d 69 6e 28 70 69 74 63 68 5b 69 5d 2c = min(pitch[i],
09b0: 20 70 69 74 63 68 5b 69 2b 31 5d 29 3b 0d 0a 09 pitch[i+1]);...
09c0: 09 09 69 6e 74 20 70 31 20 3d 20 6d 61 78 28 70 ..int p1 = max(p
09d0: 69 74 63 68 5b 69 5d 2c 20 70 69 74 63 68 5b 69 itch[i], pitch[i
09e0: 2b 31 5d 29 3b 0d 0a 09 09 09 69 66 28 70 30 20 +1]);.....if(p0
09f0: 21 3d 20 70 31 29 0d 0a 09 09 09 09 65 64 67 65 != p1)......edge
0a00: 5b 6d 61 6b 65 5f 70 61 69 72 28 70 30 2c 70 31 [make_pair(p0,p1
0a10: 29 5d 20 2b 2b 3b 0d 0a 09 09 7d 0d 0a 0d 0a 09 )] ++;....}.....
0a20: 09 4d 61 78 46 6c 6f 77 3c 3e 2a 20 6d 66 20 3d .MaxFlow<>* mf =
0a30: 20 6e 65 77 20 4d 61 78 46 6c 6f 77 3c 3e 3b 0d new MaxFlow<>;.
0a40: 0a 09 09 63 6f 6e 73 74 20 4c 4c 20 49 4e 46 20 ...const LL INF
0a50: 3d 20 30 78 33 66 66 66 66 66 66 66 3b 0d 0a 0d = 0x3fffffff;...
0a60: 0a 09 09 65 6e 75 6d 7b 42 4f 42 2c 20 4c 45 46 ...enum{BOB, LEF
0a70: 54 2c 20 52 49 47 48 54 2c 20 41 4c 49 43 45 7d T, RIGHT, ALICE}
0a80: 3b 0d 0a 09 09 70 61 69 72 3c 69 6e 74 2c 69 6e ;....pair<int,in
0a90: 74 3e 20 42 6f 62 28 42 4f 42 2c 30 29 3b 0d 0a t> Bob(BOB,0);..
0aa0: 09 09 70 61 69 72 3c 69 6e 74 2c 69 6e 74 3e 20 ..pair<int,int>
0ab0: 41 6c 69 63 65 28 41 4c 49 43 45 2c 30 29 3b 0d Alice(ALICE,0);.
0ac0: 0a 09 09 66 6f 72 28 69 6e 74 20 76 3d 31 3b 20 ...for(int v=1;
0ad0: 76 3c 3d 4e 3b 20 2b 2b 76 29 20 7b 0d 0a 09 09 v<=N; ++v) {....
0ae0: 09 69 66 28 76 3e 68 69 67 68 29 0d 0a 09 09 09 .if(v>high).....
0af0: 09 6d 66 2d 3e 61 64 64 45 64 67 65 28 42 6f 62 .mf->addEdge(Bob
0b00: 2c 20 6d 61 6b 65 5f 70 61 69 72 28 4c 45 46 54 , make_pair(LEFT
0b10: 2c 76 29 2c 20 49 4e 46 29 3b 0d 0a 09 09 09 6d ,v), INF);.....m
0b20: 66 2d 3e 61 64 64 45 64 67 65 28 6d 61 6b 65 5f f->addEdge(make_
0b30: 70 61 69 72 28 4c 45 46 54 2c 76 29 2c 20 6d 61 pair(LEFT,v), ma
0b40: 6b 65 5f 70 61 69 72 28 52 49 47 48 54 2c 76 29 ke_pair(RIGHT,v)
0b50: 2c 20 49 4e 46 29 3b 0d 0a 09 09 09 69 66 28 76 , INF);.....if(v
0b60: 3c 6c 6f 77 29 0d 0a 09 09 09 09 6d 66 2d 3e 61 <low)......mf->a
0b70: 64 64 45 64 67 65 28 6d 61 6b 65 5f 70 61 69 72 ddEdge(make_pair
0b80: 28 52 49 47 48 54 2c 76 29 2c 20 41 6c 69 63 65 (RIGHT,v), Alice
0b90: 2c 20 49 4e 46 29 3b 0d 0a 09 09 7d 0d 0a 09 09 , INF);....}....
0ba0: 66 6f 72 28 61 75 74 6f 20 65 3a 20 65 64 67 65 for(auto e: edge
0bb0: 29 20 7b 0d 0a 09 09 09 69 6e 74 20 76 20 3d 20 ) {.....int v =
0bc0: 65 2e 66 69 72 73 74 2e 66 69 72 73 74 3b 0d 0a e.first.first;..
0bd0: 09 09 09 69 6e 74 20 75 20 3d 20 65 2e 66 69 72 ...int u = e.fir
0be0: 73 74 2e 73 65 63 6f 6e 64 3b 0d 0a 09 09 09 69 st.second;.....i
0bf0: 6e 74 20 63 20 3d 20 65 2e 73 65 63 6f 6e 64 3b nt c = e.second;
0c00: 0d 0a 09 09 09 6d 66 2d 3e 61 64 64 45 64 67 65 .....mf->addEdge
0c10: 28 6d 61 6b 65 5f 70 61 69 72 28 52 49 47 48 54 (make_pair(RIGHT
0c20: 2c 76 29 2c 20 6d 61 6b 65 5f 70 61 69 72 28 4c ,v), make_pair(L
0c30: 45 46 54 2c 75 29 2c 20 63 29 3b 0d 0a 09 09 09 EFT,u), c);.....
0c40: 6d 66 2d 3e 61 64 64 45 64 67 65 28 6d 61 6b 65 mf->addEdge(make
0c50: 5f 70 61 69 72 28 52 49 47 48 54 2c 75 29 2c 20 _pair(RIGHT,u),
0c60: 6d 61 6b 65 5f 70 61 69 72 28 4c 45 46 54 2c 76 make_pair(LEFT,v
0c70: 29 2c 20 63 29 3b 0d 0a 09 09 7d 0d 0a 0d 0a 09 ), c);....}.....
0c80: 09 4c 4c 20 61 6e 73 20 3d 20 6d 66 2d 3e 63 61 .LL ans = mf->ca
0c90: 6c 63 28 42 6f 62 2c 20 41 6c 69 63 65 29 3b 0d lc(Bob, Alice);.
0ca0: 0a 09 09 64 65 6c 65 74 65 20 6d 66 3b 0d 0a 09 ...delete mf;...
0cb0: 09 72 65 74 75 72 6e 20 69 6e 74 28 61 6e 73 29 .return int(ans)
0cc0: 3b 0d 0a 09 7d 0d 0a 7d 3b 0d 0a 0d 0a 2f 2f 20 ;...}..};....//
0cd0: 42 45 47 49 4e 20 43 55 54 20 48 45 52 45 0d 0a BEGIN CUT HERE..
0ce0: 23 69 6e 63 6c 75 64 65 20 3c 63 74 69 6d 65 3e #include <ctime>
0cf0: 0d 0a 64 6f 75 62 6c 65 20 73 74 61 72 74 5f 74 ..double start_t
0d00: 69 6d 65 3b 20 73 74 72 69 6e 67 20 74 69 6d 65 ime; string time
0d10: 72 28 29 0d 0a 20 7b 20 6f 73 74 72 69 6e 67 73 r().. { ostrings
0d20: 74 72 65 61 6d 20 6f 73 3b 20 6f 73 20 3c 3c 20 tream os; os <<
0d30: 22 20 28 22 20 3c 3c 20 69 6e 74 28 28 63 6c 6f " (" << int((clo
0d40: 63 6b 28 29 2d 73 74 61 72 74 5f 74 69 6d 65 29 ck()-start_time)
0d50: 2f 43 4c 4f 43 4b 53 5f 50 45 52 5f 53 45 43 2a /CLOCKS_PER_SEC*
0d60: 31 30 30 30 29 20 3c 3c 20 22 20 6d 73 65 63 29 1000) << " msec)
0d70: 22 3b 20 72 65 74 75 72 6e 20 6f 73 2e 73 74 72 "; return os.str
0d80: 28 29 3b 20 7d 0d 0a 74 65 6d 70 6c 61 74 65 3c (); }..template<
0d90: 74 79 70 65 6e 61 6d 65 20 54 3e 20 6f 73 74 72 typename T> ostr
0da0: 65 61 6d 26 20 6f 70 65 72 61 74 6f 72 3c 3c 28 eam& operator<<(
0db0: 6f 73 74 72 65 61 6d 26 20 6f 73 2c 20 63 6f 6e ostream& os, con
0dc0: 73 74 20 76 65 63 74 6f 72 3c 54 3e 26 20 76 29 st vector<T>& v)
0dd0: 0d 0a 20 7b 20 6f 73 20 3c 3c 20 22 7b 20 22 3b .. { os << "{ ";
0de0: 0d 0a 20 20 20 66 6f 72 28 74 79 70 65 6e 61 6d .. for(typenam
0df0: 65 20 76 65 63 74 6f 72 3c 54 3e 3a 3a 63 6f 6e e vector<T>::con
0e00: 73 74 5f 69 74 65 72 61 74 6f 72 20 69 74 3d 76 st_iterator it=v
0e10: 2e 62 65 67 69 6e 28 29 3b 20 69 74 21 3d 76 2e .begin(); it!=v.
0e20: 65 6e 64 28 29 3b 20 2b 2b 69 74 29 0d 0a 20 20 end(); ++it)..
0e30: 20 6f 73 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 2a os << '\"' << *
0e40: 69 74 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 28 69 it << '\"' << (i
0e50: 74 2b 31 3d 3d 76 2e 65 6e 64 28 29 20 3f 20 22 t+1==v.end() ? "
0e60: 22 20 3a 20 22 2c 20 22 29 3b 20 6f 73 20 3c 3c " : ", "); os <<
0e70: 20 22 20 7d 22 3b 20 72 65 74 75 72 6e 20 6f 73 " }"; return os
0e80: 3b 20 7d 0d 0a 76 6f 69 64 20 76 65 72 69 66 79 ; }..void verify
0e90: 5f 63 61 73 65 28 63 6f 6e 73 74 20 69 6e 74 26 _case(const int&
0ea0: 20 45 78 70 65 63 74 65 64 2c 20 63 6f 6e 73 74 Expected, const
0eb0: 20 69 6e 74 26 20 52 65 63 65 69 76 65 64 29 20 int& Received)
0ec0: 7b 0d 0a 20 62 6f 6f 6c 20 6f 6b 20 3d 20 28 45 {.. bool ok = (E
0ed0: 78 70 65 63 74 65 64 20 3d 3d 20 52 65 63 65 69 xpected == Recei
0ee0: 76 65 64 29 3b 0d 0a 20 69 66 28 6f 6b 29 20 63 ved);.. if(ok) c
0ef0: 65 72 72 20 3c 3c 20 22 50 41 53 53 45 44 22 20 err << "PASSED"
0f00: 3c 3c 20 74 69 6d 65 72 28 29 20 3c 3c 20 65 6e << timer() << en
0f10: 64 6c 3b 20 20 65 6c 73 65 20 7b 20 63 65 72 72 dl; else { cerr
0f20: 20 3c 3c 20 22 46 41 49 4c 45 44 22 20 3c 3c 20 << "FAILED" <<
0f30: 74 69 6d 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b timer() << endl;
0f40: 0d 0a 20 63 65 72 72 20 3c 3c 20 22 5c 74 6f 3a .. cerr << "\to:
0f50: 20 5c 22 22 20 3c 3c 20 45 78 70 65 63 74 65 64 \"" << Expected
0f60: 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 65 6e 64 6c << '\"' << endl
0f70: 20 3c 3c 20 22 5c 74 78 3a 20 5c 22 22 20 3c 3c << "\tx: \"" <<
0f80: 20 52 65 63 65 69 76 65 64 20 3c 3c 20 27 5c 22 Received << '\"
0f90: 27 20 3c 3c 20 65 6e 64 6c 3b 20 7d 20 7d 0d 0a ' << endl; } }..
0fa0: 23 64 65 66 69 6e 65 20 43 41 53 45 28 4e 29 20 #define CASE(N)
0fb0: 7b 63 65 72 72 20 3c 3c 20 22 54 65 73 74 20 43 {cerr << "Test C
0fc0: 61 73 65 20 23 22 20 3c 3c 20 4e 20 3c 3c 20 22 ase #" << N << "
0fd0: 2e 2e 2e 22 20 3c 3c 20 66 6c 75 73 68 3b 20 73 ..." << flush; s
0fe0: 74 61 72 74 5f 74 69 6d 65 3d 63 6c 6f 63 6b 28 tart_time=clock(
0ff0: 29 3b 0d 0a 23 64 65 66 69 6e 65 20 45 4e 44 09 );..#define END.
1000: 20 76 65 72 69 66 79 5f 63 61 73 65 28 5f 2c 20 verify_case(_,
1010: 53 69 6e 67 69 6e 67 28 29 2e 73 6f 6c 76 65 28 Singing().solve(
1020: 4e 2c 20 6c 6f 77 2c 20 68 69 67 68 2c 20 70 69 N, low, high, pi
1030: 74 63 68 29 29 3b 7d 0d 0a 69 6e 74 20 6d 61 69 tch));}..int mai
1040: 6e 28 29 7b 0d 0a 0d 0a 43 41 53 45 28 30 29 0d n(){....CASE(0).
1050: 0a 09 69 6e 74 20 4e 20 3d 20 33 3b 20 0d 0a 09 ..int N = 3; ...
1060: 69 6e 74 20 6c 6f 77 20 3d 20 32 3b 20 0d 0a 09 int low = 2; ...
1070: 69 6e 74 20 68 69 67 68 20 3d 20 32 3b 20 0d 0a int high = 2; ..
1080: 09 69 6e 74 20 70 69 74 63 68 5f 5b 5d 20 3d 20 .int pitch_[] =
1090: 7b 31 2c 32 2c 33 2c 32 2c 31 2c 32 7d 3b 0d 0a {1,2,3,2,1,2};..
10a0: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 . vector <int>
10b0: 70 69 74 63 68 28 70 69 74 63 68 5f 2c 20 70 69 pitch(pitch_, pi
10c0: 74 63 68 5f 2b 73 69 7a 65 6f 66 28 70 69 74 63 tch_+sizeof(pitc
10d0: 68 5f 29 2f 73 69 7a 65 6f 66 28 2a 70 69 74 63 h_)/sizeof(*pitc
10e0: 68 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 5f 20 3d h_)); ...int _ =
10f0: 20 32 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 2; ..END..CASE(
1100: 31 29 0d 0a 09 69 6e 74 20 4e 20 3d 20 31 30 3b 1)...int N = 10;
1110: 20 0d 0a 09 69 6e 74 20 6c 6f 77 20 3d 20 33 3b ...int low = 3;
1120: 20 0d 0a 09 69 6e 74 20 68 69 67 68 20 3d 20 37 ...int high = 7
1130: 3b 20 0d 0a 09 69 6e 74 20 70 69 74 63 68 5f 5b ; ...int pitch_[
1140: 5d 20 3d 20 7b 34 2c 34 2c 35 2c 35 2c 36 2c 35 ] = {4,4,5,5,6,5
1150: 2c 33 2c 36 7d 3b 0d 0a 09 20 20 76 65 63 74 6f ,3,6};... vecto
1160: 72 20 3c 69 6e 74 3e 20 70 69 74 63 68 28 70 69 r <int> pitch(pi
1170: 74 63 68 5f 2c 20 70 69 74 63 68 5f 2b 73 69 7a tch_, pitch_+siz
1180: 65 6f 66 28 70 69 74 63 68 5f 29 2f 73 69 7a 65 eof(pitch_)/size
1190: 6f 66 28 2a 70 69 74 63 68 5f 29 29 3b 20 0d 0a of(*pitch_)); ..
11a0: 09 69 6e 74 20 5f 20 3d 20 30 3b 20 0d 0a 45 4e .int _ = 0; ..EN
11b0: 44 0d 0a 43 41 53 45 28 32 29 0d 0a 09 69 6e 74 D..CASE(2)...int
11c0: 20 4e 20 3d 20 36 3b 20 0d 0a 09 69 6e 74 20 6c N = 6; ...int l
11d0: 6f 77 20 3d 20 32 3b 20 0d 0a 09 69 6e 74 20 68 ow = 2; ...int h
11e0: 69 67 68 20 3d 20 35 3b 20 0d 0a 09 69 6e 74 20 igh = 5; ...int
11f0: 70 69 74 63 68 5f 5b 5d 20 3d 20 7b 35 2c 33 2c pitch_[] = {5,3,
1200: 31 2c 36 2c 34 2c 32 7d 3b 0d 0a 09 20 20 76 65 1,6,4,2};... ve
1210: 63 74 6f 72 20 3c 69 6e 74 3e 20 70 69 74 63 68 ctor <int> pitch
1220: 28 70 69 74 63 68 5f 2c 20 70 69 74 63 68 5f 2b (pitch_, pitch_+
1230: 73 69 7a 65 6f 66 28 70 69 74 63 68 5f 29 2f 73 sizeof(pitch_)/s
1240: 69 7a 65 6f 66 28 2a 70 69 74 63 68 5f 29 29 3b izeof(*pitch_));
1250: 20 0d 0a 09 69 6e 74 20 5f 20 3d 20 31 3b 20 0d ...int _ = 1; .
1260: 0a 45 4e 44 0d 0a 43 41 53 45 28 33 29 0d 0a 09 .END..CASE(3)...
1270: 69 6e 74 20 4e 20 3d 20 31 30 3b 20 0d 0a 09 69 int N = 10; ...i
1280: 6e 74 20 6c 6f 77 20 3d 20 34 3b 20 0d 0a 09 69 nt low = 4; ...i
1290: 6e 74 20 68 69 67 68 20 3d 20 35 3b 20 0d 0a 09 nt high = 5; ...
12a0: 69 6e 74 20 70 69 74 63 68 5f 5b 5d 20 3d 20 7b int pitch_[] = {
12b0: 31 2c 34 2c 33 2c 35 2c 32 2c 35 2c 37 2c 35 2c 1,4,3,5,2,5,7,5,
12c0: 39 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 9};... vector <
12d0: 69 6e 74 3e 20 70 69 74 63 68 28 70 69 74 63 68 int> pitch(pitch
12e0: 5f 2c 20 70 69 74 63 68 5f 2b 73 69 7a 65 6f 66 _, pitch_+sizeof
12f0: 28 70 69 74 63 68 5f 29 2f 73 69 7a 65 6f 66 28 (pitch_)/sizeof(
1300: 2a 70 69 74 63 68 5f 29 29 3b 20 0d 0a 09 69 6e *pitch_)); ...in
1310: 74 20 5f 20 3d 20 33 3b 20 0d 0a 45 4e 44 0d 0a t _ = 3; ..END..
1320: 43 41 53 45 28 34 29 0d 0a 09 69 6e 74 20 4e 20 CASE(4)...int N
1330: 3d 20 31 30 30 3b 20 0d 0a 09 69 6e 74 20 6c 6f = 100; ...int lo
1340: 77 20 3d 20 32 30 3b 20 0d 0a 09 69 6e 74 20 68 w = 20; ...int h
1350: 69 67 68 20 3d 20 38 30 3b 20 0d 0a 09 69 6e 74 igh = 80; ...int
1360: 20 70 69 74 63 68 5f 5b 5d 20 3d 20 7b 32 2c 32 pitch_[] = {2,2
1370: 37 2c 33 2c 35 33 2c 35 33 2c 35 32 2c 35 32 2c 7,3,53,53,52,52,
1380: 36 30 2c 38 35 2c 38 39 2c 31 30 30 2c 35 33 2c 60,85,89,100,53,
1390: 36 30 2c 32 2c 33 2c 35 33 2c 31 30 30 2c 38 39 60,2,3,53,100,89
13a0: 2c 34 30 2c 34 32 2c 32 2c 35 33 2c 32 2c 38 35 ,40,42,2,53,2,85
13b0: 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 };... vector <i
13c0: 6e 74 3e 20 70 69 74 63 68 28 70 69 74 63 68 5f nt> pitch(pitch_
13d0: 2c 20 70 69 74 63 68 5f 2b 73 69 7a 65 6f 66 28 , pitch_+sizeof(
13e0: 70 69 74 63 68 5f 29 2f 73 69 7a 65 6f 66 28 2a pitch_)/sizeof(*
13f0: 70 69 74 63 68 5f 29 29 3b 20 0d 0a 09 69 6e 74 pitch_)); ...int
1400: 20 5f 20 3d 20 35 3b 20 0d 0a 45 4e 44 0d 0a 2f _ = 5; ..END../
1410: 2a 0d 0a 43 41 53 45 28 35 29 0d 0a 09 69 6e 74 *..CASE(5)...int
1420: 20 4e 20 3d 20 3b 20 0d 0a 09 69 6e 74 20 6c 6f N = ; ...int lo
1430: 77 20 3d 20 3b 20 0d 0a 09 69 6e 74 20 68 69 67 w = ; ...int hig
1440: 68 20 3d 20 3b 20 0d 0a 09 69 6e 74 20 70 69 74 h = ; ...int pit
1450: 63 68 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 65 ch_[] = ;... ve
1460: 63 74 6f 72 20 3c 69 6e 74 3e 20 70 69 74 63 68 ctor <int> pitch
1470: 28 70 69 74 63 68 5f 2c 20 70 69 74 63 68 5f 2b (pitch_, pitch_+
1480: 73 69 7a 65 6f 66 28 70 69 74 63 68 5f 29 2f 73 sizeof(pitch_)/s
1490: 69 7a 65 6f 66 28 2a 70 69 74 63 68 5f 29 29 3b izeof(*pitch_));
14a0: 20 0d 0a 09 69 6e 74 20 5f 20 3d 20 3b 20 0d 0a ...int _ = ; ..
14b0: 45 4e 44 0d 0a 43 41 53 45 28 36 29 0d 0a 09 69 END..CASE(6)...i
14c0: 6e 74 20 4e 20 3d 20 3b 20 0d 0a 09 69 6e 74 20 nt N = ; ...int
14d0: 6c 6f 77 20 3d 20 3b 20 0d 0a 09 69 6e 74 20 68 low = ; ...int h
14e0: 69 67 68 20 3d 20 3b 20 0d 0a 09 69 6e 74 20 70 igh = ; ...int p
14f0: 69 74 63 68 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 itch_[] = ;...
1500: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 70 69 74 vector <int> pit
1510: 63 68 28 70 69 74 63 68 5f 2c 20 70 69 74 63 68 ch(pitch_, pitch
1520: 5f 2b 73 69 7a 65 6f 66 28 70 69 74 63 68 5f 29 _+sizeof(pitch_)
1530: 2f 73 69 7a 65 6f 66 28 2a 70 69 74 63 68 5f 29 /sizeof(*pitch_)
1540: 29 3b 20 0d 0a 09 69 6e 74 20 5f 20 3d 20 3b 20 ); ...int _ = ;
1550: 0d 0a 45 4e 44 0d 0a 2a 2f 0d 0a 7d 0d 0a 2f 2f ..END..*/..}..//
1560: 20 45 4e 44 20 43 55 54 20 48 45 52 45 0d 0a END CUT HERE..