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 6c 6f 6e 67 L;..typedef long
0170: 20 64 6f 75 62 6c 65 20 4c 44 3b 0d 0a 74 79 70 double LD;..typ
0180: 65 64 65 66 20 63 6f 6d 70 6c 65 78 3c 4c 44 3e edef complex<LD>
0190: 20 43 4d 50 3b 0d 0a 0d 0a 74 65 6d 70 6c 61 74 CMP;....templat
01a0: 65 3c 74 79 70 65 6e 61 6d 65 20 56 65 72 74 2c e<typename Vert,
01b0: 20 74 79 70 65 6e 61 6d 65 20 46 6c 6f 77 2c 20 typename Flow,
01c0: 69 6e 74 20 4e 56 3d 31 36 3e 0d 0a 63 6c 61 73 int NV=16>..clas
01d0: 73 20 4d 61 78 46 6c 6f 77 0d 0a 7b 0d 0a 09 76 s MaxFlow..{...v
01e0: 65 63 74 6f 72 3c 69 6e 74 3e 20 47 5b 4e 56 5d ector<int> G[NV]
01f0: 3b 0d 0a 09 46 6c 6f 77 20 46 5b 4e 56 5d 5b 4e ;...Flow F[NV][N
0200: 56 5d 3b 0d 0a 0d 0a 70 75 62 6c 69 63 3a 0d 0a V];....public:..
0210: 09 76 6f 69 64 20 61 64 64 45 64 67 65 28 20 56 .void addEdge( V
0220: 65 72 74 20 73 2c 20 56 65 72 74 20 74 2c 20 46 ert s, Vert t, F
0230: 6c 6f 77 20 66 20 29 0d 0a 09 7b 0d 0a 09 09 47 low f )...{....G
0240: 5b 73 5d 2e 70 75 73 68 5f 62 61 63 6b 28 74 29 [s].push_back(t)
0250: 3b 0d 0a 09 09 47 5b 74 5d 2e 70 75 73 68 5f 62 ;....G[t].push_b
0260: 61 63 6b 28 73 29 3b 0d 0a 09 09 46 5b 73 5d 5b ack(s);....F[s][
0270: 74 5d 20 3d 20 66 3b 0d 0a 09 09 46 5b 74 5d 5b t] = f;....F[t][
0280: 73 5d 20 3d 20 30 3b 0d 0a 09 7d 0d 0a 0d 0a 09 s] = 0;...}.....
0290: 46 6c 6f 77 20 63 61 6c 63 28 20 56 65 72 74 20 Flow calc( Vert
02a0: 53 2c 20 56 65 72 74 20 44 20 29 0d 0a 09 7b 0d S, Vert D )...{.
02b0: 0a 09 09 66 6f 72 28 20 46 6c 6f 77 20 74 6f 74 ...for( Flow tot
02c0: 61 6c 3d 30 20 3b 3b 20 29 20 7b 0d 0a 09 09 09 al=0 ;; ) {.....
02d0: 2f 2f 20 44 6f 20 42 46 53 20 61 6e 64 20 63 6f // Do BFS and co
02e0: 6d 70 75 74 65 20 74 68 65 20 6c 65 76 65 6c 20 mpute the level
02f0: 66 6f 72 20 65 61 63 68 20 6e 6f 64 65 2e 0d 0a for each node...
0300: 09 09 09 69 6e 74 20 4c 56 5b 4e 56 5d 20 3d 20 ...int LV[NV] =
0310: 7b 30 7d 3b 0d 0a 09 09 09 76 65 63 74 6f 72 3c {0};.....vector<
0320: 69 6e 74 3e 20 51 28 31 2c 20 53 29 3b 0d 0a 09 int> Q(1, S);...
0330: 09 09 66 6f 72 28 69 6e 74 20 6c 76 3d 31 3b 20 ..for(int lv=1;
0340: 21 51 2e 65 6d 70 74 79 28 29 3b 20 2b 2b 6c 76 !Q.empty(); ++lv
0350: 29 20 7b 0d 0a 09 09 09 09 76 65 63 74 6f 72 3c ) {......vector<
0360: 69 6e 74 3e 20 51 32 3b 0d 0a 09 09 09 09 66 6f int> Q2;......fo
0370: 72 28 73 69 7a 65 5f 74 20 69 3d 30 3b 20 69 21 r(size_t i=0; i!
0380: 3d 51 2e 73 69 7a 65 28 29 3b 20 2b 2b 69 29 20 =Q.size(); ++i)
0390: 7b 0d 0a 09 09 09 09 09 63 6f 6e 73 74 20 76 65 {.......const ve
03a0: 63 74 6f 72 3c 69 6e 74 3e 26 20 6e 65 20 3d 20 ctor<int>& ne =
03b0: 47 5b 51 5b 69 5d 5d 3b 0d 0a 09 09 09 09 09 66 G[Q[i]];.......f
03c0: 6f 72 28 73 69 7a 65 5f 74 20 6a 3d 30 3b 20 6a or(size_t j=0; j
03d0: 21 3d 6e 65 2e 73 69 7a 65 28 29 3b 20 2b 2b 6a !=ne.size(); ++j
03e0: 29 0d 0a 09 09 09 09 09 09 69 66 28 20 46 5b 51 )........if( F[Q
03f0: 5b 69 5d 5d 5b 6e 65 5b 6a 5d 5d 20 26 26 20 21 [i]][ne[j]] && !
0400: 4c 56 5b 6e 65 5b 6a 5d 5d 20 26 26 20 6e 65 5b LV[ne[j]] && ne[
0410: 6a 5d 21 3d 53 20 29 0d 0a 09 09 09 09 09 09 09 j]!=S ).........
0420: 4c 56 5b 6e 65 5b 6a 5d 5d 3d 6c 76 2c 20 51 32 LV[ne[j]]=lv, Q2
0430: 2e 70 75 73 68 5f 62 61 63 6b 28 6e 65 5b 6a 5d .push_back(ne[j]
0440: 29 3b 0d 0a 09 09 09 09 7d 0d 0a 09 09 09 09 51 );......}......Q
0450: 2e 73 77 61 70 28 51 32 29 3b 0d 0a 09 09 09 7d .swap(Q2);.....}
0460: 0d 0a 0d 0a 09 09 09 2f 2f 20 44 65 73 74 69 6e .......// Destin
0470: 61 74 69 6f 6e 20 69 73 20 6e 6f 77 20 75 6e 72 ation is now unr
0480: 65 61 63 68 61 62 6c 65 2e 20 44 6f 6e 65 2e 0d eachable. Done..
0490: 0a 09 09 09 69 66 28 20 21 4c 56 5b 44 5d 20 29 ....if( !LV[D] )
04a0: 0d 0a 09 09 09 09 72 65 74 75 72 6e 20 74 6f 74 ......return tot
04b0: 61 6c 3b 0d 0a 0d 0a 09 09 09 2f 2f 20 49 74 65 al;.......// Ite
04c0: 72 61 74 69 6e 67 20 44 46 53 2e 0d 0a 09 09 09 rating DFS......
04d0: 62 6f 6f 6c 20 62 6c 6f 63 6b 65 64 5b 4e 56 5d bool blocked[NV]
04e0: 20 3d 20 7b 7d 3b 0d 0a 09 09 09 74 6f 74 61 6c = {};.....total
04f0: 20 2b 3d 20 64 69 6e 69 63 5f 64 66 73 28 20 53 += dinic_dfs( S
0500: 2c 20 44 2c 20 4c 56 2c 20 30 78 37 66 66 66 66 , D, LV, 0x7ffff
0510: 66 66 66 2c 20 62 6c 6f 63 6b 65 64 20 29 3b 0d fff, blocked );.
0520: 0a 09 09 7d 0d 0a 09 7d 0d 0a 0d 0a 70 72 69 76 ...}...}....priv
0530: 61 74 65 3a 0d 0a 09 46 6c 6f 77 20 64 69 6e 69 ate:...Flow dini
0540: 63 5f 64 66 73 28 20 69 6e 74 20 76 2c 20 69 6e c_dfs( int v, in
0550: 74 20 44 2c 20 69 6e 74 20 4c 56 5b 5d 2c 20 46 t D, int LV[], F
0560: 6c 6f 77 20 66 6c 6f 77 5f 69 6e 2c 20 62 6f 6f low flow_in, boo
0570: 6c 20 62 6c 6f 63 6b 65 64 5b 5d 20 29 0d 0a 09 l blocked[] )...
0580: 7b 0d 0a 09 09 46 6c 6f 77 20 66 6c 6f 77 5f 6f {....Flow flow_o
0590: 75 74 20 3d 20 30 3b 0d 0a 09 09 66 6f 72 28 73 ut = 0;....for(s
05a0: 69 7a 65 5f 74 20 69 3d 30 3b 20 69 21 3d 47 5b ize_t i=0; i!=G[
05b0: 76 5d 2e 73 69 7a 65 28 29 3b 20 2b 2b 69 29 20 v].size(); ++i)
05c0: 7b 0d 0a 09 09 09 69 6e 74 20 75 20 3d 20 47 5b {.....int u = G[
05d0: 76 5d 5b 69 5d 3b 0d 0a 09 09 09 69 66 28 20 4c v][i];.....if( L
05e0: 56 5b 76 5d 2b 31 3d 3d 4c 56 5b 75 5d 20 26 26 V[v]+1==LV[u] &&
05f0: 20 46 5b 76 5d 5b 75 5d 20 29 20 7b 0d 0a 09 09 F[v][u] ) {....
0600: 09 09 46 6c 6f 77 20 66 20 3d 20 6d 69 6e 28 66 ..Flow f = min(f
0610: 6c 6f 77 5f 69 6e 2d 66 6c 6f 77 5f 6f 75 74 2c low_in-flow_out,
0620: 20 46 5b 76 5d 5b 75 5d 29 3b 0d 0a 09 09 09 09 F[v][u]);......
0630: 69 66 28 20 75 3d 3d 44 20 7c 7c 20 21 62 6c 6f if( u==D || !blo
0640: 63 6b 65 64 5b 75 5d 20 26 26 20 28 66 3d 64 69 cked[u] && (f=di
0650: 6e 69 63 5f 64 66 73 28 75 2c 44 2c 4c 56 2c 66 nic_dfs(u,D,LV,f
0660: 2c 62 6c 6f 63 6b 65 64 29 29 3e 30 20 29 20 7b ,blocked))>0 ) {
0670: 0d 0a 09 09 09 09 09 46 5b 76 5d 5b 75 5d 20 20 .......F[v][u]
0680: 2d 3d 20 66 3b 0d 0a 09 09 09 09 09 46 5b 75 5d -= f;.......F[u]
0690: 5b 76 5d 20 20 2b 3d 20 66 3b 0d 0a 09 09 09 09 [v] += f;......
06a0: 09 66 6c 6f 77 5f 6f 75 74 20 2b 3d 20 66 3b 0d .flow_out += f;.
06b0: 0a 09 09 09 09 09 69 66 28 20 66 6c 6f 77 5f 69 ......if( flow_i
06c0: 6e 20 3d 3d 20 66 6c 6f 77 5f 6f 75 74 20 29 20 n == flow_out )
06d0: 72 65 74 75 72 6e 20 66 6c 6f 77 5f 6f 75 74 3b return flow_out;
06e0: 0d 0a 09 09 09 09 7d 0d 0a 09 09 09 7d 0d 0a 09 ......}.....}...
06f0: 09 7d 0d 0a 09 09 62 6c 6f 63 6b 65 64 5b 76 5d .}....blocked[v]
0700: 20 3d 20 28 66 6c 6f 77 5f 6f 75 74 3d 3d 30 29 = (flow_out==0)
0710: 3b 0d 0a 09 09 72 65 74 75 72 6e 20 66 6c 6f 77 ;....return flow
0720: 5f 6f 75 74 3b 0d 0a 09 7d 0d 0a 7d 3b 0d 0a 0d _out;...}..};...
0730: 0a 63 6c 61 73 73 20 43 6f 73 6d 69 63 42 6c 6f .class CosmicBlo
0740: 63 6b 73 20 7b 20 70 75 62 6c 69 63 3a 0d 0a 09 cks { public:...
0750: 76 65 63 74 6f 72 3c 69 6e 74 3e 20 62 6c 6f 63 vector<int> bloc
0760: 6b 54 79 70 65 73 3b 0d 0a 09 69 6e 74 20 4e 2c kTypes;...int N,
0770: 20 6d 69 6e 57 61 79 73 2c 20 6d 61 78 57 61 79 minWays, maxWay
0780: 73 3b 0d 0a 0d 0a 09 69 6e 74 20 67 65 74 4e 75 s;.....int getNu
0790: 6d 4f 72 64 65 72 73 28 76 65 63 74 6f 72 20 3c mOrders(vector <
07a0: 69 6e 74 3e 20 62 6c 6f 63 6b 54 79 70 65 73 5f int> blockTypes_
07b0: 2c 20 69 6e 74 20 6d 69 6e 57 61 79 73 5f 2c 20 , int minWays_,
07c0: 69 6e 74 20 6d 61 78 57 61 79 73 5f 29 0d 0a 09 int maxWays_)...
07d0: 7b 0d 0a 09 09 4e 20 20 20 20 20 20 20 20 20 20 {....N
07e0: 3d 20 62 6c 6f 63 6b 54 79 70 65 73 5f 2e 73 69 = blockTypes_.si
07f0: 7a 65 28 29 3b 0d 0a 09 09 62 6c 6f 63 6b 54 79 ze();....blockTy
0800: 70 65 73 20 3d 20 62 6c 6f 63 6b 54 79 70 65 73 pes = blockTypes
0810: 5f 3b 0d 0a 09 09 6d 69 6e 57 61 79 73 20 20 20 _;....minWays
0820: 20 3d 20 6d 69 6e 57 61 79 73 5f 3b 0d 0a 09 09 = minWays_;....
0830: 6d 61 78 57 61 79 73 20 20 20 20 3d 20 6d 61 78 maxWays = max
0840: 57 61 79 73 5f 3b 0d 0a 09 09 72 65 74 75 72 6e Ways_;....return
0850: 20 74 72 79 5f 61 6c 6c 5f 6c 61 74 69 74 75 64 try_all_latitud
0860: 65 28 29 3b 0d 0a 09 7d 0d 0a 0d 0a 09 69 6e 74 e();...}.....int
0870: 20 74 72 79 5f 61 6c 6c 5f 6c 61 74 69 74 75 64 try_all_latitud
0880: 65 28 29 0d 0a 09 7b 0d 0a 09 09 76 65 63 74 6f e()...{....vecto
0890: 72 3c 69 6e 74 3e 20 6c 61 74 3b 0d 0a 09 09 72 r<int> lat;....r
08a0: 65 74 75 72 6e 20 72 65 63 5f 61 73 73 69 67 6e eturn rec_assign
08b0: 5f 6c 61 74 69 74 75 64 65 28 26 6c 61 74 2c 20 _latitude(&lat,
08c0: 28 31 3c 3c 4e 29 2d 31 2c 20 30 78 37 66 66 66 (1<<N)-1, 0x7fff
08d0: 66 66 66 66 29 3b 0d 0a 09 7d 0d 0a 0d 0a 09 69 ffff);...}.....i
08e0: 6e 74 20 72 65 63 5f 61 73 73 69 67 6e 5f 6c 61 nt rec_assign_la
08f0: 74 69 74 75 64 65 28 76 65 63 74 6f 72 3c 69 6e titude(vector<in
0900: 74 3e 2a 20 6c 61 74 2c 20 69 6e 74 20 6d 61 73 t>* lat, int mas
0910: 6b 2c 20 69 6e 74 20 70 72 65 76 5f 6c 61 74 5f k, int prev_lat_
0920: 62 6c 6f 63 6b 73 29 0d 0a 09 7b 0d 0a 09 09 69 blocks)...{....i
0930: 66 28 20 6d 61 73 6b 20 3d 3d 20 30 20 29 0d 0a f( mask == 0 )..
0940: 09 09 09 72 65 74 75 72 6e 20 74 72 79 5f 61 6c ...return try_al
0950: 6c 5f 74 6f 75 63 68 69 6e 67 5f 72 65 6c 65 61 l_touching_relea
0960: 74 69 6f 6e 28 2a 6c 61 74 29 3b 0d 0a 0d 0a 09 tion(*lat);.....
0970: 09 69 6e 74 20 74 6f 74 61 6c 20 3d 20 30 3b 0d .int total = 0;.
0980: 0a 09 09 66 6f 72 28 69 6e 74 20 6d 3d 31 3b 20 ...for(int m=1;
0990: 6d 3c 3d 6d 61 73 6b 3b 20 2b 2b 6d 29 20 69 66 m<=mask; ++m) if
09a0: 28 20 28 6d 26 6d 61 73 6b 29 3d 3d 6d 20 29 20 ( (m&mask)==m )
09b0: 7b 0d 0a 09 09 09 69 6e 74 20 73 75 6d 20 3d 20 {.....int sum =
09c0: 30 3b 0d 0a 09 09 09 66 6f 72 28 69 6e 74 20 69 0;.....for(int i
09d0: 3d 30 3b 20 28 31 3c 3c 69 29 3c 3d 6d 3b 20 2b =0; (1<<i)<=m; +
09e0: 2b 69 29 20 69 66 28 31 3c 3c 69 20 26 20 6d 29 +i) if(1<<i & m)
09f0: 0d 0a 09 09 09 09 73 75 6d 2b 3d 62 6c 6f 63 6b ......sum+=block
0a00: 54 79 70 65 73 5b 69 5d 3b 0d 0a 09 09 09 69 66 Types[i];.....if
0a10: 28 20 70 72 65 76 5f 6c 61 74 5f 62 6c 6f 63 6b ( prev_lat_block
0a20: 73 20 3e 3d 20 73 75 6d 20 29 20 7b 0d 0a 09 09 s >= sum ) {....
0a30: 09 09 6c 61 74 2d 3e 70 75 73 68 5f 62 61 63 6b ..lat->push_back
0a40: 28 6d 29 3b 0d 0a 09 09 09 09 74 6f 74 61 6c 20 (m);......total
0a50: 2b 3d 20 72 65 63 5f 61 73 73 69 67 6e 5f 6c 61 += rec_assign_la
0a60: 74 69 74 75 64 65 28 6c 61 74 2c 20 6d 61 73 6b titude(lat, mask
0a70: 26 7e 6d 2c 20 73 75 6d 29 3b 0d 0a 09 09 09 09 &~m, sum);......
0a80: 6c 61 74 2d 3e 70 6f 70 5f 62 61 63 6b 28 29 3b lat->pop_back();
0a90: 0d 0a 09 09 09 7d 0d 0a 09 09 7d 0d 0a 09 09 72 .....}....}....r
0aa0: 65 74 75 72 6e 20 74 6f 74 61 6c 3b 0d 0a 09 7d eturn total;...}
0ab0: 0d 0a 0d 0a 09 69 6e 74 20 74 72 79 5f 61 6c 6c .....int try_all
0ac0: 5f 74 6f 75 63 68 69 6e 67 5f 72 65 6c 65 61 74 _touching_releat
0ad0: 69 6f 6e 28 63 6f 6e 73 74 20 76 65 63 74 6f 72 ion(const vector
0ae0: 3c 69 6e 74 3e 26 20 6c 61 74 29 0d 0a 09 7b 0d <int>& lat)...{.
0af0: 0a 09 09 69 6e 74 20 74 6f 74 61 6c 20 3d 20 30 ...int total = 0
0b00: 3b 0d 0a 0d 0a 09 09 76 65 63 74 6f 72 3c 69 6e ;......vector<in
0b10: 74 3e 20 66 73 74 20 3d 20 65 6c 65 6d 65 6e 74 t> fst = element
0b20: 73 28 6c 61 74 5b 30 5d 29 3b 0d 0a 09 09 76 65 s(lat[0]);....ve
0b30: 63 74 6f 72 3c 20 70 61 69 72 3c 69 6e 74 2c 69 ctor< pair<int,i
0b40: 6e 74 3e 20 3e 20 61 64 6a 20 3d 20 61 64 6a 61 nt> > adj = adja
0b50: 63 65 6e 74 5f 62 6c 6f 63 6b 73 28 6c 61 74 29 cent_blocks(lat)
0b60: 3b 0d 0a 09 09 66 6f 72 28 69 6e 74 20 6d 3d 30 ;....for(int m=0
0b70: 3b 20 6d 3c 28 31 3c 3c 61 64 6a 2e 73 69 7a 65 ; m<(1<<adj.size
0b80: 28 29 29 3b 20 2b 2b 6d 29 0d 0a 09 09 7b 0d 0a ()); ++m)....{..
0b90: 09 09 09 76 65 63 74 6f 72 3c 20 76 65 63 74 6f ...vector< vecto
0ba0: 72 3c 69 6e 74 3e 20 3e 20 75 70 28 4e 29 3b 0d r<int> > up(N);.
0bb0: 0a 09 09 09 75 70 2e 70 75 73 68 5f 62 61 63 6b ....up.push_back
0bc0: 28 66 73 74 29 3b 0d 0a 09 09 09 66 6f 72 28 69 (fst);.....for(i
0bd0: 6e 74 20 69 3d 30 3b 20 28 31 3c 3c 69 29 3c 3d nt i=0; (1<<i)<=
0be0: 6d 3b 20 2b 2b 69 29 20 69 66 28 31 3c 3c 69 20 m; ++i) if(1<<i
0bf0: 26 20 6d 29 0d 0a 09 09 09 09 75 70 5b 61 64 6a & m)......up[adj
0c00: 5b 69 5d 2e 66 69 72 73 74 5d 2e 70 75 73 68 5f [i].first].push_
0c10: 62 61 63 6b 28 61 64 6a 5b 69 5d 2e 73 65 63 6f back(adj[i].seco
0c20: 6e 64 29 3b 0d 0a 0d 0a 09 09 09 69 66 28 20 21 nd);.......if( !
0c30: 70 6f 73 73 69 62 6c 65 28 75 70 29 20 29 0d 0a possible(up) )..
0c40: 09 09 09 09 63 6f 6e 74 69 6e 75 65 3b 0d 0a 09 ....continue;...
0c50: 09 09 69 6e 74 20 77 61 79 73 20 3d 20 63 61 6c ..int ways = cal
0c60: 63 5f 77 61 79 73 28 75 70 29 3b 0d 0a 09 09 09 c_ways(up);.....
0c70: 69 66 28 20 6d 69 6e 57 61 79 73 3c 3d 77 61 79 if( minWays<=way
0c80: 73 20 26 26 20 77 61 79 73 3c 3d 6d 61 78 57 61 s && ways<=maxWa
0c90: 79 73 20 29 0d 0a 09 09 09 09 2b 2b 74 6f 74 61 ys )......++tota
0ca0: 6c 3b 0d 0a 09 09 7d 0d 0a 09 09 72 65 74 75 72 l;....}....retur
0cb0: 6e 20 74 6f 74 61 6c 3b 0d 0a 09 7d 0d 0a 0d 0a n total;...}....
0cc0: 09 76 65 63 74 6f 72 3c 69 6e 74 3e 20 65 6c 65 .vector<int> ele
0cd0: 6d 65 6e 74 73 28 69 6e 74 20 6d 61 73 6b 29 0d ments(int mask).
0ce0: 0a 09 7b 0d 0a 09 09 76 65 63 74 6f 72 3c 69 6e ..{....vector<in
0cf0: 74 3e 20 72 65 73 75 6c 74 3b 0d 0a 09 09 66 6f t> result;....fo
0d00: 72 28 69 6e 74 20 69 3d 30 3b 20 28 31 3c 3c 69 r(int i=0; (1<<i
0d10: 29 3c 3d 6d 61 73 6b 3b 20 2b 2b 69 29 20 69 66 )<=mask; ++i) if
0d20: 28 31 3c 3c 69 20 26 20 6d 61 73 6b 29 0d 0a 09 (1<<i & mask)...
0d30: 09 09 72 65 73 75 6c 74 2e 70 75 73 68 5f 62 61 ..result.push_ba
0d40: 63 6b 28 69 29 3b 0d 0a 09 09 72 65 74 75 72 6e ck(i);....return
0d50: 20 72 65 73 75 6c 74 3b 0d 0a 09 7d 0d 0a 0d 0a result;...}....
0d60: 09 76 65 63 74 6f 72 3c 20 70 61 69 72 3c 69 6e .vector< pair<in
0d70: 74 2c 69 6e 74 3e 20 3e 20 61 64 6a 61 63 65 6e t,int> > adjacen
0d80: 74 5f 62 6c 6f 63 6b 73 28 63 6f 6e 73 74 20 76 t_blocks(const v
0d90: 65 63 74 6f 72 3c 69 6e 74 3e 26 20 6c 61 74 29 ector<int>& lat)
0da0: 0d 0a 09 7b 0d 0a 09 09 76 65 63 74 6f 72 3c 20 ...{....vector<
0db0: 70 61 69 72 3c 69 6e 74 2c 69 6e 74 3e 20 3e 20 pair<int,int> >
0dc0: 72 65 73 75 6c 74 3b 0d 0a 09 09 66 6f 72 28 69 result;....for(i
0dd0: 6e 74 20 69 3d 31 3b 20 69 3c 6c 61 74 2e 73 69 nt i=1; i<lat.si
0de0: 7a 65 28 29 3b 20 2b 2b 69 29 0d 0a 09 09 09 66 ze(); ++i).....f
0df0: 6f 72 28 69 6e 74 20 61 3d 30 3b 20 28 31 3c 3c or(int a=0; (1<<
0e00: 61 29 3c 3d 6c 61 74 5b 69 2d 31 5d 3b 20 2b 2b a)<=lat[i-1]; ++
0e10: 61 29 20 69 66 28 31 3c 3c 61 20 26 20 6c 61 74 a) if(1<<a & lat
0e20: 5b 69 2d 31 5d 29 0d 0a 09 09 09 66 6f 72 28 69 [i-1]).....for(i
0e30: 6e 74 20 62 3d 30 3b 20 28 31 3c 3c 62 29 3c 3d nt b=0; (1<<b)<=
0e40: 6c 61 74 5b 69 5d 3b 20 20 20 2b 2b 62 29 20 69 lat[i]; ++b) i
0e50: 66 28 31 3c 3c 62 20 26 20 6c 61 74 5b 69 5d 29 f(1<<b & lat[i])
0e60: 0d 0a 09 09 09 09 72 65 73 75 6c 74 2e 70 75 73 ......result.pus
0e70: 68 5f 62 61 63 6b 28 6d 61 6b 65 5f 70 61 69 72 h_back(make_pair
0e80: 28 61 2c 62 29 29 3b 0d 0a 09 09 72 65 74 75 72 (a,b));....retur
0e90: 6e 20 72 65 73 75 6c 74 3b 0d 0a 09 7d 0d 0a 0d n result;...}...
0ea0: 0a 09 69 6e 74 20 63 61 6c 63 5f 77 61 79 73 28 ..int calc_ways(
0eb0: 63 6f 6e 73 74 20 76 65 63 74 6f 72 3c 20 76 65 const vector< ve
0ec0: 63 74 6f 72 3c 69 6e 74 3e 20 3e 26 20 75 70 29 ctor<int> >& up)
0ed0: 0d 0a 09 7b 0d 0a 09 09 76 65 63 74 6f 72 3c 69 ...{....vector<i
0ee0: 6e 74 3e 20 67 6f 76 65 72 6e 65 64 28 4e 29 3b nt> governed(N);
0ef0: 0d 0a 09 09 66 6f 72 28 69 6e 74 20 69 3d 30 3b ....for(int i=0;
0f00: 20 69 3c 75 70 2e 73 69 7a 65 28 29 3b 20 2b 2b i<up.size(); ++
0f10: 69 29 0d 0a 09 09 09 66 6f 72 28 69 6e 74 20 6b i).....for(int k
0f20: 3d 30 3b 20 6b 3c 75 70 5b 69 5d 2e 73 69 7a 65 =0; k<up[i].size
0f30: 28 29 3b 20 2b 2b 6b 29 0d 0a 09 09 09 09 67 6f (); ++k)......go
0f40: 76 65 72 6e 65 64 5b 75 70 5b 69 5d 5b 6b 5d 5d verned[up[i][k]]
0f50: 20 2b 2b 3b 0d 0a 0d 0a 09 09 2f 2f 20 74 72 69 ++;......// tri
0f60: 76 69 61 6c 20 63 68 65 63 6b 20 32 0d 0a 09 09 vial check 2....
0f70: 66 6f 72 28 69 6e 74 20 69 3d 30 3b 20 69 3c 67 for(int i=0; i<g
0f80: 6f 76 65 72 6e 65 64 2e 73 69 7a 65 28 29 3b 20 overned.size();
0f90: 2b 2b 69 29 0d 0a 09 09 09 69 66 28 67 6f 76 65 ++i).....if(gove
0fa0: 72 6e 65 64 5b 69 5d 20 3d 3d 20 30 29 0d 0a 09 rned[i] == 0)...
0fb0: 09 09 09 72 65 74 75 72 6e 20 30 3b 0d 0a 0d 0a ...return 0;....
0fc0: 09 09 76 65 63 74 6f 72 3c 69 6e 74 3e 20 6d 65 ..vector<int> me
0fd0: 6d 6f 28 31 3c 3c 4e 2c 20 2d 31 29 3b 0d 0a 09 mo(1<<N, -1);...
0fe0: 09 72 65 74 75 72 6e 20 63 61 6c 63 5f 77 61 79 .return calc_way
0ff0: 73 5f 72 65 63 28 4e 2c 20 28 31 3c 3c 4e 29 2d s_rec(N, (1<<N)-
1000: 31 2c 20 75 70 2c 20 67 6f 76 65 72 6e 65 64 2c 1, up, governed,
1010: 20 6d 65 6d 6f 29 3b 0d 0a 09 7d 0d 0a 0d 0a 09 memo);...}.....
1020: 69 6e 74 20 63 61 6c 63 5f 77 61 79 73 5f 72 65 int calc_ways_re
1030: 63 28 69 6e 74 20 69 2c 20 69 6e 74 20 6d 61 73 c(int i, int mas
1040: 6b 2c 20 63 6f 6e 73 74 20 76 65 63 74 6f 72 3c k, const vector<
1050: 20 76 65 63 74 6f 72 3c 69 6e 74 3e 20 3e 26 20 vector<int> >&
1060: 75 70 2c 20 76 65 63 74 6f 72 3c 69 6e 74 3e 26 up, vector<int>&
1070: 20 67 6f 76 65 72 6e 65 64 2c 0d 0a 09 09 76 65 governed,....ve
1080: 63 74 6f 72 3c 69 6e 74 3e 26 20 6d 65 6d 6f 29 ctor<int>& memo)
1090: 0d 0a 09 7b 0d 0a 09 09 69 66 28 6d 61 73 6b 20 ...{....if(mask
10a0: 3d 3d 20 30 29 0d 0a 09 09 09 72 65 74 75 72 6e == 0).....return
10b0: 20 31 3b 0d 0a 09 09 69 66 28 6d 65 6d 6f 5b 6d 1;....if(memo[m
10c0: 61 73 6b 5d 20 21 3d 20 2d 31 29 0d 0a 09 09 09 ask] != -1).....
10d0: 72 65 74 75 72 6e 20 6d 65 6d 6f 5b 6d 61 73 6b return memo[mask
10e0: 5d 3b 0d 0a 0d 0a 09 09 69 6e 74 20 73 75 6d 20 ];......int sum
10f0: 3d 20 30 3b 0d 0a 0d 0a 09 09 66 6f 72 28 69 6e = 0;......for(in
1100: 74 20 6b 3d 30 3b 20 6b 3c 75 70 5b 69 5d 2e 73 t k=0; k<up[i].s
1110: 69 7a 65 28 29 3b 20 2b 2b 6b 29 20 67 6f 76 65 ize(); ++k) gove
1120: 72 6e 65 64 5b 75 70 5b 69 5d 5b 6b 5d 5d 2d 2d rned[up[i][k]]--
1130: 3b 0d 0a 09 09 7b 0d 0a 09 09 09 66 6f 72 28 69 ;....{.....for(i
1140: 6e 74 20 6b 3d 30 3b 20 6b 3c 67 6f 76 65 72 6e nt k=0; k<govern
1150: 65 64 2e 73 69 7a 65 28 29 3b 20 2b 2b 6b 29 0d ed.size(); ++k).
1160: 0a 09 09 09 09 69 66 28 21 67 6f 76 65 72 6e 65 .....if(!governe
1170: 64 5b 6b 5d 20 26 26 20 28 6d 61 73 6b 20 26 20 d[k] && (mask &
1180: 31 3c 3c 6b 29 29 0d 0a 09 09 09 09 09 73 75 6d 1<<k)).......sum
1190: 20 2b 3d 20 63 61 6c 63 5f 77 61 79 73 5f 72 65 += calc_ways_re
11a0: 63 28 6b 2c 20 6d 61 73 6b 20 26 7e 20 28 31 3c c(k, mask &~ (1<
11b0: 3c 6b 29 2c 20 75 70 2c 20 67 6f 76 65 72 6e 65 <k), up, governe
11c0: 64 2c 20 6d 65 6d 6f 29 3b 0d 0a 09 09 7d 0d 0a d, memo);....}..
11d0: 09 09 66 6f 72 28 69 6e 74 20 6b 3d 30 3b 20 6b ..for(int k=0; k
11e0: 3c 75 70 5b 69 5d 2e 73 69 7a 65 28 29 3b 20 2b <up[i].size(); +
11f0: 2b 6b 29 20 67 6f 76 65 72 6e 65 64 5b 75 70 5b +k) governed[up[
1200: 69 5d 5b 6b 5d 5d 2b 2b 3b 0d 0a 0d 0a 09 09 72 i][k]]++;......r
1210: 65 74 75 72 6e 20 6d 65 6d 6f 5b 6d 61 73 6b 5d eturn memo[mask]
1220: 20 3d 20 73 75 6d 3b 0d 0a 09 7d 0d 0a 0d 0a 09 = sum;...}.....
1230: 62 6f 6f 6c 20 70 6f 73 73 69 62 6c 65 28 63 6f bool possible(co
1240: 6e 73 74 20 76 65 63 74 6f 72 3c 20 76 65 63 74 nst vector< vect
1250: 6f 72 3c 69 6e 74 3e 20 3e 26 20 75 70 29 0d 0a or<int> >& up)..
1260: 09 7b 0d 0a 09 09 73 74 61 74 69 63 20 63 6f 6e .{....static con
1270: 73 74 20 69 6e 74 20 49 4e 46 20 3d 20 30 78 33 st int INF = 0x3
1280: 66 66 66 66 66 66 66 3b 0d 0a 0d 0a 09 09 69 6e fffffff;......in
1290: 74 20 74 6f 74 61 6c 20 3d 20 61 63 63 75 6d 75 t total = accumu
12a0: 6c 61 74 65 28 62 6c 6f 63 6b 54 79 70 65 73 2e late(blockTypes.
12b0: 62 65 67 69 6e 28 29 2c 20 62 6c 6f 63 6b 54 79 begin(), blockTy
12c0: 70 65 73 2e 65 6e 64 28 29 2c 20 30 29 3b 0d 0a pes.end(), 0);..
12d0: 0d 0a 09 09 4d 61 78 46 6c 6f 77 3c 69 6e 74 2c ....MaxFlow<int,
12e0: 20 69 6e 74 3e 20 6d 66 3b 0d 0a 09 09 76 65 63 int> mf;....vec
12f0: 74 6f 72 3c 69 6e 74 3e 20 6f 75 74 20 3d 20 62 tor<int> out = b
1300: 6c 6f 63 6b 54 79 70 65 73 2c 20 69 6e 20 3d 20 lockTypes, in =
1310: 62 6c 6f 63 6b 54 79 70 65 73 3b 0d 0a 09 09 66 blockTypes;....f
1320: 6f 72 28 69 6e 74 20 76 3d 30 3b 20 76 3c 3d 4e or(int v=0; v<=N
1330: 3b 20 2b 2b 76 29 20 7b 0d 0a 09 09 09 66 6f 72 ; ++v) {.....for
1340: 28 69 6e 74 20 69 3d 30 3b 20 69 3c 75 70 5b 76 (int i=0; i<up[v
1350: 5d 2e 73 69 7a 65 28 29 3b 20 2b 2b 69 29 20 7b ].size(); ++i) {
1360: 0d 0a 09 09 09 09 69 6e 74 20 75 20 3d 20 75 70 ......int u = up
1370: 5b 76 5d 5b 69 5d 3b 0d 0a 09 09 09 09 6d 66 2e [v][i];......mf.
1380: 61 64 64 45 64 67 65 28 76 2c 20 4e 2b 31 2b 75 addEdge(v, N+1+u
1390: 2c 20 49 4e 46 29 3b 0d 0a 09 09 09 09 69 66 28 , INF);......if(
13a0: 76 3c 4e 29 20 7b 69 6e 5b 76 5d 20 2d 3d 20 31 v<N) {in[v] -= 1
13b0: 3b 20 69 66 28 69 6e 5b 76 5d 3c 30 29 20 72 65 ; if(in[v]<0) re
13c0: 74 75 72 6e 20 66 61 6c 73 65 3b 20 7d 0d 0a 09 turn false; }...
13d0: 09 09 09 6f 75 74 5b 75 5d 20 2d 3d 20 31 3b 20 ...out[u] -= 1;
13e0: 69 66 28 6f 75 74 5b 75 5d 3c 30 29 20 72 65 74 if(out[u]<0) ret
13f0: 75 72 6e 20 66 61 6c 73 65 3b 0d 0a 09 09 09 09 urn false;......
1400: 74 6f 74 61 6c 20 2d 3d 20 31 3b 0d 0a 09 09 09 total -= 1;.....
1410: 7d 0d 0a 09 09 7d 0d 0a 09 09 66 6f 72 28 69 6e }....}....for(in
1420: 74 20 76 3d 30 3b 20 76 3c 3d 4e 3b 20 2b 2b 76 t v=0; v<=N; ++v
1430: 29 0d 0a 09 09 09 69 66 28 76 20 3d 3d 20 4e 29 ).....if(v == N)
1440: 20 7b 0d 0a 09 09 09 09 6d 66 2e 61 64 64 45 64 {......mf.addEd
1450: 67 65 28 4e 2b 31 2b 4e 2c 20 76 2c 20 49 4e 46 ge(N+1+N, v, INF
1460: 29 3b 0d 0a 09 09 09 7d 20 65 6c 73 65 20 7b 0d );.....} else {.
1470: 0a 09 09 09 09 6d 66 2e 61 64 64 45 64 67 65 28 .....mf.addEdge(
1480: 4e 2b 31 2b 4e 2c 20 76 2c 20 69 6e 5b 76 5d 29 N+1+N, v, in[v])
1490: 3b 0d 0a 09 09 09 09 6d 66 2e 61 64 64 45 64 67 ;......mf.addEdg
14a0: 65 28 4e 2b 31 2b 76 2c 20 4e 2b 31 2b 4e 2b 31 e(N+1+v, N+1+N+1
14b0: 2c 20 6f 75 74 5b 76 5d 29 3b 0d 0a 09 09 09 7d , out[v]);.....}
14c0: 0d 0a 09 09 72 65 74 75 72 6e 20 6d 66 2e 63 61 ....return mf.ca
14d0: 6c 63 28 4e 2b 31 2b 4e 2c 20 4e 2b 31 2b 4e 2b lc(N+1+N, N+1+N+
14e0: 31 29 20 3d 3d 20 74 6f 74 61 6c 3b 0d 0a 09 7d 1) == total;...}
14f0: 0d 0a 7d 3b 0d 0a 0d 0a 2f 2f 20 42 45 47 49 4e ..};....// BEGIN
1500: 20 43 55 54 20 48 45 52 45 0d 0a 23 69 6e 63 6c CUT HERE..#incl
1510: 75 64 65 20 3c 63 74 69 6d 65 3e 0d 0a 64 6f 75 ude <ctime>..dou
1520: 62 6c 65 20 73 74 61 72 74 5f 74 69 6d 65 3b 20 ble start_time;
1530: 73 74 72 69 6e 67 20 74 69 6d 65 72 28 29 0d 0a string timer()..
1540: 20 7b 20 6f 73 74 72 69 6e 67 73 74 72 65 61 6d { ostringstream
1550: 20 6f 73 3b 20 6f 73 20 3c 3c 20 22 20 28 22 20 os; os << " ("
1560: 3c 3c 20 69 6e 74 28 28 63 6c 6f 63 6b 28 29 2d << int((clock()-
1570: 73 74 61 72 74 5f 74 69 6d 65 29 2f 43 4c 4f 43 start_time)/CLOC
1580: 4b 53 5f 50 45 52 5f 53 45 43 2a 31 30 30 30 29 KS_PER_SEC*1000)
1590: 20 3c 3c 20 22 20 6d 73 65 63 29 22 3b 20 72 65 << " msec)"; re
15a0: 74 75 72 6e 20 6f 73 2e 73 74 72 28 29 3b 20 7d turn os.str(); }
15b0: 0d 0a 74 65 6d 70 6c 61 74 65 3c 74 79 70 65 6e ..template<typen
15c0: 61 6d 65 20 54 3e 20 6f 73 74 72 65 61 6d 26 20 ame T> ostream&
15d0: 6f 70 65 72 61 74 6f 72 3c 3c 28 6f 73 74 72 65 operator<<(ostre
15e0: 61 6d 26 20 6f 73 2c 20 63 6f 6e 73 74 20 76 65 am& os, const ve
15f0: 63 74 6f 72 3c 54 3e 26 20 76 29 0d 0a 20 7b 20 ctor<T>& v).. {
1600: 6f 73 20 3c 3c 20 22 7b 20 22 3b 0d 0a 20 20 20 os << "{ ";..
1610: 66 6f 72 28 74 79 70 65 6e 61 6d 65 20 76 65 63 for(typename vec
1620: 74 6f 72 3c 54 3e 3a 3a 63 6f 6e 73 74 5f 69 74 tor<T>::const_it
1630: 65 72 61 74 6f 72 20 69 74 3d 76 2e 62 65 67 69 erator it=v.begi
1640: 6e 28 29 3b 20 69 74 21 3d 76 2e 65 6e 64 28 29 n(); it!=v.end()
1650: 3b 20 2b 2b 69 74 29 0d 0a 20 20 20 6f 73 20 3c ; ++it).. os <
1660: 3c 20 27 5c 22 27 20 3c 3c 20 2a 69 74 20 3c 3c < '\"' << *it <<
1670: 20 27 5c 22 27 20 3c 3c 20 28 69 74 2b 31 3d 3d '\"' << (it+1==
1680: 76 2e 65 6e 64 28 29 20 3f 20 22 22 20 3a 20 22 v.end() ? "" : "
1690: 2c 20 22 29 3b 20 6f 73 20 3c 3c 20 22 20 7d 22 , "); os << " }"
16a0: 3b 20 72 65 74 75 72 6e 20 6f 73 3b 20 7d 0d 0a ; return os; }..
16b0: 76 6f 69 64 20 76 65 72 69 66 79 5f 63 61 73 65 void verify_case
16c0: 28 63 6f 6e 73 74 20 69 6e 74 26 20 45 78 70 65 (const int& Expe
16d0: 63 74 65 64 2c 20 63 6f 6e 73 74 20 69 6e 74 26 cted, const int&
16e0: 20 52 65 63 65 69 76 65 64 29 20 7b 0d 0a 20 62 Received) {.. b
16f0: 6f 6f 6c 20 6f 6b 20 3d 20 28 45 78 70 65 63 74 ool ok = (Expect
1700: 65 64 20 3d 3d 20 52 65 63 65 69 76 65 64 29 3b ed == Received);
1710: 0d 0a 20 69 66 28 6f 6b 29 20 63 65 72 72 20 3c .. if(ok) cerr <
1720: 3c 20 22 50 41 53 53 45 44 22 20 3c 3c 20 74 69 < "PASSED" << ti
1730: 6d 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b 20 20 mer() << endl;
1740: 65 6c 73 65 20 7b 20 63 65 72 72 20 3c 3c 20 22 else { cerr << "
1750: 46 41 49 4c 45 44 22 20 3c 3c 20 74 69 6d 65 72 FAILED" << timer
1760: 28 29 20 3c 3c 20 65 6e 64 6c 3b 0d 0a 20 63 65 () << endl;.. ce
1770: 72 72 20 3c 3c 20 22 5c 74 6f 3a 20 5c 22 22 20 rr << "\to: \""
1780: 3c 3c 20 45 78 70 65 63 74 65 64 20 3c 3c 20 27 << Expected << '
1790: 5c 22 27 20 3c 3c 20 65 6e 64 6c 20 3c 3c 20 22 \"' << endl << "
17a0: 5c 74 78 3a 20 5c 22 22 20 3c 3c 20 52 65 63 65 \tx: \"" << Rece
17b0: 69 76 65 64 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 ived << '\"' <<
17c0: 65 6e 64 6c 3b 20 7d 20 7d 0d 0a 23 64 65 66 69 endl; } }..#defi
17d0: 6e 65 20 43 41 53 45 28 4e 29 20 7b 63 65 72 72 ne CASE(N) {cerr
17e0: 20 3c 3c 20 22 54 65 73 74 20 43 61 73 65 20 23 << "Test Case #
17f0: 22 20 3c 3c 20 4e 20 3c 3c 20 22 2e 2e 2e 22 20 " << N << "..."
1800: 3c 3c 20 66 6c 75 73 68 3b 20 73 74 61 72 74 5f << flush; start_
1810: 74 69 6d 65 3d 63 6c 6f 63 6b 28 29 3b 0d 0a 23 time=clock();..#
1820: 64 65 66 69 6e 65 20 45 4e 44 09 20 76 65 72 69 define END. veri
1830: 66 79 5f 63 61 73 65 28 5f 2c 20 43 6f 73 6d 69 fy_case(_, Cosmi
1840: 63 42 6c 6f 63 6b 73 28 29 2e 67 65 74 4e 75 6d cBlocks().getNum
1850: 4f 72 64 65 72 73 28 62 6c 6f 63 6b 54 79 70 65 Orders(blockType
1860: 73 2c 20 6d 69 6e 57 61 79 73 2c 20 6d 61 78 57 s, minWays, maxW
1870: 61 79 73 29 29 3b 7d 0d 0a 69 6e 74 20 6d 61 69 ays));}..int mai
1880: 6e 28 29 7b 0d 0a 0d 0a 43 41 53 45 28 30 29 0d n(){....CASE(0).
1890: 0a 09 69 6e 74 20 62 6c 6f 63 6b 54 79 70 65 73 ..int blockTypes
18a0: 5f 5b 5d 20 3d 20 7b 32 2c 32 2c 32 7d 3b 0d 0a _[] = {2,2,2};..
18b0: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 . vector <int>
18c0: 62 6c 6f 63 6b 54 79 70 65 73 28 62 6c 6f 63 6b blockTypes(block
18d0: 54 79 70 65 73 5f 2c 20 62 6c 6f 63 6b 54 79 70 Types_, blockTyp
18e0: 65 73 5f 2b 73 69 7a 65 6f 66 28 62 6c 6f 63 6b es_+sizeof(block
18f0: 54 79 70 65 73 5f 29 2f 73 69 7a 65 6f 66 28 2a Types_)/sizeof(*
1900: 62 6c 6f 63 6b 54 79 70 65 73 5f 29 29 3b 20 0d blockTypes_)); .
1910: 0a 09 69 6e 74 20 6d 69 6e 57 61 79 73 20 3d 20 ..int minWays =
1920: 31 3b 20 0d 0a 09 69 6e 74 20 6d 61 78 57 61 79 1; ...int maxWay
1930: 73 20 3d 20 31 3b 20 0d 0a 09 69 6e 74 20 5f 20 s = 1; ...int _
1940: 3d 20 36 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 = 6; ..END..CASE
1950: 28 31 29 0d 0a 09 69 6e 74 20 62 6c 6f 63 6b 54 (1)...int blockT
1960: 79 70 65 73 5f 5b 5d 20 3d 20 7b 31 2c 31 2c 31 ypes_[] = {1,1,1
1970: 2c 31 2c 31 2c 31 7d 3b 0d 0a 09 20 20 76 65 63 ,1,1,1};... vec
1980: 74 6f 72 20 3c 69 6e 74 3e 20 62 6c 6f 63 6b 54 tor <int> blockT
1990: 79 70 65 73 28 62 6c 6f 63 6b 54 79 70 65 73 5f ypes(blockTypes_
19a0: 2c 20 62 6c 6f 63 6b 54 79 70 65 73 5f 2b 73 69 , blockTypes_+si
19b0: 7a 65 6f 66 28 62 6c 6f 63 6b 54 79 70 65 73 5f zeof(blockTypes_
19c0: 29 2f 73 69 7a 65 6f 66 28 2a 62 6c 6f 63 6b 54 )/sizeof(*blockT
19d0: 79 70 65 73 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 ypes_)); ...int
19e0: 6d 69 6e 57 61 79 73 20 3d 20 37 32 30 3b 20 0d minWays = 720; .
19f0: 0a 09 69 6e 74 20 6d 61 78 57 61 79 73 20 3d 20 ..int maxWays =
1a00: 37 32 30 3b 20 0d 0a 09 69 6e 74 20 5f 20 3d 20 720; ...int _ =
1a10: 31 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 32 1; ..END..CASE(2
1a20: 29 0d 0a 09 69 6e 74 20 62 6c 6f 63 6b 54 79 70 )...int blockTyp
1a30: 65 73 5f 5b 5d 20 3d 20 7b 32 2c 32 7d 3b 0d 0a es_[] = {2,2};..
1a40: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 . vector <int>
1a50: 62 6c 6f 63 6b 54 79 70 65 73 28 62 6c 6f 63 6b blockTypes(block
1a60: 54 79 70 65 73 5f 2c 20 62 6c 6f 63 6b 54 79 70 Types_, blockTyp
1a70: 65 73 5f 2b 73 69 7a 65 6f 66 28 62 6c 6f 63 6b es_+sizeof(block
1a80: 54 79 70 65 73 5f 29 2f 73 69 7a 65 6f 66 28 2a Types_)/sizeof(*
1a90: 62 6c 6f 63 6b 54 79 70 65 73 5f 29 29 3b 20 0d blockTypes_)); .
1aa0: 0a 09 69 6e 74 20 6d 69 6e 57 61 79 73 20 3d 20 ..int minWays =
1ab0: 31 3b 20 0d 0a 09 69 6e 74 20 6d 61 78 57 61 79 1; ...int maxWay
1ac0: 73 20 3d 20 32 3b 20 0d 0a 09 69 6e 74 20 5f 20 s = 2; ...int _
1ad0: 3d 20 33 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 = 3; ..END..CASE
1ae0: 28 33 29 0d 0a 09 69 6e 74 20 62 6c 6f 63 6b 54 (3)...int blockT
1af0: 79 70 65 73 5f 5b 5d 20 3d 20 7b 31 2c 32 7d 3b ypes_[] = {1,2};
1b00: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 ... vector <int
1b10: 3e 20 62 6c 6f 63 6b 54 79 70 65 73 28 62 6c 6f > blockTypes(blo
1b20: 63 6b 54 79 70 65 73 5f 2c 20 62 6c 6f 63 6b 54 ckTypes_, blockT
1b30: 79 70 65 73 5f 2b 73 69 7a 65 6f 66 28 62 6c 6f ypes_+sizeof(blo
1b40: 63 6b 54 79 70 65 73 5f 29 2f 73 69 7a 65 6f 66 ckTypes_)/sizeof
1b50: 28 2a 62 6c 6f 63 6b 54 79 70 65 73 5f 29 29 3b (*blockTypes_));
1b60: 20 0d 0a 09 69 6e 74 20 6d 69 6e 57 61 79 73 20 ...int minWays
1b70: 3d 20 31 3b 20 0d 0a 09 69 6e 74 20 6d 61 78 57 = 1; ...int maxW
1b80: 61 79 73 20 3d 20 32 3b 20 0d 0a 09 69 6e 74 20 ays = 2; ...int
1b90: 5f 20 3d 20 32 3b 20 0d 0a 45 4e 44 0d 0a 43 41 _ = 2; ..END..CA
1ba0: 53 45 28 34 29 0d 0a 09 69 6e 74 20 62 6c 6f 63 SE(4)...int bloc
1bb0: 6b 54 79 70 65 73 5f 5b 5d 20 3d 20 7b 31 7d 3b kTypes_[] = {1};
1bc0: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 ... vector <int
1bd0: 3e 20 62 6c 6f 63 6b 54 79 70 65 73 28 62 6c 6f > blockTypes(blo
1be0: 63 6b 54 79 70 65 73 5f 2c 20 62 6c 6f 63 6b 54 ckTypes_, blockT
1bf0: 79 70 65 73 5f 2b 73 69 7a 65 6f 66 28 62 6c 6f ypes_+sizeof(blo
1c00: 63 6b 54 79 70 65 73 5f 29 2f 73 69 7a 65 6f 66 ckTypes_)/sizeof
1c10: 28 2a 62 6c 6f 63 6b 54 79 70 65 73 5f 29 29 3b (*blockTypes_));
1c20: 20 0d 0a 09 69 6e 74 20 6d 69 6e 57 61 79 73 20 ...int minWays
1c30: 3d 20 31 3b 20 0d 0a 09 69 6e 74 20 6d 61 78 57 = 1; ...int maxW
1c40: 61 79 73 20 3d 20 31 3b 20 0d 0a 09 69 6e 74 20 ays = 1; ...int
1c50: 5f 20 3d 20 31 3b 20 0d 0a 45 4e 44 0d 0a 43 41 _ = 1; ..END..CA
1c60: 53 45 28 35 29 0d 0a 09 69 6e 74 20 62 6c 6f 63 SE(5)...int bloc
1c70: 6b 54 79 70 65 73 5f 5b 5d 20 3d 20 7b 31 2c 32 kTypes_[] = {1,2
1c80: 2c 34 2c 38 7d 3b 0d 0a 09 20 20 76 65 63 74 6f ,4,8};... vecto
1c90: 72 20 3c 69 6e 74 3e 20 62 6c 6f 63 6b 54 79 70 r <int> blockTyp
1ca0: 65 73 28 62 6c 6f 63 6b 54 79 70 65 73 5f 2c 20 es(blockTypes_,
1cb0: 62 6c 6f 63 6b 54 79 70 65 73 5f 2b 73 69 7a 65 blockTypes_+size
1cc0: 6f 66 28 62 6c 6f 63 6b 54 79 70 65 73 5f 29 2f of(blockTypes_)/
1cd0: 73 69 7a 65 6f 66 28 2a 62 6c 6f 63 6b 54 79 70 sizeof(*blockTyp
1ce0: 65 73 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 6d 69 es_)); ...int mi
1cf0: 6e 57 61 79 73 20 3d 20 35 3b 20 0d 0a 09 69 6e nWays = 5; ...in
1d00: 74 20 6d 61 78 57 61 79 73 20 3d 20 33 30 3b 20 t maxWays = 30;
1d10: 0d 0a 09 69 6e 74 20 5f 20 3d 20 32 37 3b 20 0d ...int _ = 27; .
1d20: 0a 45 4e 44 0d 0a 43 41 53 45 28 36 29 0d 0a 09 .END..CASE(6)...
1d30: 69 6e 74 20 62 6c 6f 63 6b 54 79 70 65 73 5f 5b int blockTypes_[
1d40: 5d 20 3d 20 7b 31 2c 32 2c 33 2c 34 2c 35 2c 36 ] = {1,2,3,4,5,6
1d50: 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 };... vector <i
1d60: 6e 74 3e 20 62 6c 6f 63 6b 54 79 70 65 73 28 62 nt> blockTypes(b
1d70: 6c 6f 63 6b 54 79 70 65 73 5f 2c 20 62 6c 6f 63 lockTypes_, bloc
1d80: 6b 54 79 70 65 73 5f 2b 73 69 7a 65 6f 66 28 62 kTypes_+sizeof(b
1d90: 6c 6f 63 6b 54 79 70 65 73 5f 29 2f 73 69 7a 65 lockTypes_)/size
1da0: 6f 66 28 2a 62 6c 6f 63 6b 54 79 70 65 73 5f 29 of(*blockTypes_)
1db0: 29 3b 20 0d 0a 09 69 6e 74 20 6d 69 6e 57 61 79 ); ...int minWay
1dc0: 73 20 3d 20 31 3b 20 0d 0a 09 69 6e 74 20 6d 61 s = 1; ...int ma
1dd0: 78 57 61 79 73 20 3d 20 37 32 30 3b 20 0d 0a 09 xWays = 720; ...
1de0: 69 6e 74 20 5f 20 3d 20 34 34 34 35 3b 20 0d 0a int _ = 4445; ..
1df0: 45 4e 44 0d 0a 43 41 53 45 28 37 29 0d 0a 09 69 END..CASE(7)...i
1e00: 6e 74 20 62 6c 6f 63 6b 54 79 70 65 73 5f 5b 5d nt blockTypes_[]
1e10: 20 3d 20 7b 37 35 30 30 2c 31 30 30 30 2c 37 35 = {7500,1000,75
1e20: 30 30 2c 31 30 30 30 2c 37 35 30 30 7d 3b 0d 0a 00,1000,7500};..
1e30: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 . vector <int>
1e40: 62 6c 6f 63 6b 54 79 70 65 73 28 62 6c 6f 63 6b blockTypes(block
1e50: 54 79 70 65 73 5f 2c 20 62 6c 6f 63 6b 54 79 70 Types_, blockTyp
1e60: 65 73 5f 2b 73 69 7a 65 6f 66 28 62 6c 6f 63 6b es_+sizeof(block
1e70: 54 79 70 65 73 5f 29 2f 73 69 7a 65 6f 66 28 2a Types_)/sizeof(*
1e80: 62 6c 6f 63 6b 54 79 70 65 73 5f 29 29 3b 20 0d blockTypes_)); .
1e90: 0a 09 69 6e 74 20 6d 69 6e 57 61 79 73 20 3d 20 ..int minWays =
1ea0: 38 3b 20 0d 0a 09 69 6e 74 20 6d 61 78 57 61 79 8; ...int maxWay
1eb0: 73 20 3d 20 38 38 3b 20 0d 0a 09 69 6e 74 20 5f s = 88; ...int _
1ec0: 20 3d 20 34 34 38 3b 20 0d 0a 45 4e 44 0d 0a 2f = 448; ..END../
1ed0: 2a 0d 0a 43 41 53 45 28 38 29 0d 0a 09 69 6e 74 *..CASE(8)...int
1ee0: 20 62 6c 6f 63 6b 54 79 70 65 73 5f 5b 5d 20 3d blockTypes_[] =
1ef0: 20 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 ;... vector <i
1f00: 6e 74 3e 20 62 6c 6f 63 6b 54 79 70 65 73 28 62 nt> blockTypes(b
1f10: 6c 6f 63 6b 54 79 70 65 73 5f 2c 20 62 6c 6f 63 lockTypes_, bloc
1f20: 6b 54 79 70 65 73 5f 2b 73 69 7a 65 6f 66 28 62 kTypes_+sizeof(b
1f30: 6c 6f 63 6b 54 79 70 65 73 5f 29 2f 73 69 7a 65 lockTypes_)/size
1f40: 6f 66 28 2a 62 6c 6f 63 6b 54 79 70 65 73 5f 29 of(*blockTypes_)
1f50: 29 3b 20 0d 0a 09 69 6e 74 20 6d 69 6e 57 61 79 ); ...int minWay
1f60: 73 20 3d 20 3b 20 0d 0a 09 69 6e 74 20 6d 61 78 s = ; ...int max
1f70: 57 61 79 73 20 3d 20 3b 20 0d 0a 09 69 6e 74 20 Ways = ; ...int
1f80: 5f 20 3d 20 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 _ = ; ..END..CAS
1f90: 45 28 39 29 0d 0a 09 69 6e 74 20 62 6c 6f 63 6b E(9)...int block
1fa0: 54 79 70 65 73 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 Types_[] = ;...
1fb0: 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 62 6c vector <int> bl
1fc0: 6f 63 6b 54 79 70 65 73 28 62 6c 6f 63 6b 54 79 ockTypes(blockTy
1fd0: 70 65 73 5f 2c 20 62 6c 6f 63 6b 54 79 70 65 73 pes_, blockTypes
1fe0: 5f 2b 73 69 7a 65 6f 66 28 62 6c 6f 63 6b 54 79 _+sizeof(blockTy
1ff0: 70 65 73 5f 29 2f 73 69 7a 65 6f 66 28 2a 62 6c pes_)/sizeof(*bl
2000: 6f 63 6b 54 79 70 65 73 5f 29 29 3b 20 0d 0a 09 ockTypes_)); ...
2010: 69 6e 74 20 6d 69 6e 57 61 79 73 20 3d 20 3b 20 int minWays = ;
2020: 0d 0a 09 69 6e 74 20 6d 61 78 57 61 79 73 20 3d ...int maxWays =
2030: 20 3b 20 0d 0a 09 69 6e 74 20 5f 20 3d 20 3b 20 ; ...int _ = ;
2040: 0d 0a 45 4e 44 0d 0a 2a 2f 0d 0a 7d 0d 0a 2f 2f ..END..*/..}..//
2050: 20 45 4e 44 20 43 55 54 20 48 45 52 45 0d 0a END CUT HERE..