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 52 6f CMP;....class Ro
01a0: 61 64 73 4f 66 4b 69 6e 67 64 6f 6d 20 7b 20 70 adsOfKingdom { p
01b0: 75 62 6c 69 63 3a 0d 0a 09 76 65 63 74 6f 72 3c ublic:...vector<
01c0: 20 76 65 63 74 6f 72 3c 64 6f 75 62 6c 65 3e 20 vector<double>
01d0: 3e 20 50 3b 0d 0a 0d 0a 09 64 6f 75 62 6c 65 20 > P;.....double
01e0: 67 65 74 50 72 6f 62 61 62 69 6c 69 74 79 28 76 getProbability(v
01f0: 65 63 74 6f 72 20 3c 73 74 72 69 6e 67 3e 20 72 ector <string> r
0200: 6f 61 64 73 29 20 0d 0a 09 7b 0d 0a 09 09 50 20 oads) ...{....P
0210: 20 20 20 20 3d 20 72 65 61 64 5f 69 6e 70 75 74 = read_input
0220: 28 72 6f 61 64 73 29 3b 0d 0a 09 09 6d 65 6d 6f (roads);....memo
0230: 20 20 3d 20 76 65 63 74 6f 72 3c 64 6f 75 62 6c = vector<doubl
0240: 65 3e 28 31 3c 3c 50 2e 73 69 7a 65 28 29 2c 20 e>(1<<P.size(),
0250: 2d 31 29 3b 0d 0a 09 09 6d 65 6d 6f 30 20 3d 20 -1);....memo0 =
0260: 76 65 63 74 6f 72 3c 64 6f 75 62 6c 65 3e 28 50 vector<double>(P
0270: 2e 73 69 7a 65 28 29 3c 3c 50 2e 73 69 7a 65 28 .size()<<P.size(
0280: 29 2c 20 2d 31 29 3b 0d 0a 09 09 6d 65 6d 6f 31 ), -1);....memo1
0290: 20 3d 20 76 65 63 74 6f 72 3c 64 6f 75 62 6c 65 = vector<double
02a0: 3e 28 50 2e 73 69 7a 65 28 29 3c 3c 50 2e 73 69 >(P.size()<<P.si
02b0: 7a 65 28 29 2c 20 2d 31 29 3b 0d 0a 09 09 72 65 ze(), -1);....re
02c0: 74 75 72 6e 20 70 72 6f 62 5f 74 72 65 65 28 20 turn prob_tree(
02d0: 28 31 3c 3c 50 2e 73 69 7a 65 28 29 29 2d 31 20 (1<<P.size())-1
02e0: 29 3b 0d 0a 09 7d 0d 0a 0d 0a 09 76 65 63 74 6f );...}.....vecto
02f0: 72 3c 20 76 65 63 74 6f 72 3c 64 6f 75 62 6c 65 r< vector<double
0300: 3e 20 3e 20 72 65 61 64 5f 69 6e 70 75 74 28 63 > > read_input(c
0310: 6f 6e 73 74 20 76 65 63 74 6f 72 3c 73 74 72 69 onst vector<stri
0320: 6e 67 3e 26 20 52 29 0d 0a 09 7b 0d 0a 09 09 69 ng>& R)...{....i
0330: 6e 74 20 4e 20 3d 20 52 2e 73 69 7a 65 28 29 3b nt N = R.size();
0340: 0d 0a 09 09 76 65 63 74 6f 72 3c 20 76 65 63 74 ....vector< vect
0350: 6f 72 3c 64 6f 75 62 6c 65 3e 20 3e 20 50 28 4e or<double> > P(N
0360: 2c 20 76 65 63 74 6f 72 3c 64 6f 75 62 6c 65 3e , vector<double>
0370: 28 4e 29 29 3b 0d 0a 09 09 66 6f 72 28 69 6e 74 (N));....for(int
0380: 20 69 3d 30 3b 20 69 3c 4e 3b 20 2b 2b 69 29 0d i=0; i<N; ++i).
0390: 0a 09 09 09 66 6f 72 28 69 6e 74 20 6a 3d 30 3b ....for(int j=0;
03a0: 20 6a 3c 4e 3b 20 2b 2b 6a 29 0d 0a 09 09 09 09 j<N; ++j)......
03b0: 50 5b 69 5d 5b 6a 5d 20 3d 20 28 52 5b 69 5d 5b P[i][j] = (R[i][
03c0: 6a 5d 20 2d 20 27 30 27 29 20 2f 20 38 2e 30 3b j] - '0') / 8.0;
03d0: 0d 0a 09 09 72 65 74 75 72 6e 20 50 3b 0d 0a 09 ....return P;...
03e0: 7d 0d 0a 0d 0a 09 76 65 63 74 6f 72 3c 64 6f 75 }.....vector<dou
03f0: 62 6c 65 3e 20 6d 65 6d 6f 3b 0d 0a 09 64 6f 75 ble> memo;...dou
0400: 62 6c 65 20 70 72 6f 62 5f 74 72 65 65 28 69 6e ble prob_tree(in
0410: 74 20 53 29 0d 0a 09 7b 0d 0a 09 09 69 66 28 20 t S)...{....if(
0420: 6d 65 6d 6f 5b 53 5d 20 3e 3d 20 30 20 29 0d 0a memo[S] >= 0 )..
0430: 09 09 09 72 65 74 75 72 6e 20 6d 65 6d 6f 5b 53 ...return memo[S
0440: 5d 3b 0d 0a 0d 0a 09 09 2f 2f 20 74 68 65 20 66 ];......// the f
0450: 69 72 73 74 20 61 6e 64 20 74 68 65 20 73 65 63 irst and the sec
0460: 6f 6e 64 20 73 6d 61 6c 6c 65 73 74 20 49 44 73 ond smallest IDs
0470: 20 69 6e 20 53 0d 0a 09 09 69 6e 74 20 76 2c 20 in S....int v,
0480: 77 3b 0d 0a 09 09 66 6f 72 28 76 3d 30 3b 20 20 w;....for(v=0;
0490: 20 28 31 3c 3c 76 29 3c 3d 53 3b 20 2b 2b 76 29 (1<<v)<=S; ++v)
04a0: 20 69 66 28 20 53 20 26 20 31 3c 3c 76 20 29 20 if( S & 1<<v )
04b0: 62 72 65 61 6b 3b 0d 0a 09 09 66 6f 72 28 77 3d break;....for(w=
04c0: 76 2b 31 3b 20 28 31 3c 3c 77 29 3c 3d 53 3b 20 v+1; (1<<w)<=S;
04d0: 2b 2b 77 29 20 69 66 28 20 53 20 26 20 31 3c 3c ++w) if( S & 1<<
04e0: 77 20 29 20 62 72 65 61 6b 3b 0d 0a 0d 0a 09 09 w ) break;......
04f0: 2f 2f 20 69 66 20 7c 53 7c 20 3c 20 32 20 74 68 // if |S| < 2 th
0500: 65 6e 20 69 74 20 61 6c 77 61 79 73 20 66 6f 72 en it always for
0510: 6d 73 20 61 20 74 72 65 65 0d 0a 09 09 69 66 28 ms a tree....if(
0520: 20 28 31 3c 3c 77 29 20 3e 20 53 20 29 0d 0a 09 (1<<w) > S )...
0530: 09 09 72 65 74 75 72 6e 20 6d 65 6d 6f 5b 53 5d ..return memo[S]
0540: 20 3d 20 31 2e 30 3b 0d 0a 0d 0a 09 09 2f 2f 20 = 1.0;......//
0550: 4c 65 74 27 73 20 63 6f 6e 73 69 64 65 72 20 76 Let's consider v
0560: 20 61 73 20 74 68 65 20 22 72 6f 6f 74 20 6e 6f as the "root no
0570: 64 65 22 20 6f 66 20 53 2c 20 61 6e 64 20 74 72 de" of S, and tr
0580: 79 20 61 6c 6c 20 70 6f 73 73 69 62 6c 65 20 22 y all possible "
0590: 73 75 62 74 72 65 65 73 22 20 54 20 63 6f 6e 74 subtrees" T cont
05a0: 61 69 6e 69 6e 67 20 77 2e 0d 0a 09 09 2f 2f 20 aining w.....//
05b0: 54 68 65 20 73 69 74 75 61 74 69 6f 6e 20 69 73 The situation is
05c0: 20 28 6f 74 68 65 72 20 6e 6f 64 65 73 29 3d 76 (other nodes)=v
05d0: 2d 28 54 20 77 68 65 72 65 20 77 20 69 6e 20 54 -(T where w in T
05e0: 29 0d 0a 09 09 64 6f 75 62 6c 65 20 70 20 3d 20 )....double p =
05f0: 30 2e 30 3b 0d 0a 09 09 66 6f 72 28 69 6e 74 20 0.0;....for(int
0600: 54 3d 53 3b 20 54 3b 20 54 3d 28 54 2d 31 29 26 T=S; T; T=(T-1)&
0610: 53 29 0d 0a 09 09 09 69 66 28 20 28 54 20 26 20 S).....if( (T &
0620: 31 3c 3c 77 29 20 26 26 20 21 28 54 20 26 20 31 1<<w) && !(T & 1
0630: 3c 3c 76 29 20 29 0d 0a 09 09 09 09 70 20 2b 3d <<v) )......p +=
0640: 20 70 72 6f 62 5f 74 72 65 65 28 54 29 20 2a 20 prob_tree(T) *
0650: 70 72 6f 62 5f 74 72 65 65 28 53 26 7e 54 29 20 prob_tree(S&~T)
0660: 2a 20 70 72 6f 62 5f 73 65 70 61 72 61 74 65 64 * prob_separated
0670: 28 54 2c 20 53 26 7e 54 26 7e 28 31 3c 3c 76 29 (T, S&~T&~(1<<v)
0680: 29 20 2a 20 70 72 6f 62 5f 6f 6e 65 28 76 2c 20 ) * prob_one(v,
0690: 54 29 3b 0d 0a 09 09 72 65 74 75 72 6e 20 6d 65 T);....return me
06a0: 6d 6f 5b 53 5d 20 3d 20 70 3b 0d 0a 09 7d 0d 0a mo[S] = p;...}..
06b0: 0d 0a 09 64 6f 75 62 6c 65 20 70 72 6f 62 5f 73 ...double prob_s
06c0: 65 70 61 72 61 74 65 64 28 69 6e 74 20 53 31 2c eparated(int S1,
06d0: 20 69 6e 74 20 53 32 29 0d 0a 09 7b 0d 0a 09 09 int S2)...{....
06e0: 64 6f 75 62 6c 65 20 70 20 3d 20 31 2e 30 3b 0d double p = 1.0;.
06f0: 0a 09 09 66 6f 72 28 69 6e 74 20 76 3d 30 3b 20 ...for(int v=0;
0700: 28 31 3c 3c 76 29 3c 3d 53 31 3b 20 2b 2b 76 29 (1<<v)<=S1; ++v)
0710: 20 69 66 28 20 53 31 20 26 20 31 3c 3c 76 20 29 if( S1 & 1<<v )
0720: 0d 0a 09 09 09 70 20 2a 3d 20 70 72 6f 62 5f 7a .....p *= prob_z
0730: 65 72 6f 28 76 2c 20 53 32 29 3b 0d 0a 09 09 72 ero(v, S2);....r
0740: 65 74 75 72 6e 20 70 3b 0d 0a 09 7d 0d 0a 0d 0a eturn p;...}....
0750: 09 76 65 63 74 6f 72 3c 64 6f 75 62 6c 65 3e 20 .vector<double>
0760: 6d 65 6d 6f 30 3b 0d 0a 09 64 6f 75 62 6c 65 20 memo0;...double
0770: 70 72 6f 62 5f 7a 65 72 6f 28 69 6e 74 20 76 2c prob_zero(int v,
0780: 20 69 6e 74 20 53 29 20 2f 2f 20 30 20 63 6f 6e int S) // 0 con
0790: 6e 65 63 74 69 6f 6e 20 62 65 74 77 65 65 6e 20 nection between
07a0: 76 20 61 6e 64 20 53 0d 0a 09 7b 0d 0a 09 09 63 v and S...{....c
07b0: 6f 6e 73 74 20 69 6e 74 20 6b 65 79 20 3d 20 76 onst int key = v
07c0: 3c 3c 50 2e 73 69 7a 65 28 29 20 7c 20 53 3b 0d <<P.size() | S;.
07d0: 0a 09 09 69 66 28 20 6d 65 6d 6f 30 5b 6b 65 79 ...if( memo0[key
07e0: 5d 20 3e 3d 20 30 20 29 20 72 65 74 75 72 6e 20 ] >= 0 ) return
07f0: 6d 65 6d 6f 30 5b 6b 65 79 5d 3b 0d 0a 0d 0a 09 memo0[key];.....
0800: 09 64 6f 75 62 6c 65 20 70 20 3d 20 31 2e 30 3b .double p = 1.0;
0810: 0d 0a 09 09 66 6f 72 28 69 6e 74 20 75 3d 30 3b ....for(int u=0;
0820: 20 28 31 3c 3c 75 29 3c 3d 53 3b 20 2b 2b 75 29 (1<<u)<=S; ++u)
0830: 20 69 66 28 20 53 20 26 20 31 3c 3c 75 20 29 0d if( S & 1<<u ).
0840: 0a 09 09 09 70 20 2a 3d 20 31 2e 30 20 2d 20 50 ....p *= 1.0 - P
0850: 5b 76 5d 5b 75 5d 3b 0d 0a 09 09 72 65 74 75 72 [v][u];....retur
0860: 6e 20 6d 65 6d 6f 30 5b 6b 65 79 5d 20 3d 20 70 n memo0[key] = p
0870: 3b 0d 0a 09 7d 0d 0a 0d 0a 09 76 65 63 74 6f 72 ;...}.....vector
0880: 3c 64 6f 75 62 6c 65 3e 20 6d 65 6d 6f 31 3b 0d <double> memo1;.
0890: 0a 09 64 6f 75 62 6c 65 20 70 72 6f 62 5f 6f 6e ..double prob_on
08a0: 65 28 69 6e 74 20 76 2c 20 69 6e 74 20 53 29 20 e(int v, int S)
08b0: 2f 2f 20 65 78 61 63 74 6c 79 20 31 20 63 6f 6e // exactly 1 con
08c0: 6e 65 63 74 69 6f 6e 20 62 65 74 77 65 65 6e 20 nection between
08d0: 76 20 61 6e 64 20 53 0d 0a 09 7b 0d 0a 09 09 63 v and S...{....c
08e0: 6f 6e 73 74 20 69 6e 74 20 6b 65 79 20 3d 20 76 onst int key = v
08f0: 3c 3c 50 2e 73 69 7a 65 28 29 20 7c 20 53 3b 0d <<P.size() | S;.
0900: 0a 09 09 69 66 28 20 6d 65 6d 6f 31 5b 6b 65 79 ...if( memo1[key
0910: 5d 20 3e 3d 20 30 20 29 20 72 65 74 75 72 6e 20 ] >= 0 ) return
0920: 6d 65 6d 6f 31 5b 6b 65 79 5d 3b 0d 0a 0d 0a 09 memo1[key];.....
0930: 09 64 6f 75 62 6c 65 20 70 20 3d 20 30 2e 30 3b .double p = 0.0;
0940: 0d 0a 09 09 66 6f 72 28 69 6e 74 20 63 3d 30 3b ....for(int c=0;
0950: 20 28 31 3c 3c 63 29 3c 3d 53 3b 20 2b 2b 63 29 (1<<c)<=S; ++c)
0960: 20 69 66 28 20 53 20 26 20 31 3c 3c 63 20 29 0d if( S & 1<<c ).
0970: 0a 09 09 7b 0d 0a 09 09 09 64 6f 75 62 6c 65 20 ...{.....double
0980: 71 20 3d 20 31 2e 30 3b 0d 0a 09 09 09 66 6f 72 q = 1.0;.....for
0990: 28 69 6e 74 20 75 3d 30 3b 20 28 31 3c 3c 75 29 (int u=0; (1<<u)
09a0: 3c 3d 53 3b 20 2b 2b 75 29 20 69 66 28 20 53 20 <=S; ++u) if( S
09b0: 26 20 31 3c 3c 75 20 29 0d 0a 09 09 09 09 71 20 & 1<<u )......q
09c0: 2a 3d 20 75 3d 3d 63 20 3f 20 50 5b 76 5d 5b 75 *= u==c ? P[v][u
09d0: 5d 20 3a 20 31 2d 50 5b 76 5d 5b 75 5d 3b 0d 0a ] : 1-P[v][u];..
09e0: 09 09 09 70 20 2b 3d 20 71 3b 0d 0a 09 09 7d 0d ...p += q;....}.
09f0: 0a 09 09 72 65 74 75 72 6e 20 6d 65 6d 6f 31 5b ...return memo1[
0a00: 6b 65 79 5d 20 3d 20 70 3b 0d 0a 09 7d 0d 0a 7d key] = p;...}..}
0a10: 3b 0d 0a 0d 0a 2f 2f 20 42 45 47 49 4e 20 43 55 ;....// BEGIN CU
0a20: 54 20 48 45 52 45 0d 0a 23 69 6e 63 6c 75 64 65 T HERE..#include
0a30: 20 3c 63 74 69 6d 65 3e 0d 0a 64 6f 75 62 6c 65 <ctime>..double
0a40: 20 73 74 61 72 74 5f 74 69 6d 65 3b 20 73 74 72 start_time; str
0a50: 69 6e 67 20 74 69 6d 65 72 28 29 0d 0a 20 7b 20 ing timer().. {
0a60: 6f 73 74 72 69 6e 67 73 74 72 65 61 6d 20 6f 73 ostringstream os
0a70: 3b 20 6f 73 20 3c 3c 20 22 20 28 22 20 3c 3c 20 ; os << " (" <<
0a80: 69 6e 74 28 28 63 6c 6f 63 6b 28 29 2d 73 74 61 int((clock()-sta
0a90: 72 74 5f 74 69 6d 65 29 2f 43 4c 4f 43 4b 53 5f rt_time)/CLOCKS_
0aa0: 50 45 52 5f 53 45 43 2a 31 30 30 30 29 20 3c 3c PER_SEC*1000) <<
0ab0: 20 22 20 6d 73 65 63 29 22 3b 20 72 65 74 75 72 " msec)"; retur
0ac0: 6e 20 6f 73 2e 73 74 72 28 29 3b 20 7d 0d 0a 74 n os.str(); }..t
0ad0: 65 6d 70 6c 61 74 65 3c 74 79 70 65 6e 61 6d 65 emplate<typename
0ae0: 20 54 3e 20 6f 73 74 72 65 61 6d 26 20 6f 70 65 T> ostream& ope
0af0: 72 61 74 6f 72 3c 3c 28 6f 73 74 72 65 61 6d 26 rator<<(ostream&
0b00: 20 6f 73 2c 20 63 6f 6e 73 74 20 76 65 63 74 6f os, const vecto
0b10: 72 3c 54 3e 26 20 76 29 0d 0a 20 7b 20 6f 73 20 r<T>& v).. { os
0b20: 3c 3c 20 22 7b 20 22 3b 0d 0a 20 20 20 66 6f 72 << "{ ";.. for
0b30: 28 74 79 70 65 6e 61 6d 65 20 76 65 63 74 6f 72 (typename vector
0b40: 3c 54 3e 3a 3a 63 6f 6e 73 74 5f 69 74 65 72 61 <T>::const_itera
0b50: 74 6f 72 20 69 74 3d 76 2e 62 65 67 69 6e 28 29 tor it=v.begin()
0b60: 3b 20 69 74 21 3d 76 2e 65 6e 64 28 29 3b 20 2b ; it!=v.end(); +
0b70: 2b 69 74 29 0d 0a 20 20 20 6f 73 20 3c 3c 20 27 +it).. os << '
0b80: 5c 22 27 20 3c 3c 20 2a 69 74 20 3c 3c 20 27 5c \"' << *it << '\
0b90: 22 27 20 3c 3c 20 28 69 74 2b 31 3d 3d 76 2e 65 "' << (it+1==v.e
0ba0: 6e 64 28 29 20 3f 20 22 22 20 3a 20 22 2c 20 22 nd() ? "" : ", "
0bb0: 29 3b 20 6f 73 20 3c 3c 20 22 20 7d 22 3b 20 72 ); os << " }"; r
0bc0: 65 74 75 72 6e 20 6f 73 3b 20 7d 0d 0a 76 6f 69 eturn os; }..voi
0bd0: 64 20 76 65 72 69 66 79 5f 63 61 73 65 28 63 6f d verify_case(co
0be0: 6e 73 74 20 64 6f 75 62 6c 65 26 20 45 78 70 65 nst double& Expe
0bf0: 63 74 65 64 2c 20 63 6f 6e 73 74 20 64 6f 75 62 cted, const doub
0c00: 6c 65 26 20 52 65 63 65 69 76 65 64 29 20 7b 0d le& Received) {.
0c10: 0a 20 62 6f 6f 6c 20 6f 6b 20 3d 20 28 61 62 73 . bool ok = (abs
0c20: 28 45 78 70 65 63 74 65 64 20 2d 20 52 65 63 65 (Expected - Rece
0c30: 69 76 65 64 29 20 3c 20 31 65 2d 39 29 3b 0d 0a ived) < 1e-9);..
0c40: 20 69 66 28 6f 6b 29 20 63 65 72 72 20 3c 3c 20 if(ok) cerr <<
0c50: 22 50 41 53 53 45 44 22 20 3c 3c 20 74 69 6d 65 "PASSED" << time
0c60: 72 28 29 20 3c 3c 20 65 6e 64 6c 3b 20 20 65 6c r() << endl; el
0c70: 73 65 20 7b 20 63 65 72 72 20 3c 3c 20 22 46 41 se { cerr << "FA
0c80: 49 4c 45 44 22 20 3c 3c 20 74 69 6d 65 72 28 29 ILED" << timer()
0c90: 20 3c 3c 20 65 6e 64 6c 3b 0d 0a 20 63 65 72 72 << endl;.. cerr
0ca0: 20 3c 3c 20 22 5c 74 6f 3a 20 5c 22 22 20 3c 3c << "\to: \"" <<
0cb0: 20 45 78 70 65 63 74 65 64 20 3c 3c 20 27 5c 22 Expected << '\"
0cc0: 27 20 3c 3c 20 65 6e 64 6c 20 3c 3c 20 22 5c 74 ' << endl << "\t
0cd0: 78 3a 20 5c 22 22 20 3c 3c 20 52 65 63 65 69 76 x: \"" << Receiv
0ce0: 65 64 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 65 6e ed << '\"' << en
0cf0: 64 6c 3b 20 7d 20 7d 0d 0a 23 64 65 66 69 6e 65 dl; } }..#define
0d00: 20 43 41 53 45 28 4e 29 20 7b 63 65 72 72 20 3c CASE(N) {cerr <
0d10: 3c 20 22 54 65 73 74 20 43 61 73 65 20 23 22 20 < "Test Case #"
0d20: 3c 3c 20 4e 20 3c 3c 20 22 2e 2e 2e 22 20 3c 3c << N << "..." <<
0d30: 20 66 6c 75 73 68 3b 20 73 74 61 72 74 5f 74 69 flush; start_ti
0d40: 6d 65 3d 63 6c 6f 63 6b 28 29 3b 0d 0a 23 64 65 me=clock();..#de
0d50: 66 69 6e 65 20 45 4e 44 09 20 76 65 72 69 66 79 fine END. verify
0d60: 5f 63 61 73 65 28 5f 2c 20 52 6f 61 64 73 4f 66 _case(_, RoadsOf
0d70: 4b 69 6e 67 64 6f 6d 28 29 2e 67 65 74 50 72 6f Kingdom().getPro
0d80: 62 61 62 69 6c 69 74 79 28 72 6f 61 64 73 29 29 bability(roads))
0d90: 3b 7d 0d 0a 69 6e 74 20 6d 61 69 6e 28 29 7b 0d ;}..int main(){.
0da0: 0a 0d 0a 43 41 53 45 28 30 29 0d 0a 09 73 74 72 ...CASE(0)...str
0db0: 69 6e 67 20 72 6f 61 64 73 5f 5b 5d 20 3d 20 7b ing roads_[] = {
0dc0: 22 30 34 22 2c 0d 0a 20 22 34 30 22 7d 3b 0d 0a "04",.. "40"};..
0dd0: 09 20 20 76 65 63 74 6f 72 20 3c 73 74 72 69 6e . vector <strin
0de0: 67 3e 20 72 6f 61 64 73 28 72 6f 61 64 73 5f 2c g> roads(roads_,
0df0: 20 72 6f 61 64 73 5f 2b 73 69 7a 65 6f 66 28 72 roads_+sizeof(r
0e00: 6f 61 64 73 5f 29 2f 73 69 7a 65 6f 66 28 2a 72 oads_)/sizeof(*r
0e10: 6f 61 64 73 5f 29 29 3b 20 0d 0a 09 64 6f 75 62 oads_)); ...doub
0e20: 6c 65 20 5f 20 3d 20 30 2e 35 3b 20 0d 0a 45 4e le _ = 0.5; ..EN
0e30: 44 0d 0a 43 41 53 45 28 31 29 0d 0a 09 73 74 72 D..CASE(1)...str
0e40: 69 6e 67 20 72 6f 61 64 73 5f 5b 5d 20 3d 20 7b ing roads_[] = {
0e50: 22 30 38 22 2c 0d 0a 20 22 38 30 22 7d 3b 0d 0a "08",.. "80"};..
0e60: 09 20 20 76 65 63 74 6f 72 20 3c 73 74 72 69 6e . vector <strin
0e70: 67 3e 20 72 6f 61 64 73 28 72 6f 61 64 73 5f 2c g> roads(roads_,
0e80: 20 72 6f 61 64 73 5f 2b 73 69 7a 65 6f 66 28 72 roads_+sizeof(r
0e90: 6f 61 64 73 5f 29 2f 73 69 7a 65 6f 66 28 2a 72 oads_)/sizeof(*r
0ea0: 6f 61 64 73 5f 29 29 3b 20 0d 0a 09 64 6f 75 62 oads_)); ...doub
0eb0: 6c 65 20 5f 20 3d 20 31 2e 30 3b 20 0d 0a 45 4e le _ = 1.0; ..EN
0ec0: 44 0d 0a 43 41 53 45 28 32 29 0d 0a 09 73 74 72 D..CASE(2)...str
0ed0: 69 6e 67 20 72 6f 61 64 73 5f 5b 5d 20 3d 20 7b ing roads_[] = {
0ee0: 22 30 30 22 2c 0d 0a 20 22 30 30 22 7d 3b 0d 0a "00",.. "00"};..
0ef0: 09 20 20 76 65 63 74 6f 72 20 3c 73 74 72 69 6e . vector <strin
0f00: 67 3e 20 72 6f 61 64 73 28 72 6f 61 64 73 5f 2c g> roads(roads_,
0f10: 20 72 6f 61 64 73 5f 2b 73 69 7a 65 6f 66 28 72 roads_+sizeof(r
0f20: 6f 61 64 73 5f 29 2f 73 69 7a 65 6f 66 28 2a 72 oads_)/sizeof(*r
0f30: 6f 61 64 73 5f 29 29 3b 20 0d 0a 09 64 6f 75 62 oads_)); ...doub
0f40: 6c 65 20 5f 20 3d 20 30 2e 30 3b 20 0d 0a 45 4e le _ = 0.0; ..EN
0f50: 44 0d 0a 43 41 53 45 28 33 29 0d 0a 09 73 74 72 D..CASE(3)...str
0f60: 69 6e 67 20 72 6f 61 64 73 5f 5b 5d 20 3d 20 7b ing roads_[] = {
0f70: 22 30 38 38 22 2c 0d 0a 20 22 38 30 38 22 2c 0d "088",.. "808",.
0f80: 0a 20 22 38 38 30 22 7d 3b 0d 0a 09 20 20 76 65 . "880"};... ve
0f90: 63 74 6f 72 20 3c 73 74 72 69 6e 67 3e 20 72 6f ctor <string> ro
0fa0: 61 64 73 28 72 6f 61 64 73 5f 2c 20 72 6f 61 64 ads(roads_, road
0fb0: 73 5f 2b 73 69 7a 65 6f 66 28 72 6f 61 64 73 5f s_+sizeof(roads_
0fc0: 29 2f 73 69 7a 65 6f 66 28 2a 72 6f 61 64 73 5f )/sizeof(*roads_
0fd0: 29 29 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 )); ...double _
0fe0: 3d 20 30 2e 30 3b 20 0d 0a 45 4e 44 0d 0a 43 41 = 0.0; ..END..CA
0ff0: 53 45 28 34 29 0d 0a 09 73 74 72 69 6e 67 20 72 SE(4)...string r
1000: 6f 61 64 73 5f 5b 5d 20 3d 20 7b 22 30 34 34 22 oads_[] = {"044"
1010: 2c 0d 0a 20 22 34 30 34 22 2c 0d 0a 20 22 34 34 ,.. "404",.. "44
1020: 30 22 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 0"};... vector
1030: 3c 73 74 72 69 6e 67 3e 20 72 6f 61 64 73 28 72 <string> roads(r
1040: 6f 61 64 73 5f 2c 20 72 6f 61 64 73 5f 2b 73 69 oads_, roads_+si
1050: 7a 65 6f 66 28 72 6f 61 64 73 5f 29 2f 73 69 7a zeof(roads_)/siz
1060: 65 6f 66 28 2a 72 6f 61 64 73 5f 29 29 3b 20 0d eof(*roads_)); .
1070: 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20 30 2e 33 ..double _ = 0.3
1080: 37 35 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 75; ..END..CASE(
1090: 35 29 0d 0a 09 73 74 72 69 6e 67 20 72 6f 61 64 5)...string road
10a0: 73 5f 5b 5d 20 3d 20 7b 22 30 37 30 31 22 2c 0d s_[] = {"0701",.
10b0: 0a 20 22 37 30 37 31 22 2c 0d 0a 20 22 30 37 30 . "7071",.. "070
10c0: 38 22 2c 0d 0a 20 22 31 31 38 30 22 7d 3b 0d 0a 8",.. "1180"};..
10d0: 09 20 20 76 65 63 74 6f 72 20 3c 73 74 72 69 6e . vector <strin
10e0: 67 3e 20 72 6f 61 64 73 28 72 6f 61 64 73 5f 2c g> roads(roads_,
10f0: 20 72 6f 61 64 73 5f 2b 73 69 7a 65 6f 66 28 72 roads_+sizeof(r
1100: 6f 61 64 73 5f 29 2f 73 69 7a 65 6f 66 28 2a 72 oads_)/sizeof(*r
1110: 6f 61 64 73 5f 29 29 3b 20 0d 0a 09 64 6f 75 62 oads_)); ...doub
1120: 6c 65 20 5f 20 3d 20 30 2e 36 32 32 33 31 34 34 le _ = 0.6223144
1130: 35 33 31 32 35 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53125; ..END..CA
1140: 53 45 28 36 29 0d 0a 09 73 74 72 69 6e 67 20 72 SE(6)...string r
1150: 6f 61 64 73 5f 5b 5d 20 3d 20 7b 22 30 36 32 32 oads_[] = {"0622
1160: 31 30 30 31 30 32 33 30 30 30 32 30 22 2c 20 22 100102300020", "
1170: 36 30 37 30 30 30 30 31 30 37 32 30 38 30 30 38 6070000107208008
1180: 22 2c 20 22 32 37 30 30 32 30 37 30 31 30 31 31 ", "270020701011
1190: 30 34 31 30 22 2c 20 22 32 30 30 30 31 38 38 30 0410", "20001880
11a0: 30 31 30 31 30 38 32 30 22 2c 20 22 31 30 32 31 01010820", "1021
11b0: 30 33 32 30 30 38 30 30 30 30 32 31 22 2c 20 22 032008000021", "
11c0: 30 30 30 38 33 30 30 30 31 33 30 30 31 30 30 30 0008300013001000
11d0: 22 2c 20 22 30 30 37 38 32 30 30 32 30 30 30 30 ", "007820020000
11e0: 30 32 30 31 22 2c 20 22 31 31 30 30 30 30 32 30 0201", "11000020
11f0: 30 37 30 30 30 30 31 31 22 2c 20 22 30 30 31 30 07000011", "0010
1200: 30 31 30 30 30 30 36 30 30 30 30 30 22 2c 20 22 010000600000", "
1210: 32 37 30 31 38 33 30 37 30 30 35 31 38 31 37 32 2701830700518172
1220: 22 2c 20 22 33 32 31 30 30 30 30 30 36 35 30 30 ", "321000006500
1230: 30 32 30 30 22 2c 20 22 30 30 31 31 30 30 30 30 0200", "00110000
1240: 30 31 30 30 32 30 38 31 22 2c 20 22 30 38 30 30 01002081", "0800
1250: 30 31 30 30 30 38 30 32 30 31 30 30 22 2c 20 22 010008020100", "
1260: 30 30 34 38 30 30 32 30 30 31 32 30 31 30 30 30 0048002001201000
1270: 22 2c 20 22 32 30 31 32 32 30 30 31 30 37 30 38 ", "201220010708
1280: 30 30 30 30 22 2c 20 22 30 38 30 30 31 30 31 31 0000", "08001011
1290: 30 32 30 31 30 30 30 30 22 7d 3b 0d 0a 09 20 20 02010000"};...
12a0: 76 65 63 74 6f 72 20 3c 73 74 72 69 6e 67 3e 20 vector <string>
12b0: 72 6f 61 64 73 28 72 6f 61 64 73 5f 2c 20 72 6f roads(roads_, ro
12c0: 61 64 73 5f 2b 73 69 7a 65 6f 66 28 72 6f 61 64 ads_+sizeof(road
12d0: 73 5f 29 2f 73 69 7a 65 6f 66 28 2a 72 6f 61 64 s_)/sizeof(*road
12e0: 73 5f 29 29 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 s_)); ...double
12f0: 5f 20 3d 20 35 2e 35 31 39 30 33 32 34 37 31 33 _ = 5.5190324713
1300: 35 38 33 34 31 45 2d 35 3b 20 0d 0a 45 4e 44 0d 58341E-5; ..END.
1310: 0a 2f 2a 0d 0a 43 41 53 45 28 37 29 0d 0a 09 73 ./*..CASE(7)...s
1320: 74 72 69 6e 67 20 72 6f 61 64 73 5f 5b 5d 20 3d tring roads_[] =
1330: 20 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 73 ;... vector <s
1340: 74 72 69 6e 67 3e 20 72 6f 61 64 73 28 72 6f 61 tring> roads(roa
1350: 64 73 5f 2c 20 72 6f 61 64 73 5f 2b 73 69 7a 65 ds_, roads_+size
1360: 6f 66 28 72 6f 61 64 73 5f 29 2f 73 69 7a 65 6f of(roads_)/sizeo
1370: 66 28 2a 72 6f 61 64 73 5f 29 29 3b 20 0d 0a 09 f(*roads_)); ...
1380: 64 6f 75 62 6c 65 20 5f 20 3d 20 3b 20 0d 0a 45 double _ = ; ..E
1390: 4e 44 0d 0a 2a 2f 0d 0a 7d 0d 0a 2f 2f 20 45 4e ND..*/..}..// EN
13a0: 44 20 43 55 54 20 48 45 52 45 0d 0a D CUT HERE..