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 63 6f 6d 70 ..#include <comp
00d0: 6c 65 78 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c lex>..#include <
00e0: 71 75 65 75 65 3e 0d 0a 23 69 6e 63 6c 75 64 65 queue>..#include
00f0: 20 3c 73 74 61 63 6b 3e 0d 0a 23 69 6e 63 6c 75 <stack>..#inclu
0100: 64 65 20 3c 63 6d 61 74 68 3e 0d 0a 23 69 6e 63 de <cmath>..#inc
0110: 6c 75 64 65 20 3c 63 61 73 73 65 72 74 3e 0d 0a lude <cassert>..
0120: 23 69 6e 63 6c 75 64 65 20 3c 63 73 74 72 69 6e #include <cstrin
0130: 67 3e 0d 0a 75 73 69 6e 67 20 6e 61 6d 65 73 70 g>..using namesp
0140: 61 63 65 20 73 74 64 3b 0d 0a 74 79 70 65 64 65 ace std;..typede
0150: 66 20 6c 6f 6e 67 20 6c 6f 6e 67 20 4c 4c 3b 0d f long long LL;.
0160: 0a 0d 0a 63 6c 61 73 73 20 54 68 72 65 65 42 75 ...class ThreeBu
0170: 73 65 73 20 7b 0d 0a 70 75 62 6c 69 63 3a 0d 0a ses {..public:..
0180: 09 64 6f 75 62 6c 65 20 67 65 74 50 72 6f 62 61 .double getProba
0190: 62 69 6c 69 74 79 28 76 65 63 74 6f 72 20 3c 69 bility(vector <i
01a0: 6e 74 3e 20 77 61 69 74 2c 20 76 65 63 74 6f 72 nt> wait, vector
01b0: 20 3c 69 6e 74 3e 20 74 72 61 76 65 6c 2c 20 69 <int> travel, i
01c0: 6e 74 20 74 69 6d 65 4c 65 66 74 29 20 0d 0a 09 nt timeLeft) ...
01d0: 7b 0d 0a 09 09 64 6f 75 62 6c 65 20 74 20 3d 20 {....double t =
01e0: 74 69 6d 65 4c 65 66 74 20 2d 20 61 63 63 75 6d timeLeft - accum
01f0: 75 6c 61 74 65 28 74 72 61 76 65 6c 2e 62 65 67 ulate(travel.beg
0200: 69 6e 28 29 2c 20 74 72 61 76 65 6c 2e 65 6e 64 in(), travel.end
0210: 28 29 2c 20 30 29 3b 0d 0a 09 09 73 6f 72 74 28 (), 0);....sort(
0220: 77 61 69 74 2e 62 65 67 69 6e 28 29 2c 20 77 61 wait.begin(), wa
0230: 69 74 2e 65 6e 64 28 29 29 3b 0d 0a 09 09 64 6f it.end());....do
0240: 75 62 6c 65 20 77 31 20 3d 20 77 61 69 74 5b 30 uble w1 = wait[0
0250: 5d 3b 0d 0a 09 09 64 6f 75 62 6c 65 20 77 32 20 ];....double w2
0260: 3d 20 77 61 69 74 5b 31 5d 3b 0d 0a 09 09 64 6f = wait[1];....do
0270: 75 62 6c 65 20 77 33 20 3d 20 77 61 69 74 5b 32 uble w3 = wait[2
0280: 5d 3b 0d 0a 09 09 2f 2f 20 5b 30 2c 77 31 5d 20 ];....// [0,w1]
0290: 2b 20 5b 30 2c 77 32 5d 20 2b 20 5b 30 2c 77 33 + [0,w2] + [0,w3
02a0: 5d 20 3c 3d 20 74 20 3f 3f 0d 0a 09 09 69 66 28 ] <= t ??....if(
02b0: 20 74 20 3e 3d 20 77 31 2b 77 32 2b 77 33 20 29 t >= w1+w2+w3 )
02c0: 20 72 65 74 75 72 6e 20 31 3b 0d 0a 09 09 69 66 return 1;....if
02d0: 28 20 74 20 3c 3d 20 30 20 29 20 20 20 20 20 20 ( t <= 0 )
02e0: 20 20 72 65 74 75 72 6e 20 30 3b 0d 0a 0d 0a 09 return 0;.....
02f0: 09 69 66 28 20 77 32 3d 3d 30 20 29 20 72 65 74 .if( w2==0 ) ret
0300: 75 72 6e 20 74 20 2f 20 77 33 3b 0d 0a 09 09 69 urn t / w3;....i
0310: 66 28 20 77 31 3d 3d 30 20 29 20 7b 0d 0a 09 09 f( w1==0 ) {....
0320: 09 64 6f 75 62 6c 65 20 76 6f 6c 20 3d 20 74 2a .double vol = t*
0330: 74 2f 32 3b 0d 0a 09 09 09 69 66 28 20 74 20 3e t/2;.....if( t >
0340: 20 77 32 20 29 20 76 6f 6c 20 2d 3d 20 28 74 2d w2 ) vol -= (t-
0350: 77 32 29 2a 28 74 2d 77 32 29 2f 32 3b 0d 0a 09 w2)*(t-w2)/2;...
0360: 09 09 69 66 28 20 74 20 3e 20 77 33 20 29 20 76 ..if( t > w3 ) v
0370: 6f 6c 20 2d 3d 20 28 74 2d 77 33 29 2a 28 74 2d ol -= (t-w3)*(t-
0380: 77 33 29 2f 32 3b 0d 0a 09 09 09 72 65 74 75 72 w3)/2;.....retur
0390: 6e 20 76 6f 6c 20 2f 20 28 77 32 2a 77 33 29 3b n vol / (w2*w3);
03a0: 0d 0a 09 09 7d 0d 0a 0d 0a 09 09 64 6f 75 62 6c ....}......doubl
03b0: 65 20 76 6f 6c 20 3d 20 74 2a 74 2a 74 2f 36 3b e vol = t*t*t/6;
03c0: 0d 0a 09 09 69 66 28 20 74 20 3e 20 77 31 20 29 ....if( t > w1 )
03d0: 20 20 20 20 76 6f 6c 20 2d 3d 20 28 74 2d 77 31 vol -= (t-w1
03e0: 29 2a 28 74 2d 77 31 29 2a 28 74 2d 77 31 29 2f )*(t-w1)*(t-w1)/
03f0: 36 3b 0d 0a 09 09 69 66 28 20 74 20 3e 20 77 32 6;....if( t > w2
0400: 20 29 20 20 20 20 76 6f 6c 20 2d 3d 20 28 74 2d ) vol -= (t-
0410: 77 32 29 2a 28 74 2d 77 32 29 2a 28 74 2d 77 32 w2)*(t-w2)*(t-w2
0420: 29 2f 36 3b 0d 0a 09 09 69 66 28 20 74 20 3e 20 )/6;....if( t >
0430: 77 33 20 29 20 20 20 20 76 6f 6c 20 2d 3d 20 28 w3 ) vol -= (
0440: 74 2d 77 33 29 2a 28 74 2d 77 33 29 2a 28 74 2d t-w3)*(t-w3)*(t-
0450: 77 33 29 2f 36 3b 0d 0a 09 09 69 66 28 20 74 20 w3)/6;....if( t
0460: 3e 20 77 31 2b 77 32 20 29 20 76 6f 6c 20 2b 3d > w1+w2 ) vol +=
0470: 20 28 74 2d 77 31 2d 77 32 29 2a 28 74 2d 77 31 (t-w1-w2)*(t-w1
0480: 2d 77 32 29 2a 28 74 2d 77 31 2d 77 32 29 2f 36 -w2)*(t-w1-w2)/6
0490: 3b 0d 0a 09 09 69 66 28 20 74 20 3e 20 77 32 2b ;....if( t > w2+
04a0: 77 33 20 29 20 76 6f 6c 20 2b 3d 20 28 74 2d 77 w3 ) vol += (t-w
04b0: 32 2d 77 33 29 2a 28 74 2d 77 32 2d 77 33 29 2a 2-w3)*(t-w2-w3)*
04c0: 28 74 2d 77 32 2d 77 33 29 2f 36 3b 0d 0a 09 09 (t-w2-w3)/6;....
04d0: 69 66 28 20 74 20 3e 20 77 33 2b 77 31 20 29 20 if( t > w3+w1 )
04e0: 76 6f 6c 20 2b 3d 20 28 74 2d 77 33 2d 77 31 29 vol += (t-w3-w1)
04f0: 2a 28 74 2d 77 33 2d 77 31 29 2a 28 74 2d 77 33 *(t-w3-w1)*(t-w3
0500: 2d 77 31 29 2f 36 3b 0d 0a 09 09 72 65 74 75 72 -w1)/6;....retur
0510: 6e 20 76 6f 6c 20 2f 20 28 77 31 2a 77 32 2a 77 n vol / (w1*w2*w
0520: 33 29 3b 0d 0a 09 7d 0d 0a 7d 3b 0d 0a 0d 0a 2f 3);...}..};..../
0530: 2f 20 42 45 47 49 4e 20 43 55 54 20 48 45 52 45 / BEGIN CUT HERE
0540: 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 63 74 69 6d ..#include <ctim
0550: 65 3e 0d 0a 64 6f 75 62 6c 65 20 73 74 61 72 74 e>..double start
0560: 5f 74 69 6d 65 3b 73 74 72 69 6e 67 20 74 69 6d _time;string tim
0570: 65 72 28 29 20 7b 20 6f 73 74 72 69 6e 67 73 74 er() { ostringst
0580: 72 65 61 6d 20 6f 73 3b 20 6f 73 20 3c 3c 20 22 ream os; os << "
0590: 20 28 22 20 3c 3c 20 69 6e 74 28 28 63 6c 6f 63 (" << int((cloc
05a0: 6b 28 29 2d 73 74 61 72 74 5f 74 69 6d 65 29 2f k()-start_time)/
05b0: 43 4c 4f 43 4b 53 5f 50 45 52 5f 53 45 43 2a 31 CLOCKS_PER_SEC*1
05c0: 30 30 30 29 20 3c 3c 20 22 20 6d 73 65 63 29 22 000) << " msec)"
05d0: 3b 20 72 65 74 75 72 6e 20 6f 73 2e 73 74 72 28 ; return os.str(
05e0: 29 3b 20 7d 0d 0a 0d 0a 74 65 6d 70 6c 61 74 65 ); }....template
05f0: 20 3c 74 79 70 65 6e 61 6d 65 20 54 3e 20 73 74 <typename T> st
0600: 72 69 6e 67 20 70 72 69 6e 74 5f 61 72 72 61 79 ring print_array
0610: 28 63 6f 6e 73 74 20 76 65 63 74 6f 72 3c 54 3e (const vector<T>
0620: 20 26 56 29 20 7b 20 6f 73 74 72 69 6e 67 73 74 &V) { ostringst
0630: 72 65 61 6d 20 6f 73 3b 20 6f 73 20 3c 3c 20 22 ream os; os << "
0640: 7b 20 22 3b 20 66 6f 72 20 28 74 79 70 65 6e 61 { "; for (typena
0650: 6d 65 20 76 65 63 74 6f 72 3c 54 3e 3a 3a 63 6f me vector<T>::co
0660: 6e 73 74 5f 69 74 65 72 61 74 6f 72 20 69 74 65 nst_iterator ite
0670: 72 20 3d 20 56 2e 62 65 67 69 6e 28 29 3b 20 69 r = V.begin(); i
0680: 74 65 72 20 21 3d 20 56 2e 65 6e 64 28 29 3b 20 ter != V.end();
0690: 2b 2b 69 74 65 72 29 20 6f 73 20 3c 3c 20 27 5c ++iter) os << '\
06a0: 22 27 20 3c 3c 20 2a 69 74 65 72 20 3c 3c 20 22 "' << *iter << "
06b0: 5c 22 2c 22 3b 20 6f 73 20 3c 3c 20 22 20 7d 22 \","; os << " }"
06c0: 3b 20 72 65 74 75 72 6e 20 6f 73 2e 73 74 72 28 ; return os.str(
06d0: 29 3b 20 7d 0d 0a 69 6e 74 20 76 65 72 69 66 79 ); }..int verify
06e0: 5f 63 61 73 65 28 63 6f 6e 73 74 20 64 6f 75 62 _case(const doub
06f0: 6c 65 20 26 45 78 70 65 63 74 65 64 2c 20 63 6f le &Expected, co
0700: 6e 73 74 20 64 6f 75 62 6c 65 20 26 52 65 63 65 nst double &Rece
0710: 69 76 65 64 29 20 7b 20 64 6f 75 62 6c 65 20 64 ived) { double d
0720: 69 66 66 20 3d 20 45 78 70 65 63 74 65 64 20 2d iff = Expected -
0730: 20 52 65 63 65 69 76 65 64 3b 20 69 66 20 28 64 Received; if (d
0740: 69 66 66 20 3c 20 30 29 20 64 69 66 66 20 3d 20 iff < 0) diff =
0750: 2d 64 69 66 66 3b 20 69 66 20 28 64 69 66 66 20 -diff; if (diff
0760: 3c 20 31 65 2d 39 29 20 63 65 72 72 20 3c 3c 20 < 1e-9) cerr <<
0770: 22 50 41 53 53 45 44 22 20 3c 3c 20 74 69 6d 65 "PASSED" << time
0780: 72 28 29 20 3c 3c 20 65 6e 64 6c 3b 20 65 6c 73 r() << endl; els
0790: 65 20 7b 20 63 65 72 72 20 3c 3c 20 22 46 41 49 e { cerr << "FAI
07a0: 4c 45 44 22 20 3c 3c 20 74 69 6d 65 72 28 29 20 LED" << timer()
07b0: 3c 3c 20 65 6e 64 6c 3b 20 63 65 72 72 20 3c 3c << endl; cerr <<
07c0: 20 22 5c 74 45 78 70 65 63 74 65 64 3a 20 5c 22 "\tExpected: \"
07d0: 22 20 3c 3c 20 45 78 70 65 63 74 65 64 20 3c 3c " << Expected <<
07e0: 20 27 5c 22 27 20 3c 3c 20 65 6e 64 6c 3b 20 63 '\"' << endl; c
07f0: 65 72 72 20 3c 3c 20 22 5c 74 52 65 63 65 69 76 err << "\tReceiv
0800: 65 64 3a 20 5c 22 22 20 3c 3c 20 52 65 63 65 69 ed: \"" << Recei
0810: 76 65 64 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 65 ved << '\"' << e
0820: 6e 64 6c 3b 20 7d 20 72 65 74 75 72 6e 20 30 3b ndl; } return 0;
0830: 7d 0d 0a 0d 0a 74 65 6d 70 6c 61 74 65 3c 69 6e }....template<in
0840: 74 20 4e 3e 20 73 74 72 75 63 74 20 43 61 73 65 t N> struct Case
0850: 5f 20 7b 20 43 61 73 65 5f 28 29 7b 73 74 61 72 _ { Case_(){star
0860: 74 5f 74 69 6d 65 3d 63 6c 6f 63 6b 28 29 3b 7d t_time=clock();}
0870: 20 7d 3b 0d 0a 63 68 61 72 20 54 65 73 74 5f 28 };..char Test_(
0880: 2e 2e 2e 29 3b 0d 0a 69 6e 74 20 54 65 73 74 5f ...);..int Test_
0890: 28 43 61 73 65 5f 3c 30 3e 29 20 7b 0d 0a 09 69 (Case_<0>) {...i
08a0: 6e 74 20 77 61 69 74 5f 5b 5d 20 3d 20 7b 30 2c nt wait_[] = {0,
08b0: 20 30 2c 20 30 7d 3b 0d 0a 09 20 20 76 65 63 74 0, 0};... vect
08c0: 6f 72 20 3c 69 6e 74 3e 20 77 61 69 74 28 77 61 or <int> wait(wa
08d0: 69 74 5f 2c 20 77 61 69 74 5f 2b 73 69 7a 65 6f it_, wait_+sizeo
08e0: 66 28 77 61 69 74 5f 29 2f 73 69 7a 65 6f 66 28 f(wait_)/sizeof(
08f0: 2a 77 61 69 74 5f 29 29 3b 20 0d 0a 09 69 6e 74 *wait_)); ...int
0900: 20 74 72 61 76 65 6c 5f 5b 5d 20 3d 20 7b 31 30 travel_[] = {10
0910: 2c 20 31 35 2c 20 31 30 7d 3b 0d 0a 09 20 20 76 , 15, 10};... v
0920: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 74 72 61 76 ector <int> trav
0930: 65 6c 28 74 72 61 76 65 6c 5f 2c 20 74 72 61 76 el(travel_, trav
0940: 65 6c 5f 2b 73 69 7a 65 6f 66 28 74 72 61 76 65 el_+sizeof(trave
0950: 6c 5f 29 2f 73 69 7a 65 6f 66 28 2a 74 72 61 76 l_)/sizeof(*trav
0960: 65 6c 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 74 69 el_)); ...int ti
0970: 6d 65 4c 65 66 74 20 3d 20 34 37 3b 20 0d 0a 09 meLeft = 47; ...
0980: 64 6f 75 62 6c 65 20 52 65 74 56 61 6c 20 3d 20 double RetVal =
0990: 31 2e 30 3b 20 0d 0a 09 72 65 74 75 72 6e 20 76 1.0; ...return v
09a0: 65 72 69 66 79 5f 63 61 73 65 28 52 65 74 56 61 erify_case(RetVa
09b0: 6c 2c 20 54 68 72 65 65 42 75 73 65 73 28 29 2e l, ThreeBuses().
09c0: 67 65 74 50 72 6f 62 61 62 69 6c 69 74 79 28 77 getProbability(w
09d0: 61 69 74 2c 20 74 72 61 76 65 6c 2c 20 74 69 6d ait, travel, tim
09e0: 65 4c 65 66 74 29 29 3b 20 7d 0d 0a 69 6e 74 20 eLeft)); }..int
09f0: 54 65 73 74 5f 28 43 61 73 65 5f 3c 31 3e 29 20 Test_(Case_<1>)
0a00: 7b 0d 0a 09 69 6e 74 20 77 61 69 74 5f 5b 5d 20 {...int wait_[]
0a10: 3d 20 7b 30 2c 20 30 2c 20 30 7d 3b 0d 0a 09 20 = {0, 0, 0};...
0a20: 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 77 61 vector <int> wa
0a30: 69 74 28 77 61 69 74 5f 2c 20 77 61 69 74 5f 2b it(wait_, wait_+
0a40: 73 69 7a 65 6f 66 28 77 61 69 74 5f 29 2f 73 69 sizeof(wait_)/si
0a50: 7a 65 6f 66 28 2a 77 61 69 74 5f 29 29 3b 20 0d zeof(*wait_)); .
0a60: 0a 09 69 6e 74 20 74 72 61 76 65 6c 5f 5b 5d 20 ..int travel_[]
0a70: 3d 20 7b 31 30 2c 20 31 35 2c 20 31 30 7d 3b 0d = {10, 15, 10};.
0a80: 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e .. vector <int>
0a90: 20 74 72 61 76 65 6c 28 74 72 61 76 65 6c 5f 2c travel(travel_,
0aa0: 20 74 72 61 76 65 6c 5f 2b 73 69 7a 65 6f 66 28 travel_+sizeof(
0ab0: 74 72 61 76 65 6c 5f 29 2f 73 69 7a 65 6f 66 28 travel_)/sizeof(
0ac0: 2a 74 72 61 76 65 6c 5f 29 29 3b 20 0d 0a 09 69 *travel_)); ...i
0ad0: 6e 74 20 74 69 6d 65 4c 65 66 74 20 3d 20 33 35 nt timeLeft = 35
0ae0: 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 52 65 74 56 ; ...double RetV
0af0: 61 6c 20 3d 20 31 2e 30 3b 20 0d 0a 09 72 65 74 al = 1.0; ...ret
0b00: 75 72 6e 20 76 65 72 69 66 79 5f 63 61 73 65 28 urn verify_case(
0b10: 52 65 74 56 61 6c 2c 20 54 68 72 65 65 42 75 73 RetVal, ThreeBus
0b20: 65 73 28 29 2e 67 65 74 50 72 6f 62 61 62 69 6c es().getProbabil
0b30: 69 74 79 28 77 61 69 74 2c 20 74 72 61 76 65 6c ity(wait, travel
0b40: 2c 20 74 69 6d 65 4c 65 66 74 29 29 3b 20 7d 0d , timeLeft)); }.
0b50: 0a 69 6e 74 20 54 65 73 74 5f 28 43 61 73 65 5f .int Test_(Case_
0b60: 3c 32 3e 29 20 7b 0d 0a 09 69 6e 74 20 77 61 69 <2>) {...int wai
0b70: 74 5f 5b 5d 20 3d 20 7b 31 2c 20 31 30 30 2c 20 t_[] = {1, 100,
0b80: 31 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 1};... vector <
0b90: 69 6e 74 3e 20 77 61 69 74 28 77 61 69 74 5f 2c int> wait(wait_,
0ba0: 20 77 61 69 74 5f 2b 73 69 7a 65 6f 66 28 77 61 wait_+sizeof(wa
0bb0: 69 74 5f 29 2f 73 69 7a 65 6f 66 28 2a 77 61 69 it_)/sizeof(*wai
0bc0: 74 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 74 72 61 t_)); ...int tra
0bd0: 76 65 6c 5f 5b 5d 20 3d 20 7b 31 30 2c 20 31 30 vel_[] = {10, 10
0be0: 2c 20 31 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f , 10};... vecto
0bf0: 72 20 3c 69 6e 74 3e 20 74 72 61 76 65 6c 28 74 r <int> travel(t
0c00: 72 61 76 65 6c 5f 2c 20 74 72 61 76 65 6c 5f 2b ravel_, travel_+
0c10: 73 69 7a 65 6f 66 28 74 72 61 76 65 6c 5f 29 2f sizeof(travel_)/
0c20: 73 69 7a 65 6f 66 28 2a 74 72 61 76 65 6c 5f 29 sizeof(*travel_)
0c30: 29 3b 20 0d 0a 09 69 6e 74 20 74 69 6d 65 4c 65 ); ...int timeLe
0c40: 66 74 20 3d 20 35 32 3b 20 0d 0a 09 64 6f 75 62 ft = 52; ...doub
0c50: 6c 65 20 52 65 74 56 61 6c 20 3d 20 30 2e 32 31 le RetVal = 0.21
0c60: 3b 20 0d 0a 09 72 65 74 75 72 6e 20 76 65 72 69 ; ...return veri
0c70: 66 79 5f 63 61 73 65 28 52 65 74 56 61 6c 2c 20 fy_case(RetVal,
0c80: 54 68 72 65 65 42 75 73 65 73 28 29 2e 67 65 74 ThreeBuses().get
0c90: 50 72 6f 62 61 62 69 6c 69 74 79 28 77 61 69 74 Probability(wait
0ca0: 2c 20 74 72 61 76 65 6c 2c 20 74 69 6d 65 4c 65 , travel, timeLe
0cb0: 66 74 29 29 3b 20 7d 0d 0a 69 6e 74 20 54 65 73 ft)); }..int Tes
0cc0: 74 5f 28 43 61 73 65 5f 3c 33 3e 29 20 7b 0d 0a t_(Case_<3>) {..
0cd0: 09 69 6e 74 20 77 61 69 74 5f 5b 5d 20 3d 20 7b .int wait_[] = {
0ce0: 31 30 30 2c 20 31 30 30 2c 20 37 30 7d 3b 0d 0a 100, 100, 70};..
0cf0: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 . vector <int>
0d00: 77 61 69 74 28 77 61 69 74 5f 2c 20 77 61 69 74 wait(wait_, wait
0d10: 5f 2b 73 69 7a 65 6f 66 28 77 61 69 74 5f 29 2f _+sizeof(wait_)/
0d20: 73 69 7a 65 6f 66 28 2a 77 61 69 74 5f 29 29 3b sizeof(*wait_));
0d30: 20 0d 0a 09 69 6e 74 20 74 72 61 76 65 6c 5f 5b ...int travel_[
0d40: 5d 20 3d 20 7b 31 2c 20 31 2c 20 31 7d 3b 0d 0a ] = {1, 1, 1};..
0d50: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 . vector <int>
0d60: 74 72 61 76 65 6c 28 74 72 61 76 65 6c 5f 2c 20 travel(travel_,
0d70: 74 72 61 76 65 6c 5f 2b 73 69 7a 65 6f 66 28 74 travel_+sizeof(t
0d80: 72 61 76 65 6c 5f 29 2f 73 69 7a 65 6f 66 28 2a ravel_)/sizeof(*
0d90: 74 72 61 76 65 6c 5f 29 29 3b 20 0d 0a 09 69 6e travel_)); ...in
0da0: 74 20 74 69 6d 65 4c 65 66 74 20 3d 20 34 37 3b t timeLeft = 47;
0db0: 20 0d 0a 09 64 6f 75 62 6c 65 20 52 65 74 56 61 ...double RetVa
0dc0: 6c 20 3d 20 30 2e 30 32 30 32 38 31 39 30 34 37 l = 0.0202819047
0dd0: 36 31 39 30 34 37 33 37 3b 20 0d 0a 09 72 65 74 61904737; ...ret
0de0: 75 72 6e 20 76 65 72 69 66 79 5f 63 61 73 65 28 urn verify_case(
0df0: 52 65 74 56 61 6c 2c 20 54 68 72 65 65 42 75 73 RetVal, ThreeBus
0e00: 65 73 28 29 2e 67 65 74 50 72 6f 62 61 62 69 6c es().getProbabil
0e10: 69 74 79 28 77 61 69 74 2c 20 74 72 61 76 65 6c ity(wait, travel
0e20: 2c 20 74 69 6d 65 4c 65 66 74 29 29 3b 20 7d 0d , timeLeft)); }.
0e30: 0a 69 6e 74 20 54 65 73 74 5f 28 43 61 73 65 5f .int Test_(Case_
0e40: 3c 34 3e 29 20 7b 0d 0a 09 69 6e 74 20 77 61 69 <4>) {...int wai
0e50: 74 5f 5b 5d 20 3d 20 7b 34 33 32 31 30 2c 20 33 t_[] = {43210, 3
0e60: 31 37 33 30 2c 20 30 7d 3b 0d 0a 09 20 20 76 65 1730, 0};... ve
0e70: 63 74 6f 72 20 3c 69 6e 74 3e 20 77 61 69 74 28 ctor <int> wait(
0e80: 77 61 69 74 5f 2c 20 77 61 69 74 5f 2b 73 69 7a wait_, wait_+siz
0e90: 65 6f 66 28 77 61 69 74 5f 29 2f 73 69 7a 65 6f eof(wait_)/sizeo
0ea0: 66 28 2a 77 61 69 74 5f 29 29 3b 20 0d 0a 09 69 f(*wait_)); ...i
0eb0: 6e 74 20 74 72 61 76 65 6c 5f 5b 5d 20 3d 20 7b nt travel_[] = {
0ec0: 31 2c 20 31 2c 20 31 7d 3b 0d 0a 09 20 20 76 65 1, 1, 1};... ve
0ed0: 63 74 6f 72 20 3c 69 6e 74 3e 20 74 72 61 76 65 ctor <int> trave
0ee0: 6c 28 74 72 61 76 65 6c 5f 2c 20 74 72 61 76 65 l(travel_, trave
0ef0: 6c 5f 2b 73 69 7a 65 6f 66 28 74 72 61 76 65 6c l_+sizeof(travel
0f00: 5f 29 2f 73 69 7a 65 6f 66 28 2a 74 72 61 76 65 _)/sizeof(*trave
0f10: 6c 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 74 69 6d l_)); ...int tim
0f20: 65 4c 65 66 74 20 3d 20 34 33 32 31 32 3b 20 0d eLeft = 43212; .
0f30: 0a 09 64 6f 75 62 6c 65 20 52 65 74 56 61 6c 20 ..double RetVal
0f40: 3d 20 30 2e 36 33 32 38 31 36 34 37 37 36 36 37 = 0.632816477667
0f50: 32 30 36 36 3b 20 0d 0a 09 72 65 74 75 72 6e 20 2066; ...return
0f60: 76 65 72 69 66 79 5f 63 61 73 65 28 52 65 74 56 verify_case(RetV
0f70: 61 6c 2c 20 54 68 72 65 65 42 75 73 65 73 28 29 al, ThreeBuses()
0f80: 2e 67 65 74 50 72 6f 62 61 62 69 6c 69 74 79 28 .getProbability(
0f90: 77 61 69 74 2c 20 74 72 61 76 65 6c 2c 20 74 69 wait, travel, ti
0fa0: 6d 65 4c 65 66 74 29 29 3b 20 7d 0d 0a 0d 0a 74 meLeft)); }....t
0fb0: 65 6d 70 6c 61 74 65 3c 69 6e 74 20 4e 3e 20 76 emplate<int N> v
0fc0: 6f 69 64 20 52 75 6e 5f 28 29 20 7b 20 63 65 72 oid Run_() { cer
0fd0: 72 20 3c 3c 20 22 54 65 73 74 20 43 61 73 65 20 r << "Test Case
0fe0: 23 22 20 3c 3c 20 4e 20 3c 3c 20 22 2e 2e 2e 22 #" << N << "..."
0ff0: 20 3c 3c 20 66 6c 75 73 68 3b 20 54 65 73 74 5f << flush; Test_
1000: 28 43 61 73 65 5f 3c 4e 3e 28 29 29 3b 20 52 75 (Case_<N>()); Ru
1010: 6e 5f 3c 73 69 7a 65 6f 66 28 54 65 73 74 5f 28 n_<sizeof(Test_(
1020: 43 61 73 65 5f 3c 4e 2b 31 3e 28 29 29 29 3d 3d Case_<N+1>()))==
1030: 31 20 3f 20 2d 31 20 3a 20 4e 2b 31 3e 28 29 3b 1 ? -1 : N+1>();
1040: 20 7d 0d 0a 74 65 6d 70 6c 61 74 65 3c 3e 20 20 }..template<>
1050: 20 20 20 20 76 6f 69 64 20 52 75 6e 5f 3c 2d 31 void Run_<-1
1060: 3e 28 29 20 7b 7d 0d 0a 69 6e 74 20 6d 61 69 6e >() {}..int main
1070: 28 29 20 7b 20 52 75 6e 5f 3c 30 3e 28 29 3b 20 () { Run_<0>();
1080: 7d 0d 0a 2f 2f 20 45 4e 44 20 43 55 54 20 48 45 }..// END CUT HE
1090: 52 45 0d 0a 0d 0a RE....