Hex Artifact Content
Not logged in

Artifact 9180662f25d62836a12c4545c257ee3a342283ae:


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 4b 74 68 50 72 6f 62  ...class KthProb
0170: 61 62 6c 65 45 6c 65 6d 65 6e 74 20 7b 0d 0a 70  ableElement {..p
0180: 75 62 6c 69 63 3a 0d 0a 20 20 64 6f 75 62 6c 65  ublic:..  double
0190: 20 70 72 6f 62 61 62 69 6c 69 74 79 28 69 6e 74   probability(int
01a0: 20 4d 2c 20 69 6e 74 20 6c 6f 77 65 72 42 6f 75   M, int lowerBou
01b0: 6e 64 2c 20 69 6e 74 20 75 70 70 65 72 42 6f 75  nd, int upperBou
01c0: 6e 64 2c 20 69 6e 74 20 4e 2c 20 69 6e 74 20 4b  nd, int N, int K
01d0: 29 20 0d 0a 20 20 7b 0d 0a 20 20 20 20 64 6f 75  ) ..  {..    dou
01e0: 62 6c 65 20 64 70 31 5b 31 30 31 5d 5b 31 30 31  ble dp1[101][101
01f0: 5d 3d 7b 30 7d 2c 20 64 70 32 5b 31 30 31 5d 5b  ]={0}, dp2[101][
0200: 31 30 31 5d 2c 20 28 2a 50 29 5b 31 30 31 5d 3d  101], (*P)[101]=
0210: 64 70 31 2c 20 28 2a 51 29 5b 31 30 31 5d 3d 64  dp1, (*Q)[101]=d
0220: 70 32 3b 0d 0a 0d 0a 20 20 20 20 50 5b 30 5d 5b  p2;....    P[0][
0230: 30 5d 20 3d 20 31 2e 30 3b 0d 0a 20 20 20 20 66  0] = 1.0;..    f
0240: 6f 72 28 69 6e 74 20 69 3d 30 3b 20 69 3c 4d 3b  or(int i=0; i<M;
0250: 20 2b 2b 69 29 0d 0a 20 20 20 20 7b 0d 0a 20 20   ++i)..    {..  
0260: 20 20 20 20 6d 65 6d 73 65 74 28 20 51 2c 20 30      memset( Q, 0
0270: 2c 20 73 69 7a 65 6f 66 28 64 6f 75 62 6c 65 29  , sizeof(double)
0280: 2a 31 30 31 2a 31 30 31 20 29 3b 0d 0a 20 20 20  *101*101 );..   
0290: 20 20 20 66 6f 72 28 69 6e 74 20 61 3d 30 3b 20     for(int a=0; 
02a0: 61 3c 3d 69 3b 20 2b 2b 61 29 0d 0a 20 20 20 20  a<=i; ++a)..    
02b0: 20 20 20 20 66 6f 72 28 69 6e 74 20 62 3d 30 3b      for(int b=0;
02c0: 20 61 2b 62 3c 3d 69 3b 20 2b 2b 62 29 0d 0a 20   a+b<=i; ++b).. 
02d0: 20 20 20 20 20 20 20 7b 0d 0a 20 20 20 20 20 20         {..      
02e0: 20 20 20 20 51 5b 61 2b 31 5d 5b 62 5d 20 2b 3d      Q[a+1][b] +=
02f0: 20 50 5b 61 5d 5b 62 5d 20 2a 20 28 4e 2d 6c 6f   P[a][b] * (N-lo
0300: 77 65 72 42 6f 75 6e 64 29 20 2f 20 28 75 70 70  werBound) / (upp
0310: 65 72 42 6f 75 6e 64 2d 6c 6f 77 65 72 42 6f 75  erBound-lowerBou
0320: 6e 64 2b 31 29 3b 0d 0a 20 20 20 20 20 20 20 20  nd+1);..        
0330: 20 20 51 5b 61 5d 5b 62 2b 31 5d 20 2b 3d 20 50    Q[a][b+1] += P
0340: 5b 61 5d 5b 62 5d 20 2a 20 31 20 2f 20 28 75 70  [a][b] * 1 / (up
0350: 70 65 72 42 6f 75 6e 64 2d 6c 6f 77 65 72 42 6f  perBound-lowerBo
0360: 75 6e 64 2b 31 29 3b 0d 0a 20 20 20 20 20 20 20  und+1);..       
0370: 20 20 20 51 5b 61 5d 5b 62 5d 20 20 20 2b 3d 20     Q[a][b]   += 
0380: 50 5b 61 5d 5b 62 5d 20 2a 20 28 75 70 70 65 72  P[a][b] * (upper
0390: 42 6f 75 6e 64 2d 4e 29 20 2f 20 28 75 70 70 65  Bound-N) / (uppe
03a0: 72 42 6f 75 6e 64 2d 6c 6f 77 65 72 42 6f 75 6e  rBound-lowerBoun
03b0: 64 2b 31 29 3b 0d 0a 20 20 20 20 20 20 20 20 7d  d+1);..        }
03c0: 0d 0a 20 20 20 20 20 20 73 77 61 70 28 50 2c 20  ..      swap(P, 
03d0: 51 29 3b 0d 0a 20 20 20 20 7d 0d 0a 0d 0a 20 20  Q);..    }....  
03e0: 20 20 64 6f 75 62 6c 65 20 70 20 3d 20 30 2e 30    double p = 0.0
03f0: 3b 0d 0a 20 20 20 20 66 6f 72 28 69 6e 74 20 61  ;..    for(int a
0400: 3d 30 3b 20 61 3c 3d 4d 3b 20 2b 2b 61 29 0d 0a  =0; a<=M; ++a)..
0410: 20 20 20 20 20 20 66 6f 72 28 69 6e 74 20 62 3d        for(int b=
0420: 30 3b 20 61 2b 62 3c 3d 4d 3b 20 2b 2b 62 29 0d  0; a+b<=M; ++b).
0430: 0a 20 20 20 20 20 20 20 20 69 66 28 20 61 3c 4b  .        if( a<K
0440: 20 26 26 20 4b 3c 3d 61 2b 62 20 29 0d 0a 20 20   && K<=a+b )..  
0450: 20 20 20 20 20 20 20 20 70 20 2b 3d 20 50 5b 61          p += P[a
0460: 5d 5b 62 5d 3b 0d 0a 20 20 20 20 72 65 74 75 72  ][b];..    retur
0470: 6e 20 70 3b 0d 0a 20 20 7d 0d 0a 7d 3b 0d 0a 0d  n p;..  }..};...
0480: 0a 2f 2f 20 42 45 47 49 4e 20 43 55 54 20 48 45  .// BEGIN CUT HE
0490: 52 45 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 63 74  RE..#include <ct
04a0: 69 6d 65 3e 0d 0a 64 6f 75 62 6c 65 20 73 74 61  ime>..double sta
04b0: 72 74 5f 74 69 6d 65 3b 73 74 72 69 6e 67 20 74  rt_time;string t
04c0: 69 6d 65 72 28 29 20 7b 20 6f 73 74 72 69 6e 67  imer() { ostring
04d0: 73 74 72 65 61 6d 20 6f 73 3b 20 6f 73 20 3c 3c  stream os; os <<
04e0: 20 22 20 28 22 20 3c 3c 20 69 6e 74 28 28 63 6c   " (" << int((cl
04f0: 6f 63 6b 28 29 2d 73 74 61 72 74 5f 74 69 6d 65  ock()-start_time
0500: 29 2f 43 4c 4f 43 4b 53 5f 50 45 52 5f 53 45 43  )/CLOCKS_PER_SEC
0510: 2a 31 30 30 30 29 20 3c 3c 20 22 20 6d 73 65 63  *1000) << " msec
0520: 29 22 3b 20 72 65 74 75 72 6e 20 6f 73 2e 73 74  )"; return os.st
0530: 72 28 29 3b 20 7d 0d 0a 0d 0a 74 65 6d 70 6c 61  r(); }....templa
0540: 74 65 20 3c 74 79 70 65 6e 61 6d 65 20 54 3e 20  te <typename T> 
0550: 73 74 72 69 6e 67 20 70 72 69 6e 74 5f 61 72 72  string print_arr
0560: 61 79 28 63 6f 6e 73 74 20 76 65 63 74 6f 72 3c  ay(const vector<
0570: 54 3e 20 26 56 29 20 7b 20 6f 73 74 72 69 6e 67  T> &V) { ostring
0580: 73 74 72 65 61 6d 20 6f 73 3b 20 6f 73 20 3c 3c  stream os; os <<
0590: 20 22 7b 20 22 3b 20 66 6f 72 20 28 74 79 70 65   "{ "; for (type
05a0: 6e 61 6d 65 20 76 65 63 74 6f 72 3c 54 3e 3a 3a  name vector<T>::
05b0: 63 6f 6e 73 74 5f 69 74 65 72 61 74 6f 72 20 69  const_iterator i
05c0: 74 65 72 20 3d 20 56 2e 62 65 67 69 6e 28 29 3b  ter = V.begin();
05d0: 20 69 74 65 72 20 21 3d 20 56 2e 65 6e 64 28 29   iter != V.end()
05e0: 3b 20 2b 2b 69 74 65 72 29 20 6f 73 20 3c 3c 20  ; ++iter) os << 
05f0: 27 5c 22 27 20 3c 3c 20 2a 69 74 65 72 20 3c 3c  '\"' << *iter <<
0600: 20 22 5c 22 2c 22 3b 20 6f 73 20 3c 3c 20 22 20   "\","; os << " 
0610: 7d 22 3b 20 72 65 74 75 72 6e 20 6f 73 2e 73 74  }"; return os.st
0620: 72 28 29 3b 20 7d 0d 0a 69 6e 74 20 76 65 72 69  r(); }..int veri
0630: 66 79 5f 63 61 73 65 28 63 6f 6e 73 74 20 64 6f  fy_case(const do
0640: 75 62 6c 65 20 26 45 78 70 65 63 74 65 64 2c 20  uble &Expected, 
0650: 63 6f 6e 73 74 20 64 6f 75 62 6c 65 20 26 52 65  const double &Re
0660: 63 65 69 76 65 64 29 20 7b 20 64 6f 75 62 6c 65  ceived) { double
0670: 20 64 69 66 66 20 3d 20 45 78 70 65 63 74 65 64   diff = Expected
0680: 20 2d 20 52 65 63 65 69 76 65 64 3b 20 69 66 20   - Received; if 
0690: 28 64 69 66 66 20 3c 20 30 29 20 64 69 66 66 20  (diff < 0) diff 
06a0: 3d 20 2d 64 69 66 66 3b 20 69 66 20 28 64 69 66  = -diff; if (dif
06b0: 66 20 3c 20 31 65 2d 39 29 20 63 65 72 72 20 3c  f < 1e-9) cerr <
06c0: 3c 20 22 50 41 53 53 45 44 22 20 3c 3c 20 74 69  < "PASSED" << ti
06d0: 6d 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b 20 65  mer() << endl; e
06e0: 6c 73 65 20 7b 20 63 65 72 72 20 3c 3c 20 22 46  lse { cerr << "F
06f0: 41 49 4c 45 44 22 20 3c 3c 20 74 69 6d 65 72 28  AILED" << timer(
0700: 29 20 3c 3c 20 65 6e 64 6c 3b 20 63 65 72 72 20  ) << endl; cerr 
0710: 3c 3c 20 22 5c 74 45 78 70 65 63 74 65 64 3a 20  << "\tExpected: 
0720: 5c 22 22 20 3c 3c 20 45 78 70 65 63 74 65 64 20  \"" << Expected 
0730: 3c 3c 20 27 5c 22 27 20 3c 3c 20 65 6e 64 6c 3b  << '\"' << endl;
0740: 20 63 65 72 72 20 3c 3c 20 22 5c 74 52 65 63 65   cerr << "\tRece
0750: 69 76 65 64 3a 20 5c 22 22 20 3c 3c 20 52 65 63  ived: \"" << Rec
0760: 65 69 76 65 64 20 3c 3c 20 27 5c 22 27 20 3c 3c  eived << '\"' <<
0770: 20 65 6e 64 6c 3b 20 7d 20 72 65 74 75 72 6e 20   endl; } return 
0780: 30 3b 7d 0d 0a 0d 0a 74 65 6d 70 6c 61 74 65 3c  0;}....template<
0790: 69 6e 74 20 4e 3e 20 73 74 72 75 63 74 20 43 61  int N> struct Ca
07a0: 73 65 5f 20 7b 20 43 61 73 65 5f 28 29 7b 73 74  se_ { Case_(){st
07b0: 61 72 74 5f 74 69 6d 65 3d 63 6c 6f 63 6b 28 29  art_time=clock()
07c0: 3b 7d 20 7d 3b 0d 0a 63 68 61 72 20 54 65 73 74  ;} };..char Test
07d0: 5f 28 2e 2e 2e 29 3b 0d 0a 69 6e 74 20 54 65 73  _(...);..int Tes
07e0: 74 5f 28 43 61 73 65 5f 3c 30 3e 29 20 7b 0d 0a  t_(Case_<0>) {..
07f0: 20 20 69 6e 74 20 4d 20 3d 20 31 3b 20 0d 0a 20    int M = 1; .. 
0800: 20 69 6e 74 20 6c 6f 77 65 72 42 6f 75 6e 64 20   int lowerBound 
0810: 3d 20 31 3b 20 0d 0a 20 20 69 6e 74 20 75 70 70  = 1; ..  int upp
0820: 65 72 42 6f 75 6e 64 20 3d 20 31 30 3b 20 0d 0a  erBound = 10; ..
0830: 20 20 69 6e 74 20 4e 20 3d 20 33 3b 20 0d 0a 20    int N = 3; .. 
0840: 20 69 6e 74 20 4b 20 3d 20 31 3b 20 0d 0a 20 20   int K = 1; ..  
0850: 64 6f 75 62 6c 65 20 52 65 74 56 61 6c 20 3d 20  double RetVal = 
0860: 30 2e 31 3b 20 0d 0a 20 20 72 65 74 75 72 6e 20  0.1; ..  return 
0870: 76 65 72 69 66 79 5f 63 61 73 65 28 52 65 74 56  verify_case(RetV
0880: 61 6c 2c 20 4b 74 68 50 72 6f 62 61 62 6c 65 45  al, KthProbableE
0890: 6c 65 6d 65 6e 74 28 29 2e 70 72 6f 62 61 62 69  lement().probabi
08a0: 6c 69 74 79 28 4d 2c 20 6c 6f 77 65 72 42 6f 75  lity(M, lowerBou
08b0: 6e 64 2c 20 75 70 70 65 72 42 6f 75 6e 64 2c 20  nd, upperBound, 
08c0: 4e 2c 20 4b 29 29 3b 20 7d 0d 0a 69 6e 74 20 54  N, K)); }..int T
08d0: 65 73 74 5f 28 43 61 73 65 5f 3c 31 3e 29 20 7b  est_(Case_<1>) {
08e0: 0d 0a 20 20 69 6e 74 20 4d 20 3d 20 33 3b 20 0d  ..  int M = 3; .
08f0: 0a 20 20 69 6e 74 20 6c 6f 77 65 72 42 6f 75 6e  .  int lowerBoun
0900: 64 20 3d 20 31 3b 20 0d 0a 20 20 69 6e 74 20 75  d = 1; ..  int u
0910: 70 70 65 72 42 6f 75 6e 64 20 3d 20 32 3b 20 0d  pperBound = 2; .
0920: 0a 20 20 69 6e 74 20 4e 20 3d 20 32 3b 20 0d 0a  .  int N = 2; ..
0930: 20 20 69 6e 74 20 4b 20 3d 20 32 3b 20 0d 0a 20    int K = 2; .. 
0940: 20 64 6f 75 62 6c 65 20 52 65 74 56 61 6c 20 3d   double RetVal =
0950: 20 30 2e 35 3b 20 0d 0a 20 20 72 65 74 75 72 6e   0.5; ..  return
0960: 20 76 65 72 69 66 79 5f 63 61 73 65 28 52 65 74   verify_case(Ret
0970: 56 61 6c 2c 20 4b 74 68 50 72 6f 62 61 62 6c 65  Val, KthProbable
0980: 45 6c 65 6d 65 6e 74 28 29 2e 70 72 6f 62 61 62  Element().probab
0990: 69 6c 69 74 79 28 4d 2c 20 6c 6f 77 65 72 42 6f  ility(M, lowerBo
09a0: 75 6e 64 2c 20 75 70 70 65 72 42 6f 75 6e 64 2c  und, upperBound,
09b0: 20 4e 2c 20 4b 29 29 3b 20 7d 0d 0a 69 6e 74 20   N, K)); }..int 
09c0: 54 65 73 74 5f 28 43 61 73 65 5f 3c 32 3e 29 20  Test_(Case_<2>) 
09d0: 7b 0d 0a 20 20 69 6e 74 20 4d 20 3d 20 33 3b 20  {..  int M = 3; 
09e0: 0d 0a 20 20 69 6e 74 20 6c 6f 77 65 72 42 6f 75  ..  int lowerBou
09f0: 6e 64 20 3d 20 31 3b 20 0d 0a 20 20 69 6e 74 20  nd = 1; ..  int 
0a00: 75 70 70 65 72 42 6f 75 6e 64 20 3d 20 33 3b 20  upperBound = 3; 
0a10: 0d 0a 20 20 69 6e 74 20 4e 20 3d 20 32 3b 20 0d  ..  int N = 2; .
0a20: 0a 20 20 69 6e 74 20 4b 20 3d 20 32 3b 20 0d 0a  .  int K = 2; ..
0a30: 20 20 64 6f 75 62 6c 65 20 52 65 74 56 61 6c 20    double RetVal 
0a40: 3d 20 30 2e 34 38 31 34 38 31 34 38 31 34 38 31  = 0.481481481481
0a50: 34 38 31 34 35 3b 20 0d 0a 20 20 72 65 74 75 72  48145; ..  retur
0a60: 6e 20 76 65 72 69 66 79 5f 63 61 73 65 28 52 65  n verify_case(Re
0a70: 74 56 61 6c 2c 20 4b 74 68 50 72 6f 62 61 62 6c  tVal, KthProbabl
0a80: 65 45 6c 65 6d 65 6e 74 28 29 2e 70 72 6f 62 61  eElement().proba
0a90: 62 69 6c 69 74 79 28 4d 2c 20 6c 6f 77 65 72 42  bility(M, lowerB
0aa0: 6f 75 6e 64 2c 20 75 70 70 65 72 42 6f 75 6e 64  ound, upperBound
0ab0: 2c 20 4e 2c 20 4b 29 29 3b 20 7d 0d 0a 69 6e 74  , N, K)); }..int
0ac0: 20 54 65 73 74 5f 28 43 61 73 65 5f 3c 33 3e 29   Test_(Case_<3>)
0ad0: 20 7b 0d 0a 20 20 69 6e 74 20 4d 20 3d 20 31 30   {..  int M = 10
0ae0: 3b 20 0d 0a 20 20 69 6e 74 20 6c 6f 77 65 72 42  ; ..  int lowerB
0af0: 6f 75 6e 64 20 3d 20 31 3b 20 0d 0a 20 20 69 6e  ound = 1; ..  in
0b00: 74 20 75 70 70 65 72 42 6f 75 6e 64 20 3d 20 31  t upperBound = 1
0b10: 30 3b 20 0d 0a 20 20 69 6e 74 20 4e 20 3d 20 31  0; ..  int N = 1
0b20: 3b 20 0d 0a 20 20 69 6e 74 20 4b 20 3d 20 31 30  ; ..  int K = 10
0b30: 3b 20 0d 0a 20 20 64 6f 75 62 6c 65 20 52 65 74  ; ..  double Ret
0b40: 56 61 6c 20 3d 20 31 2e 30 30 30 30 30 30 30 30  Val = 1.00000000
0b50: 30 30 30 30 30 30 30 33 45 2d 31 30 3b 20 0d 0a  00000003E-10; ..
0b60: 20 20 72 65 74 75 72 6e 20 76 65 72 69 66 79 5f    return verify_
0b70: 63 61 73 65 28 52 65 74 56 61 6c 2c 20 4b 74 68  case(RetVal, Kth
0b80: 50 72 6f 62 61 62 6c 65 45 6c 65 6d 65 6e 74 28  ProbableElement(
0b90: 29 2e 70 72 6f 62 61 62 69 6c 69 74 79 28 4d 2c  ).probability(M,
0ba0: 20 6c 6f 77 65 72 42 6f 75 6e 64 2c 20 75 70 70   lowerBound, upp
0bb0: 65 72 42 6f 75 6e 64 2c 20 4e 2c 20 4b 29 29 3b  erBound, N, K));
0bc0: 20 7d 0d 0a 69 6e 74 20 54 65 73 74 5f 28 43 61   }..int Test_(Ca
0bd0: 73 65 5f 3c 34 3e 29 20 7b 0d 0a 20 20 69 6e 74  se_<4>) {..  int
0be0: 20 4d 20 3d 20 34 3b 20 0d 0a 20 20 69 6e 74 20   M = 4; ..  int 
0bf0: 6c 6f 77 65 72 42 6f 75 6e 64 20 3d 20 36 31 3b  lowerBound = 61;
0c00: 20 0d 0a 20 20 69 6e 74 20 75 70 70 65 72 42 6f   ..  int upperBo
0c10: 75 6e 64 20 3d 20 36 35 3b 20 0d 0a 20 20 69 6e  und = 65; ..  in
0c20: 74 20 4e 20 3d 20 36 32 3b 20 0d 0a 20 20 69 6e  t N = 62; ..  in
0c30: 74 20 4b 20 3d 20 33 3b 20 0d 0a 20 20 64 6f 75  t K = 3; ..  dou
0c40: 62 6c 65 20 52 65 74 56 61 6c 20 3d 20 30 2e 31  ble RetVal = 0.1
0c50: 35 32 30 30 30 30 30 30 30 30 30 30 30 30 30 32  5200000000000002
0c60: 3b 20 0d 0a 20 20 72 65 74 75 72 6e 20 76 65 72  ; ..  return ver
0c70: 69 66 79 5f 63 61 73 65 28 52 65 74 56 61 6c 2c  ify_case(RetVal,
0c80: 20 4b 74 68 50 72 6f 62 61 62 6c 65 45 6c 65 6d   KthProbableElem
0c90: 65 6e 74 28 29 2e 70 72 6f 62 61 62 69 6c 69 74  ent().probabilit
0ca0: 79 28 4d 2c 20 6c 6f 77 65 72 42 6f 75 6e 64 2c  y(M, lowerBound,
0cb0: 20 75 70 70 65 72 42 6f 75 6e 64 2c 20 4e 2c 20   upperBound, N, 
0cc0: 4b 29 29 3b 20 7d 0d 0a 0d 0a 74 65 6d 70 6c 61  K)); }....templa
0cd0: 74 65 3c 69 6e 74 20 4e 3e 20 76 6f 69 64 20 52  te<int N> void R
0ce0: 75 6e 5f 28 29 20 7b 20 63 65 72 72 20 3c 3c 20  un_() { cerr << 
0cf0: 22 54 65 73 74 20 43 61 73 65 20 23 22 20 3c 3c  "Test Case #" <<
0d00: 20 4e 20 3c 3c 20 22 2e 2e 2e 22 20 3c 3c 20 66   N << "..." << f
0d10: 6c 75 73 68 3b 20 54 65 73 74 5f 28 43 61 73 65  lush; Test_(Case
0d20: 5f 3c 4e 3e 28 29 29 3b 20 52 75 6e 5f 3c 73 69  _<N>()); Run_<si
0d30: 7a 65 6f 66 28 54 65 73 74 5f 28 43 61 73 65 5f  zeof(Test_(Case_
0d40: 3c 4e 2b 31 3e 28 29 29 29 3d 3d 31 20 3f 20 2d  <N+1>()))==1 ? -
0d50: 31 20 3a 20 4e 2b 31 3e 28 29 3b 20 7d 0d 0a 74  1 : N+1>(); }..t
0d60: 65 6d 70 6c 61 74 65 3c 3e 20 20 20 20 20 20 76  emplate<>      v
0d70: 6f 69 64 20 52 75 6e 5f 3c 2d 31 3e 28 29 20 7b  oid Run_<-1>() {
0d80: 7d 0d 0a 69 6e 74 20 6d 61 69 6e 28 29 20 7b 20  }..int main() { 
0d90: 52 75 6e 5f 3c 30 3e 28 29 3b 20 7d 0d 0a 2f 2f  Run_<0>(); }..//
0da0: 20 45 4e 44 20 43 55 54 20 48 45 52 45 0d 0a 0d   END CUT HERE...
0db0: 0a                                               .