Hex Artifact Content
Not logged in

Artifact 3680215ae4f12899056694d2a64d0b272945fadf:


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 73 74 61 74 69 63 20 63  CMP;....static c
01a0: 6f 6e 73 74 20 69 6e 74 20 4d 4f 44 56 41 4c 20  onst int MODVAL 
01b0: 3d 20 31 30 30 30 30 30 30 30 30 37 3b 20 2f 2f  = 1000000007; //
01c0: 20 6d 75 73 74 20 62 65 20 70 72 69 6d 65 20 66   must be prime f
01d0: 6f 72 20 6f 70 2f 0d 0a 73 74 72 75 63 74 20 6d  or op/..struct m
01e0: 69 6e 74 0d 0a 7b 0d 0a 09 69 6e 74 20 76 61 6c  int..{...int val
01f0: 3b 0d 0a 09 6d 69 6e 74 28 29 3a 76 61 6c 28 30  ;...mint():val(0
0200: 29 7b 7d 0d 0a 09 6d 69 6e 74 28 69 6e 74 20 20  ){}...mint(int  
0210: 20 20 78 29 3a 76 61 6c 28 78 25 4d 4f 44 56 41    x):val(x%MODVA
0220: 4c 29 20 7b 7d 0d 0a 09 6d 69 6e 74 28 73 69 7a  L) {}...mint(siz
0230: 65 5f 74 20 78 29 3a 76 61 6c 28 78 25 4d 4f 44  e_t x):val(x%MOD
0240: 56 41 4c 29 20 7b 7d 0d 0a 09 6d 69 6e 74 28 4c  VAL) {}...mint(L
0250: 4c 20 20 20 20 20 78 29 3a 76 61 6c 28 78 25 4d  L     x):val(x%M
0260: 4f 44 56 41 4c 29 20 7b 7d 0d 0a 7d 3b 0d 0a 6d  ODVAL) {}..};..m
0270: 69 6e 74 26 20 6f 70 65 72 61 74 6f 72 2b 3d 28  int& operator+=(
0280: 6d 69 6e 74 26 20 78 2c 20 6d 69 6e 74 20 79 29  mint& x, mint y)
0290: 20 7b 20 72 65 74 75 72 6e 20 78 20 3d 20 78 2e   { return x = x.
02a0: 76 61 6c 2b 79 2e 76 61 6c 3b 20 7d 0d 0a 6d 69  val+y.val; }..mi
02b0: 6e 74 26 20 6f 70 65 72 61 74 6f 72 2d 3d 28 6d  nt& operator-=(m
02c0: 69 6e 74 26 20 78 2c 20 6d 69 6e 74 20 79 29 20  int& x, mint y) 
02d0: 7b 20 72 65 74 75 72 6e 20 78 20 3d 20 78 2e 76  { return x = x.v
02e0: 61 6c 2d 79 2e 76 61 6c 2b 4d 4f 44 56 41 4c 3b  al-y.val+MODVAL;
02f0: 20 7d 0d 0a 6d 69 6e 74 26 20 6f 70 65 72 61 74   }..mint& operat
0300: 6f 72 2a 3d 28 6d 69 6e 74 26 20 78 2c 20 6d 69  or*=(mint& x, mi
0310: 6e 74 20 79 29 20 7b 20 72 65 74 75 72 6e 20 78  nt y) { return x
0320: 20 3d 20 4c 4c 28 78 2e 76 61 6c 29 2a 79 2e 76   = LL(x.val)*y.v
0330: 61 6c 3b 20 7d 0d 0a 6d 69 6e 74 20 50 4f 57 28  al; }..mint POW(
0340: 6d 69 6e 74 20 78 2c 20 4c 4c 20 65 29 20 7b 20  mint x, LL e) { 
0350: 6d 69 6e 74 20 76 3d 31 3b 20 66 6f 72 28 3b 65  mint v=1; for(;e
0360: 3b 78 2a 3d 78 2c 65 3e 3e 3d 31 29 20 69 66 28  ;x*=x,e>>=1) if(
0370: 65 26 31 29 20 76 2a 3d 78 3b 20 72 65 74 75 72  e&1) v*=x; retur
0380: 6e 20 76 3b 20 7d 0d 0a 6d 69 6e 74 26 20 6f 70  n v; }..mint& op
0390: 65 72 61 74 6f 72 2f 3d 28 6d 69 6e 74 26 20 78  erator/=(mint& x
03a0: 2c 20 6d 69 6e 74 20 79 29 20 7b 20 72 65 74 75  , mint y) { retu
03b0: 72 6e 20 78 20 2a 3d 20 50 4f 57 28 79 2c 20 4d  rn x *= POW(y, M
03c0: 4f 44 56 41 4c 2d 32 29 3b 20 7d 0d 0a 6d 69 6e  ODVAL-2); }..min
03d0: 74 20 6f 70 65 72 61 74 6f 72 2b 28 6d 69 6e 74  t operator+(mint
03e0: 20 78 2c 20 6d 69 6e 74 20 79 29 20 7b 20 72 65   x, mint y) { re
03f0: 74 75 72 6e 20 78 2b 3d 79 3b 20 7d 0d 0a 6d 69  turn x+=y; }..mi
0400: 6e 74 20 6f 70 65 72 61 74 6f 72 2d 28 6d 69 6e  nt operator-(min
0410: 74 20 78 2c 20 6d 69 6e 74 20 79 29 20 7b 20 72  t x, mint y) { r
0420: 65 74 75 72 6e 20 78 2d 3d 79 3b 20 7d 0d 0a 6d  eturn x-=y; }..m
0430: 69 6e 74 20 6f 70 65 72 61 74 6f 72 2a 28 6d 69  int operator*(mi
0440: 6e 74 20 78 2c 20 6d 69 6e 74 20 79 29 20 7b 20  nt x, mint y) { 
0450: 72 65 74 75 72 6e 20 78 2a 3d 79 3b 20 7d 0d 0a  return x*=y; }..
0460: 6d 69 6e 74 20 6f 70 65 72 61 74 6f 72 2f 28 6d  mint operator/(m
0470: 69 6e 74 20 78 2c 20 6d 69 6e 74 20 79 29 20 7b  int x, mint y) {
0480: 20 72 65 74 75 72 6e 20 78 2f 3d 79 3b 20 7d 0d   return x/=y; }.
0490: 0a 76 65 63 74 6f 72 3c 6d 69 6e 74 3e 20 46 41  .vector<mint> FA
04a0: 43 5f 28 31 2c 31 29 3b 0d 0a 6d 69 6e 74 20 46  C_(1,1);..mint F
04b0: 41 43 28 4c 4c 20 6e 29 20 7b 20 77 68 69 6c 65  AC(LL n) { while
04c0: 28 20 46 41 43 5f 2e 73 69 7a 65 28 29 3c 3d 6e  ( FAC_.size()<=n
04d0: 20 29 20 46 41 43 5f 2e 70 75 73 68 5f 62 61 63   ) FAC_.push_bac
04e0: 6b 28 20 46 41 43 5f 2e 62 61 63 6b 28 29 2a 46  k( FAC_.back()*F
04f0: 41 43 5f 2e 73 69 7a 65 28 29 20 29 3b 20 72 65  AC_.size() ); re
0500: 74 75 72 6e 20 46 41 43 5f 5b 6e 5d 3b 20 7d 0d  turn FAC_[n]; }.
0510: 0a 6d 69 6e 74 20 43 28 4c 4c 20 6e 2c 20 4c 4c  .mint C(LL n, LL
0520: 20 6b 29 20 7b 20 72 65 74 75 72 6e 20 6b 3c 30   k) { return k<0
0530: 20 7c 7c 20 6e 3c 6b 20 3f 20 30 20 3a 20 46 41   || n<k ? 0 : FA
0540: 43 28 6e 29 20 2f 20 28 46 41 43 28 6b 29 20 2a  C(n) / (FAC(k) *
0550: 20 46 41 43 28 6e 2d 6b 29 29 3b 20 7d 0d 0a 0d   FAC(n-k)); }...
0560: 0a 74 65 6d 70 6c 61 74 65 3c 74 79 70 65 6e 61  .template<typena
0570: 6d 65 20 54 3e 0d 0a 73 74 72 75 63 74 20 44 50  me T>..struct DP
0580: 34 0d 0a 7b 0d 0a 09 69 6e 74 20 4e 31 2c 20 4e  4..{...int N1, N
0590: 32 2c 20 4e 33 2c 20 4e 34 3b 0d 0a 09 76 65 63  2, N3, N4;...vec
05a0: 74 6f 72 3c 54 3e 20 64 61 74 61 3b 0d 0a 09 44  tor<T> data;...D
05b0: 50 34 28 69 6e 74 20 4e 31 2c 20 69 6e 74 20 4e  P4(int N1, int N
05c0: 32 2c 20 69 6e 74 20 4e 33 2c 20 69 6e 74 20 4e  2, int N3, int N
05d0: 34 2c 20 63 6f 6e 73 74 20 54 26 20 74 20 3d 20  4, const T& t = 
05e0: 54 28 29 29 0d 0a 09 09 3a 20 4e 31 28 4e 31 29  T())....: N1(N1)
05f0: 2c 20 4e 32 28 4e 32 29 2c 20 4e 33 28 4e 33 29  , N2(N2), N3(N3)
0600: 2c 20 4e 34 28 4e 34 29 2c 20 64 61 74 61 28 4e  , N4(N4), data(N
0610: 31 2a 4e 32 2a 4e 33 2a 4e 34 2c 20 74 29 20 7b  1*N2*N3*N4, t) {
0620: 20 61 73 73 65 72 74 28 64 61 74 61 2e 73 69 7a   assert(data.siz
0630: 65 28 29 2a 73 69 7a 65 6f 66 28 54 29 3c 28 31  e()*sizeof(T)<(1
0640: 3c 3c 32 36 29 29 3b 20 7d 0d 0a 09 54 26 20 6f  <<26)); }...T& o
0650: 70 65 72 61 74 6f 72 28 29 28 69 6e 74 20 69 31  perator()(int i1
0660: 2c 20 69 6e 74 20 69 32 2c 20 69 6e 74 20 69 33  , int i2, int i3
0670: 2c 20 69 6e 74 20 69 34 29 0d 0a 09 09 7b 20 72  , int i4)....{ r
0680: 65 74 75 72 6e 20 64 61 74 61 5b 20 28 28 28 69  eturn data[ (((i
0690: 31 2a 4e 32 29 2b 69 32 29 2a 4e 33 2b 69 33 29  1*N2)+i2)*N3+i3)
06a0: 2a 4e 34 2b 69 34 20 5d 3b 20 7d 0d 0a 09 76 6f  *N4+i4 ]; }...vo
06b0: 69 64 20 73 77 61 70 28 44 50 34 26 20 72 68 73  id swap(DP4& rhs
06c0: 29 0d 0a 09 09 7b 20 64 61 74 61 2e 73 77 61 70  )....{ data.swap
06d0: 28 72 68 73 2e 64 61 74 61 29 3b 20 7d 0d 0a 7d  (rhs.data); }..}
06e0: 3b 0d 0a 0d 0a 63 6c 61 73 73 20 53 52 4d 49 6e  ;....class SRMIn
06f0: 74 65 72 6d 69 73 73 69 6f 6e 50 68 61 73 65 20  termissionPhase 
0700: 7b 20 70 75 62 6c 69 63 3a 0d 0a 09 69 6e 74 20  { public:...int 
0710: 63 6f 75 6e 74 57 61 79 73 28 76 65 63 74 6f 72  countWays(vector
0720: 20 3c 69 6e 74 3e 20 70 6f 69 6e 74 73 2c 20 76   <int> points, v
0730: 65 63 74 6f 72 20 3c 73 74 72 69 6e 67 3e 20 64  ector <string> d
0740: 65 73 63 72 69 70 74 69 6f 6e 29 0d 0a 09 7b 0d  escription)...{.
0750: 0a 09 09 63 6f 6e 73 74 20 69 6e 74 20 54 4f 54  ...const int TOT
0760: 41 4c 20 3d 20 61 63 63 75 6d 75 6c 61 74 65 28  AL = accumulate(
0770: 70 6f 69 6e 74 73 2e 62 65 67 69 6e 28 29 2c 20  points.begin(), 
0780: 70 6f 69 6e 74 73 2e 65 6e 64 28 29 2c 20 30 29  points.end(), 0)
0790: 3b 0d 0a 0d 0a 09 09 44 50 34 3c 6d 69 6e 74 3e  ;......DP4<mint>
07a0: 20 70 61 74 6e 75 6d 28 32 2c 32 2c 32 2c 54 4f   patnum(2,2,2,TO
07b0: 54 41 4c 2b 31 29 3b 0d 0a 09 09 69 6e 74 20 61  TAL+1);....int a
07c0: 5b 33 5d 3b 0d 0a 09 09 66 6f 72 28 61 5b 30 5d  [3];....for(a[0]
07d0: 3d 30 3b 20 61 5b 30 5d 3c 32 3b 20 2b 2b 61 5b  =0; a[0]<2; ++a[
07e0: 30 5d 29 0d 0a 09 09 66 6f 72 28 61 5b 31 5d 3d  0])....for(a[1]=
07f0: 30 3b 20 61 5b 31 5d 3c 32 3b 20 2b 2b 61 5b 31  0; a[1]<2; ++a[1
0800: 5d 29 0d 0a 09 09 66 6f 72 28 61 5b 32 5d 3d 30  ])....for(a[2]=0
0810: 3b 20 61 5b 32 5d 3c 32 3b 20 2b 2b 61 5b 32 5d  ; a[2]<2; ++a[2]
0820: 29 0d 0a 09 09 7b 0d 0a 09 09 09 70 61 74 6e 75  )....{.....patnu
0830: 6d 28 61 5b 30 5d 2c 61 5b 31 5d 2c 61 5b 32 5d  m(a[0],a[1],a[2]
0840: 2c 30 29 20 3d 20 31 3b 0d 0a 09 09 09 66 6f 72  ,0) = 1;.....for
0850: 28 69 6e 74 20 69 3d 30 3b 20 69 3c 33 3b 20 2b  (int i=0; i<3; +
0860: 2b 69 29 20 69 66 28 61 5b 69 5d 29 20 7b 0d 0a  +i) if(a[i]) {..
0870: 09 09 09 09 66 6f 72 28 69 6e 74 20 70 3d 31 3b  ....for(int p=1;
0880: 20 70 3c 3d 54 4f 54 41 4c 3b 20 2b 2b 70 29 0d   p<=TOTAL; ++p).
0890: 0a 09 09 09 09 09 70 61 74 6e 75 6d 28 61 5b 30  ......patnum(a[0
08a0: 5d 2c 61 5b 31 5d 2c 61 5b 32 5d 2c 70 29 20 2b  ],a[1],a[2],p) +
08b0: 3d 20 70 61 74 6e 75 6d 28 61 5b 30 5d 2c 61 5b  = patnum(a[0],a[
08c0: 31 5d 2c 61 5b 32 5d 2c 70 2d 31 29 3b 0d 0a 09  1],a[2],p-1);...
08d0: 09 09 09 66 6f 72 28 69 6e 74 20 70 3d 54 4f 54  ...for(int p=TOT
08e0: 41 4c 3b 20 70 3e 3d 70 6f 69 6e 74 73 5b 69 5d  AL; p>=points[i]
08f0: 2b 31 3b 20 2d 2d 70 29 0d 0a 09 09 09 09 09 70  +1; --p).......p
0900: 61 74 6e 75 6d 28 61 5b 30 5d 2c 61 5b 31 5d 2c  atnum(a[0],a[1],
0910: 61 5b 32 5d 2c 70 29 20 3d 20 70 61 74 6e 75 6d  a[2],p) = patnum
0920: 28 61 5b 30 5d 2c 61 5b 31 5d 2c 61 5b 32 5d 2c  (a[0],a[1],a[2],
0930: 70 2d 31 29 20 2d 20 70 61 74 6e 75 6d 28 61 5b  p-1) - patnum(a[
0940: 30 5d 2c 61 5b 31 5d 2c 61 5b 32 5d 2c 70 2d 70  0],a[1],a[2],p-p
0950: 6f 69 6e 74 73 5b 69 5d 2d 31 29 3b 0d 0a 09 09  oints[i]-1);....
0960: 09 09 66 6f 72 28 69 6e 74 20 70 3d 70 6f 69 6e  ..for(int p=poin
0970: 74 73 5b 69 5d 3b 20 70 3e 3d 30 3b 20 2d 2d 70  ts[i]; p>=0; --p
0980: 29 0d 0a 09 09 09 09 09 70 61 74 6e 75 6d 28 61  ).......patnum(a
0990: 5b 30 5d 2c 61 5b 31 5d 2c 61 5b 32 5d 2c 70 29  [0],a[1],a[2],p)
09a0: 20 3d 20 70 3e 30 20 3f 20 70 61 74 6e 75 6d 28   = p>0 ? patnum(
09b0: 61 5b 30 5d 2c 61 5b 31 5d 2c 61 5b 32 5d 2c 70  a[0],a[1],a[2],p
09c0: 2d 31 29 20 3a 20 30 3b 0d 0a 09 09 09 7d 0d 0a  -1) : 0;.....}..
09d0: 09 09 7d 0d 0a 0d 0a 09 09 76 65 63 74 6f 72 3c  ..}......vector<
09e0: 6d 69 6e 74 3e 20 64 70 28 54 4f 54 41 4c 2b 32  mint> dp(TOTAL+2
09f0: 29 3b 0d 0a 09 09 64 70 5b 30 5d 20 3d 20 31 3b  );....dp[0] = 1;
0a00: 20 2f 2f 20 2d 31 0d 0a 09 09 66 6f 72 28 69 6e   // -1....for(in
0a10: 74 20 69 3d 64 65 73 63 72 69 70 74 69 6f 6e 2e  t i=description.
0a20: 73 69 7a 65 28 29 2d 31 3b 20 69 3e 3d 30 3b 20  size()-1; i>=0; 
0a30: 2d 2d 69 29 0d 0a 09 09 7b 0d 0a 09 09 09 70 61  --i)....{.....pa
0a40: 72 74 69 61 6c 5f 73 75 6d 28 64 70 2e 62 65 67  rtial_sum(dp.beg
0a50: 69 6e 28 29 2c 20 64 70 2e 65 6e 64 28 29 2c 20  in(), dp.end(), 
0a60: 64 70 2e 62 65 67 69 6e 28 29 29 3b 0d 0a 09 09  dp.begin());....
0a70: 09 66 6f 72 28 69 6e 74 20 73 3d 54 4f 54 41 4c  .for(int s=TOTAL
0a80: 3b 20 73 3e 3d 30 3b 20 2d 2d 73 29 0d 0a 09 09  ; s>=0; --s)....
0a90: 09 09 64 70 5b 73 2b 31 5d 20 3d 20 64 70 5b 73  ..dp[s+1] = dp[s
0aa0: 5d 20 2a 20 70 61 74 6e 75 6d 28 0d 0a 09 09 09  ] * patnum(.....
0ab0: 09 09 64 65 73 63 72 69 70 74 69 6f 6e 5b 69 5d  ..description[i]
0ac0: 5b 30 5d 3d 3d 27 59 27 2c 0d 0a 09 09 09 09 09  [0]=='Y',.......
0ad0: 64 65 73 63 72 69 70 74 69 6f 6e 5b 69 5d 5b 31  description[i][1
0ae0: 5d 3d 3d 27 59 27 2c 0d 0a 09 09 09 09 09 64 65  ]=='Y',.......de
0af0: 73 63 72 69 70 74 69 6f 6e 5b 69 5d 5b 32 5d 3d  scription[i][2]=
0b00: 3d 27 59 27 2c 0d 0a 09 09 09 09 09 73 0d 0a 09  ='Y',.......s...
0b10: 09 09 09 29 3b 0d 0a 09 09 09 64 70 5b 30 5d 20  ...);.....dp[0] 
0b20: 3d 20 30 3b 0d 0a 09 09 7d 0d 0a 09 09 72 65 74  = 0;....}....ret
0b30: 75 72 6e 20 61 63 63 75 6d 75 6c 61 74 65 28 64  urn accumulate(d
0b40: 70 2e 62 65 67 69 6e 28 29 2b 31 2c 20 64 70 2e  p.begin()+1, dp.
0b50: 65 6e 64 28 29 2c 20 6d 69 6e 74 28 30 29 29 2e  end(), mint(0)).
0b60: 76 61 6c 3b 0d 0a 09 7d 0d 0a 7d 3b 0d 0a 0d 0a  val;...}..};....
0b70: 2f 2f 20 42 45 47 49 4e 20 43 55 54 20 48 45 52  // BEGIN CUT HER
0b80: 45 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 63 74 69  E..#include <cti
0b90: 6d 65 3e 0d 0a 64 6f 75 62 6c 65 20 73 74 61 72  me>..double star
0ba0: 74 5f 74 69 6d 65 3b 20 73 74 72 69 6e 67 20 74  t_time; string t
0bb0: 69 6d 65 72 28 29 0d 0a 20 7b 20 6f 73 74 72 69  imer().. { ostri
0bc0: 6e 67 73 74 72 65 61 6d 20 6f 73 3b 20 6f 73 20  ngstream os; os 
0bd0: 3c 3c 20 22 20 28 22 20 3c 3c 20 69 6e 74 28 28  << " (" << int((
0be0: 63 6c 6f 63 6b 28 29 2d 73 74 61 72 74 5f 74 69  clock()-start_ti
0bf0: 6d 65 29 2f 43 4c 4f 43 4b 53 5f 50 45 52 5f 53  me)/CLOCKS_PER_S
0c00: 45 43 2a 31 30 30 30 29 20 3c 3c 20 22 20 6d 73  EC*1000) << " ms
0c10: 65 63 29 22 3b 20 72 65 74 75 72 6e 20 6f 73 2e  ec)"; return os.
0c20: 73 74 72 28 29 3b 20 7d 0d 0a 74 65 6d 70 6c 61  str(); }..templa
0c30: 74 65 3c 74 79 70 65 6e 61 6d 65 20 54 3e 20 6f  te<typename T> o
0c40: 73 74 72 65 61 6d 26 20 6f 70 65 72 61 74 6f 72  stream& operator
0c50: 3c 3c 28 6f 73 74 72 65 61 6d 26 20 6f 73 2c 20  <<(ostream& os, 
0c60: 63 6f 6e 73 74 20 76 65 63 74 6f 72 3c 54 3e 26  const vector<T>&
0c70: 20 76 29 0d 0a 20 7b 20 6f 73 20 3c 3c 20 22 7b   v).. { os << "{
0c80: 20 22 3b 0d 0a 20 20 20 66 6f 72 28 74 79 70 65   ";..   for(type
0c90: 6e 61 6d 65 20 76 65 63 74 6f 72 3c 54 3e 3a 3a  name vector<T>::
0ca0: 63 6f 6e 73 74 5f 69 74 65 72 61 74 6f 72 20 69  const_iterator i
0cb0: 74 3d 76 2e 62 65 67 69 6e 28 29 3b 20 69 74 21  t=v.begin(); it!
0cc0: 3d 76 2e 65 6e 64 28 29 3b 20 2b 2b 69 74 29 0d  =v.end(); ++it).
0cd0: 0a 20 20 20 6f 73 20 3c 3c 20 27 5c 22 27 20 3c  .   os << '\"' <
0ce0: 3c 20 2a 69 74 20 3c 3c 20 27 5c 22 27 20 3c 3c  < *it << '\"' <<
0cf0: 20 28 69 74 2b 31 3d 3d 76 2e 65 6e 64 28 29 20   (it+1==v.end() 
0d00: 3f 20 22 22 20 3a 20 22 2c 20 22 29 3b 20 6f 73  ? "" : ", "); os
0d10: 20 3c 3c 20 22 20 7d 22 3b 20 72 65 74 75 72 6e   << " }"; return
0d20: 20 6f 73 3b 20 7d 0d 0a 76 6f 69 64 20 76 65 72   os; }..void ver
0d30: 69 66 79 5f 63 61 73 65 28 63 6f 6e 73 74 20 69  ify_case(const i
0d40: 6e 74 26 20 45 78 70 65 63 74 65 64 2c 20 63 6f  nt& Expected, co
0d50: 6e 73 74 20 69 6e 74 26 20 52 65 63 65 69 76 65  nst int& Receive
0d60: 64 29 20 7b 0d 0a 20 62 6f 6f 6c 20 6f 6b 20 3d  d) {.. bool ok =
0d70: 20 28 45 78 70 65 63 74 65 64 20 3d 3d 20 52 65   (Expected == Re
0d80: 63 65 69 76 65 64 29 3b 0d 0a 20 69 66 28 6f 6b  ceived);.. if(ok
0d90: 29 20 63 65 72 72 20 3c 3c 20 22 50 41 53 53 45  ) cerr << "PASSE
0da0: 44 22 20 3c 3c 20 74 69 6d 65 72 28 29 20 3c 3c  D" << timer() <<
0db0: 20 65 6e 64 6c 3b 20 20 65 6c 73 65 20 7b 20 63   endl;  else { c
0dc0: 65 72 72 20 3c 3c 20 22 46 41 49 4c 45 44 22 20  err << "FAILED" 
0dd0: 3c 3c 20 74 69 6d 65 72 28 29 20 3c 3c 20 65 6e  << timer() << en
0de0: 64 6c 3b 0d 0a 20 63 65 72 72 20 3c 3c 20 22 5c  dl;.. cerr << "\
0df0: 74 6f 3a 20 5c 22 22 20 3c 3c 20 45 78 70 65 63  to: \"" << Expec
0e00: 74 65 64 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 65  ted << '\"' << e
0e10: 6e 64 6c 20 3c 3c 20 22 5c 74 78 3a 20 5c 22 22  ndl << "\tx: \""
0e20: 20 3c 3c 20 52 65 63 65 69 76 65 64 20 3c 3c 20   << Received << 
0e30: 27 5c 22 27 20 3c 3c 20 65 6e 64 6c 3b 20 7d 20  '\"' << endl; } 
0e40: 7d 0d 0a 23 64 65 66 69 6e 65 20 43 41 53 45 28  }..#define CASE(
0e50: 4e 29 20 7b 63 65 72 72 20 3c 3c 20 22 54 65 73  N) {cerr << "Tes
0e60: 74 20 43 61 73 65 20 23 22 20 3c 3c 20 4e 20 3c  t Case #" << N <
0e70: 3c 20 22 2e 2e 2e 22 20 3c 3c 20 66 6c 75 73 68  < "..." << flush
0e80: 3b 20 73 74 61 72 74 5f 74 69 6d 65 3d 63 6c 6f  ; start_time=clo
0e90: 63 6b 28 29 3b 0d 0a 23 64 65 66 69 6e 65 20 45  ck();..#define E
0ea0: 4e 44 09 20 76 65 72 69 66 79 5f 63 61 73 65 28  ND. verify_case(
0eb0: 5f 2c 20 53 52 4d 49 6e 74 65 72 6d 69 73 73 69  _, SRMIntermissi
0ec0: 6f 6e 50 68 61 73 65 28 29 2e 63 6f 75 6e 74 57  onPhase().countW
0ed0: 61 79 73 28 70 6f 69 6e 74 73 2c 20 64 65 73 63  ays(points, desc
0ee0: 72 69 70 74 69 6f 6e 29 29 3b 7d 0d 0a 69 6e 74  ription));}..int
0ef0: 20 6d 61 69 6e 28 29 7b 0d 0a 0d 0a 43 41 53 45   main(){....CASE
0f00: 28 34 29 0d 0a 69 6e 74 20 70 6f 69 6e 74 73 5f  (4)..int points_
0f10: 5b 5d 20 3d 20 7b 31 2c 32 2c 33 7d 3b 0d 0a 09  [] = {1,2,3};...
0f20: 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e 20 70    vector <int> p
0f30: 6f 69 6e 74 73 28 70 6f 69 6e 74 73 5f 2c 20 70  oints(points_, p
0f40: 6f 69 6e 74 73 5f 2b 73 69 7a 65 6f 66 28 70 6f  oints_+sizeof(po
0f50: 69 6e 74 73 5f 29 2f 73 69 7a 65 6f 66 28 2a 70  ints_)/sizeof(*p
0f60: 6f 69 6e 74 73 5f 29 29 3b 20 0d 0a 09 20 20 73  oints_)); ...  s
0f70: 74 72 69 6e 67 20 64 65 73 63 72 69 70 74 69 6f  tring descriptio
0f80: 6e 5f 5b 5d 20 3d 20 7b 22 59 4e 4e 22 2c 22 4e  n_[] = {"YNN","N
0f90: 4e 4e 22 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72  NN"};...  vector
0fa0: 20 3c 73 74 72 69 6e 67 3e 20 64 65 73 63 72 69   <string> descri
0fb0: 70 74 69 6f 6e 28 64 65 73 63 72 69 70 74 69 6f  ption(descriptio
0fc0: 6e 5f 2c 20 64 65 73 63 72 69 70 74 69 6f 6e 5f  n_, description_
0fd0: 2b 73 69 7a 65 6f 66 28 64 65 73 63 72 69 70 74  +sizeof(descript
0fe0: 69 6f 6e 5f 29 2f 73 69 7a 65 6f 66 28 2a 64 65  ion_)/sizeof(*de
0ff0: 73 63 72 69 70 74 69 6f 6e 5f 29 29 3b 20 0d 0a  scription_)); ..
1000: 09 69 6e 74 20 5f 20 3d 20 31 3b 20 0d 0a 45 4e  .int _ = 1; ..EN
1010: 44 0d 0a 43 41 53 45 28 30 29 0d 0a 09 69 6e 74  D..CASE(0)...int
1020: 20 70 6f 69 6e 74 73 5f 5b 5d 20 3d 20 7b 32 35   points_[] = {25
1030: 30 30 30 2c 20 35 30 30 30 30 2c 20 31 30 30 30  000, 50000, 1000
1040: 30 30 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20  00};...  vector 
1050: 3c 69 6e 74 3e 20 70 6f 69 6e 74 73 28 70 6f 69  <int> points(poi
1060: 6e 74 73 5f 2c 20 70 6f 69 6e 74 73 5f 2b 73 69  nts_, points_+si
1070: 7a 65 6f 66 28 70 6f 69 6e 74 73 5f 29 2f 73 69  zeof(points_)/si
1080: 7a 65 6f 66 28 2a 70 6f 69 6e 74 73 5f 29 29 3b  zeof(*points_));
1090: 20 0d 0a 09 73 74 72 69 6e 67 20 64 65 73 63 72   ...string descr
10a0: 69 70 74 69 6f 6e 5f 5b 5d 20 3d 20 7b 22 59 4e  iption_[] = {"YN
10b0: 4e 22 2c 0d 0a 20 22 4e 4e 4e 22 7d 3b 0d 0a 09  N",.. "NNN"};...
10c0: 20 20 76 65 63 74 6f 72 20 3c 73 74 72 69 6e 67    vector <string
10d0: 3e 20 64 65 73 63 72 69 70 74 69 6f 6e 28 64 65  > description(de
10e0: 73 63 72 69 70 74 69 6f 6e 5f 2c 20 64 65 73 63  scription_, desc
10f0: 72 69 70 74 69 6f 6e 5f 2b 73 69 7a 65 6f 66 28  ription_+sizeof(
1100: 64 65 73 63 72 69 70 74 69 6f 6e 5f 29 2f 73 69  description_)/si
1110: 7a 65 6f 66 28 2a 64 65 73 63 72 69 70 74 69 6f  zeof(*descriptio
1120: 6e 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 5f 20 3d  n_)); ...int _ =
1130: 20 32 35 30 30 30 3b 20 0d 0a 45 4e 44 0d 0a 43   25000; ..END..C
1140: 41 53 45 28 31 29 0d 0a 09 69 6e 74 20 70 6f 69  ASE(1)...int poi
1150: 6e 74 73 5f 5b 5d 20 3d 20 7b 33 30 30 30 30 2c  nts_[] = {30000,
1160: 20 36 30 30 30 30 2c 20 39 30 30 30 30 7d 3b 0d   60000, 90000};.
1170: 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e  ..  vector <int>
1180: 20 70 6f 69 6e 74 73 28 70 6f 69 6e 74 73 5f 2c   points(points_,
1190: 20 70 6f 69 6e 74 73 5f 2b 73 69 7a 65 6f 66 28   points_+sizeof(
11a0: 70 6f 69 6e 74 73 5f 29 2f 73 69 7a 65 6f 66 28  points_)/sizeof(
11b0: 2a 70 6f 69 6e 74 73 5f 29 29 3b 20 0d 0a 09 73  *points_)); ...s
11c0: 74 72 69 6e 67 20 64 65 73 63 72 69 70 74 69 6f  tring descriptio
11d0: 6e 5f 5b 5d 20 3d 20 7b 22 4e 59 4e 22 2c 0d 0a  n_[] = {"NYN",..
11e0: 20 22 4e 59 4e 22 7d 3b 0d 0a 09 20 20 76 65 63   "NYN"};...  vec
11f0: 74 6f 72 20 3c 73 74 72 69 6e 67 3e 20 64 65 73  tor <string> des
1200: 63 72 69 70 74 69 6f 6e 28 64 65 73 63 72 69 70  cription(descrip
1210: 74 69 6f 6e 5f 2c 20 64 65 73 63 72 69 70 74 69  tion_, descripti
1220: 6f 6e 5f 2b 73 69 7a 65 6f 66 28 64 65 73 63 72  on_+sizeof(descr
1230: 69 70 74 69 6f 6e 5f 29 2f 73 69 7a 65 6f 66 28  iption_)/sizeof(
1240: 2a 64 65 73 63 72 69 70 74 69 6f 6e 5f 29 29 3b  *description_));
1250: 20 0d 0a 09 69 6e 74 20 5f 20 3d 20 37 39 39 39   ...int _ = 7999
1260: 36 39 39 39 33 3b 20 0d 0a 45 4e 44 0d 0a 43 41  69993; ..END..CA
1270: 53 45 28 32 29 0d 0a 09 69 6e 74 20 70 6f 69 6e  SE(2)...int poin
1280: 74 73 5f 5b 5d 20 3d 20 7b 32 35 30 30 30 2c 20  ts_[] = {25000, 
1290: 34 35 30 30 30 2c 20 31 31 30 30 30 30 7d 3b 0d  45000, 110000};.
12a0: 0a 09 20 20 76 65 63 74 6f 72 20 3c 69 6e 74 3e  ..  vector <int>
12b0: 20 70 6f 69 6e 74 73 28 70 6f 69 6e 74 73 5f 2c   points(points_,
12c0: 20 70 6f 69 6e 74 73 5f 2b 73 69 7a 65 6f 66 28   points_+sizeof(
12d0: 70 6f 69 6e 74 73 5f 29 2f 73 69 7a 65 6f 66 28  points_)/sizeof(
12e0: 2a 70 6f 69 6e 74 73 5f 29 29 3b 20 0d 0a 09 73  *points_)); ...s
12f0: 74 72 69 6e 67 20 64 65 73 63 72 69 70 74 69 6f  tring descriptio
1300: 6e 5f 5b 5d 20 3d 20 7b 22 4e 4e 4e 22 2c 0d 0a  n_[] = {"NNN",..
1310: 20 22 59 59 59 22 7d 3b 0d 0a 09 20 20 76 65 63   "YYY"};...  vec
1320: 74 6f 72 20 3c 73 74 72 69 6e 67 3e 20 64 65 73  tor <string> des
1330: 63 72 69 70 74 69 6f 6e 28 64 65 73 63 72 69 70  cription(descrip
1340: 74 69 6f 6e 5f 2c 20 64 65 73 63 72 69 70 74 69  tion_, descripti
1350: 6f 6e 5f 2b 73 69 7a 65 6f 66 28 64 65 73 63 72  on_+sizeof(descr
1360: 69 70 74 69 6f 6e 5f 29 2f 73 69 7a 65 6f 66 28  iption_)/sizeof(
1370: 2a 64 65 73 63 72 69 70 74 69 6f 6e 5f 29 29 3b  *description_));
1380: 20 0d 0a 09 69 6e 74 20 5f 20 3d 20 30 3b 20 0d   ...int _ = 0; .
1390: 0a 45 4e 44 0d 0a 43 41 53 45 28 33 29 0d 0a 09  .END..CASE(3)...
13a0: 69 6e 74 20 70 6f 69 6e 74 73 5f 5b 5d 20 3d 20  int points_[] = 
13b0: 7b 32 35 36 30 30 2c 20 35 31 32 30 30 2c 20 31  {25600, 51200, 1
13c0: 30 32 34 30 30 7d 3b 0d 0a 09 20 20 76 65 63 74  02400};...  vect
13d0: 6f 72 20 3c 69 6e 74 3e 20 70 6f 69 6e 74 73 28  or <int> points(
13e0: 70 6f 69 6e 74 73 5f 2c 20 70 6f 69 6e 74 73 5f  points_, points_
13f0: 2b 73 69 7a 65 6f 66 28 70 6f 69 6e 74 73 5f 29  +sizeof(points_)
1400: 2f 73 69 7a 65 6f 66 28 2a 70 6f 69 6e 74 73 5f  /sizeof(*points_
1410: 29 29 3b 20 0d 0a 09 73 74 72 69 6e 67 20 64 65  )); ...string de
1420: 73 63 72 69 70 74 69 6f 6e 5f 5b 5d 20 3d 20 7b  scription_[] = {
1430: 22 4e 59 59 22 2c 0d 0a 20 22 59 4e 59 22 2c 0d  "NYY",.. "YNY",.
1440: 0a 20 22 59 59 59 22 2c 0d 0a 20 22 59 4e 4e 22  . "YYY",.. "YNN"
1450: 2c 0d 0a 20 22 59 59 4e 22 2c 0d 0a 20 22 4e 4e  ,.. "YYN",.. "NN
1460: 59 22 2c 0d 0a 20 22 4e 59 4e 22 2c 0d 0a 20 22  Y",.. "NYN",.. "
1470: 4e 4e 4e 22 7d 3b 0d 0a 09 20 20 76 65 63 74 6f  NNN"};...  vecto
1480: 72 20 3c 73 74 72 69 6e 67 3e 20 64 65 73 63 72  r <string> descr
1490: 69 70 74 69 6f 6e 28 64 65 73 63 72 69 70 74 69  iption(descripti
14a0: 6f 6e 5f 2c 20 64 65 73 63 72 69 70 74 69 6f 6e  on_, description
14b0: 5f 2b 73 69 7a 65 6f 66 28 64 65 73 63 72 69 70  _+sizeof(descrip
14c0: 74 69 6f 6e 5f 29 2f 73 69 7a 65 6f 66 28 2a 64  tion_)/sizeof(*d
14d0: 65 73 63 72 69 70 74 69 6f 6e 5f 29 29 3b 20 0d  escription_)); .
14e0: 0a 09 69 6e 74 20 5f 20 3d 20 38 36 37 35 36 30  ..int _ = 867560
14f0: 38 30 35 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45  805; ..END..CASE
1500: 28 35 29 0d 0a 69 6e 74 20 70 6f 69 6e 74 73 5f  (5)..int points_
1510: 5b 5d 20 3d 20 7b 33 30 30 30 30 2c 36 30 30 30  [] = {30000,6000
1520: 30 2c 31 31 30 30 30 30 7d 3b 0d 0a 09 20 20 76  0,110000};...  v
1530: 65 63 74 6f 72 20 3c 69 6e 74 3e 20 70 6f 69 6e  ector <int> poin
1540: 74 73 28 70 6f 69 6e 74 73 5f 2c 20 70 6f 69 6e  ts(points_, poin
1550: 74 73 5f 2b 73 69 7a 65 6f 66 28 70 6f 69 6e 74  ts_+sizeof(point
1560: 73 5f 29 2f 73 69 7a 65 6f 66 28 2a 70 6f 69 6e  s_)/sizeof(*poin
1570: 74 73 5f 29 29 3b 20 0d 0a 09 73 74 72 69 6e 67  ts_)); ...string
1580: 20 64 65 73 63 72 69 70 74 69 6f 6e 5f 5b 5d 20   description_[] 
1590: 3d 20 7b 0d 0a 09 09 22 59 59 59 22 2c 0d 0a 09  = {...."YYY",...
15a0: 09 22 59 59 59 22 2c 0d 0a 09 09 22 59 59 59 22  ."YYY",...."YYY"
15b0: 2c 0d 0a 09 09 22 59 59 59 22 2c 0d 0a 09 09 22  ,...."YYY",...."
15c0: 59 59 59 22 2c 0d 0a 09 09 22 59 59 59 22 2c 0d  YYY",...."YYY",.
15d0: 0a 09 09 22 59 59 59 22 2c 0d 0a 09 09 22 59 59  ..."YYY",...."YY
15e0: 59 22 2c 0d 0a 09 09 22 59 59 59 22 2c 0d 0a 09  Y",...."YYY",...
15f0: 09 22 59 59 59 22 2c 0d 0a 09 09 22 59 59 59 22  ."YYY",...."YYY"
1600: 2c 0d 0a 09 09 22 59 59 59 22 2c 0d 0a 09 09 22  ,...."YYY",...."
1610: 59 59 59 22 2c 0d 0a 09 09 22 59 59 59 22 2c 0d  YYY",...."YYY",.
1620: 0a 09 09 22 59 59 59 22 2c 0d 0a 09 09 22 59 59  ..."YYY",...."YY
1630: 59 22 2c 0d 0a 09 09 22 59 59 59 22 2c 0d 0a 09  Y",...."YYY",...
1640: 09 22 59 59 59 22 2c 0d 0a 09 09 22 59 59 59 22  ."YYY",...."YYY"
1650: 2c 0d 0a 09 09 22 59 59 59 22 2c 0d 0a 09 7d 3b  ,...."YYY",...};
1660: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 73 74 72  ...  vector <str
1670: 69 6e 67 3e 20 64 65 73 63 72 69 70 74 69 6f 6e  ing> description
1680: 28 64 65 73 63 72 69 70 74 69 6f 6e 5f 2c 20 64  (description_, d
1690: 65 73 63 72 69 70 74 69 6f 6e 5f 2b 73 69 7a 65  escription_+size
16a0: 6f 66 28 64 65 73 63 72 69 70 74 69 6f 6e 5f 29  of(description_)
16b0: 2f 73 69 7a 65 6f 66 28 2a 64 65 73 63 72 69 70  /sizeof(*descrip
16c0: 74 69 6f 6e 5f 29 29 3b 20 0d 0a 09 69 6e 74 20  tion_)); ...int 
16d0: 5f 20 3d 20 2d 31 3b 20 0d 0a 45 4e 44 0d 0a 7d  _ = -1; ..END..}
16e0: 0d 0a 2f 2f 20 45 4e 44 20 43 55 54 20 48 45 52  ..// END CUT HER
16f0: 45 0d 0a                                         E..