Hex Artifact Content
Not logged in

Artifact 48b50023e1108913a83b4121a32f69c6092b5bcc:


0000: 23 69 6e 63 6c 75 64 65 20 3c 69 6f 73 74 72 65  #include <iostre
0010: 61 6d 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 73  am>..#include <s
0020: 73 74 72 65 61 6d 3e 0d 0a 23 69 6e 63 6c 75 64  stream>..#includ
0030: 65 20 3c 69 6f 6d 61 6e 69 70 3e 0d 0a 23 69 6e  e <iomanip>..#in
0040: 63 6c 75 64 65 20 3c 76 65 63 74 6f 72 3e 0d 0a  clude <vector>..
0050: 23 69 6e 63 6c 75 64 65 20 3c 73 74 72 69 6e 67  #include <string
0060: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 6d 61 70  >..#include <map
0070: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 73 65 74  >..#include <set
0080: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 61 6c 67  >..#include <alg
0090: 6f 72 69 74 68 6d 3e 0d 0a 23 69 6e 63 6c 75 64  orithm>..#includ
00a0: 65 20 3c 6e 75 6d 65 72 69 63 3e 0d 0a 23 69 6e  e <numeric>..#in
00b0: 63 6c 75 64 65 20 3c 69 74 65 72 61 74 6f 72 3e  clude <iterator>
00c0: 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 66 75 6e 63  ..#include <func
00d0: 74 69 6f 6e 61 6c 3e 0d 0a 23 69 6e 63 6c 75 64  tional>..#includ
00e0: 65 20 3c 63 6f 6d 70 6c 65 78 3e 0d 0a 23 69 6e  e <complex>..#in
00f0: 63 6c 75 64 65 20 3c 71 75 65 75 65 3e 0d 0a 23  clude <queue>..#
0100: 69 6e 63 6c 75 64 65 20 3c 73 74 61 63 6b 3e 0d  include <stack>.
0110: 0a 23 69 6e 63 6c 75 64 65 20 3c 63 6d 61 74 68  .#include <cmath
0120: 3e 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 63 61 73  >..#include <cas
0130: 73 65 72 74 3e 0d 0a 23 69 6e 63 6c 75 64 65 20  sert>..#include 
0140: 3c 74 75 70 6c 65 3e 0d 0a 75 73 69 6e 67 20 6e  <tuple>..using n
0150: 61 6d 65 73 70 61 63 65 20 73 74 64 3b 0d 0a 74  amespace std;..t
0160: 79 70 65 64 65 66 20 6c 6f 6e 67 20 6c 6f 6e 67  ypedef long long
0170: 20 4c 4c 3b 0d 0a 74 79 70 65 64 65 66 20 63 6f   LL;..typedef co
0180: 6d 70 6c 65 78 3c 64 6f 75 62 6c 65 3e 20 43 4d  mplex<double> CM
0190: 50 3b 0d 0a 0d 0a 63 6c 61 73 73 20 54 65 73 74  P;....class Test
01a0: 50 72 6f 63 74 6f 72 69 6e 67 20 7b 20 70 75 62  Proctoring { pub
01b0: 6c 69 63 3a 0d 0a 09 64 6f 75 62 6c 65 20 65 78  lic:...double ex
01c0: 70 65 63 74 65 64 54 69 6d 65 28 76 65 63 74 6f  pectedTime(vecto
01d0: 72 20 3c 69 6e 74 3e 20 70 2c 20 76 65 63 74 6f  r <int> p, vecto
01e0: 72 20 3c 69 6e 74 3e 20 71 29 0d 0a 09 7b 0d 0a  r <int> q)...{..
01f0: 09 09 63 6f 6e 73 74 20 69 6e 74 20 4e 20 3d 20  ..const int N = 
0200: 70 2e 73 69 7a 65 28 29 3b 0d 0a 09 09 76 65 63  p.size();....vec
0210: 74 6f 72 3c 64 6f 75 62 6c 65 3e 20 64 70 28 31  tor<double> dp(1
0220: 3c 3c 4e 29 2c 20 64 70 32 28 28 4e 2b 31 29 2a  <<N), dp2((N+1)*
0230: 28 31 3c 3c 4e 29 29 3b 0d 0a 09 09 66 6f 72 20  (1<<N));....for 
0240: 28 69 6e 74 20 6d 61 73 6b 20 3d 20 31 3b 20 6d  (int mask = 1; m
0250: 61 73 6b 20 3c 20 28 31 3c 3c 4e 29 3b 20 2b 2b  ask < (1<<N); ++
0260: 6d 61 73 6b 29 20 7b 0d 0a 09 09 09 66 6f 72 20  mask) {.....for 
0270: 28 69 6e 74 20 69 20 3d 20 4e 3b 20 69 20 3e 3d  (int i = N; i >=
0280: 20 30 3b 20 2d 2d 69 29 0d 0a 09 09 09 09 69 66   0; --i)......if
0290: 20 28 69 20 3d 3d 20 4e 29 0d 0a 09 09 09 09 09   (i == N).......
02a0: 64 70 32 5b 28 69 20 3c 3c 20 4e 29 20 7c 20 6d  dp2[(i << N) | m
02b0: 61 73 6b 5d 20 3d 20 64 70 5b 6d 61 73 6b 5d 3b  ask] = dp[mask];
02c0: 0d 0a 09 09 09 09 65 6c 73 65 20 69 66 20 28 28  ......else if ((
02d0: 6d 61 73 6b 26 28 31 20 3c 3c 20 69 29 29 20 3d  mask&(1 << i)) =
02e0: 3d 20 30 29 0d 0a 09 09 09 09 09 64 70 32 5b 28  = 0).......dp2[(
02f0: 69 20 3c 3c 20 4e 29 20 7c 20 6d 61 73 6b 5d 20  i << N) | mask] 
0300: 3d 20 64 70 32 5b 28 28 69 20 2b 20 31 29 20 3c  = dp2[((i + 1) <
0310: 3c 20 4e 29 20 7c 20 6d 61 73 6b 5d 3b 0d 0a 09  < N) | mask];...
0320: 09 09 09 65 6c 73 65 20 7b 0d 0a 09 09 09 09 09  ...else {.......
0330: 64 6f 75 62 6c 65 20 65 31 20 3d 20 64 70 32 5b  double e1 = dp2[
0340: 28 28 69 20 2b 20 31 29 20 3c 3c 20 4e 29 20 7c  ((i + 1) << N) |
0350: 20 6d 61 73 6b 20 26 7e 28 31 20 3c 3c 20 69 29   mask &~(1 << i)
0360: 5d 3b 0d 0a 09 09 09 09 09 64 6f 75 62 6c 65 20  ];.......double 
0370: 65 32 20 3d 20 64 70 32 5b 28 28 69 20 2b 20 31  e2 = dp2[((i + 1
0380: 29 20 3c 3c 20 4e 29 20 7c 20 6d 61 73 6b 5d 3b  ) << N) | mask];
0390: 0d 0a 09 09 09 09 09 64 70 32 5b 28 69 20 3c 3c  .......dp2[(i <<
03a0: 20 4e 29 20 7c 20 6d 61 73 6b 5d 20 3d 20 28 70   N) | mask] = (p
03b0: 5b 69 5d 20 2a 20 65 31 20 2b 20 28 71 5b 69 5d  [i] * e1 + (q[i]
03c0: 20 2d 20 70 5b 69 5d 29 2a 65 32 29 20 2f 20 71   - p[i])*e2) / q
03d0: 5b 69 5d 3b 0d 0a 09 09 09 09 7d 0d 0a 0d 0a 09  [i];......}.....
03e0: 09 09 64 6f 75 62 6c 65 20 70 70 20 3d 20 31 2e  ..double pp = 1.
03f0: 30 3b 0d 0a 09 09 09 66 6f 72 20 28 69 6e 74 20  0;.....for (int 
0400: 69 20 3d 20 30 3b 20 69 3c 4e 3b 20 2b 2b 69 29  i = 0; i<N; ++i)
0410: 20 69 66 20 28 6d 61 73 6b 26 28 31 20 3c 3c 20   if (mask&(1 << 
0420: 69 29 29 0d 0a 09 09 09 09 70 70 20 3d 20 70 70  i))......pp = pp
0430: 20 2a 20 28 71 5b 69 5d 20 2d 20 70 5b 69 5d 29   * (q[i] - p[i])
0440: 20 2f 20 71 5b 69 5d 3b 0d 0a 09 09 09 64 70 5b   / q[i];.....dp[
0450: 6d 61 73 6b 5d 20 3d 20 28 64 70 32 5b 6d 61 73  mask] = (dp2[mas
0460: 6b 5d 2b 31 29 20 2f 20 28 31 20 2d 20 70 70 29  k]+1) / (1 - pp)
0470: 3b 0d 0a 0d 0a 09 09 09 2f 2f 20 52 65 72 75 6e  ;.......// Rerun
0480: 0d 0a 09 09 09 66 6f 72 20 28 69 6e 74 20 69 20  .....for (int i 
0490: 3d 20 4e 3b 20 69 20 3e 3d 20 30 3b 20 2d 2d 69  = N; i >= 0; --i
04a0: 29 0d 0a 09 09 09 09 69 66 20 28 69 20 3d 3d 20  )......if (i == 
04b0: 4e 29 0d 0a 09 09 09 09 09 64 70 32 5b 28 69 20  N).......dp2[(i 
04c0: 3c 3c 20 4e 29 20 7c 20 6d 61 73 6b 5d 20 3d 20  << N) | mask] = 
04d0: 64 70 5b 6d 61 73 6b 5d 3b 0d 0a 09 09 09 09 65  dp[mask];......e
04e0: 6c 73 65 20 69 66 20 28 28 6d 61 73 6b 26 28 31  lse if ((mask&(1
04f0: 20 3c 3c 20 69 29 29 20 3d 3d 20 30 29 0d 0a 09   << i)) == 0)...
0500: 09 09 09 09 64 70 32 5b 28 69 20 3c 3c 20 4e 29  ....dp2[(i << N)
0510: 20 7c 20 6d 61 73 6b 5d 20 3d 20 64 70 32 5b 28   | mask] = dp2[(
0520: 28 69 20 2b 20 31 29 20 3c 3c 20 4e 29 20 7c 20  (i + 1) << N) | 
0530: 6d 61 73 6b 5d 3b 0d 0a 09 09 09 09 65 6c 73 65  mask];......else
0540: 20 7b 0d 0a 09 09 09 09 09 64 6f 75 62 6c 65 20   {.......double 
0550: 65 31 20 3d 20 64 70 32 5b 28 28 69 20 2b 20 31  e1 = dp2[((i + 1
0560: 29 20 3c 3c 20 4e 29 20 7c 20 6d 61 73 6b 20 26  ) << N) | mask &
0570: 7e 28 31 20 3c 3c 20 69 29 5d 3b 0d 0a 09 09 09  ~(1 << i)];.....
0580: 09 09 64 6f 75 62 6c 65 20 65 32 20 3d 20 64 70  ..double e2 = dp
0590: 32 5b 28 28 69 20 2b 20 31 29 20 3c 3c 20 4e 29  2[((i + 1) << N)
05a0: 20 7c 20 6d 61 73 6b 5d 3b 0d 0a 09 09 09 09 09   | mask];.......
05b0: 64 70 32 5b 28 69 20 3c 3c 20 4e 29 20 7c 20 6d  dp2[(i << N) | m
05c0: 61 73 6b 5d 20 3d 20 28 70 5b 69 5d 20 2a 20 65  ask] = (p[i] * e
05d0: 31 20 2b 20 28 71 5b 69 5d 20 2d 20 70 5b 69 5d  1 + (q[i] - p[i]
05e0: 29 2a 65 32 29 20 2f 20 71 5b 69 5d 3b 0d 0a 09  )*e2) / q[i];...
05f0: 09 09 09 7d 0d 0a 09 09 7d 0d 0a 0d 0a 09 09 72  ...}....}......r
0600: 65 74 75 72 6e 20 64 70 5b 28 31 3c 3c 4e 29 20  eturn dp[(1<<N) 
0610: 2d 20 31 5d 3b 0d 0a 09 7d 0d 0a 7d 3b 0d 0a 0d  - 1];...}..};...
0620: 0a 2f 2f 20 42 45 47 49 4e 20 43 55 54 20 48 45  .// BEGIN CUT HE
0630: 52 45 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 63 74  RE..#include <ct
0640: 69 6d 65 3e 0d 0a 64 6f 75 62 6c 65 20 73 74 61  ime>..double sta
0650: 72 74 5f 74 69 6d 65 3b 20 73 74 72 69 6e 67 20  rt_time; string 
0660: 74 69 6d 65 72 28 29 0d 0a 20 7b 20 6f 73 74 72  timer().. { ostr
0670: 69 6e 67 73 74 72 65 61 6d 20 6f 73 3b 20 6f 73  ingstream os; os
0680: 20 3c 3c 20 22 20 28 22 20 3c 3c 20 69 6e 74 28   << " (" << int(
0690: 28 63 6c 6f 63 6b 28 29 2d 73 74 61 72 74 5f 74  (clock()-start_t
06a0: 69 6d 65 29 2f 43 4c 4f 43 4b 53 5f 50 45 52 5f  ime)/CLOCKS_PER_
06b0: 53 45 43 2a 31 30 30 30 29 20 3c 3c 20 22 20 6d  SEC*1000) << " m
06c0: 73 65 63 29 22 3b 20 72 65 74 75 72 6e 20 6f 73  sec)"; return os
06d0: 2e 73 74 72 28 29 3b 20 7d 0d 0a 74 65 6d 70 6c  .str(); }..templ
06e0: 61 74 65 3c 74 79 70 65 6e 61 6d 65 20 54 3e 20  ate<typename T> 
06f0: 6f 73 74 72 65 61 6d 26 20 6f 70 65 72 61 74 6f  ostream& operato
0700: 72 3c 3c 28 6f 73 74 72 65 61 6d 26 20 6f 73 2c  r<<(ostream& os,
0710: 20 63 6f 6e 73 74 20 76 65 63 74 6f 72 3c 54 3e   const vector<T>
0720: 26 20 76 29 0d 0a 20 7b 20 6f 73 20 3c 3c 20 22  & v).. { os << "
0730: 7b 20 22 3b 0d 0a 20 20 20 66 6f 72 28 74 79 70  { ";..   for(typ
0740: 65 6e 61 6d 65 20 76 65 63 74 6f 72 3c 54 3e 3a  ename vector<T>:
0750: 3a 63 6f 6e 73 74 5f 69 74 65 72 61 74 6f 72 20  :const_iterator 
0760: 69 74 3d 76 2e 62 65 67 69 6e 28 29 3b 20 69 74  it=v.begin(); it
0770: 21 3d 76 2e 65 6e 64 28 29 3b 20 2b 2b 69 74 29  !=v.end(); ++it)
0780: 0d 0a 20 20 20 6f 73 20 3c 3c 20 27 5c 22 27 20  ..   os << '\"' 
0790: 3c 3c 20 2a 69 74 20 3c 3c 20 27 5c 22 27 20 3c  << *it << '\"' <
07a0: 3c 20 28 69 74 2b 31 3d 3d 76 2e 65 6e 64 28 29  < (it+1==v.end()
07b0: 20 3f 20 22 22 20 3a 20 22 2c 20 22 29 3b 20 6f   ? "" : ", "); o
07c0: 73 20 3c 3c 20 22 20 7d 22 3b 20 72 65 74 75 72  s << " }"; retur
07d0: 6e 20 6f 73 3b 20 7d 0d 0a 76 6f 69 64 20 76 65  n os; }..void ve
07e0: 72 69 66 79 5f 63 61 73 65 28 63 6f 6e 73 74 20  rify_case(const 
07f0: 64 6f 75 62 6c 65 26 20 45 78 70 65 63 74 65 64  double& Expected
0800: 2c 20 63 6f 6e 73 74 20 64 6f 75 62 6c 65 26 20  , const double& 
0810: 52 65 63 65 69 76 65 64 29 20 7b 0d 0a 20 62 6f  Received) {.. bo
0820: 6f 6c 20 6f 6b 20 3d 20 28 61 62 73 28 45 78 70  ol ok = (abs(Exp
0830: 65 63 74 65 64 20 2d 20 52 65 63 65 69 76 65 64  ected - Received
0840: 29 2f 61 62 73 28 45 78 70 65 63 74 65 64 29 20  )/abs(Expected) 
0850: 3c 20 31 65 2d 36 29 3b 0d 0a 20 69 66 28 6f 6b  < 1e-6);.. if(ok
0860: 29 20 63 65 72 72 20 3c 3c 20 22 50 41 53 53 45  ) cerr << "PASSE
0870: 44 22 20 3c 3c 20 74 69 6d 65 72 28 29 20 3c 3c  D" << timer() <<
0880: 20 65 6e 64 6c 3b 20 20 65 6c 73 65 20 7b 20 63   endl;  else { c
0890: 65 72 72 20 3c 3c 20 22 46 41 49 4c 45 44 22 20  err << "FAILED" 
08a0: 3c 3c 20 74 69 6d 65 72 28 29 20 3c 3c 20 65 6e  << timer() << en
08b0: 64 6c 3b 0d 0a 20 63 65 72 72 20 3c 3c 20 22 5c  dl;.. cerr << "\
08c0: 74 6f 3a 20 5c 22 22 20 3c 3c 20 45 78 70 65 63  to: \"" << Expec
08d0: 74 65 64 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 65  ted << '\"' << e
08e0: 6e 64 6c 20 3c 3c 20 22 5c 74 78 3a 20 5c 22 22  ndl << "\tx: \""
08f0: 20 3c 3c 20 52 65 63 65 69 76 65 64 20 3c 3c 20   << Received << 
0900: 27 5c 22 27 20 3c 3c 20 65 6e 64 6c 3b 20 7d 20  '\"' << endl; } 
0910: 7d 0d 0a 23 64 65 66 69 6e 65 20 43 41 53 45 28  }..#define CASE(
0920: 4e 29 20 7b 63 65 72 72 20 3c 3c 20 22 54 65 73  N) {cerr << "Tes
0930: 74 20 43 61 73 65 20 23 22 20 3c 3c 20 4e 20 3c  t Case #" << N <
0940: 3c 20 22 2e 2e 2e 22 20 3c 3c 20 66 6c 75 73 68  < "..." << flush
0950: 3b 20 73 74 61 72 74 5f 74 69 6d 65 3d 63 6c 6f  ; start_time=clo
0960: 63 6b 28 29 3b 0d 0a 23 64 65 66 69 6e 65 20 45  ck();..#define E
0970: 4e 44 09 20 76 65 72 69 66 79 5f 63 61 73 65 28  ND. verify_case(
0980: 5f 2c 20 54 65 73 74 50 72 6f 63 74 6f 72 69 6e  _, TestProctorin
0990: 67 28 29 2e 65 78 70 65 63 74 65 64 54 69 6d 65  g().expectedTime
09a0: 28 70 2c 20 71 29 29 3b 7d 0d 0a 69 6e 74 20 6d  (p, q));}..int m
09b0: 61 69 6e 28 29 7b 0d 0a 43 41 53 45 28 30 29 0d  ain(){..CASE(0).
09c0: 0a 09 69 6e 74 20 70 5f 5b 5d 20 3d 20 7b 32 7d  ..int p_[] = {2}
09d0: 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e  ;...  vector <in
09e0: 74 3e 20 70 28 70 5f 2c 20 70 5f 2b 73 69 7a 65  t> p(p_, p_+size
09f0: 6f 66 28 70 5f 29 2f 73 69 7a 65 6f 66 28 2a 70  of(p_)/sizeof(*p
0a00: 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 71 5f 5b 5d  _)); ...int q_[]
0a10: 20 3d 20 7b 34 7d 3b 0d 0a 09 20 20 76 65 63 74   = {4};...  vect
0a20: 6f 72 20 3c 69 6e 74 3e 20 71 28 71 5f 2c 20 71  or <int> q(q_, q
0a30: 5f 2b 73 69 7a 65 6f 66 28 71 5f 29 2f 73 69 7a  _+sizeof(q_)/siz
0a40: 65 6f 66 28 2a 71 5f 29 29 3b 20 0d 0a 09 64 6f  eof(*q_)); ...do
0a50: 75 62 6c 65 20 5f 20 3d 20 32 2e 30 3b 20 0d 0a  uble _ = 2.0; ..
0a60: 45 4e 44 0d 0a 43 41 53 45 28 31 29 0d 0a 09 69  END..CASE(1)...i
0a70: 6e 74 20 70 5f 5b 5d 20 3d 20 7b 31 2c 32 7d 3b  nt p_[] = {1,2};
0a80: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74  ...  vector <int
0a90: 3e 20 70 28 70 5f 2c 20 70 5f 2b 73 69 7a 65 6f  > p(p_, p_+sizeo
0aa0: 66 28 70 5f 29 2f 73 69 7a 65 6f 66 28 2a 70 5f  f(p_)/sizeof(*p_
0ab0: 29 29 3b 20 0d 0a 09 69 6e 74 20 71 5f 5b 5d 20  )); ...int q_[] 
0ac0: 3d 20 7b 32 2c 34 7d 3b 0d 0a 09 20 20 76 65 63  = {2,4};...  vec
0ad0: 74 6f 72 20 3c 69 6e 74 3e 20 71 28 71 5f 2c 20  tor <int> q(q_, 
0ae0: 71 5f 2b 73 69 7a 65 6f 66 28 71 5f 29 2f 73 69  q_+sizeof(q_)/si
0af0: 7a 65 6f 66 28 2a 71 5f 29 29 3b 20 0d 0a 09 64  zeof(*q_)); ...d
0b00: 6f 75 62 6c 65 20 5f 20 3d 20 32 2e 36 36 36 36  ouble _ = 2.6666
0b10: 36 36 36 36 36 36 36 36 36 36 36 35 3b 20 0d 0a  666666666665; ..
0b20: 45 4e 44 0d 0a 43 41 53 45 28 32 29 0d 0a 09 69  END..CASE(2)...i
0b30: 6e 74 20 70 5f 5b 5d 20 3d 20 7b 33 2c 31 2c 32  nt p_[] = {3,1,2
0b40: 2c 34 2c 32 2c 35 7d 3b 0d 0a 09 20 20 76 65 63  ,4,2,5};...  vec
0b50: 74 6f 72 20 3c 69 6e 74 3e 20 70 28 70 5f 2c 20  tor <int> p(p_, 
0b60: 70 5f 2b 73 69 7a 65 6f 66 28 70 5f 29 2f 73 69  p_+sizeof(p_)/si
0b70: 7a 65 6f 66 28 2a 70 5f 29 29 3b 20 0d 0a 09 69  zeof(*p_)); ...i
0b80: 6e 74 20 71 5f 5b 5d 20 3d 20 7b 33 2c 31 2c 32  nt q_[] = {3,1,2
0b90: 2c 34 2c 32 2c 35 7d 3b 0d 0a 09 20 20 76 65 63  ,4,2,5};...  vec
0ba0: 74 6f 72 20 3c 69 6e 74 3e 20 71 28 71 5f 2c 20  tor <int> q(q_, 
0bb0: 71 5f 2b 73 69 7a 65 6f 66 28 71 5f 29 2f 73 69  q_+sizeof(q_)/si
0bc0: 7a 65 6f 66 28 2a 71 5f 29 29 3b 20 0d 0a 09 64  zeof(*q_)); ...d
0bd0: 6f 75 62 6c 65 20 5f 20 3d 20 31 2e 30 3b 20 0d  ouble _ = 1.0; .
0be0: 0a 45 4e 44 0d 0a 43 41 53 45 28 33 29 0d 0a 09  .END..CASE(3)...
0bf0: 69 6e 74 20 70 5f 5b 5d 20 3d 20 7b 31 2c 31 2c  int p_[] = {1,1,
0c00: 31 2c 31 2c 31 2c 31 2c 31 2c 31 7d 3b 0d 0a 09  1,1,1,1,1,1};...
0c10: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 70    vector <int> p
0c20: 28 70 5f 2c 20 70 5f 2b 73 69 7a 65 6f 66 28 70  (p_, p_+sizeof(p
0c30: 5f 29 2f 73 69 7a 65 6f 66 28 2a 70 5f 29 29 3b  _)/sizeof(*p_));
0c40: 20 0d 0a 09 69 6e 74 20 71 5f 5b 5d 20 3d 20 7b   ...int q_[] = {
0c50: 31 2c 32 2c 33 2c 34 2c 35 2c 36 2c 37 2c 38 7d  1,2,3,4,5,6,7,8}
0c60: 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e  ;...  vector <in
0c70: 74 3e 20 71 28 71 5f 2c 20 71 5f 2b 73 69 7a 65  t> q(q_, q_+size
0c80: 6f 66 28 71 5f 29 2f 73 69 7a 65 6f 66 28 2a 71  of(q_)/sizeof(*q
0c90: 5f 29 29 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 5f  _)); ...double _
0ca0: 20 3d 20 31 33 2e 36 30 34 38 33 34 36 36 35 31   = 13.6048346651
0cb0: 32 30 39 39 31 3b 20 0d 0a 45 4e 44 0d 0a 43 41  20991; ..END..CA
0cc0: 53 45 28 34 29 0d 0a 09 69 6e 74 20 70 5f 5b 5d  SE(4)...int p_[]
0cd0: 20 3d 20 7b 32 2c 33 2c 35 2c 37 2c 31 31 2c 31   = {2,3,5,7,11,1
0ce0: 33 2c 31 37 7d 3b 0d 0a 09 20 20 76 65 63 74 6f  3,17};...  vecto
0cf0: 72 20 3c 69 6e 74 3e 20 70 28 70 5f 2c 20 70 5f  r <int> p(p_, p_
0d00: 2b 73 69 7a 65 6f 66 28 70 5f 29 2f 73 69 7a 65  +sizeof(p_)/size
0d10: 6f 66 28 2a 70 5f 29 29 3b 20 0d 0a 09 69 6e 74  of(*p_)); ...int
0d20: 20 71 5f 5b 5d 20 3d 20 7b 33 2c 35 2c 37 2c 31   q_[] = {3,5,7,1
0d30: 31 2c 31 33 2c 31 37 2c 31 39 7d 3b 0d 0a 09 20  1,13,17,19};... 
0d40: 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 71 28   vector <int> q(
0d50: 71 5f 2c 20 71 5f 2b 73 69 7a 65 6f 66 28 71 5f  q_, q_+sizeof(q_
0d60: 29 2f 73 69 7a 65 6f 66 28 2a 71 5f 29 29 3b 20  )/sizeof(*q_)); 
0d70: 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20 32 2e  ...double _ = 2.
0d80: 36 32 39 39 34 34 35 30 36 35 39 32 34 32 37 36  6299445065924276
0d90: 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 35 29  ; ..END..CASE(5)
0da0: 0d 0a 09 69 6e 74 20 70 5f 5b 5d 20 3d 20 7b 31  ...int p_[] = {1
0db0: 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31  ,1,1,1,1,1,1,1,1
0dc0: 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31 2c 31  ,1,1,1,1,1,1,1,1
0dd0: 2c 31 2c 31 2c 31 7d 3b 0d 0a 09 20 20 76 65 63  ,1,1,1};...  vec
0de0: 74 6f 72 20 3c 69 6e 74 3e 20 70 28 70 5f 2c 20  tor <int> p(p_, 
0df0: 70 5f 2b 73 69 7a 65 6f 66 28 70 5f 29 2f 73 69  p_+sizeof(p_)/si
0e00: 7a 65 6f 66 28 2a 70 5f 29 29 3b 20 0d 0a 09 69  zeof(*p_)); ...i
0e10: 6e 74 20 71 5f 5b 5d 20 3d 20 7b 31 30 30 30 30  nt q_[] = {10000
0e20: 30 30 30 30 30 2c 31 30 30 30 30 30 30 30 30 30  00000,1000000000
0e30: 2c 31 30 30 30 30 30 30 30 30 30 2c 31 30 30 30  ,1000000000,1000
0e40: 30 30 30 30 30 30 2c 31 30 30 30 30 30 30 30 30  000000,100000000
0e50: 30 2c 0d 0a 31 30 30 30 30 30 30 30 30 30 2c 31  0,..1000000000,1
0e60: 30 30 30 30 30 30 30 30 30 2c 31 30 30 30 30 30  000000000,100000
0e70: 30 30 30 30 2c 31 30 30 30 30 30 30 30 30 30 2c  0000,1000000000,
0e80: 31 30 30 30 30 30 30 30 30 30 2c 0d 0a 31 30 30  1000000000,..100
0e90: 30 30 30 30 30 30 30 2c 31 30 30 30 30 30 30 30  0000000,10000000
0ea0: 30 30 2c 31 30 30 30 30 30 30 30 30 30 2c 31 30  00,1000000000,10
0eb0: 30 30 30 30 30 30 30 30 2c 31 30 30 30 30 30 30  00000000,1000000
0ec0: 30 30 30 2c 0d 0a 31 30 30 30 30 30 30 30 30 30  000,..1000000000
0ed0: 2c 31 30 30 30 30 30 30 30 30 30 2c 31 30 30 30  ,1000000000,1000
0ee0: 30 30 30 30 30 30 2c 31 30 30 30 30 30 30 30 30  000000,100000000
0ef0: 30 2c 31 30 30 30 30 30 30 30 30 30 7d 3b 0d 0a  0,1000000000};..
0f00: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20  .  vector <int> 
0f10: 71 28 71 5f 2c 20 71 5f 2b 73 69 7a 65 6f 66 28  q(q_, q_+sizeof(
0f20: 71 5f 29 2f 73 69 7a 65 6f 66 28 2a 71 5f 29 29  q_)/sizeof(*q_))
0f30: 3b 20 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20  ; ...double _ = 
0f40: 33 2e 35 39 37 37 34 30 39 32 34 34 39 31 36 33  3.59774092449163
0f50: 38 45 39 3b 20 0d 0a 45 4e 44 0d 0a 2f 2a 0d 0a  8E9; ..END../*..
0f60: 43 41 53 45 28 36 29 0d 0a 09 69 6e 74 20 70 5f  CASE(6)...int p_
0f70: 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 65 63 74 6f  [] = ;...  vecto
0f80: 72 20 3c 69 6e 74 3e 20 70 28 70 5f 2c 20 70 5f  r <int> p(p_, p_
0f90: 2b 73 69 7a 65 6f 66 28 70 5f 29 2f 73 69 7a 65  +sizeof(p_)/size
0fa0: 6f 66 28 2a 70 5f 29 29 3b 20 0d 0a 09 69 6e 74  of(*p_)); ...int
0fb0: 20 71 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 65   q_[] = ;...  ve
0fc0: 63 74 6f 72 20 3c 69 6e 74 3e 20 71 28 71 5f 2c  ctor <int> q(q_,
0fd0: 20 71 5f 2b 73 69 7a 65 6f 66 28 71 5f 29 2f 73   q_+sizeof(q_)/s
0fe0: 69 7a 65 6f 66 28 2a 71 5f 29 29 3b 20 0d 0a 09  izeof(*q_)); ...
0ff0: 64 6f 75 62 6c 65 20 5f 20 3d 20 3b 20 0d 0a 45  double _ = ; ..E
1000: 4e 44 0d 0a 43 41 53 45 28 37 29 0d 0a 09 69 6e  ND..CASE(7)...in
1010: 74 20 70 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76  t p_[] = ;...  v
1020: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 70 28 70 5f  ector <int> p(p_
1030: 2c 20 70 5f 2b 73 69 7a 65 6f 66 28 70 5f 29 2f  , p_+sizeof(p_)/
1040: 73 69 7a 65 6f 66 28 2a 70 5f 29 29 3b 20 0d 0a  sizeof(*p_)); ..
1050: 09 69 6e 74 20 71 5f 5b 5d 20 3d 20 3b 0d 0a 09  .int q_[] = ;...
1060: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 71    vector <int> q
1070: 28 71 5f 2c 20 71 5f 2b 73 69 7a 65 6f 66 28 71  (q_, q_+sizeof(q
1080: 5f 29 2f 73 69 7a 65 6f 66 28 2a 71 5f 29 29 3b  _)/sizeof(*q_));
1090: 20 0d 0a 09 64 6f 75 62 6c 65 20 5f 20 3d 20 3b   ...double _ = ;
10a0: 20 0d 0a 45 4e 44 0d 0a 2a 2f 0d 0a 7d 0d 0a 2f   ..END..*/..}../
10b0: 2f 20 45 4e 44 20 43 55 54 20 48 45 52 45 0d 0a  / END CUT HERE..