Hex Artifact Content
Not logged in

Artifact 2941affd3a3de0674bdc74f1a528f90f6f5a149d:


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 68 65 55  P;....class TheU
01a0: 6c 74 69 6d 61 74 65 50 61 63 6b 61 67 65 73 20  ltimatePackages 
01b0: 7b 20 70 75 62 6c 69 63 3a 0d 0a 09 76 65 63 74  { public:...vect
01c0: 6f 72 20 3c 69 6e 74 3e 20 63 6f 75 6e 74 28 69  or <int> count(i
01d0: 6e 74 20 4e 2c 20 69 6e 74 20 44 2c 20 76 65 63  nt N, int D, vec
01e0: 74 6f 72 20 3c 69 6e 74 3e 20 58 70 72 65 66 69  tor <int> Xprefi
01f0: 78 2c 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20  x, vector <int> 
0200: 59 70 72 65 66 69 78 2c 20 69 6e 74 20 4c 2c 20  Yprefix, int L, 
0210: 69 6e 74 20 73 65 65 64 29 0d 0a 09 7b 0d 0a 09  int seed)...{...
0220: 09 76 65 63 74 6f 72 3c 69 6e 74 3e 20 58 28 44  .vector<int> X(D
0230: 29 2c 20 59 28 44 29 3b 0d 0a 09 09 69 6e 74 20  ), Y(D);....int 
0240: 58 4c 20 3d 20 58 70 72 65 66 69 78 2e 73 69 7a  XL = Xprefix.siz
0250: 65 28 29 3b 0d 0a 0d 0a 09 09 66 6f 72 20 28 69  e();......for (i
0260: 6e 74 20 69 20 3d 20 30 3b 20 69 3c 58 4c 3b 20  nt i = 0; i<XL; 
0270: 2b 2b 69 29 20 7b 0d 0a 09 09 09 58 5b 69 5d 20  ++i) {.....X[i] 
0280: 3d 20 58 70 72 65 66 69 78 5b 69 5d 3b 0d 0a 09  = Xprefix[i];...
0290: 09 09 59 5b 69 5d 20 3d 20 59 70 72 65 66 69 78  ..Y[i] = Yprefix
02a0: 5b 69 5d 3b 0d 0a 09 09 7d 0d 0a 0d 0a 09 09 6c  [i];....}......l
02b0: 6f 6e 67 20 6c 6f 6e 67 20 73 74 61 74 65 20 3d  ong long state =
02c0: 20 73 65 65 64 3b 0d 0a 09 09 66 6f 72 20 28 69   seed;....for (i
02d0: 6e 74 20 69 20 3d 20 58 4c 3b 20 69 3c 44 3b 20  nt i = XL; i<D; 
02e0: 2b 2b 69 29 20 7b 0d 0a 09 09 09 73 74 61 74 65  ++i) {.....state
02f0: 20 3d 20 28 73 74 61 74 65 20 2a 20 31 31 30 33   = (state * 1103
0300: 35 31 35 32 34 35 20 2b 20 31 32 33 34 35 29 20  515245 + 12345) 
0310: 25 20 28 31 4c 4c 20 3c 3c 20 33 31 29 3b 0d 0a  % (1LL << 31);..
0320: 09 09 09 69 6e 74 20 65 6c 65 6e 20 3d 20 31 20  ...int elen = 1 
0330: 2b 20 73 74 61 74 65 25 4c 3b 0d 0a 09 09 09 73  + state%L;.....s
0340: 74 61 74 65 20 3d 20 28 73 74 61 74 65 20 2a 20  tate = (state * 
0350: 31 31 30 33 35 31 35 32 34 35 20 2b 20 31 32 33  1103515245 + 123
0360: 34 35 29 20 25 20 28 31 4c 4c 20 3c 3c 20 33 31  45) % (1LL << 31
0370: 29 3b 0d 0a 09 09 09 59 5b 69 5d 20 3d 20 73 74  );.....Y[i] = st
0380: 61 74 65 20 25 20 28 4e 20 2d 20 65 6c 65 6e 29  ate % (N - elen)
0390: 3b 0d 0a 09 09 09 58 5b 69 5d 20 3d 20 59 5b 69  ;.....X[i] = Y[i
03a0: 5d 20 2b 20 65 6c 65 6e 3b 0d 0a 09 09 7d 0d 0a  ] + elen;....}..
03b0: 0d 0a 09 09 76 65 63 74 6f 72 3c 69 6e 74 3e 20  ....vector<int> 
03c0: 6d 61 78 64 65 70 28 4e 2c 20 2d 31 29 3b 0d 0a  maxdep(N, -1);..
03d0: 09 09 66 6f 72 20 28 69 6e 74 20 69 20 3d 20 30  ..for (int i = 0
03e0: 3b 20 69 20 3c 20 44 3b 20 2b 2b 69 29 0d 0a 09  ; i < D; ++i)...
03f0: 09 09 6d 61 78 64 65 70 5b 58 5b 69 5d 5d 20 3d  ..maxdep[X[i]] =
0400: 20 6d 61 78 28 6d 61 78 64 65 70 5b 58 5b 69 5d   max(maxdep[X[i]
0410: 5d 2c 20 59 5b 69 5d 29 3b 0d 0a 0d 0a 09 09 76  ], Y[i]);......v
0420: 65 63 74 6f 72 3c 69 6e 74 3e 20 64 65 70 65 64  ector<int> deped
0430: 5f 62 79 5f 61 6c 6c 5f 6c 61 72 67 65 72 3b 0d  _by_all_larger;.
0440: 0a 09 09 69 6e 74 20 63 75 72 20 3d 20 4e 20 2d  ...int cur = N -
0450: 20 31 3b 0d 0a 09 09 66 6f 72 20 28 69 6e 74 20   1;....for (int 
0460: 76 20 3d 20 4e 20 2d 20 31 3b 20 76 20 3e 3d 20  v = N - 1; v >= 
0470: 30 3b 20 2d 2d 76 29 20 7b 0d 0a 09 09 09 69 66  0; --v) {.....if
0480: 20 28 76 20 3d 3d 20 63 75 72 29 0d 0a 09 09 09   (v == cur).....
0490: 09 64 65 70 65 64 5f 62 79 5f 61 6c 6c 5f 6c 61  .deped_by_all_la
04a0: 72 67 65 72 2e 70 75 73 68 5f 62 61 63 6b 28 76  rger.push_back(v
04b0: 29 3b 0d 0a 09 09 09 69 66 20 28 6d 61 78 64 65  );.....if (maxde
04c0: 70 5b 76 5d 20 3c 20 63 75 72 29 0d 0a 09 09 09  p[v] < cur).....
04d0: 09 63 75 72 20 3d 20 6d 61 78 64 65 70 5b 76 5d  .cur = maxdep[v]
04e0: 3b 0d 0a 09 09 7d 0d 0a 0d 0a 09 09 76 65 63 74  ;....}......vect
04f0: 6f 72 3c 69 6e 74 3e 20 6d 69 6e 64 65 70 28 4e  or<int> mindep(N
0500: 2c 20 4e 29 3b 0d 0a 09 09 66 6f 72 20 28 69 6e  , N);....for (in
0510: 74 20 69 20 3d 20 30 3b 20 69 20 3c 20 44 3b 20  t i = 0; i < D; 
0520: 2b 2b 69 29 0d 0a 09 09 09 6d 69 6e 64 65 70 5b  ++i).....mindep[
0530: 59 5b 69 5d 5d 20 3d 20 6d 69 6e 28 6d 69 6e 64  Y[i]] = min(mind
0540: 65 70 5b 59 5b 69 5d 5d 2c 20 58 5b 69 5d 29 3b  ep[Y[i]], X[i]);
0550: 0d 0a 0d 0a 09 09 73 65 74 3c 69 6e 74 3e 20 64  ......set<int> d
0560: 65 70 65 64 5f 62 79 5f 61 6c 6c 5f 73 6d 61 6c  eped_by_all_smal
0570: 6c 65 72 3b 0d 0a 09 09 63 75 72 20 3d 20 30 3b  ler;....cur = 0;
0580: 0d 0a 09 09 66 6f 72 20 28 69 6e 74 20 76 20 3d  ....for (int v =
0590: 20 30 3b 20 76 20 3c 3d 20 4e 2d 31 3b 20 2b 2b   0; v <= N-1; ++
05a0: 76 29 20 7b 0d 0a 09 09 09 69 66 20 28 76 20 3d  v) {.....if (v =
05b0: 3d 20 63 75 72 29 0d 0a 09 09 09 09 64 65 70 65  = cur)......depe
05c0: 64 5f 62 79 5f 61 6c 6c 5f 73 6d 61 6c 6c 65 72  d_by_all_smaller
05d0: 2e 69 6e 73 65 72 74 28 76 29 3b 0d 0a 09 09 09  .insert(v);.....
05e0: 69 66 20 28 63 75 72 20 3c 20 6d 69 6e 64 65 70  if (cur < mindep
05f0: 5b 76 5d 29 0d 0a 09 09 09 09 63 75 72 20 3d 20  [v])......cur = 
0600: 6d 69 6e 64 65 70 5b 76 5d 3b 0d 0a 09 09 7d 0d  mindep[v];....}.
0610: 0a 0d 0a 09 09 69 6e 74 20 73 75 6d 20 3d 20 30  .....int sum = 0
0620: 3b 0d 0a 09 09 69 6e 74 20 63 6e 74 20 3d 20 30  ;....int cnt = 0
0630: 3b 0d 0a 09 09 66 6f 72 20 28 69 6e 74 20 76 20  ;....for (int v 
0640: 3a 20 64 65 70 65 64 5f 62 79 5f 61 6c 6c 5f 6c  : deped_by_all_l
0650: 61 72 67 65 72 29 0d 0a 09 09 09 69 66 20 28 64  arger).....if (d
0660: 65 70 65 64 5f 62 79 5f 61 6c 6c 5f 73 6d 61 6c  eped_by_all_smal
0670: 6c 65 72 2e 63 6f 75 6e 74 28 76 29 29 20 7b 0d  ler.count(v)) {.
0680: 0a 09 09 09 09 73 75 6d 20 2b 3d 20 76 3b 0d 0a  .....sum += v;..
0690: 09 09 09 09 63 6e 74 20 2b 3d 20 31 3b 0d 0a 09  ....cnt += 1;...
06a0: 09 09 7d 0d 0a 0d 0a 09 09 76 65 63 74 6f 72 3c  ..}......vector<
06b0: 69 6e 74 3e 20 61 6e 73 20 3d 20 7b 20 63 6e 74  int> ans = { cnt
06c0: 2c 20 73 75 6d 20 7d 3b 0d 0a 09 09 72 65 74 75  , sum };....retu
06d0: 72 6e 20 61 6e 73 3b 0d 0a 09 7d 0d 0a 7d 3b 0d  rn ans;...}..};.
06e0: 0a 0d 0a 2f 2f 20 42 45 47 49 4e 20 43 55 54 20  ...// BEGIN CUT 
06f0: 48 45 52 45 0d 0a 23 69 6e 63 6c 75 64 65 20 3c  HERE..#include <
0700: 63 74 69 6d 65 3e 0d 0a 64 6f 75 62 6c 65 20 73  ctime>..double s
0710: 74 61 72 74 5f 74 69 6d 65 3b 20 73 74 72 69 6e  tart_time; strin
0720: 67 20 74 69 6d 65 72 28 29 0d 0a 20 7b 20 6f 73  g timer().. { os
0730: 74 72 69 6e 67 73 74 72 65 61 6d 20 6f 73 3b 20  tringstream os; 
0740: 6f 73 20 3c 3c 20 22 20 28 22 20 3c 3c 20 69 6e  os << " (" << in
0750: 74 28 28 63 6c 6f 63 6b 28 29 2d 73 74 61 72 74  t((clock()-start
0760: 5f 74 69 6d 65 29 2f 43 4c 4f 43 4b 53 5f 50 45  _time)/CLOCKS_PE
0770: 52 5f 53 45 43 2a 31 30 30 30 29 20 3c 3c 20 22  R_SEC*1000) << "
0780: 20 6d 73 65 63 29 22 3b 20 72 65 74 75 72 6e 20   msec)"; return 
0790: 6f 73 2e 73 74 72 28 29 3b 20 7d 0d 0a 74 65 6d  os.str(); }..tem
07a0: 70 6c 61 74 65 3c 74 79 70 65 6e 61 6d 65 20 54  plate<typename T
07b0: 3e 20 6f 73 74 72 65 61 6d 26 20 6f 70 65 72 61  > ostream& opera
07c0: 74 6f 72 3c 3c 28 6f 73 74 72 65 61 6d 26 20 6f  tor<<(ostream& o
07d0: 73 2c 20 63 6f 6e 73 74 20 76 65 63 74 6f 72 3c  s, const vector<
07e0: 54 3e 26 20 76 29 0d 0a 20 7b 20 6f 73 20 3c 3c  T>& v).. { os <<
07f0: 20 22 7b 20 22 3b 0d 0a 20 20 20 66 6f 72 28 74   "{ ";..   for(t
0800: 79 70 65 6e 61 6d 65 20 76 65 63 74 6f 72 3c 54  ypename vector<T
0810: 3e 3a 3a 63 6f 6e 73 74 5f 69 74 65 72 61 74 6f  >::const_iterato
0820: 72 20 69 74 3d 76 2e 62 65 67 69 6e 28 29 3b 20  r it=v.begin(); 
0830: 69 74 21 3d 76 2e 65 6e 64 28 29 3b 20 2b 2b 69  it!=v.end(); ++i
0840: 74 29 0d 0a 20 20 20 6f 73 20 3c 3c 20 27 5c 22  t)..   os << '\"
0850: 27 20 3c 3c 20 2a 69 74 20 3c 3c 20 27 5c 22 27  ' << *it << '\"'
0860: 20 3c 3c 20 28 69 74 2b 31 3d 3d 76 2e 65 6e 64   << (it+1==v.end
0870: 28 29 20 3f 20 22 22 20 3a 20 22 2c 20 22 29 3b  () ? "" : ", ");
0880: 20 6f 73 20 3c 3c 20 22 20 7d 22 3b 20 72 65 74   os << " }"; ret
0890: 75 72 6e 20 6f 73 3b 20 7d 0d 0a 76 6f 69 64 20  urn os; }..void 
08a0: 76 65 72 69 66 79 5f 63 61 73 65 28 63 6f 6e 73  verify_case(cons
08b0: 74 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 26 20  t vector <int>& 
08c0: 45 78 70 65 63 74 65 64 2c 20 63 6f 6e 73 74 20  Expected, const 
08d0: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 26 20 52 65  vector <int>& Re
08e0: 63 65 69 76 65 64 29 20 7b 0d 0a 20 62 6f 6f 6c  ceived) {.. bool
08f0: 20 6f 6b 20 3d 20 28 45 78 70 65 63 74 65 64 20   ok = (Expected 
0900: 3d 3d 20 52 65 63 65 69 76 65 64 29 3b 0d 0a 20  == Received);.. 
0910: 69 66 28 6f 6b 29 20 63 65 72 72 20 3c 3c 20 22  if(ok) cerr << "
0920: 50 41 53 53 45 44 22 20 3c 3c 20 74 69 6d 65 72  PASSED" << timer
0930: 28 29 20 3c 3c 20 65 6e 64 6c 3b 20 20 65 6c 73  () << endl;  els
0940: 65 20 7b 20 63 65 72 72 20 3c 3c 20 22 46 41 49  e { cerr << "FAI
0950: 4c 45 44 22 20 3c 3c 20 74 69 6d 65 72 28 29 20  LED" << timer() 
0960: 3c 3c 20 65 6e 64 6c 3b 0d 0a 20 63 65 72 72 20  << endl;.. cerr 
0970: 3c 3c 20 22 5c 74 6f 3a 20 22 20 3c 3c 20 45 78  << "\to: " << Ex
0980: 70 65 63 74 65 64 20 3c 3c 20 65 6e 64 6c 20 3c  pected << endl <
0990: 3c 20 22 5c 74 78 3a 20 22 20 3c 3c 20 52 65 63  < "\tx: " << Rec
09a0: 65 69 76 65 64 20 3c 3c 20 65 6e 64 6c 3b 20 7d  eived << endl; }
09b0: 20 7d 0d 0a 23 64 65 66 69 6e 65 20 43 41 53 45   }..#define CASE
09c0: 28 4e 29 20 7b 63 65 72 72 20 3c 3c 20 22 54 65  (N) {cerr << "Te
09d0: 73 74 20 43 61 73 65 20 23 22 20 3c 3c 20 4e 20  st Case #" << N 
09e0: 3c 3c 20 22 2e 2e 2e 22 20 3c 3c 20 66 6c 75 73  << "..." << flus
09f0: 68 3b 20 73 74 61 72 74 5f 74 69 6d 65 3d 63 6c  h; start_time=cl
0a00: 6f 63 6b 28 29 3b 0d 0a 23 64 65 66 69 6e 65 20  ock();..#define 
0a10: 45 4e 44 09 20 76 65 72 69 66 79 5f 63 61 73 65  END. verify_case
0a20: 28 5f 2c 20 54 68 65 55 6c 74 69 6d 61 74 65 50  (_, TheUltimateP
0a30: 61 63 6b 61 67 65 73 28 29 2e 63 6f 75 6e 74 28  ackages().count(
0a40: 4e 2c 20 44 2c 20 58 70 72 65 66 69 78 2c 20 59  N, D, Xprefix, Y
0a50: 70 72 65 66 69 78 2c 20 4c 2c 20 73 65 65 64 29  prefix, L, seed)
0a60: 29 3b 7d 0d 0a 69 6e 74 20 6d 61 69 6e 28 29 7b  );}..int main(){
0a70: 0d 0a 0d 0a 43 41 53 45 28 30 29 0d 0a 09 69 6e  ....CASE(0)...in
0a80: 74 20 4e 20 3d 20 35 3b 20 0d 0a 09 69 6e 74 20  t N = 5; ...int 
0a90: 44 20 3d 20 34 3b 20 0d 0a 09 69 6e 74 20 58 70  D = 4; ...int Xp
0aa0: 72 65 66 69 78 5f 5b 5d 20 3d 20 7b 31 2c 20 33  refix_[] = {1, 3
0ab0: 2c 20 32 2c 20 34 7d 3b 0d 0a 09 20 20 76 65 63  , 2, 4};...  vec
0ac0: 74 6f 72 20 3c 69 6e 74 3e 20 58 70 72 65 66 69  tor <int> Xprefi
0ad0: 78 28 58 70 72 65 66 69 78 5f 2c 20 58 70 72 65  x(Xprefix_, Xpre
0ae0: 66 69 78 5f 2b 73 69 7a 65 6f 66 28 58 70 72 65  fix_+sizeof(Xpre
0af0: 66 69 78 5f 29 2f 73 69 7a 65 6f 66 28 2a 58 70  fix_)/sizeof(*Xp
0b00: 72 65 66 69 78 5f 29 29 3b 20 0d 0a 09 69 6e 74  refix_)); ...int
0b10: 20 59 70 72 65 66 69 78 5f 5b 5d 20 3d 20 7b 30   Yprefix_[] = {0
0b20: 2c 20 32 2c 20 31 2c 20 33 7d 3b 0d 0a 09 20 20  , 2, 1, 3};...  
0b30: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 59 70 72  vector <int> Ypr
0b40: 65 66 69 78 28 59 70 72 65 66 69 78 5f 2c 20 59  efix(Yprefix_, Y
0b50: 70 72 65 66 69 78 5f 2b 73 69 7a 65 6f 66 28 59  prefix_+sizeof(Y
0b60: 70 72 65 66 69 78 5f 29 2f 73 69 7a 65 6f 66 28  prefix_)/sizeof(
0b70: 2a 59 70 72 65 66 69 78 5f 29 29 3b 20 0d 0a 09  *Yprefix_)); ...
0b80: 69 6e 74 20 4c 20 3d 20 31 3b 20 0d 0a 09 69 6e  int L = 1; ...in
0b90: 74 20 73 65 65 64 20 3d 20 34 37 3b 20 0d 0a 09  t seed = 47; ...
0ba0: 69 6e 74 20 5f 5f 5b 5d 20 3d 20 7b 35 2c 20 31  int __[] = {5, 1
0bb0: 30 20 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20  0 };...  vector 
0bc0: 3c 69 6e 74 3e 20 5f 28 5f 5f 2c 20 5f 5f 2b 73  <int> _(__, __+s
0bd0: 69 7a 65 6f 66 28 5f 5f 29 2f 73 69 7a 65 6f 66  izeof(__)/sizeof
0be0: 28 2a 5f 5f 29 29 3b 20 0d 0a 45 4e 44 0d 0a 43  (*__)); ..END..C
0bf0: 41 53 45 28 31 29 0d 0a 09 69 6e 74 20 4e 20 3d  ASE(1)...int N =
0c00: 20 35 3b 20 0d 0a 09 69 6e 74 20 44 20 3d 20 36   5; ...int D = 6
0c10: 3b 20 0d 0a 09 69 6e 74 20 58 70 72 65 66 69 78  ; ...int Xprefix
0c20: 5f 5b 5d 20 3d 20 7b 31 2c 20 32 2c 20 33 2c 20  _[] = {1, 2, 3, 
0c30: 34 2c 20 34 2c 20 34 7d 3b 0d 0a 09 20 20 76 65  4, 4, 4};...  ve
0c40: 63 74 6f 72 20 3c 69 6e 74 3e 20 58 70 72 65 66  ctor <int> Xpref
0c50: 69 78 28 58 70 72 65 66 69 78 5f 2c 20 58 70 72  ix(Xprefix_, Xpr
0c60: 65 66 69 78 5f 2b 73 69 7a 65 6f 66 28 58 70 72  efix_+sizeof(Xpr
0c70: 65 66 69 78 5f 29 2f 73 69 7a 65 6f 66 28 2a 58  efix_)/sizeof(*X
0c80: 70 72 65 66 69 78 5f 29 29 3b 20 0d 0a 09 69 6e  prefix_)); ...in
0c90: 74 20 59 70 72 65 66 69 78 5f 5b 5d 20 3d 20 7b  t Yprefix_[] = {
0ca0: 30 2c 20 30 2c 20 30 2c 20 31 2c 20 33 2c 20 32  0, 0, 0, 1, 3, 2
0cb0: 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69  };...  vector <i
0cc0: 6e 74 3e 20 59 70 72 65 66 69 78 28 59 70 72 65  nt> Yprefix(Ypre
0cd0: 66 69 78 5f 2c 20 59 70 72 65 66 69 78 5f 2b 73  fix_, Yprefix_+s
0ce0: 69 7a 65 6f 66 28 59 70 72 65 66 69 78 5f 29 2f  izeof(Yprefix_)/
0cf0: 73 69 7a 65 6f 66 28 2a 59 70 72 65 66 69 78 5f  sizeof(*Yprefix_
0d00: 29 29 3b 20 0d 0a 09 69 6e 74 20 4c 20 3d 20 31  )); ...int L = 1
0d10: 3b 20 0d 0a 09 69 6e 74 20 73 65 65 64 20 3d 20  ; ...int seed = 
0d20: 36 34 3b 20 0d 0a 09 69 6e 74 20 5f 5f 5b 5d 20  64; ...int __[] 
0d30: 3d 20 7b 32 2c 20 34 20 7d 3b 0d 0a 09 20 20 76  = {2, 4 };...  v
0d40: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 5f 28 5f 5f  ector <int> _(__
0d50: 2c 20 5f 5f 2b 73 69 7a 65 6f 66 28 5f 5f 29 2f  , __+sizeof(__)/
0d60: 73 69 7a 65 6f 66 28 2a 5f 5f 29 29 3b 20 0d 0a  sizeof(*__)); ..
0d70: 45 4e 44 0d 0a 43 41 53 45 28 32 29 0d 0a 09 69  END..CASE(2)...i
0d80: 6e 74 20 4e 20 3d 20 35 3b 20 0d 0a 09 69 6e 74  nt N = 5; ...int
0d90: 20 44 20 3d 20 34 3b 20 0d 0a 09 69 6e 74 20 58   D = 4; ...int X
0da0: 70 72 65 66 69 78 5f 5b 5d 20 3d 20 7b 32 2c 20  prefix_[] = {2, 
0db0: 32 2c 20 33 2c 20 34 7d 3b 0d 0a 09 20 20 76 65  2, 3, 4};...  ve
0dc0: 63 74 6f 72 20 3c 69 6e 74 3e 20 58 70 72 65 66  ctor <int> Xpref
0dd0: 69 78 28 58 70 72 65 66 69 78 5f 2c 20 58 70 72  ix(Xprefix_, Xpr
0de0: 65 66 69 78 5f 2b 73 69 7a 65 6f 66 28 58 70 72  efix_+sizeof(Xpr
0df0: 65 66 69 78 5f 29 2f 73 69 7a 65 6f 66 28 2a 58  efix_)/sizeof(*X
0e00: 70 72 65 66 69 78 5f 29 29 3b 20 0d 0a 09 69 6e  prefix_)); ...in
0e10: 74 20 59 70 72 65 66 69 78 5f 5b 5d 20 3d 20 7b  t Yprefix_[] = {
0e20: 30 2c 20 31 2c 20 32 2c 20 32 7d 3b 0d 0a 09 20  0, 1, 2, 2};... 
0e30: 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 59 70   vector <int> Yp
0e40: 72 65 66 69 78 28 59 70 72 65 66 69 78 5f 2c 20  refix(Yprefix_, 
0e50: 59 70 72 65 66 69 78 5f 2b 73 69 7a 65 6f 66 28  Yprefix_+sizeof(
0e60: 59 70 72 65 66 69 78 5f 29 2f 73 69 7a 65 6f 66  Yprefix_)/sizeof
0e70: 28 2a 59 70 72 65 66 69 78 5f 29 29 3b 20 0d 0a  (*Yprefix_)); ..
0e80: 09 69 6e 74 20 4c 20 3d 20 31 3b 20 0d 0a 09 69  .int L = 1; ...i
0e90: 6e 74 20 73 65 65 64 20 3d 20 33 32 35 33 32 3b  nt seed = 32532;
0ea0: 20 0d 0a 09 69 6e 74 20 5f 5f 5b 5d 20 3d 20 7b   ...int __[] = {
0eb0: 31 2c 20 32 20 7d 3b 0d 0a 09 20 20 76 65 63 74  1, 2 };...  vect
0ec0: 6f 72 20 3c 69 6e 74 3e 20 5f 28 5f 5f 2c 20 5f  or <int> _(__, _
0ed0: 5f 2b 73 69 7a 65 6f 66 28 5f 5f 29 2f 73 69 7a  _+sizeof(__)/siz
0ee0: 65 6f 66 28 2a 5f 5f 29 29 3b 20 0d 0a 45 4e 44  eof(*__)); ..END
0ef0: 0d 0a 43 41 53 45 28 33 29 0d 0a 09 69 6e 74 20  ..CASE(3)...int 
0f00: 4e 20 3d 20 34 3b 20 0d 0a 09 69 6e 74 20 44 20  N = 4; ...int D 
0f10: 3d 20 33 3b 20 0d 0a 09 69 6e 74 20 58 70 72 65  = 3; ...int Xpre
0f20: 66 69 78 5f 5b 5d 20 3d 20 7b 33 2c 20 33 2c 20  fix_[] = {3, 3, 
0f30: 32 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c  2};...  vector <
0f40: 69 6e 74 3e 20 58 70 72 65 66 69 78 28 58 70 72  int> Xprefix(Xpr
0f50: 65 66 69 78 5f 2c 20 58 70 72 65 66 69 78 5f 2b  efix_, Xprefix_+
0f60: 73 69 7a 65 6f 66 28 58 70 72 65 66 69 78 5f 29  sizeof(Xprefix_)
0f70: 2f 73 69 7a 65 6f 66 28 2a 58 70 72 65 66 69 78  /sizeof(*Xprefix
0f80: 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 59 70 72 65  _)); ...int Ypre
0f90: 66 69 78 5f 5b 5d 20 3d 20 7b 30 2c 20 32 2c 20  fix_[] = {0, 2, 
0fa0: 31 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c  1};...  vector <
0fb0: 69 6e 74 3e 20 59 70 72 65 66 69 78 28 59 70 72  int> Yprefix(Ypr
0fc0: 65 66 69 78 5f 2c 20 59 70 72 65 66 69 78 5f 2b  efix_, Yprefix_+
0fd0: 73 69 7a 65 6f 66 28 59 70 72 65 66 69 78 5f 29  sizeof(Yprefix_)
0fe0: 2f 73 69 7a 65 6f 66 28 2a 59 70 72 65 66 69 78  /sizeof(*Yprefix
0ff0: 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 4c 20 3d 20  _)); ...int L = 
1000: 31 3b 20 0d 0a 09 69 6e 74 20 73 65 65 64 20 3d  1; ...int seed =
1010: 20 32 35 35 35 3b 20 0d 0a 09 69 6e 74 20 5f 5f   2555; ...int __
1020: 5b 5d 20 3d 20 7b 31 2c 20 33 20 7d 3b 0d 0a 09  [] = {1, 3 };...
1030: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 5f    vector <int> _
1040: 28 5f 5f 2c 20 5f 5f 2b 73 69 7a 65 6f 66 28 5f  (__, __+sizeof(_
1050: 5f 29 2f 73 69 7a 65 6f 66 28 2a 5f 5f 29 29 3b  _)/sizeof(*__));
1060: 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 34 29 0d   ..END..CASE(4).
1070: 0a 09 69 6e 74 20 4e 20 3d 20 37 3b 20 0d 0a 09  ..int N = 7; ...
1080: 69 6e 74 20 44 20 3d 20 30 3b 20 0d 0a 09 76 65  int D = 0; ...ve
1090: 63 74 6f 72 20 3c 69 6e 74 3e 20 58 70 72 65 66  ctor <int> Xpref
10a0: 69 78 3b 20 0d 0a 09 76 65 63 74 6f 72 20 3c 69  ix; ...vector <i
10b0: 6e 74 3e 20 59 70 72 65 66 69 78 3b 20 0d 0a 09  nt> Yprefix; ...
10c0: 69 6e 74 20 4c 20 3d 20 31 3b 20 0d 0a 09 69 6e  int L = 1; ...in
10d0: 74 20 73 65 65 64 20 3d 20 30 3b 20 0d 0a 09 69  t seed = 0; ...i
10e0: 6e 74 20 5f 5f 5b 5d 20 3d 20 7b 30 2c 20 30 20  nt __[] = {0, 0 
10f0: 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69  };...  vector <i
1100: 6e 74 3e 20 5f 28 5f 5f 2c 20 5f 5f 2b 73 69 7a  nt> _(__, __+siz
1110: 65 6f 66 28 5f 5f 29 2f 73 69 7a 65 6f 66 28 2a  eof(__)/sizeof(*
1120: 5f 5f 29 29 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53  __)); ..END..CAS
1130: 45 28 35 29 0d 0a 09 69 6e 74 20 4e 20 3d 20 32  E(5)...int N = 2
1140: 3b 20 0d 0a 09 69 6e 74 20 44 20 3d 20 34 3b 20  ; ...int D = 4; 
1150: 0d 0a 09 69 6e 74 20 58 70 72 65 66 69 78 5f 5b  ...int Xprefix_[
1160: 5d 20 3d 20 7b 31 2c 20 31 2c 20 31 2c 20 31 7d  ] = {1, 1, 1, 1}
1170: 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e  ;...  vector <in
1180: 74 3e 20 58 70 72 65 66 69 78 28 58 70 72 65 66  t> Xprefix(Xpref
1190: 69 78 5f 2c 20 58 70 72 65 66 69 78 5f 2b 73 69  ix_, Xprefix_+si
11a0: 7a 65 6f 66 28 58 70 72 65 66 69 78 5f 29 2f 73  zeof(Xprefix_)/s
11b0: 69 7a 65 6f 66 28 2a 58 70 72 65 66 69 78 5f 29  izeof(*Xprefix_)
11c0: 29 3b 20 0d 0a 09 69 6e 74 20 59 70 72 65 66 69  ); ...int Yprefi
11d0: 78 5f 5b 5d 20 3d 20 7b 30 2c 20 30 2c 20 30 2c  x_[] = {0, 0, 0,
11e0: 20 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20   0};...  vector 
11f0: 3c 69 6e 74 3e 20 59 70 72 65 66 69 78 28 59 70  <int> Yprefix(Yp
1200: 72 65 66 69 78 5f 2c 20 59 70 72 65 66 69 78 5f  refix_, Yprefix_
1210: 2b 73 69 7a 65 6f 66 28 59 70 72 65 66 69 78 5f  +sizeof(Yprefix_
1220: 29 2f 73 69 7a 65 6f 66 28 2a 59 70 72 65 66 69  )/sizeof(*Yprefi
1230: 78 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 4c 20 3d  x_)); ...int L =
1240: 20 31 3b 20 0d 0a 09 69 6e 74 20 73 65 65 64 20   1; ...int seed 
1250: 3d 20 30 3b 20 0d 0a 09 69 6e 74 20 5f 5f 5b 5d  = 0; ...int __[]
1260: 20 3d 20 7b 32 2c 20 31 20 7d 3b 0d 0a 09 20 20   = {2, 1 };...  
1270: 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 5f 28 5f  vector <int> _(_
1280: 5f 2c 20 5f 5f 2b 73 69 7a 65 6f 66 28 5f 5f 29  _, __+sizeof(__)
1290: 2f 73 69 7a 65 6f 66 28 2a 5f 5f 29 29 3b 20 0d  /sizeof(*__)); .
12a0: 0a 45 4e 44 0d 0a 43 41 53 45 28 36 29 0d 0a 09  .END..CASE(6)...
12b0: 69 6e 74 20 4e 20 3d 20 37 3b 20 0d 0a 09 69 6e  int N = 7; ...in
12c0: 74 20 44 20 3d 20 31 32 3b 20 0d 0a 09 69 6e 74  t D = 12; ...int
12d0: 20 58 70 72 65 66 69 78 5f 5b 5d 20 3d 20 7b 32   Xprefix_[] = {2
12e0: 2c 20 33 2c 20 34 7d 3b 0d 0a 09 20 20 76 65 63  , 3, 4};...  vec
12f0: 74 6f 72 20 3c 69 6e 74 3e 20 58 70 72 65 66 69  tor <int> Xprefi
1300: 78 28 58 70 72 65 66 69 78 5f 2c 20 58 70 72 65  x(Xprefix_, Xpre
1310: 66 69 78 5f 2b 73 69 7a 65 6f 66 28 58 70 72 65  fix_+sizeof(Xpre
1320: 66 69 78 5f 29 2f 73 69 7a 65 6f 66 28 2a 58 70  fix_)/sizeof(*Xp
1330: 72 65 66 69 78 5f 29 29 3b 20 0d 0a 09 69 6e 74  refix_)); ...int
1340: 20 59 70 72 65 66 69 78 5f 5b 5d 20 3d 20 7b 31   Yprefix_[] = {1
1350: 2c 20 32 2c 20 33 7d 3b 0d 0a 09 20 20 76 65 63  , 2, 3};...  vec
1360: 74 6f 72 20 3c 69 6e 74 3e 20 59 70 72 65 66 69  tor <int> Yprefi
1370: 78 28 59 70 72 65 66 69 78 5f 2c 20 59 70 72 65  x(Yprefix_, Ypre
1380: 66 69 78 5f 2b 73 69 7a 65 6f 66 28 59 70 72 65  fix_+sizeof(Ypre
1390: 66 69 78 5f 29 2f 73 69 7a 65 6f 66 28 2a 59 70  fix_)/sizeof(*Yp
13a0: 72 65 66 69 78 5f 29 29 3b 20 0d 0a 09 69 6e 74  refix_)); ...int
13b0: 20 4c 20 3d 20 35 3b 20 0d 0a 09 69 6e 74 20 73   L = 5; ...int s
13c0: 65 65 64 20 3d 20 34 37 34 37 3b 20 0d 0a 09 69  eed = 4747; ...i
13d0: 6e 74 20 5f 5f 5b 5d 20 3d 20 7b 30 2c 20 30 20  nt __[] = {0, 0 
13e0: 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69  };...  vector <i
13f0: 6e 74 3e 20 5f 28 5f 5f 2c 20 5f 5f 2b 73 69 7a  nt> _(__, __+siz
1400: 65 6f 66 28 5f 5f 29 2f 73 69 7a 65 6f 66 28 2a  eof(__)/sizeof(*
1410: 5f 5f 29 29 3b 20 0d 0a 45 4e 44 0d 0a 2f 2a 0d  __)); ..END../*.
1420: 0a 43 41 53 45 28 37 29 0d 0a 09 69 6e 74 20 4e  .CASE(7)...int N
1430: 20 3d 20 3b 20 0d 0a 09 69 6e 74 20 44 20 3d 20   = ; ...int D = 
1440: 3b 20 0d 0a 09 69 6e 74 20 58 70 72 65 66 69 78  ; ...int Xprefix
1450: 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 65 63 74  _[] = ;...  vect
1460: 6f 72 20 3c 69 6e 74 3e 20 58 70 72 65 66 69 78  or <int> Xprefix
1470: 28 58 70 72 65 66 69 78 5f 2c 20 58 70 72 65 66  (Xprefix_, Xpref
1480: 69 78 5f 2b 73 69 7a 65 6f 66 28 58 70 72 65 66  ix_+sizeof(Xpref
1490: 69 78 5f 29 2f 73 69 7a 65 6f 66 28 2a 58 70 72  ix_)/sizeof(*Xpr
14a0: 65 66 69 78 5f 29 29 3b 20 0d 0a 09 69 6e 74 20  efix_)); ...int 
14b0: 59 70 72 65 66 69 78 5f 5b 5d 20 3d 20 3b 0d 0a  Yprefix_[] = ;..
14c0: 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20  .  vector <int> 
14d0: 59 70 72 65 66 69 78 28 59 70 72 65 66 69 78 5f  Yprefix(Yprefix_
14e0: 2c 20 59 70 72 65 66 69 78 5f 2b 73 69 7a 65 6f  , Yprefix_+sizeo
14f0: 66 28 59 70 72 65 66 69 78 5f 29 2f 73 69 7a 65  f(Yprefix_)/size
1500: 6f 66 28 2a 59 70 72 65 66 69 78 5f 29 29 3b 20  of(*Yprefix_)); 
1510: 0d 0a 09 69 6e 74 20 4c 20 3d 20 3b 20 0d 0a 09  ...int L = ; ...
1520: 69 6e 74 20 73 65 65 64 20 3d 20 3b 20 0d 0a 09  int seed = ; ...
1530: 69 6e 74 20 5f 5f 5b 5d 20 3d 20 3b 0d 0a 09 20  int __[] = ;... 
1540: 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 5f 28   vector <int> _(
1550: 5f 5f 2c 20 5f 5f 2b 73 69 7a 65 6f 66 28 5f 5f  __, __+sizeof(__
1560: 29 2f 73 69 7a 65 6f 66 28 2a 5f 5f 29 29 3b 20  )/sizeof(*__)); 
1570: 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 38 29 0d 0a  ..END..CASE(8)..
1580: 09 69 6e 74 20 4e 20 3d 20 3b 20 0d 0a 09 69 6e  .int N = ; ...in
1590: 74 20 44 20 3d 20 3b 20 0d 0a 09 69 6e 74 20 58  t D = ; ...int X
15a0: 70 72 65 66 69 78 5f 5b 5d 20 3d 20 3b 0d 0a 09  prefix_[] = ;...
15b0: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 58    vector <int> X
15c0: 70 72 65 66 69 78 28 58 70 72 65 66 69 78 5f 2c  prefix(Xprefix_,
15d0: 20 58 70 72 65 66 69 78 5f 2b 73 69 7a 65 6f 66   Xprefix_+sizeof
15e0: 28 58 70 72 65 66 69 78 5f 29 2f 73 69 7a 65 6f  (Xprefix_)/sizeo
15f0: 66 28 2a 58 70 72 65 66 69 78 5f 29 29 3b 20 0d  f(*Xprefix_)); .
1600: 0a 09 69 6e 74 20 59 70 72 65 66 69 78 5f 5b 5d  ..int Yprefix_[]
1610: 20 3d 20 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20   = ;...  vector 
1620: 3c 69 6e 74 3e 20 59 70 72 65 66 69 78 28 59 70  <int> Yprefix(Yp
1630: 72 65 66 69 78 5f 2c 20 59 70 72 65 66 69 78 5f  refix_, Yprefix_
1640: 2b 73 69 7a 65 6f 66 28 59 70 72 65 66 69 78 5f  +sizeof(Yprefix_
1650: 29 2f 73 69 7a 65 6f 66 28 2a 59 70 72 65 66 69  )/sizeof(*Yprefi
1660: 78 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 4c 20 3d  x_)); ...int L =
1670: 20 3b 20 0d 0a 09 69 6e 74 20 73 65 65 64 20 3d   ; ...int seed =
1680: 20 3b 20 0d 0a 09 69 6e 74 20 5f 5f 5b 5d 20 3d   ; ...int __[] =
1690: 20 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 69   ;...  vector <i
16a0: 6e 74 3e 20 5f 28 5f 5f 2c 20 5f 5f 2b 73 69 7a  nt> _(__, __+siz
16b0: 65 6f 66 28 5f 5f 29 2f 73 69 7a 65 6f 66 28 2a  eof(__)/sizeof(*
16c0: 5f 5f 29 29 3b 20 0d 0a 45 4e 44 0d 0a 2a 2f 0d  __)); ..END..*/.
16d0: 0a 7d 0d 0a 2f 2f 20 45 4e 44 20 43 55 54 20 48  .}..// END CUT H
16e0: 45 52 45 0d 0a                                   ERE..