Hex Artifact Content
Not logged in

Artifact 527ef29abc31b65231b9849c635227471eb24494:


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 73 74 61 74 69 63 20 63 6f 6e  P;....static con
01a0: 73 74 20 75 6e 73 69 67 6e 65 64 20 4d 4f 44 56  st unsigned MODV
01b0: 41 4c 20 3d 20 31 30 30 30 30 30 30 30 30 37 3b  AL = 1000000007;
01c0: 0d 0a 73 74 72 75 63 74 20 6d 69 6e 74 0d 0a 7b  ..struct mint..{
01d0: 0d 0a 09 75 6e 73 69 67 6e 65 64 20 76 61 6c 3b  ...unsigned val;
01e0: 0d 0a 09 6d 69 6e 74 28 29 3a 76 61 6c 28 30 29  ...mint():val(0)
01f0: 7b 7d 0d 0a 09 6d 69 6e 74 28 69 6e 74 20 20 20  {}...mint(int   
0200: 20 20 20 78 29 3a 76 61 6c 28 78 25 4d 4f 44 56     x):val(x%MODV
0210: 41 4c 29 20 7b 7d 0d 0a 09 6d 69 6e 74 28 75 6e  AL) {}...mint(un
0220: 73 69 67 6e 65 64 20 78 29 3a 76 61 6c 28 78 25  signed x):val(x%
0230: 4d 4f 44 56 41 4c 29 20 7b 7d 0d 0a 09 6d 69 6e  MODVAL) {}...min
0240: 74 28 4c 4c 20 20 20 20 20 20 20 78 29 3a 76 61  t(LL       x):va
0250: 6c 28 78 25 4d 4f 44 56 41 4c 29 20 7b 7d 0d 0a  l(x%MODVAL) {}..
0260: 7d 3b 0d 0a 6d 69 6e 74 26 20 6f 70 65 72 61 74  };..mint& operat
0270: 6f 72 2b 3d 28 6d 69 6e 74 26 20 78 2c 20 6d 69  or+=(mint& x, mi
0280: 6e 74 20 79 29 20 7b 20 72 65 74 75 72 6e 20 78  nt y) { return x
0290: 20 3d 20 78 2e 76 61 6c 2b 79 2e 76 61 6c 3b 20   = x.val+y.val; 
02a0: 7d 0d 0a 6d 69 6e 74 26 20 6f 70 65 72 61 74 6f  }..mint& operato
02b0: 72 2d 3d 28 6d 69 6e 74 26 20 78 2c 20 6d 69 6e  r-=(mint& x, min
02c0: 74 20 79 29 20 7b 20 72 65 74 75 72 6e 20 78 20  t y) { return x 
02d0: 3d 20 78 2e 76 61 6c 2d 79 2e 76 61 6c 2b 4d 4f  = x.val-y.val+MO
02e0: 44 56 41 4c 3b 20 7d 0d 0a 6d 69 6e 74 26 20 6f  DVAL; }..mint& o
02f0: 70 65 72 61 74 6f 72 2a 3d 28 6d 69 6e 74 26 20  perator*=(mint& 
0300: 78 2c 20 6d 69 6e 74 20 79 29 20 7b 20 72 65 74  x, mint y) { ret
0310: 75 72 6e 20 78 20 3d 20 4c 4c 28 78 2e 76 61 6c  urn x = LL(x.val
0320: 29 2a 79 2e 76 61 6c 3b 20 7d 0d 0a 6d 69 6e 74  )*y.val; }..mint
0330: 20 6f 70 65 72 61 74 6f 72 2b 28 6d 69 6e 74 20   operator+(mint 
0340: 78 2c 20 6d 69 6e 74 20 79 29 20 7b 20 72 65 74  x, mint y) { ret
0350: 75 72 6e 20 78 2b 3d 79 3b 20 7d 0d 0a 6d 69 6e  urn x+=y; }..min
0360: 74 20 6f 70 65 72 61 74 6f 72 2d 28 6d 69 6e 74  t operator-(mint
0370: 20 78 2c 20 6d 69 6e 74 20 79 29 20 7b 20 72 65   x, mint y) { re
0380: 74 75 72 6e 20 78 2d 3d 79 3b 20 7d 0d 0a 6d 69  turn x-=y; }..mi
0390: 6e 74 20 6f 70 65 72 61 74 6f 72 2a 28 6d 69 6e  nt operator*(min
03a0: 74 20 78 2c 20 6d 69 6e 74 20 79 29 20 7b 20 72  t x, mint y) { r
03b0: 65 74 75 72 6e 20 78 2a 3d 79 3b 20 7d 0d 0a 0d  eturn x*=y; }...
03c0: 0a 76 65 63 74 6f 72 3c 20 76 65 63 74 6f 72 3c  .vector< vector<
03d0: 6d 69 6e 74 3e 20 3e 20 43 50 5f 3b 0d 0a 6d 69  mint> > CP_;..mi
03e0: 6e 74 20 43 28 69 6e 74 20 6e 2c 20 69 6e 74 20  nt C(int n, int 
03f0: 6b 29 20 7b 0d 0a 09 77 68 69 6c 65 28 20 43 50  k) {...while( CP
0400: 5f 2e 73 69 7a 65 28 29 20 3c 3d 20 6e 20 29 20  _.size() <= n ) 
0410: 7b 0d 0a 09 09 69 6e 74 20 6e 6e 20 3d 20 43 50  {....int nn = CP
0420: 5f 2e 73 69 7a 65 28 29 3b 0d 0a 09 09 43 50 5f  _.size();....CP_
0430: 2e 70 75 73 68 5f 62 61 63 6b 28 76 65 63 74 6f  .push_back(vecto
0440: 72 3c 6d 69 6e 74 3e 28 6e 6e 2b 31 2c 31 29 29  r<mint>(nn+1,1))
0450: 3b 0d 0a 09 09 66 6f 72 28 69 6e 74 20 6b 6b 3d  ;....for(int kk=
0460: 31 3b 20 6b 6b 3c 6e 6e 3b 20 2b 2b 6b 6b 29 0d  1; kk<nn; ++kk).
0470: 0a 09 09 09 43 50 5f 5b 6e 6e 5d 5b 6b 6b 5d 20  ....CP_[nn][kk] 
0480: 3d 20 43 50 5f 5b 6e 6e 2d 31 5d 5b 6b 6b 2d 31  = CP_[nn-1][kk-1
0490: 5d 20 2b 20 43 50 5f 5b 6e 6e 2d 31 5d 5b 6b 6b  ] + CP_[nn-1][kk
04a0: 5d 3b 0d 0a 09 7d 0d 0a 09 72 65 74 75 72 6e 20  ];...}...return 
04b0: 6b 3c 30 20 7c 7c 20 6e 3c 6b 20 3f 20 30 20 3a  k<0 || n<k ? 0 :
04c0: 20 43 50 5f 5b 6e 5d 5b 6b 5d 3b 0d 0a 7d 0d 0a   CP_[n][k];..}..
04d0: 0d 0a 63 6c 61 73 73 20 53 75 6d 4f 76 65 72 50  ..class SumOverP
04e0: 65 72 6d 75 74 61 74 69 6f 6e 73 20 7b 20 70 75  ermutations { pu
04f0: 62 6c 69 63 3a 0d 0a 09 69 6e 74 20 4e 3b 0d 0a  blic:...int N;..
0500: 09 76 65 63 74 6f 72 3c 69 6e 74 3e 20 6d 65 6d  .vector<int> mem
0510: 6f 4f 4f 2c 20 6d 65 6d 6f 43 4f 2c 20 6d 65 6d  oOO, memoCO, mem
0520: 6f 43 43 3b 0d 0a 0d 0a 09 69 6e 74 20 66 69 6e  oCC;.....int fin
0530: 64 53 75 6d 28 69 6e 74 20 6e 29 0d 0a 09 7b 0d  dSum(int n)...{.
0540: 0a 09 09 4e 20 3d 20 6e 3b 0d 0a 09 09 6d 65 6d  ...N = n;....mem
0550: 6f 4f 4f 2e 61 73 73 69 67 6e 28 4e 2b 31 2c 20  oOO.assign(N+1, 
0560: 2d 31 29 3b 0d 0a 09 09 6d 65 6d 6f 43 4f 2e 61  -1);....memoCO.a
0570: 73 73 69 67 6e 28 4e 2b 31 2c 20 2d 31 29 3b 0d  ssign(N+1, -1);.
0580: 0a 09 09 6d 65 6d 6f 43 43 2e 61 73 73 69 67 6e  ...memoCC.assign
0590: 28 4e 2b 31 2c 20 2d 31 29 3b 0d 0a 09 09 72 65  (N+1, -1);....re
05a0: 74 75 72 6e 20 72 65 63 4f 4f 28 6e 29 3b 0d 0a  turn recOO(n);..
05b0: 09 7d 0d 0a 0d 0a 09 69 6e 74 20 72 65 63 4f 4f  .}.....int recOO
05c0: 28 69 6e 74 20 6b 29 20 7b 0d 0a 09 09 69 66 28  (int k) {....if(
05d0: 6b 20 3d 3d 20 30 29 20 72 65 74 75 72 6e 20 31  k == 0) return 1
05e0: 3b 0d 0a 09 09 69 66 28 6b 20 3d 3d 20 31 29 20  ;....if(k == 1) 
05f0: 72 65 74 75 72 6e 20 4e 3b 0d 0a 09 09 69 66 28  return N;....if(
0600: 6d 65 6d 6f 4f 4f 5b 6b 5d 20 21 3d 20 2d 31 29  memoOO[k] != -1)
0610: 20 72 65 74 75 72 6e 20 6d 65 6d 6f 4f 4f 5b 6b   return memoOO[k
0620: 5d 3b 0d 0a 0d 0a 09 09 6d 69 6e 74 20 74 6f 74  ];......mint tot
0630: 61 6c 20 3d 20 30 3b 0d 0a 09 09 74 6f 74 61 6c  al = 0;....total
0640: 20 2b 3d 20 6d 69 6e 74 28 4e 29 20 2a 20 72 65   += mint(N) * re
0650: 63 43 4f 28 6b 2d 31 29 3b 20 2f 2f 20 6c 65 66  cCO(k-1); // lef
0660: 74 0d 0a 09 09 74 6f 74 61 6c 20 2b 3d 20 6d 69  t....total += mi
0670: 6e 74 28 4e 29 20 2a 20 72 65 63 43 4f 28 6b 2d  nt(N) * recCO(k-
0680: 31 29 3b 20 2f 2f 20 72 69 67 68 74 0d 0a 09 09  1); // right....
0690: 66 6f 72 28 69 6e 74 20 69 3d 31 3b 20 69 2b 31  for(int i=1; i+1
06a0: 3c 6b 3b 20 2b 2b 69 29 20 2f 2f 20 6d 69 64 0d  <k; ++i) // mid.
06b0: 0a 09 09 09 74 6f 74 61 6c 20 2b 3d 20 6d 69 6e  ....total += min
06c0: 74 28 4e 29 20 2a 20 72 65 63 43 4f 28 69 29 20  t(N) * recCO(i) 
06d0: 2a 20 72 65 63 43 4f 28 6b 2d 31 2d 69 29 20 2a  * recCO(k-1-i) *
06e0: 20 43 28 6b 2d 31 2c 20 69 29 3b 0d 0a 09 09 72   C(k-1, i);....r
06f0: 65 74 75 72 6e 20 6d 65 6d 6f 4f 4f 5b 6b 5d 20  eturn memoOO[k] 
0700: 3d 20 74 6f 74 61 6c 2e 76 61 6c 3b 0d 0a 09 7d  = total.val;...}
0710: 0d 0a 0d 0a 09 69 6e 74 20 72 65 63 43 4f 28 69  .....int recCO(i
0720: 6e 74 20 6b 29 20 7b 0d 0a 09 09 69 66 28 6b 20  nt k) {....if(k 
0730: 3d 3d 20 30 29 20 72 65 74 75 72 6e 20 31 3b 0d  == 0) return 1;.
0740: 0a 09 09 69 66 28 6b 20 3d 3d 20 31 29 20 72 65  ...if(k == 1) re
0750: 74 75 72 6e 20 4e 2d 31 3b 0d 0a 09 09 69 66 28  turn N-1;....if(
0760: 6d 65 6d 6f 43 4f 5b 6b 5d 20 21 3d 20 2d 31 29  memoCO[k] != -1)
0770: 20 72 65 74 75 72 6e 20 6d 65 6d 6f 43 4f 5b 6b   return memoCO[k
0780: 5d 3b 0d 0a 0d 0a 09 09 6d 69 6e 74 20 74 6f 74  ];......mint tot
0790: 61 6c 20 3d 20 30 3b 0d 0a 09 09 74 6f 74 61 6c  al = 0;....total
07a0: 20 2b 3d 20 6d 69 6e 74 28 4e 2d 31 29 20 2a 20   += mint(N-1) * 
07b0: 72 65 63 43 4f 28 6b 2d 31 29 3b 20 2f 2f 20 6c  recCO(k-1); // l
07c0: 65 66 74 0d 0a 09 09 74 6f 74 61 6c 20 2b 3d 20  eft....total += 
07d0: 6d 69 6e 74 28 4e 29 20 2a 20 72 65 63 43 43 28  mint(N) * recCC(
07e0: 6b 2d 31 29 3b 20 2f 2f 20 72 69 67 68 74 0d 0a  k-1); // right..
07f0: 09 09 66 6f 72 28 69 6e 74 20 69 3d 31 3b 20 69  ..for(int i=1; i
0800: 2b 31 3c 6b 3b 20 2b 2b 69 29 20 2f 2f 20 6d 69  +1<k; ++i) // mi
0810: 64 0d 0a 09 09 09 74 6f 74 61 6c 20 2b 3d 20 6d  d.....total += m
0820: 69 6e 74 28 4e 29 20 2a 20 72 65 63 43 43 28 69  int(N) * recCC(i
0830: 29 20 2a 20 72 65 63 43 4f 28 6b 2d 31 2d 69 29  ) * recCO(k-1-i)
0840: 20 2a 20 43 28 6b 2d 31 2c 20 69 29 3b 0d 0a 09   * C(k-1, i);...
0850: 09 72 65 74 75 72 6e 20 6d 65 6d 6f 43 4f 5b 6b  .return memoCO[k
0860: 5d 20 3d 20 74 6f 74 61 6c 2e 76 61 6c 3b 0d 0a  ] = total.val;..
0870: 09 7d 0d 0a 0d 0a 09 6d 69 6e 74 20 72 65 63 43  .}.....mint recC
0880: 43 28 69 6e 74 20 6b 29 20 7b 0d 0a 09 09 69 66  C(int k) {....if
0890: 28 6b 20 3d 3d 20 30 29 20 72 65 74 75 72 6e 20  (k == 0) return 
08a0: 31 3b 0d 0a 09 09 69 66 28 6b 20 3d 3d 20 31 29  1;....if(k == 1)
08b0: 20 72 65 74 75 72 6e 20 4e 2d 32 3b 0d 0a 09 09   return N-2;....
08c0: 69 66 28 6d 65 6d 6f 43 43 5b 6b 5d 20 21 3d 20  if(memoCC[k] != 
08d0: 2d 31 29 20 72 65 74 75 72 6e 20 6d 65 6d 6f 43  -1) return memoC
08e0: 43 5b 6b 5d 3b 0d 0a 0d 0a 09 09 6d 69 6e 74 20  C[k];......mint 
08f0: 74 6f 74 61 6c 20 3d 20 30 3b 0d 0a 09 09 74 6f  total = 0;....to
0900: 74 61 6c 20 2b 3d 20 6d 69 6e 74 28 4e 2d 31 29  tal += mint(N-1)
0910: 20 2a 20 72 65 63 43 43 28 6b 2d 31 29 3b 20 2f   * recCC(k-1); /
0920: 2f 20 6c 65 66 74 0d 0a 09 09 74 6f 74 61 6c 20  / left....total 
0930: 2b 3d 20 6d 69 6e 74 28 4e 2d 31 29 20 2a 20 72  += mint(N-1) * r
0940: 65 63 43 43 28 6b 2d 31 29 3b 20 2f 2f 20 72 69  ecCC(k-1); // ri
0950: 67 68 74 0d 0a 09 09 66 6f 72 28 69 6e 74 20 69  ght....for(int i
0960: 3d 31 3b 20 69 2b 31 3c 6b 3b 20 2b 2b 69 29 20  =1; i+1<k; ++i) 
0970: 2f 2f 20 6d 69 64 0d 0a 09 09 09 74 6f 74 61 6c  // mid.....total
0980: 20 2b 3d 20 6d 69 6e 74 28 4e 29 20 2a 20 72 65   += mint(N) * re
0990: 63 43 43 28 69 29 20 2a 20 72 65 63 43 43 28 6b  cCC(i) * recCC(k
09a0: 2d 31 2d 69 29 20 2a 20 43 28 6b 2d 31 2c 20 69  -1-i) * C(k-1, i
09b0: 29 3b 0d 0a 09 09 72 65 74 75 72 6e 20 6d 65 6d  );....return mem
09c0: 6f 43 43 5b 6b 5d 20 3d 20 74 6f 74 61 6c 2e 76  oCC[k] = total.v
09d0: 61 6c 3b 0d 0a 09 7d 0d 0a 7d 3b 0d 0a 0d 0a 2f  al;...}..};..../
09e0: 2f 20 42 45 47 49 4e 20 43 55 54 20 48 45 52 45  / BEGIN CUT HERE
09f0: 0d 0a 23 69 6e 63 6c 75 64 65 20 3c 63 74 69 6d  ..#include <ctim
0a00: 65 3e 0d 0a 64 6f 75 62 6c 65 20 73 74 61 72 74  e>..double start
0a10: 5f 74 69 6d 65 3b 20 73 74 72 69 6e 67 20 74 69  _time; string ti
0a20: 6d 65 72 28 29 0d 0a 20 7b 20 6f 73 74 72 69 6e  mer().. { ostrin
0a30: 67 73 74 72 65 61 6d 20 6f 73 3b 20 6f 73 20 3c  gstream os; os <
0a40: 3c 20 22 20 28 22 20 3c 3c 20 69 6e 74 28 28 63  < " (" << int((c
0a50: 6c 6f 63 6b 28 29 2d 73 74 61 72 74 5f 74 69 6d  lock()-start_tim
0a60: 65 29 2f 43 4c 4f 43 4b 53 5f 50 45 52 5f 53 45  e)/CLOCKS_PER_SE
0a70: 43 2a 31 30 30 30 29 20 3c 3c 20 22 20 6d 73 65  C*1000) << " mse
0a80: 63 29 22 3b 20 72 65 74 75 72 6e 20 6f 73 2e 73  c)"; return os.s
0a90: 74 72 28 29 3b 20 7d 0d 0a 74 65 6d 70 6c 61 74  tr(); }..templat
0aa0: 65 3c 74 79 70 65 6e 61 6d 65 20 54 3e 20 6f 73  e<typename T> os
0ab0: 74 72 65 61 6d 26 20 6f 70 65 72 61 74 6f 72 3c  tream& operator<
0ac0: 3c 28 6f 73 74 72 65 61 6d 26 20 6f 73 2c 20 63  <(ostream& os, c
0ad0: 6f 6e 73 74 20 76 65 63 74 6f 72 3c 54 3e 26 20  onst vector<T>& 
0ae0: 76 29 0d 0a 20 7b 20 6f 73 20 3c 3c 20 22 7b 20  v).. { os << "{ 
0af0: 22 3b 0d 0a 20 20 20 66 6f 72 28 74 79 70 65 6e  ";..   for(typen
0b00: 61 6d 65 20 76 65 63 74 6f 72 3c 54 3e 3a 3a 63  ame vector<T>::c
0b10: 6f 6e 73 74 5f 69 74 65 72 61 74 6f 72 20 69 74  onst_iterator it
0b20: 3d 76 2e 62 65 67 69 6e 28 29 3b 20 69 74 21 3d  =v.begin(); it!=
0b30: 76 2e 65 6e 64 28 29 3b 20 2b 2b 69 74 29 0d 0a  v.end(); ++it)..
0b40: 20 20 20 6f 73 20 3c 3c 20 27 5c 22 27 20 3c 3c     os << '\"' <<
0b50: 20 2a 69 74 20 3c 3c 20 27 5c 22 27 20 3c 3c 20   *it << '\"' << 
0b60: 28 69 74 2b 31 3d 3d 76 2e 65 6e 64 28 29 20 3f  (it+1==v.end() ?
0b70: 20 22 22 20 3a 20 22 2c 20 22 29 3b 20 6f 73 20   "" : ", "); os 
0b80: 3c 3c 20 22 20 7d 22 3b 20 72 65 74 75 72 6e 20  << " }"; return 
0b90: 6f 73 3b 20 7d 0d 0a 76 6f 69 64 20 76 65 72 69  os; }..void veri
0ba0: 66 79 5f 63 61 73 65 28 63 6f 6e 73 74 20 69 6e  fy_case(const in
0bb0: 74 26 20 45 78 70 65 63 74 65 64 2c 20 63 6f 6e  t& Expected, con
0bc0: 73 74 20 69 6e 74 26 20 52 65 63 65 69 76 65 64  st int& Received
0bd0: 29 20 7b 0d 0a 20 62 6f 6f 6c 20 6f 6b 20 3d 20  ) {.. bool ok = 
0be0: 28 45 78 70 65 63 74 65 64 20 3d 3d 20 52 65 63  (Expected == Rec
0bf0: 65 69 76 65 64 29 3b 0d 0a 20 69 66 28 6f 6b 29  eived);.. if(ok)
0c00: 20 63 65 72 72 20 3c 3c 20 22 50 41 53 53 45 44   cerr << "PASSED
0c10: 22 20 3c 3c 20 74 69 6d 65 72 28 29 20 3c 3c 20  " << timer() << 
0c20: 65 6e 64 6c 3b 20 20 65 6c 73 65 20 7b 20 63 65  endl;  else { ce
0c30: 72 72 20 3c 3c 20 22 46 41 49 4c 45 44 22 20 3c  rr << "FAILED" <
0c40: 3c 20 74 69 6d 65 72 28 29 20 3c 3c 20 65 6e 64  < timer() << end
0c50: 6c 3b 0d 0a 20 63 65 72 72 20 3c 3c 20 22 5c 74  l;.. cerr << "\t
0c60: 6f 3a 20 5c 22 22 20 3c 3c 20 45 78 70 65 63 74  o: \"" << Expect
0c70: 65 64 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 65 6e  ed << '\"' << en
0c80: 64 6c 20 3c 3c 20 22 5c 74 78 3a 20 5c 22 22 20  dl << "\tx: \"" 
0c90: 3c 3c 20 52 65 63 65 69 76 65 64 20 3c 3c 20 27  << Received << '
0ca0: 5c 22 27 20 3c 3c 20 65 6e 64 6c 3b 20 7d 20 7d  \"' << endl; } }
0cb0: 0d 0a 23 64 65 66 69 6e 65 20 43 41 53 45 28 4e  ..#define CASE(N
0cc0: 29 20 7b 63 65 72 72 20 3c 3c 20 22 54 65 73 74  ) {cerr << "Test
0cd0: 20 43 61 73 65 20 23 22 20 3c 3c 20 4e 20 3c 3c   Case #" << N <<
0ce0: 20 22 2e 2e 2e 22 20 3c 3c 20 66 6c 75 73 68 3b   "..." << flush;
0cf0: 20 73 74 61 72 74 5f 74 69 6d 65 3d 63 6c 6f 63   start_time=cloc
0d00: 6b 28 29 3b 0d 0a 23 64 65 66 69 6e 65 20 45 4e  k();..#define EN
0d10: 44 09 20 76 65 72 69 66 79 5f 63 61 73 65 28 5f  D. verify_case(_
0d20: 2c 20 53 75 6d 4f 76 65 72 50 65 72 6d 75 74 61  , SumOverPermuta
0d30: 74 69 6f 6e 73 28 29 2e 66 69 6e 64 53 75 6d 28  tions().findSum(
0d40: 6e 29 29 3b 7d 0d 0a 69 6e 74 20 6d 61 69 6e 28  n));}..int main(
0d50: 29 7b 0d 0a 0d 0a 43 41 53 45 28 30 29 0d 0a 09  ){....CASE(0)...
0d60: 69 6e 74 20 6e 20 3d 20 32 3b 20 0d 0a 09 69 6e  int n = 2; ...in
0d70: 74 20 5f 20 3d 20 34 3b 20 0d 0a 45 4e 44 0d 0a  t _ = 4; ..END..
0d80: 43 41 53 45 28 31 29 0d 0a 09 69 6e 74 20 6e 20  CASE(1)...int n 
0d90: 3d 20 33 3b 20 0d 0a 09 69 6e 74 20 5f 20 3d 20  = 3; ...int _ = 
0da0: 36 36 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28  66; ..END..CASE(
0db0: 32 29 0d 0a 09 69 6e 74 20 6e 20 3d 20 31 30 3b  2)...int n = 10;
0dc0: 20 0d 0a 09 69 6e 74 20 5f 20 3d 20 35 38 33 31   ...int _ = 5831
0dd0: 30 31 31 34 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53  0114; ..END..CAS
0de0: 45 28 33 29 0d 0a 09 69 6e 74 20 6e 20 3d 20 33  E(3)...int n = 3
0df0: 39 30 30 3b 20 0d 0a 09 69 6e 74 20 5f 20 3d 20  900; ...int _ = 
0e00: 39 34 30 35 36 30 38 31 34 3b 20 0d 0a 45 4e 44  940560814; ..END
0e10: 0d 0a 2f 2a 0d 0a 43 41 53 45 28 34 29 0d 0a 09  ../*..CASE(4)...
0e20: 69 6e 74 20 6e 20 3d 20 3b 20 0d 0a 09 69 6e 74  int n = ; ...int
0e30: 20 5f 20 3d 20 3b 20 0d 0a 45 4e 44 0d 0a 43 41   _ = ; ..END..CA
0e40: 53 45 28 35 29 0d 0a 09 69 6e 74 20 6e 20 3d 20  SE(5)...int n = 
0e50: 3b 20 0d 0a 09 69 6e 74 20 5f 20 3d 20 3b 20 0d  ; ...int _ = ; .
0e60: 0a 45 4e 44 0d 0a 2a 2f 0d 0a 7d 0d 0a 2f 2f 20  .END..*/..}..// 
0e70: 45 4e 44 20 43 55 54 20 48 45 52 45 0d 0a        END CUT HERE..