Hex Artifact Content
Not logged in

Artifact c220b48258d1e6335347484fe8dd5d23331c239f:


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 23 69 6e 63 6c 75 64 65 20 3c 63 73  g>..#include <cs
0140: 74 64 6c 69 62 3e 0d 0a 75 73 69 6e 67 20 6e 61  tdlib>..using na
0150: 6d 65 73 70 61 63 65 20 73 74 64 3b 0d 0a 74 79  mespace std;..ty
0160: 70 65 64 65 66 20 6c 6f 6e 67 20 6c 6f 6e 67 20  pedef long long 
0170: 4c 4c 3b 0d 0a 0d 0a 74 65 6d 70 6c 61 74 65 3c  LL;....template<
0180: 74 79 70 65 6e 61 6d 65 20 54 3e 0d 0a 73 74 72  typename T>..str
0190: 75 63 74 20 52 4d 51 0d 0a 7b 0d 0a 09 76 65 63  uct RMQ..{...vec
01a0: 74 6f 72 3c 20 76 65 63 74 6f 72 3c 69 6e 74 3e  tor< vector<int>
01b0: 20 3e 20 72 6d 3b 0d 0a 09 54 2a 20 20 64 3b 0d   > rm;...T*  d;.
01c0: 0a 09 69 6e 74 20 6e 3b 0d 0a 0d 0a 09 52 4d 51  ..int n;.....RMQ
01d0: 28 20 54 2a 20 64 2c 20 69 6e 74 20 6e 20 29 20  ( T* d, int n ) 
01e0: 3a 20 64 28 64 29 2c 20 6e 28 6e 29 20 7b 0d 0a  : d(d), n(n) {..
01f0: 09 09 72 6d 2e 70 75 73 68 5f 62 61 63 6b 28 20  ..rm.push_back( 
0200: 76 65 63 74 6f 72 3c 69 6e 74 3e 28 6e 29 20 29  vector<int>(n) )
0210: 3b 0d 0a 09 09 66 6f 72 28 69 6e 74 20 78 3d 30  ;....for(int x=0
0220: 3b 20 78 3c 6e 3b 20 2b 2b 78 29 0d 0a 09 09 09  ; x<n; ++x).....
0230: 72 6d 5b 30 5d 5b 78 5d 20 3d 20 78 3b 0d 0a 09  rm[0][x] = x;...
0240: 09 66 6f 72 28 69 6e 74 20 6b 3d 31 3b 20 28 31  .for(int k=1; (1
0250: 3c 3c 6b 29 3c 3d 6e 3b 20 2b 2b 6b 29 20 7b 0d  <<k)<=n; ++k) {.
0260: 0a 09 09 09 72 6d 2e 70 75 73 68 5f 62 61 63 6b  ....rm.push_back
0270: 28 20 72 6d 5b 6b 2d 31 5d 20 29 3b 0d 0a 09 09  ( rm[k-1] );....
0280: 09 66 6f 72 28 69 6e 74 20 78 3d 30 3b 20 78 2b  .for(int x=0; x+
0290: 28 31 3c 3c 6b 2d 31 29 3c 6e 3b 20 2b 2b 78 29  (1<<k-1)<n; ++x)
02a0: 0d 0a 09 09 09 09 69 66 28 20 64 5b 72 6d 5b 6b  ......if( d[rm[k
02b0: 5d 5b 78 5d 5d 20 3e 20 64 5b 72 6d 5b 6b 2d 31  ][x]] > d[rm[k-1
02c0: 5d 5b 78 20 2b 20 28 31 3c 3c 6b 2d 31 29 5d 5d  ][x + (1<<k-1)]]
02d0: 20 29 0d 0a 09 09 09 09 09 72 6d 5b 6b 5d 5b 78   ).......rm[k][x
02e0: 5d 20 3d 20 72 6d 5b 6b 2d 31 5d 5b 78 20 2b 20  ] = rm[k-1][x + 
02f0: 28 31 3c 3c 6b 2d 31 29 5d 3b 0d 0a 09 09 7d 0d  (1<<k-1)];....}.
0300: 0a 09 7d 0d 0a 0d 0a 09 69 6e 74 20 6f 70 65 72  ..}.....int oper
0310: 61 74 6f 72 28 29 28 69 6e 74 20 4c 2c 20 69 6e  ator()(int L, in
0320: 74 20 52 29 20 63 6f 6e 73 74 20 7b 0d 0a 09 09  t R) const {....
0330: 69 6e 74 20 6b 3d 30 3b 0d 0a 09 09 66 6f 72 28  int k=0;....for(
0340: 3b 20 4c 2b 28 31 3c 3c 6b 29 20 3c 20 52 2d 28  ; L+(1<<k) < R-(
0350: 31 3c 3c 6b 29 2b 31 3b 20 2b 2b 6b 29 20 7b 7d  1<<k)+1; ++k) {}
0360: 0d 0a 09 09 4c 4c 20 69 20 3d 20 72 6d 5b 6b 5d  ....LL i = rm[k]
0370: 5b 4c 5d 3b 0d 0a 09 09 4c 4c 20 6a 20 3d 20 72  [L];....LL j = r
0380: 6d 5b 6b 5d 5b 52 2d 28 31 3c 3c 6b 29 2b 31 5d  m[k][R-(1<<k)+1]
0390: 3b 0d 0a 09 09 72 65 74 75 72 6e 20 28 64 5b 69  ;....return (d[i
03a0: 5d 3c 3d 64 5b 6a 5d 20 3f 20 69 20 3a 20 6a 29  ]<=d[j] ? i : j)
03b0: 3b 0d 0a 09 7d 0d 0a 0d 0a 09 76 65 63 74 6f 72  ;...}.....vector
03c0: 3c 69 6e 74 3e 20 61 6c 6c 28 69 6e 74 20 4c 2c  <int> all(int L,
03d0: 20 69 6e 74 20 52 29 20 63 6f 6e 73 74 20 7b 0d   int R) const {.
03e0: 0a 09 09 76 65 63 74 6f 72 3c 69 6e 74 3e 20 61  ...vector<int> a
03f0: 6e 73 3b 0d 0a 09 09 69 6e 74 20 6d 69 6e 56 61  ns;....int minVa
0400: 6c 75 65 20 3d 20 64 5b 28 2a 74 68 69 73 29 28  lue = d[(*this)(
0410: 4c 2c 20 52 29 5d 3b 0d 0a 09 09 77 68 69 6c 65  L, R)];....while
0420: 28 20 4c 20 3c 3d 20 52 20 29 20 7b 0d 0a 09 09  ( L <= R ) {....
0430: 09 69 6e 74 20 43 20 3d 20 28 2a 74 68 69 73 29  .int C = (*this)
0440: 28 4c 2c 20 52 29 3b 0d 0a 09 09 09 69 66 28 20  (L, R);.....if( 
0450: 6d 69 6e 56 61 6c 75 65 20 3c 20 64 5b 43 5d 20  minValue < d[C] 
0460: 29 0d 0a 09 09 09 09 62 72 65 61 6b 3b 0d 0a 09  )......break;...
0470: 09 09 61 6e 73 2e 70 75 73 68 5f 62 61 63 6b 28  ..ans.push_back(
0480: 43 29 3b 0d 0a 09 09 09 4c 20 3d 20 43 2b 31 3b  C);.....L = C+1;
0490: 0d 0a 09 09 7d 0d 0a 09 09 72 65 74 75 72 6e 20  ....}....return 
04a0: 61 6e 73 3b 0d 0a 09 7d 0d 0a 7d 3b 0d 0a 0d 0a  ans;...}..};....
04b0: 69 6e 74 20 42 5b 32 30 34 38 5d 5b 32 30 34 38  int B[2048][2048
04c0: 5d 3b 0d 0a 69 6e 74 20 55 5b 32 30 34 38 5d 5b  ];..int U[2048][
04d0: 32 30 34 38 5d 3b 0d 0a 69 6e 74 20 44 5b 32 30  2048];..int D[20
04e0: 34 38 5d 5b 32 30 34 38 5d 3b 0d 0a 0d 0a 63 6c  48][2048];....cl
04f0: 61 73 73 20 50 65 72 66 65 63 74 52 65 63 74 61  ass PerfectRecta
0500: 6e 67 6c 65 73 0d 0a 7b 0d 0a 70 75 62 6c 69 63  ngles..{..public
0510: 3a 0d 0a 09 69 6e 74 20 6e 75 6d 62 65 72 4f 66  :...int numberOf
0520: 52 65 63 74 61 6e 67 6c 65 73 28 69 6e 74 20 4e  Rectangles(int N
0530: 2c 20 69 6e 74 20 4d 2c 20 69 6e 74 20 58 30 2c  , int M, int X0,
0540: 20 69 6e 74 20 41 2c 20 69 6e 74 20 5f 42 2c 20   int A, int _B, 
0550: 69 6e 74 20 59 30 2c 20 69 6e 74 20 43 2c 20 69  int Y0, int C, i
0560: 6e 74 20 5f 44 29 20 0d 0a 09 7b 0d 0a 09 09 2f  nt _D) ...{..../
0570: 2f 20 69 6e 70 75 74 0d 0a 09 09 69 6e 74 20 78  / input....int x
0580: 3d 58 30 25 4e 2c 20 79 3d 59 30 25 4e 3b 0d 0a  =X0%N, y=Y0%N;..
0590: 09 09 6d 65 6d 73 65 74 28 42 2c 20 30 2c 20 73  ..memset(B, 0, s
05a0: 69 7a 65 6f 66 28 42 29 29 3b 0d 0a 09 09 66 6f  izeof(B));....fo
05b0: 72 28 69 6e 74 20 69 3d 30 3b 20 69 3c 4d 3b 20  r(int i=0; i<M; 
05c0: 2b 2b 69 29 20 7b 0d 0a 09 09 09 42 5b 78 5d 5b  ++i) {.....B[x][
05d0: 79 5d 20 3d 20 31 3b 0d 0a 09 09 09 78 20 3d 20  y] = 1;.....x = 
05e0: 28 78 2a 41 2b 5f 42 29 25 4e 3b 0d 0a 09 09 09  (x*A+_B)%N;.....
05f0: 79 20 3d 20 28 79 2a 43 2b 5f 44 29 25 4e 3b 0d  y = (y*C+_D)%N;.
0600: 0a 09 09 7d 0d 0a 0d 0a 09 09 2f 2f 20 75 70 0d  ...}......// up.
0610: 0a 09 09 66 6f 72 28 69 6e 74 20 79 3d 30 3b 20  ...for(int y=0; 
0620: 79 3c 4e 3b 20 2b 2b 79 29 0d 0a 09 09 09 66 6f  y<N; ++y).....fo
0630: 72 28 69 6e 74 20 78 3d 30 3b 20 78 3c 4e 3b 20  r(int x=0; x<N; 
0640: 2b 2b 78 29 0d 0a 09 09 09 09 55 5b 79 5d 5b 78  ++x)......U[y][x
0650: 5d 20 3d 20 28 42 5b 79 5d 5b 78 5d 20 3f 20 30  ] = (B[y][x] ? 0
0660: 20 3a 20 79 3d 3d 30 20 3f 20 31 20 3a 20 55 5b   : y==0 ? 1 : U[
0670: 79 2d 31 5d 5b 78 5d 2b 31 29 3b 0d 0a 0d 0a 09  y-1][x]+1);.....
0680: 09 2f 2f 20 64 6f 77 6e 0d 0a 09 09 66 6f 72 28  .// down....for(
0690: 69 6e 74 20 79 3d 4e 2d 31 3b 20 79 3e 3d 30 3b  int y=N-1; y>=0;
06a0: 20 2d 2d 79 29 0d 0a 09 09 09 66 6f 72 28 69 6e   --y).....for(in
06b0: 74 20 78 3d 30 3b 20 78 3c 4e 3b 20 2b 2b 78 29  t x=0; x<N; ++x)
06c0: 0d 0a 09 09 09 09 44 5b 79 5d 5b 78 5d 20 3d 20  ......D[y][x] = 
06d0: 28 42 5b 79 5d 5b 78 5d 20 3f 20 30 20 3a 20 79  (B[y][x] ? 0 : y
06e0: 3d 3d 4e 2d 31 20 3f 20 31 20 3a 20 44 5b 79 2b  ==N-1 ? 1 : D[y+
06f0: 31 5d 5b 78 5d 2b 31 29 3b 0d 0a 0d 0a 09 09 2f  1][x]+1);....../
0700: 2f 20 73 6f 6c 76 65 0d 0a 09 09 69 6e 74 20 63  / solve....int c
0710: 6e 74 20 3d 20 30 3b 0d 0a 09 09 66 6f 72 28 69  nt = 0;....for(i
0720: 6e 74 20 79 3d 30 3b 20 79 3c 4e 3b 20 2b 2b 79  nt y=0; y<N; ++y
0730: 29 20 7b 0d 0a 09 09 09 52 4d 51 3c 69 6e 74 3e  ) {.....RMQ<int>
0740: 20 72 6d 55 28 55 5b 79 5d 2c 20 4e 29 3b 0d 0a   rmU(U[y], N);..
0750: 09 09 09 52 4d 51 3c 69 6e 74 3e 20 72 6d 44 28  ...RMQ<int> rmD(
0760: 44 5b 79 5d 2c 20 4e 29 3b 0d 0a 0d 0a 09 09 09  D[y], N);.......
0770: 73 74 61 63 6b 3c 20 70 61 69 72 3c 69 6e 74 2c  stack< pair<int,
0780: 69 6e 74 3e 20 3e 20 53 3b 0d 0a 09 09 09 53 2e  int> > S;.....S.
0790: 70 75 73 68 28 20 6d 61 6b 65 5f 70 61 69 72 28  push( make_pair(
07a0: 30 2c 20 4e 2d 31 29 20 29 3b 0d 0a 09 09 09 77  0, N-1) );.....w
07b0: 68 69 6c 65 28 20 21 53 2e 65 6d 70 74 79 28 29  hile( !S.empty()
07c0: 20 29 20 7b 0d 0a 09 09 09 09 69 6e 74 20 4c 20   ) {......int L 
07d0: 3d 20 53 2e 74 6f 70 28 29 2e 66 69 72 73 74 3b  = S.top().first;
07e0: 0d 0a 09 09 09 09 69 6e 74 20 52 20 3d 20 53 2e  ......int R = S.
07f0: 74 6f 70 28 29 2e 73 65 63 6f 6e 64 3b 0d 0a 09  top().second;...
0800: 09 09 09 53 2e 70 6f 70 28 29 3b 0d 0a 0d 0a 09  ...S.pop();.....
0810: 09 09 09 69 6e 74 20 64 20 3d 20 72 6d 44 28 4c  ...int d = rmD(L
0820: 2c 52 29 3b 0d 0a 09 09 09 09 69 66 28 20 44 5b  ,R);......if( D[
0830: 79 5d 5b 64 5d 20 3d 3d 20 31 20 29 0d 0a 09 09  y][d] == 1 )....
0840: 09 09 09 2b 2b 63 6e 74 3b 0d 0a 09 09 09 09 69  ...++cnt;......i
0850: 66 28 20 44 5b 79 5d 5b 64 5d 20 3c 3d 20 31 20  f( D[y][d] <= 1 
0860: 29 20 7b 0d 0a 09 09 09 09 09 76 65 63 74 6f 72  ) {.......vector
0870: 3c 69 6e 74 3e 20 6d 20 3d 20 72 6d 55 2e 61 6c  <int> m = rmU.al
0880: 6c 28 4c 2c 20 52 29 3b 0d 0a 09 09 09 09 09 6d  l(L, R);.......m
0890: 2e 70 75 73 68 5f 62 61 63 6b 28 52 2b 31 29 3b  .push_back(R+1);
08a0: 0d 0a 09 09 09 09 09 66 6f 72 28 69 6e 74 20 69  .......for(int i
08b0: 3d 30 3b 20 69 3c 6d 2e 73 69 7a 65 28 29 3b 20  =0; i<m.size(); 
08c0: 4c 3d 6d 5b 69 2b 2b 5d 2b 31 29 0d 0a 09 09 09  L=m[i++]+1).....
08d0: 09 09 09 69 66 28 20 4c 20 3c 3d 20 6d 5b 69 5d  ...if( L <= m[i]
08e0: 2d 31 20 29 0d 0a 09 09 09 09 09 09 09 53 2e 70  -1 ).........S.p
08f0: 75 73 68 28 20 6d 61 6b 65 5f 70 61 69 72 28 4c  ush( make_pair(L
0900: 2c 20 6d 5b 69 5d 2d 31 29 20 29 3b 0d 0a 09 09  , m[i]-1) );....
0910: 09 09 7d 0d 0a 09 09 09 7d 0d 0a 09 09 7d 0d 0a  ..}.....}....}..
0920: 09 09 72 65 74 75 72 6e 20 63 6e 74 3b 0d 0a 09  ..return cnt;...
0930: 7d 0d 0a 7d 3b 0d 0a 0d 0a 2f 2f 20 42 45 47 49  }..};....// BEGI
0940: 4e 20 43 55 54 20 48 45 52 45 0d 0a 23 69 6e 63  N CUT HERE..#inc
0950: 6c 75 64 65 20 3c 63 74 69 6d 65 3e 0d 0a 64 6f  lude <ctime>..do
0960: 75 62 6c 65 20 73 74 61 72 74 5f 74 69 6d 65 3b  uble start_time;
0970: 73 74 72 69 6e 67 20 74 69 6d 65 72 28 29 20 7b  string timer() {
0980: 20 6f 73 74 72 69 6e 67 73 74 72 65 61 6d 20 6f   ostringstream o
0990: 73 3b 20 6f 73 20 3c 3c 20 22 20 28 22 20 3c 3c  s; os << " (" <<
09a0: 20 69 6e 74 28 28 63 6c 6f 63 6b 28 29 2d 73 74   int((clock()-st
09b0: 61 72 74 5f 74 69 6d 65 29 2f 43 4c 4f 43 4b 53  art_time)/CLOCKS
09c0: 5f 50 45 52 5f 53 45 43 2a 31 30 30 30 29 20 3c  _PER_SEC*1000) <
09d0: 3c 20 22 20 6d 73 65 63 29 22 3b 20 72 65 74 75  < " msec)"; retu
09e0: 72 6e 20 6f 73 2e 73 74 72 28 29 3b 20 7d 0d 0a  rn os.str(); }..
09f0: 0d 0a 74 65 6d 70 6c 61 74 65 20 3c 74 79 70 65  ..template <type
0a00: 6e 61 6d 65 20 54 3e 20 73 74 72 69 6e 67 20 70  name T> string p
0a10: 72 69 6e 74 5f 61 72 72 61 79 28 63 6f 6e 73 74  rint_array(const
0a20: 20 76 65 63 74 6f 72 3c 54 3e 20 26 56 29 20 7b   vector<T> &V) {
0a30: 20 6f 73 74 72 69 6e 67 73 74 72 65 61 6d 20 6f   ostringstream o
0a40: 73 3b 20 6f 73 20 3c 3c 20 22 7b 20 22 3b 20 66  s; os << "{ "; f
0a50: 6f 72 20 28 74 79 70 65 6e 61 6d 65 20 76 65 63  or (typename vec
0a60: 74 6f 72 3c 54 3e 3a 3a 63 6f 6e 73 74 5f 69 74  tor<T>::const_it
0a70: 65 72 61 74 6f 72 20 69 74 65 72 20 3d 20 56 2e  erator iter = V.
0a80: 62 65 67 69 6e 28 29 3b 20 69 74 65 72 20 21 3d  begin(); iter !=
0a90: 20 56 2e 65 6e 64 28 29 3b 20 2b 2b 69 74 65 72   V.end(); ++iter
0aa0: 29 20 6f 73 20 3c 3c 20 27 5c 22 27 20 3c 3c 20  ) os << '\"' << 
0ab0: 2a 69 74 65 72 20 3c 3c 20 22 5c 22 2c 22 3b 20  *iter << "\","; 
0ac0: 6f 73 20 3c 3c 20 22 20 7d 22 3b 20 72 65 74 75  os << " }"; retu
0ad0: 72 6e 20 6f 73 2e 73 74 72 28 29 3b 20 7d 0d 0a  rn os.str(); }..
0ae0: 69 6e 74 20 76 65 72 69 66 79 5f 63 61 73 65 28  int verify_case(
0af0: 63 6f 6e 73 74 20 69 6e 74 20 26 45 78 70 65 63  const int &Expec
0b00: 74 65 64 2c 20 63 6f 6e 73 74 20 69 6e 74 20 26  ted, const int &
0b10: 52 65 63 65 69 76 65 64 29 20 7b 20 69 66 20 28  Received) { if (
0b20: 45 78 70 65 63 74 65 64 20 3d 3d 20 52 65 63 65  Expected == Rece
0b30: 69 76 65 64 29 20 63 65 72 72 20 3c 3c 20 22 50  ived) cerr << "P
0b40: 41 53 53 45 44 22 20 3c 3c 20 74 69 6d 65 72 28  ASSED" << timer(
0b50: 29 20 3c 3c 20 65 6e 64 6c 3b 20 65 6c 73 65 20  ) << endl; else 
0b60: 7b 20 63 65 72 72 20 3c 3c 20 22 46 41 49 4c 45  { cerr << "FAILE
0b70: 44 22 20 3c 3c 20 74 69 6d 65 72 28 29 20 3c 3c  D" << timer() <<
0b80: 20 65 6e 64 6c 3b 20 63 65 72 72 20 3c 3c 20 22   endl; cerr << "
0b90: 5c 74 45 78 70 65 63 74 65 64 3a 20 5c 22 22 20  \tExpected: \"" 
0ba0: 3c 3c 20 45 78 70 65 63 74 65 64 20 3c 3c 20 27  << Expected << '
0bb0: 5c 22 27 20 3c 3c 20 65 6e 64 6c 3b 20 63 65 72  \"' << endl; cer
0bc0: 72 20 3c 3c 20 22 5c 74 52 65 63 65 69 76 65 64  r << "\tReceived
0bd0: 3a 20 5c 22 22 20 3c 3c 20 52 65 63 65 69 76 65  : \"" << Receive
0be0: 64 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 65 6e 64  d << '\"' << end
0bf0: 6c 3b 20 7d 20 72 65 74 75 72 6e 20 30 3b 7d 0d  l; } return 0;}.
0c00: 0a 0d 0a 74 65 6d 70 6c 61 74 65 3c 69 6e 74 20  ...template<int 
0c10: 4e 3e 20 73 74 72 75 63 74 20 43 61 73 65 5f 20  N> struct Case_ 
0c20: 7b 20 43 61 73 65 5f 28 29 7b 73 74 61 72 74 5f  { Case_(){start_
0c30: 74 69 6d 65 3d 63 6c 6f 63 6b 28 29 3b 7d 20 7d  time=clock();} }
0c40: 3b 0d 0a 63 68 61 72 20 54 65 73 74 5f 28 2e 2e  ;..char Test_(..
0c50: 2e 29 3b 0d 0a 69 6e 74 20 54 65 73 74 5f 28 43  .);..int Test_(C
0c60: 61 73 65 5f 3c 30 3e 29 20 7b 0d 0a 09 69 6e 74  ase_<0>) {...int
0c70: 20 4e 20 3d 20 35 3b 20 0d 0a 09 69 6e 74 20 4d   N = 5; ...int M
0c80: 20 3d 20 31 3b 20 0d 0a 09 69 6e 74 20 58 30 20   = 1; ...int X0 
0c90: 3d 20 32 3b 20 0d 0a 09 69 6e 74 20 41 20 3d 20  = 2; ...int A = 
0ca0: 30 3b 20 0d 0a 09 69 6e 74 20 42 20 3d 20 30 3b  0; ...int B = 0;
0cb0: 20 0d 0a 09 69 6e 74 20 59 30 20 3d 20 32 3b 20   ...int Y0 = 2; 
0cc0: 0d 0a 09 69 6e 74 20 43 20 3d 20 30 3b 20 0d 0a  ...int C = 0; ..
0cd0: 09 69 6e 74 20 44 20 3d 20 30 3b 20 0d 0a 09 69  .int D = 0; ...i
0ce0: 6e 74 20 52 65 74 56 61 6c 20 3d 20 34 3b 20 0d  nt RetVal = 4; .
0cf0: 0a 09 72 65 74 75 72 6e 20 76 65 72 69 66 79 5f  ..return verify_
0d00: 63 61 73 65 28 52 65 74 56 61 6c 2c 20 50 65 72  case(RetVal, Per
0d10: 66 65 63 74 52 65 63 74 61 6e 67 6c 65 73 28 29  fectRectangles()
0d20: 2e 6e 75 6d 62 65 72 4f 66 52 65 63 74 61 6e 67  .numberOfRectang
0d30: 6c 65 73 28 4e 2c 20 4d 2c 20 58 30 2c 20 41 2c  les(N, M, X0, A,
0d40: 20 42 2c 20 59 30 2c 20 43 2c 20 44 29 29 3b 20   B, Y0, C, D)); 
0d50: 7d 0d 0a 69 6e 74 20 54 65 73 74 5f 28 43 61 73  }..int Test_(Cas
0d60: 65 5f 3c 31 3e 29 20 7b 0d 0a 09 69 6e 74 20 4e  e_<1>) {...int N
0d70: 20 3d 20 34 3b 20 0d 0a 09 69 6e 74 20 4d 20 3d   = 4; ...int M =
0d80: 20 34 3b 20 0d 0a 09 69 6e 74 20 58 30 20 3d 20   4; ...int X0 = 
0d90: 30 3b 20 0d 0a 09 69 6e 74 20 41 20 3d 20 31 3b  0; ...int A = 1;
0da0: 20 0d 0a 09 69 6e 74 20 42 20 3d 20 31 3b 20 0d   ...int B = 1; .
0db0: 0a 09 69 6e 74 20 59 30 20 3d 20 30 3b 20 0d 0a  ..int Y0 = 0; ..
0dc0: 09 69 6e 74 20 43 20 3d 20 31 3b 20 0d 0a 09 69  .int C = 1; ...i
0dd0: 6e 74 20 44 20 3d 20 31 3b 20 0d 0a 09 69 6e 74  nt D = 1; ...int
0de0: 20 52 65 74 56 61 6c 20 3d 20 36 3b 20 0d 0a 09   RetVal = 6; ...
0df0: 72 65 74 75 72 6e 20 76 65 72 69 66 79 5f 63 61  return verify_ca
0e00: 73 65 28 52 65 74 56 61 6c 2c 20 50 65 72 66 65  se(RetVal, Perfe
0e10: 63 74 52 65 63 74 61 6e 67 6c 65 73 28 29 2e 6e  ctRectangles().n
0e20: 75 6d 62 65 72 4f 66 52 65 63 74 61 6e 67 6c 65  umberOfRectangle
0e30: 73 28 4e 2c 20 4d 2c 20 58 30 2c 20 41 2c 20 42  s(N, M, X0, A, B
0e40: 2c 20 59 30 2c 20 43 2c 20 44 29 29 3b 20 7d 0d  , Y0, C, D)); }.
0e50: 0a 69 6e 74 20 54 65 73 74 5f 28 43 61 73 65 5f  .int Test_(Case_
0e60: 3c 32 3e 29 20 7b 0d 0a 09 69 6e 74 20 4e 20 3d  <2>) {...int N =
0e70: 20 31 3b 20 0d 0a 09 69 6e 74 20 4d 20 3d 20 31   1; ...int M = 1
0e80: 30 30 30 30 30 30 3b 20 0d 0a 09 69 6e 74 20 58  000000; ...int X
0e90: 30 20 3d 20 31 3b 20 0d 0a 09 69 6e 74 20 41 20  0 = 1; ...int A 
0ea0: 3d 20 32 3b 20 0d 0a 09 69 6e 74 20 42 20 3d 20  = 2; ...int B = 
0eb0: 33 3b 20 0d 0a 09 69 6e 74 20 59 30 20 3d 20 31  3; ...int Y0 = 1
0ec0: 3b 20 0d 0a 09 69 6e 74 20 43 20 3d 20 34 3b 20  ; ...int C = 4; 
0ed0: 0d 0a 09 69 6e 74 20 44 20 3d 20 37 3b 20 0d 0a  ...int D = 7; ..
0ee0: 09 69 6e 74 20 52 65 74 56 61 6c 20 3d 20 30 3b  .int RetVal = 0;
0ef0: 20 0d 0a 09 72 65 74 75 72 6e 20 76 65 72 69 66   ...return verif
0f00: 79 5f 63 61 73 65 28 52 65 74 56 61 6c 2c 20 50  y_case(RetVal, P
0f10: 65 72 66 65 63 74 52 65 63 74 61 6e 67 6c 65 73  erfectRectangles
0f20: 28 29 2e 6e 75 6d 62 65 72 4f 66 52 65 63 74 61  ().numberOfRecta
0f30: 6e 67 6c 65 73 28 4e 2c 20 4d 2c 20 58 30 2c 20  ngles(N, M, X0, 
0f40: 41 2c 20 42 2c 20 59 30 2c 20 43 2c 20 44 29 29  A, B, Y0, C, D))
0f50: 3b 20 7d 0d 0a 69 6e 74 20 54 65 73 74 5f 28 43  ; }..int Test_(C
0f60: 61 73 65 5f 3c 33 3e 29 20 7b 0d 0a 09 69 6e 74  ase_<3>) {...int
0f70: 20 4e 20 3d 20 31 30 3b 20 0d 0a 09 69 6e 74 20   N = 10; ...int 
0f80: 4d 20 3d 20 32 30 3b 20 0d 0a 09 69 6e 74 20 58  M = 20; ...int X
0f90: 30 20 3d 20 34 3b 20 0d 0a 09 69 6e 74 20 41 20  0 = 4; ...int A 
0fa0: 3d 20 37 36 3b 20 0d 0a 09 69 6e 74 20 42 20 3d  = 76; ...int B =
0fb0: 20 32 3b 20 0d 0a 09 69 6e 74 20 59 30 20 3d 20   2; ...int Y0 = 
0fc0: 36 3b 20 0d 0a 09 69 6e 74 20 43 20 3d 20 32 3b  6; ...int C = 2;
0fd0: 20 0d 0a 09 69 6e 74 20 44 20 3d 20 34 33 3b 20   ...int D = 43; 
0fe0: 0d 0a 09 69 6e 74 20 52 65 74 56 61 6c 20 3d 20  ...int RetVal = 
0ff0: 31 32 3b 20 0d 0a 09 72 65 74 75 72 6e 20 76 65  12; ...return ve
1000: 72 69 66 79 5f 63 61 73 65 28 52 65 74 56 61 6c  rify_case(RetVal
1010: 2c 20 50 65 72 66 65 63 74 52 65 63 74 61 6e 67  , PerfectRectang
1020: 6c 65 73 28 29 2e 6e 75 6d 62 65 72 4f 66 52 65  les().numberOfRe
1030: 63 74 61 6e 67 6c 65 73 28 4e 2c 20 4d 2c 20 58  ctangles(N, M, X
1040: 30 2c 20 41 2c 20 42 2c 20 59 30 2c 20 43 2c 20  0, A, B, Y0, C, 
1050: 44 29 29 3b 20 7d 0d 0a 69 6e 74 20 54 65 73 74  D)); }..int Test
1060: 5f 28 43 61 73 65 5f 3c 34 3e 29 20 7b 0d 0a 09  _(Case_<4>) {...
1070: 69 6e 74 20 4e 20 3d 20 31 39 39 39 3b 20 0d 0a  int N = 1999; ..
1080: 09 69 6e 74 20 4d 20 3d 20 34 30 30 30 30 30 30  .int M = 4000000
1090: 3b 20 0d 0a 09 69 6e 74 20 58 30 20 3d 20 31 30  ; ...int X0 = 10
10a0: 30 30 3b 20 0d 0a 09 69 6e 74 20 41 20 3d 20 31  00; ...int A = 1
10b0: 3b 20 0d 0a 09 69 6e 74 20 42 20 3d 20 31 3b 20  ; ...int B = 1; 
10c0: 0d 0a 09 69 6e 74 20 59 30 20 3d 20 30 3b 20 0d  ...int Y0 = 0; .
10d0: 0a 09 69 6e 74 20 43 20 3d 20 31 3b 20 0d 0a 09  ..int C = 1; ...
10e0: 69 6e 74 20 44 20 3d 20 31 3b 20 0d 0a 09 69 6e  int D = 1; ...in
10f0: 74 20 52 65 74 56 61 6c 20 3d 20 31 30 30 32 39  t RetVal = 10029
1100: 39 35 3b 20 0d 0a 09 72 65 74 75 72 6e 20 76 65  95; ...return ve
1110: 72 69 66 79 5f 63 61 73 65 28 52 65 74 56 61 6c  rify_case(RetVal
1120: 2c 20 50 65 72 66 65 63 74 52 65 63 74 61 6e 67  , PerfectRectang
1130: 6c 65 73 28 29 2e 6e 75 6d 62 65 72 4f 66 52 65  les().numberOfRe
1140: 63 74 61 6e 67 6c 65 73 28 4e 2c 20 4d 2c 20 58  ctangles(N, M, X
1150: 30 2c 20 41 2c 20 42 2c 20 59 30 2c 20 43 2c 20  0, A, B, Y0, C, 
1160: 44 29 29 3b 20 7d 0d 0a 69 6e 74 20 54 65 73 74  D)); }..int Test
1170: 5f 28 43 61 73 65 5f 3c 35 3e 29 20 7b 0d 0a 09  _(Case_<5>) {...
1180: 69 6e 74 20 4e 20 3d 20 32 30 30 30 3b 20 0d 0a  int N = 2000; ..
1190: 09 69 6e 74 20 4d 20 3d 20 34 30 30 30 30 30 30  .int M = 4000000
11a0: 3b 20 0d 0a 09 69 6e 74 20 58 30 20 3d 20 31 35  ; ...int X0 = 15
11b0: 30 30 3b 20 0d 0a 09 69 6e 74 20 41 20 3d 20 31  00; ...int A = 1
11c0: 3b 20 0d 0a 09 69 6e 74 20 42 20 3d 20 31 3b 20  ; ...int B = 1; 
11d0: 0d 0a 09 69 6e 74 20 59 30 20 3d 20 35 30 30 3b  ...int Y0 = 500;
11e0: 20 0d 0a 09 69 6e 74 20 43 20 3d 20 31 3b 20 0d   ...int C = 1; .
11f0: 0a 09 69 6e 74 20 44 20 3d 20 31 3b 20 0d 0a 09  ..int D = 1; ...
1200: 69 6e 74 20 52 65 74 56 61 6c 20 3d 20 31 30 30  int RetVal = 100
1210: 33 39 39 37 3b 20 0d 0a 09 72 65 74 75 72 6e 20  3997; ...return 
1220: 76 65 72 69 66 79 5f 63 61 73 65 28 52 65 74 56  verify_case(RetV
1230: 61 6c 2c 20 50 65 72 66 65 63 74 52 65 63 74 61  al, PerfectRecta
1240: 6e 67 6c 65 73 28 29 2e 6e 75 6d 62 65 72 4f 66  ngles().numberOf
1250: 52 65 63 74 61 6e 67 6c 65 73 28 4e 2c 20 4d 2c  Rectangles(N, M,
1260: 20 58 30 2c 20 41 2c 20 42 2c 20 59 30 2c 20 43   X0, A, B, Y0, C
1270: 2c 20 44 29 29 3b 20 7d 0d 0a 74 65 6d 70 6c 61  , D)); }..templa
1280: 74 65 3c 69 6e 74 20 4e 3e 20 76 6f 69 64 20 52  te<int N> void R
1290: 75 6e 5f 28 29 20 7b 20 63 65 72 72 20 3c 3c 20  un_() { cerr << 
12a0: 22 54 65 73 74 20 43 61 73 65 20 23 22 20 3c 3c  "Test Case #" <<
12b0: 20 4e 20 3c 3c 20 22 2e 2e 2e 22 20 3c 3c 20 66   N << "..." << f
12c0: 6c 75 73 68 3b 20 54 65 73 74 5f 28 43 61 73 65  lush; Test_(Case
12d0: 5f 3c 4e 3e 28 29 29 3b 20 52 75 6e 5f 3c 73 69  _<N>()); Run_<si
12e0: 7a 65 6f 66 28 54 65 73 74 5f 28 43 61 73 65 5f  zeof(Test_(Case_
12f0: 3c 4e 2b 31 3e 28 29 29 29 3d 3d 31 20 3f 20 2d  <N+1>()))==1 ? -
1300: 31 20 3a 20 4e 2b 31 3e 28 29 3b 20 7d 0d 0a 74  1 : N+1>(); }..t
1310: 65 6d 70 6c 61 74 65 3c 3e 20 76 6f 69 64 20 52  emplate<> void R
1320: 75 6e 5f 3c 2d 31 3e 28 29 20 7b 7d 0d 0a 69 6e  un_<-1>() {}..in
1330: 74 20 6d 61 69 6e 28 29 20 7b 20 52 75 6e 5f 3c  t main() { Run_<
1340: 30 3e 28 29 3b 20 7d 0d 0a 2f 2f 20 45 4e 44 20  0>(); }..// END 
1350: 43 55 54 20 48 45 52 45 0d 0a 0d 0a              CUT HERE....