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 75 73 69 6e 67 20 6e 61 6d sert>..using nam
0140: 65 73 70 61 63 65 20 73 74 64 3b 0d 0a 74 79 70 espace std;..typ
0150: 65 64 65 66 20 6c 6f 6e 67 20 6c 6f 6e 67 20 4c edef long long L
0160: 4c 3b 0d 0a 74 79 70 65 64 65 66 20 6c 6f 6e 67 L;..typedef long
0170: 20 64 6f 75 62 6c 65 20 4c 44 3b 0d 0a 74 79 70 double LD;..typ
0180: 65 64 65 66 20 63 6f 6d 70 6c 65 78 3c 4c 44 3e edef complex<LD>
0190: 20 43 4d 50 3b 0d 0a 0d 0a 73 74 61 74 69 63 20 CMP;....static
01a0: 63 6f 6e 73 74 20 75 6e 73 69 67 6e 65 64 20 4d const unsigned M
01b0: 4f 44 56 41 4c 20 3d 20 31 30 30 30 30 30 30 30 ODVAL = 10000000
01c0: 30 39 3b 0d 0a 73 74 72 75 63 74 20 6d 69 6e 74 09;..struct mint
01d0: 0d 0a 7b 0d 0a 09 75 6e 73 69 67 6e 65 64 20 76 ..{...unsigned v
01e0: 61 6c 3b 0d 0a 09 6d 69 6e 74 28 29 3a 76 61 6c al;...mint():val
01f0: 28 30 29 7b 7d 0d 0a 09 6d 69 6e 74 28 69 6e 74 (0){}...mint(int
0200: 20 20 20 20 20 20 78 29 3a 76 61 6c 28 78 25 4d x):val(x%M
0210: 4f 44 56 41 4c 29 20 7b 7d 0d 0a 09 6d 69 6e 74 ODVAL) {}...mint
0220: 28 75 6e 73 69 67 6e 65 64 20 78 29 3a 76 61 6c (unsigned x):val
0230: 28 78 25 4d 4f 44 56 41 4c 29 20 7b 7d 0d 0a 09 (x%MODVAL) {}...
0240: 6d 69 6e 74 28 4c 4c 20 20 20 20 20 20 20 78 29 mint(LL x)
0250: 3a 76 61 6c 28 78 25 4d 4f 44 56 41 4c 29 20 7b :val(x%MODVAL) {
0260: 7d 0d 0a 7d 3b 0d 0a 6d 69 6e 74 26 20 6f 70 65 }..};..mint& ope
0270: 72 61 74 6f 72 2b 3d 28 6d 69 6e 74 26 20 78 2c rator+=(mint& x,
0280: 20 6d 69 6e 74 20 79 29 20 7b 20 72 65 74 75 72 mint y) { retur
0290: 6e 20 78 20 3d 20 78 2e 76 61 6c 2b 79 2e 76 61 n x = x.val+y.va
02a0: 6c 3b 20 7d 0d 0a 6d 69 6e 74 26 20 6f 70 65 72 l; }..mint& oper
02b0: 61 74 6f 72 2d 3d 28 6d 69 6e 74 26 20 78 2c 20 ator-=(mint& x,
02c0: 6d 69 6e 74 20 79 29 20 7b 20 72 65 74 75 72 6e mint y) { return
02d0: 20 78 20 3d 20 78 2e 76 61 6c 2d 79 2e 76 61 6c x = x.val-y.val
02e0: 2b 4d 4f 44 56 41 4c 3b 20 7d 0d 0a 6d 69 6e 74 +MODVAL; }..mint
02f0: 26 20 6f 70 65 72 61 74 6f 72 2a 3d 28 6d 69 6e & operator*=(min
0300: 74 26 20 78 2c 20 6d 69 6e 74 20 79 29 20 7b 20 t& x, mint y) {
0310: 72 65 74 75 72 6e 20 78 20 3d 20 4c 4c 28 78 2e return x = LL(x.
0320: 76 61 6c 29 2a 79 2e 76 61 6c 3b 20 7d 0d 0a 6d val)*y.val; }..m
0330: 69 6e 74 20 6f 70 65 72 61 74 6f 72 2b 28 6d 69 int operator+(mi
0340: 6e 74 20 78 2c 20 6d 69 6e 74 20 79 29 20 7b 20 nt x, mint y) {
0350: 72 65 74 75 72 6e 20 78 2b 3d 79 3b 20 7d 0d 0a return x+=y; }..
0360: 6d 69 6e 74 20 6f 70 65 72 61 74 6f 72 2d 28 6d mint operator-(m
0370: 69 6e 74 20 78 2c 20 6d 69 6e 74 20 79 29 20 7b int x, mint y) {
0380: 20 72 65 74 75 72 6e 20 78 2d 3d 79 3b 20 7d 0d return x-=y; }.
0390: 0a 6d 69 6e 74 20 6f 70 65 72 61 74 6f 72 2a 28 .mint operator*(
03a0: 6d 69 6e 74 20 78 2c 20 6d 69 6e 74 20 79 29 20 mint x, mint y)
03b0: 7b 20 72 65 74 75 72 6e 20 78 2a 3d 79 3b 20 7d { return x*=y; }
03c0: 0d 0a 0d 0a 74 65 6d 70 6c 61 74 65 3c 74 79 70 ....template<typ
03d0: 65 6e 61 6d 65 20 54 3e 0d 0a 73 74 72 75 63 74 ename T>..struct
03e0: 20 44 50 34 78 0d 0a 7b 0d 0a 09 69 6e 74 20 4e DP4x..{...int N
03f0: 31 2c 20 4e 32 2c 20 4e 33 2c 20 4e 34 3b 0d 0a 1, N2, N3, N4;..
0400: 09 76 65 63 74 6f 72 3c 54 3e 20 64 61 74 61 3b .vector<T> data;
0410: 0d 0a 09 44 50 34 78 28 69 6e 74 2c 20 69 6e 74 ...DP4x(int, int
0420: 20 4e 32 2c 20 69 6e 74 20 4e 33 2c 20 69 6e 74 N2, int N3, int
0430: 20 4e 34 2c 20 63 6f 6e 73 74 20 54 26 20 74 20 N4, const T& t
0440: 3d 20 54 28 29 29 0d 0a 09 09 3a 20 4e 31 28 32 = T())....: N1(2
0450: 29 2c 20 4e 32 28 4e 32 29 2c 20 4e 33 28 4e 33 ), N2(N2), N3(N3
0460: 29 2c 20 4e 34 28 4e 34 29 2c 20 64 61 74 61 28 ), N4(N4), data(
0470: 4e 31 2a 4e 32 2a 4e 33 2a 4e 34 2c 20 74 29 20 N1*N2*N3*N4, t)
0480: 7b 20 61 73 73 65 72 74 28 64 61 74 61 2e 73 69 { assert(data.si
0490: 7a 65 28 29 2a 73 69 7a 65 6f 66 28 54 29 3c 28 ze()*sizeof(T)<(
04a0: 31 3c 3c 32 36 29 29 3b 20 7d 0d 0a 09 54 26 20 1<<26)); }...T&
04b0: 6f 70 65 72 61 74 6f 72 28 29 28 69 6e 74 20 69 operator()(int i
04c0: 31 2c 20 69 6e 74 20 69 32 2c 20 69 6e 74 20 69 1, int i2, int i
04d0: 33 2c 20 69 6e 74 20 69 34 29 0d 0a 09 09 7b 20 3, int i4)....{
04e0: 69 31 26 3d 31 3b 20 72 65 74 75 72 6e 20 64 61 i1&=1; return da
04f0: 74 61 5b 20 28 28 28 69 31 2a 4e 32 29 2b 69 32 ta[ (((i1*N2)+i2
0500: 29 2a 4e 33 2b 69 33 29 2a 4e 34 2b 69 34 20 5d )*N3+i3)*N4+i4 ]
0510: 3b 20 7d 0d 0a 09 76 6f 69 64 20 73 77 61 70 28 ; }...void swap(
0520: 44 50 34 78 26 20 72 68 73 29 0d 0a 09 09 7b 20 DP4x& rhs)....{
0530: 64 61 74 61 2e 73 77 61 70 28 72 68 73 2e 64 61 data.swap(rhs.da
0540: 74 61 29 3b 20 7d 0d 0a 7d 3b 0d 0a 0d 0a 63 6c ta); }..};....cl
0550: 61 73 73 20 54 68 65 45 78 70 65 72 69 6d 65 6e ass TheExperimen
0560: 74 20 7b 20 70 75 62 6c 69 63 3a 0d 0a 09 69 6e t { public:...in
0570: 74 20 63 6f 75 6e 74 50 6c 61 63 65 6d 65 6e 74 t countPlacement
0580: 73 28 76 65 63 74 6f 72 20 3c 73 74 72 69 6e 67 s(vector <string
0590: 3e 20 69 6e 74 65 6e 73 69 74 79 2c 20 69 6e 74 > intensity, int
05a0: 20 4d 2c 20 69 6e 74 20 4c 2c 20 69 6e 74 20 41 M, int L, int A
05b0: 2c 20 69 6e 74 20 42 29 0d 0a 09 7b 0d 0a 09 09 , int B)...{....
05c0: 73 74 72 69 6e 67 20 53 20 3d 20 61 63 63 75 6d string S = accum
05d0: 75 6c 61 74 65 28 69 6e 74 65 6e 73 69 74 79 2e ulate(intensity.
05e0: 62 65 67 69 6e 28 29 2c 20 69 6e 74 65 6e 73 69 begin(), intensi
05f0: 74 79 2e 65 6e 64 28 29 2c 20 73 74 72 69 6e 67 ty.end(), string
0600: 28 29 29 3b 0d 0a 09 09 69 6e 74 20 4e 20 3d 20 ());....int N =
0610: 53 2e 73 69 7a 65 28 29 3b 0d 0a 0d 0a 09 09 76 S.size();......v
0620: 65 63 74 6f 72 3c 69 6e 74 3e 20 72 61 6e 67 65 ector<int> range
0630: 5f 4c 2c 20 72 61 6e 67 65 5f 52 3b 20 2f 2f 20 _L, range_R; //
0640: 5b 4c 2c 52 29 0d 0a 09 09 7b 0d 0a 09 09 09 76 [L,R)....{.....v
0650: 65 63 74 6f 72 3c 69 6e 74 3e 20 50 53 28 4e 2b ector<int> PS(N+
0660: 31 2c 20 30 29 3b 0d 0a 09 09 09 66 6f 72 28 69 1, 0);.....for(i
0670: 6e 74 20 69 3d 30 3b 20 69 3c 4e 3b 20 2b 2b 69 nt i=0; i<N; ++i
0680: 29 0d 0a 09 09 09 09 50 53 5b 69 2b 31 5d 20 3d )......PS[i+1] =
0690: 20 50 53 5b 69 5d 20 2b 20 28 53 5b 69 5d 2d 27 PS[i] + (S[i]-'
06a0: 30 27 29 3b 0d 0a 09 09 09 66 6f 72 28 69 6e 74 0');.....for(int
06b0: 20 69 3d 30 3b 20 69 3c 4e 3b 20 2b 2b 69 29 0d i=0; i<N; ++i).
06c0: 0a 09 09 09 7b 0d 0a 09 09 09 09 69 6e 74 20 4c ....{......int L
06d0: 4c 20 3d 20 50 53 5b 69 5d 2b 41 3b 0d 0a 09 09 L = PS[i]+A;....
06e0: 09 09 69 6e 74 20 52 52 20 3d 20 50 53 5b 69 5d ..int RR = PS[i]
06f0: 2b 42 2b 31 3b 0d 0a 09 09 09 09 72 61 6e 67 65 +B+1;......range
0700: 5f 4c 2e 70 75 73 68 5f 62 61 63 6b 28 6c 6f 77 _L.push_back(low
0710: 65 72 5f 62 6f 75 6e 64 28 50 53 2e 62 65 67 69 er_bound(PS.begi
0720: 6e 28 29 2c 20 50 53 2e 65 6e 64 28 29 2c 20 4c n(), PS.end(), L
0730: 4c 29 20 2d 20 50 53 2e 62 65 67 69 6e 28 29 29 L) - PS.begin())
0740: 3b 0d 0a 09 09 09 09 72 61 6e 67 65 5f 52 2e 70 ;......range_R.p
0750: 75 73 68 5f 62 61 63 6b 28 6c 6f 77 65 72 5f 62 ush_back(lower_b
0760: 6f 75 6e 64 28 50 53 2e 62 65 67 69 6e 28 29 2c ound(PS.begin(),
0770: 20 50 53 2e 65 6e 64 28 29 2c 20 52 52 29 20 2d PS.end(), RR) -
0780: 20 50 53 2e 62 65 67 69 6e 28 29 29 3b 0d 0a 09 PS.begin());...
0790: 09 09 09 2f 2f 63 65 72 72 3c 3c 69 3c 3c 22 20 ...//cerr<<i<<"
07a0: 5b 22 3c 3c 72 61 6e 67 65 5f 4c 2e 62 61 63 6b ["<<range_L.back
07b0: 28 29 3c 3c 22 2d 22 3c 3c 72 61 6e 67 65 5f 52 ()<<"-"<<range_R
07c0: 2e 62 61 63 6b 28 29 3c 3c 22 29 22 3c 3c 22 20 .back()<<")"<<"
07d0: 22 3c 3c 4c 4c 3c 3c 22 2c 22 3c 3c 52 52 3c 3c "<<LL<<","<<RR<<
07e0: 65 6e 64 6c 3b 0d 0a 09 09 09 7d 0d 0a 09 09 7d endl;.....}....}
07f0: 0d 0a 0d 0a 09 09 44 50 34 78 3c 6d 69 6e 74 3e ......DP4x<mint>
0800: 20 64 70 28 4d 2b 31 2c 20 4e 2b 31 2c 20 4d 2c dp(M+1, N+1, M,
0810: 20 32 29 3b 0d 0a 09 09 64 70 28 30 2c 20 30 2c 2);....dp(0, 0,
0820: 20 30 2c 20 30 29 20 3d 20 31 3b 0d 0a 09 09 66 0, 0) = 1;....f
0830: 6f 72 28 69 6e 74 20 69 3d 30 3b 20 69 3c 4d 3b or(int i=0; i<M;
0840: 20 2b 2b 69 29 0d 0a 09 09 7b 0d 0a 09 09 09 2f ++i)....{...../
0850: 2f 20 43 6c 65 61 72 0d 0a 09 09 09 66 6f 72 28 / Clear.....for(
0860: 69 6e 74 20 78 3d 30 3b 20 78 3c 3d 4e 3b 20 2b int x=0; x<=N; +
0870: 2b 78 29 0d 0a 09 09 09 66 6f 72 28 69 6e 74 20 +x).....for(int
0880: 72 61 6e 6b 3d 30 3b 20 72 61 6e 6b 3c 4d 3b 20 rank=0; rank<M;
0890: 2b 2b 72 61 6e 6b 29 20 7b 0d 0a 09 09 09 66 6f ++rank) {.....fo
08a0: 72 28 69 6e 74 20 68 61 6d 69 3d 30 3b 20 68 61 r(int hami=0; ha
08b0: 6d 69 3c 32 3b 20 2b 2b 68 61 6d 69 29 20 7b 0d mi<2; ++hami) {.
08c0: 0a 09 09 09 09 69 66 28 64 70 28 69 2c 78 2c 72 .....if(dp(i,x,r
08d0: 61 6e 6b 2c 68 61 6d 69 29 2e 76 61 6c 29 0d 0a ank,hami).val)..
08e0: 09 09 09 09 63 65 72 72 3c 3c 69 3c 3c 22 20 22 ....cerr<<i<<" "
08f0: 3c 3c 78 3c 3c 22 20 22 3c 3c 72 61 6e 6b 3c 3c <<x<<" "<<rank<<
0900: 22 20 22 3c 3c 68 61 6d 69 3c 3c 22 20 3a 20 22 " "<<hami<<" : "
0910: 3c 3c 20 64 70 28 69 2c 78 2c 72 61 6e 6b 2c 68 << dp(i,x,rank,h
0920: 61 6d 69 29 2e 76 61 6c 3c 3c 65 6e 64 6c 3b 0d ami).val<<endl;.
0930: 0a 09 09 09 09 64 70 28 69 2b 31 2c 20 78 2c 20 .....dp(i+1, x,
0940: 72 61 6e 6b 2c 20 68 61 6d 69 29 20 3d 20 30 3b rank, hami) = 0;
0950: 0d 0a 09 09 09 7d 0d 0a 0d 0a 09 09 09 2f 2f 20 .....}.......//
0960: 44 50 0d 0a 09 09 09 66 6f 72 28 69 6e 74 20 78 DP.....for(int x
0970: 3d 30 3b 20 78 3c 4e 3b 20 2b 2b 78 29 0d 0a 09 =0; x<N; ++x)...
0980: 09 09 66 6f 72 28 69 6e 74 20 72 61 6e 6b 3d 30 ..for(int rank=0
0990: 3b 20 72 61 6e 6b 3c 3d 69 3b 20 2b 2b 72 61 6e ; rank<=i; ++ran
09a0: 6b 29 0d 0a 09 09 09 66 6f 72 28 69 6e 74 20 68 k).....for(int h
09b0: 61 6d 69 3d 30 3b 20 68 61 6d 69 3c 32 3b 20 2b ami=0; hami<2; +
09c0: 2b 68 61 6d 69 29 0d 0a 09 09 09 7b 0d 0a 09 09 +hami).....{....
09d0: 09 09 6d 69 6e 74 20 76 20 3d 20 64 70 28 69 2c ..mint v = dp(i,
09e0: 78 2c 72 61 6e 6b 2c 68 61 6d 69 29 3b 0d 0a 09 x,rank,hami);...
09f0: 09 09 09 69 66 28 76 2e 76 61 6c 20 21 3d 20 30 ...if(v.val != 0
0a00: 29 0d 0a 09 09 09 09 09 66 6f 72 28 69 6e 74 20 ).......for(int
0a10: 79 3d 72 61 6e 67 65 5f 4c 5b 78 5d 3b 20 79 3c y=range_L[x]; y<
0a20: 72 61 6e 67 65 5f 52 5b 78 5d 3b 20 2b 2b 79 29 range_R[x]; ++y)
0a30: 0d 0a 09 09 09 09 09 7b 0d 0a 09 09 09 09 09 09 .......{........
0a40: 69 66 28 79 2d 78 20 3c 20 4c 29 20 7b 0d 0a 09 if(y-x < L) {...
0a50: 09 09 09 09 09 09 69 66 28 68 61 6d 69 29 20 7b ......if(hami) {
0a60: 0d 0a 09 09 09 09 09 09 09 09 66 6f 72 28 69 6e ..........for(in
0a70: 74 20 72 3d 30 3b 20 72 3c 3d 72 61 6e 6b 3b 20 t r=0; r<=rank;
0a80: 2b 2b 72 29 0d 0a 09 09 09 09 09 09 09 09 09 64 ++r)...........d
0a90: 70 28 69 2b 31 2c 20 79 2c 20 61 29 20 2b 3d 20 p(i+1, y, a) +=
0aa0: 76 3b 0d 0a 09 09 09 09 09 09 09 7d 20 65 6c 73 v;.........} els
0ab0: 65 20 7b 0d 0a 09 09 09 09 09 09 09 7d 0d 0a 09 e {.........}...
0ac0: 09 09 09 09 09 7d 20 65 6c 73 65 20 69 66 28 79 .....} else if(y
0ad0: 2d 78 20 3d 3d 20 4c 29 20 7b 0d 0a 09 09 09 09 -x == L) {......
0ae0: 09 09 09 69 66 28 68 61 6d 69 29 0d 0a 09 09 09 ...if(hami).....
0af0: 09 09 09 09 09 64 70 28 69 2b 31 2c 20 79 2c 20 .....dp(i+1, y,
0b00: 69 2c 20 30 29 20 2b 3d 20 76 20 2a 20 28 72 61 i, 0) += v * (ra
0b10: 6e 6b 2b 31 29 3b 0d 0a 09 09 09 09 09 09 09 65 nk+1);.........e
0b20: 6c 73 65 0d 0a 09 09 09 09 09 09 09 09 64 70 28 lse..........dp(
0b30: 69 2b 31 2c 20 79 2c 20 69 2c 20 30 29 20 2b 3d i+1, y, i, 0) +=
0b40: 20 76 20 2a 20 28 72 61 6e 6b 2b 31 29 3b 0d 0a v * (rank+1);..
0b50: 09 09 09 09 09 09 7d 20 65 6c 73 65 20 7b 0d 0a ......} else {..
0b60: 09 09 09 09 09 09 09 62 72 65 61 6b 3b 0d 0a 09 .......break;...
0b70: 09 09 09 09 09 7d 0d 0a 09 09 09 09 09 7d 0d 0a .....}.......}..
0b80: 09 09 09 7d 0d 0a 09 09 7d 0d 0a 09 09 6d 69 6e ...}....}....min
0b90: 74 20 61 63 63 20 3d 20 30 3b 0d 0a 09 09 66 6f t acc = 0;....fo
0ba0: 72 28 69 6e 74 20 72 61 6e 6b 3d 30 3b 20 72 61 r(int rank=0; ra
0bb0: 6e 6b 3c 4d 3b 20 2b 2b 72 61 6e 6b 29 0d 0a 09 nk<M; ++rank)...
0bc0: 09 09 61 63 63 20 2b 3d 20 64 70 28 4d 2c 20 4e ..acc += dp(M, N
0bd0: 2c 20 72 61 6e 6b 2c 20 30 29 3b 0d 0a 09 09 72 , rank, 0);....r
0be0: 65 74 75 72 6e 20 61 63 63 2e 76 61 6c 3b 0d 0a eturn acc.val;..
0bf0: 09 7d 0d 0a 7d 3b 0d 0a 0d 0a 2f 2f 20 42 45 47 .}..};....// BEG
0c00: 49 4e 20 43 55 54 20 48 45 52 45 0d 0a 23 69 6e IN CUT HERE..#in
0c10: 63 6c 75 64 65 20 3c 63 74 69 6d 65 3e 0d 0a 64 clude <ctime>..d
0c20: 6f 75 62 6c 65 20 73 74 61 72 74 5f 74 69 6d 65 ouble start_time
0c30: 3b 20 73 74 72 69 6e 67 20 74 69 6d 65 72 28 29 ; string timer()
0c40: 0d 0a 20 7b 20 6f 73 74 72 69 6e 67 73 74 72 65 .. { ostringstre
0c50: 61 6d 20 6f 73 3b 20 6f 73 20 3c 3c 20 22 20 28 am os; os << " (
0c60: 22 20 3c 3c 20 69 6e 74 28 28 63 6c 6f 63 6b 28 " << int((clock(
0c70: 29 2d 73 74 61 72 74 5f 74 69 6d 65 29 2f 43 4c )-start_time)/CL
0c80: 4f 43 4b 53 5f 50 45 52 5f 53 45 43 2a 31 30 30 OCKS_PER_SEC*100
0c90: 30 29 20 3c 3c 20 22 20 6d 73 65 63 29 22 3b 20 0) << " msec)";
0ca0: 72 65 74 75 72 6e 20 6f 73 2e 73 74 72 28 29 3b return os.str();
0cb0: 20 7d 0d 0a 74 65 6d 70 6c 61 74 65 3c 74 79 70 }..template<typ
0cc0: 65 6e 61 6d 65 20 54 3e 20 6f 73 74 72 65 61 6d ename T> ostream
0cd0: 26 20 6f 70 65 72 61 74 6f 72 3c 3c 28 6f 73 74 & operator<<(ost
0ce0: 72 65 61 6d 26 20 6f 73 2c 20 63 6f 6e 73 74 20 ream& os, const
0cf0: 76 65 63 74 6f 72 3c 54 3e 26 20 76 29 0d 0a 20 vector<T>& v)..
0d00: 7b 20 6f 73 20 3c 3c 20 22 7b 20 22 3b 0d 0a 20 { os << "{ ";..
0d10: 20 20 66 6f 72 28 74 79 70 65 6e 61 6d 65 20 76 for(typename v
0d20: 65 63 74 6f 72 3c 54 3e 3a 3a 63 6f 6e 73 74 5f ector<T>::const_
0d30: 69 74 65 72 61 74 6f 72 20 69 74 3d 76 2e 62 65 iterator it=v.be
0d40: 67 69 6e 28 29 3b 20 69 74 21 3d 76 2e 65 6e 64 gin(); it!=v.end
0d50: 28 29 3b 20 2b 2b 69 74 29 0d 0a 20 20 20 6f 73 (); ++it).. os
0d60: 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 2a 69 74 20 << '\"' << *it
0d70: 3c 3c 20 27 5c 22 27 20 3c 3c 20 28 69 74 2b 31 << '\"' << (it+1
0d80: 3d 3d 76 2e 65 6e 64 28 29 20 3f 20 22 22 20 3a ==v.end() ? "" :
0d90: 20 22 2c 20 22 29 3b 20 6f 73 20 3c 3c 20 22 20 ", "); os << "
0da0: 7d 22 3b 20 72 65 74 75 72 6e 20 6f 73 3b 20 7d }"; return os; }
0db0: 0d 0a 76 6f 69 64 20 76 65 72 69 66 79 5f 63 61 ..void verify_ca
0dc0: 73 65 28 63 6f 6e 73 74 20 69 6e 74 26 20 45 78 se(const int& Ex
0dd0: 70 65 63 74 65 64 2c 20 63 6f 6e 73 74 20 69 6e pected, const in
0de0: 74 26 20 52 65 63 65 69 76 65 64 29 20 7b 0d 0a t& Received) {..
0df0: 20 62 6f 6f 6c 20 6f 6b 20 3d 20 28 45 78 70 65 bool ok = (Expe
0e00: 63 74 65 64 20 3d 3d 20 52 65 63 65 69 76 65 64 cted == Received
0e10: 29 3b 0d 0a 20 69 66 28 6f 6b 29 20 63 65 72 72 );.. if(ok) cerr
0e20: 20 3c 3c 20 22 50 41 53 53 45 44 22 20 3c 3c 20 << "PASSED" <<
0e30: 74 69 6d 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b timer() << endl;
0e40: 20 20 65 6c 73 65 20 7b 20 63 65 72 72 20 3c 3c else { cerr <<
0e50: 20 22 46 41 49 4c 45 44 22 20 3c 3c 20 74 69 6d "FAILED" << tim
0e60: 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b 0d 0a 20 er() << endl;..
0e70: 63 65 72 72 20 3c 3c 20 22 5c 74 6f 3a 20 5c 22 cerr << "\to: \"
0e80: 22 20 3c 3c 20 45 78 70 65 63 74 65 64 20 3c 3c " << Expected <<
0e90: 20 27 5c 22 27 20 3c 3c 20 65 6e 64 6c 20 3c 3c '\"' << endl <<
0ea0: 20 22 5c 74 78 3a 20 5c 22 22 20 3c 3c 20 52 65 "\tx: \"" << Re
0eb0: 63 65 69 76 65 64 20 3c 3c 20 27 5c 22 27 20 3c ceived << '\"' <
0ec0: 3c 20 65 6e 64 6c 3b 20 7d 20 7d 0d 0a 23 64 65 < endl; } }..#de
0ed0: 66 69 6e 65 20 43 41 53 45 28 4e 29 20 7b 63 65 fine CASE(N) {ce
0ee0: 72 72 20 3c 3c 20 22 54 65 73 74 20 43 61 73 65 rr << "Test Case
0ef0: 20 23 22 20 3c 3c 20 4e 20 3c 3c 20 22 2e 2e 2e #" << N << "...
0f00: 22 20 3c 3c 20 66 6c 75 73 68 3b 20 73 74 61 72 " << flush; star
0f10: 74 5f 74 69 6d 65 3d 63 6c 6f 63 6b 28 29 3b 0d t_time=clock();.
0f20: 0a 23 64 65 66 69 6e 65 20 45 4e 44 09 20 76 65 .#define END. ve
0f30: 72 69 66 79 5f 63 61 73 65 28 5f 2c 20 54 68 65 rify_case(_, The
0f40: 45 78 70 65 72 69 6d 65 6e 74 28 29 2e 63 6f 75 Experiment().cou
0f50: 6e 74 50 6c 61 63 65 6d 65 6e 74 73 28 69 6e 74 ntPlacements(int
0f60: 65 6e 73 69 74 79 2c 20 4d 2c 20 4c 2c 20 41 2c ensity, M, L, A,
0f70: 20 42 29 29 3b 7d 0d 0a 69 6e 74 20 6d 61 69 6e B));}..int main
0f80: 28 29 7b 0d 0a 0d 0a 43 41 53 45 28 30 29 0d 0a (){....CASE(0)..
0f90: 09 73 74 72 69 6e 67 20 69 6e 74 65 6e 73 69 74 .string intensit
0fa0: 79 5f 5b 5d 20 3d 20 7b 22 33 34 31 31 35 36 22 y_[] = {"341156"
0fb0: 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 73 };... vector <s
0fc0: 74 72 69 6e 67 3e 20 69 6e 74 65 6e 73 69 74 79 tring> intensity
0fd0: 28 69 6e 74 65 6e 73 69 74 79 5f 2c 20 69 6e 74 (intensity_, int
0fe0: 65 6e 73 69 74 79 5f 2b 73 69 7a 65 6f 66 28 69 ensity_+sizeof(i
0ff0: 6e 74 65 6e 73 69 74 79 5f 29 2f 73 69 7a 65 6f ntensity_)/sizeo
1000: 66 28 2a 69 6e 74 65 6e 73 69 74 79 5f 29 29 3b f(*intensity_));
1010: 20 0d 0a 09 69 6e 74 20 4d 20 3d 20 33 3b 20 0d ...int M = 3; .
1020: 0a 09 69 6e 74 20 4c 20 3d 20 33 3b 20 0d 0a 09 ..int L = 3; ...
1030: 69 6e 74 20 41 20 3d 20 36 3b 20 0d 0a 09 69 6e int A = 6; ...in
1040: 74 20 42 20 3d 20 31 30 3b 20 0d 0a 09 69 6e 74 t B = 10; ...int
1050: 20 5f 20 3d 20 32 3b 20 0d 0a 45 4e 44 0d 0a 43 _ = 2; ..END..C
1060: 41 53 45 28 31 29 0d 0a 09 73 74 72 69 6e 67 20 ASE(1)...string
1070: 69 6e 74 65 6e 73 69 74 79 5f 5b 5d 20 3d 20 7b intensity_[] = {
1080: 22 39 39 39 31 31 32 39 39 39 22 7d 3b 0d 0a 09 "999112999"};...
1090: 20 20 76 65 63 74 6f 72 20 3c 73 74 72 69 6e 67 vector <string
10a0: 3e 20 69 6e 74 65 6e 73 69 74 79 28 69 6e 74 65 > intensity(inte
10b0: 6e 73 69 74 79 5f 2c 20 69 6e 74 65 6e 73 69 74 nsity_, intensit
10c0: 79 5f 2b 73 69 7a 65 6f 66 28 69 6e 74 65 6e 73 y_+sizeof(intens
10d0: 69 74 79 5f 29 2f 73 69 7a 65 6f 66 28 2a 69 6e ity_)/sizeof(*in
10e0: 74 65 6e 73 69 74 79 5f 29 29 3b 20 0d 0a 09 69 tensity_)); ...i
10f0: 6e 74 20 4d 20 3d 20 32 3b 20 0d 0a 09 69 6e 74 nt M = 2; ...int
1100: 20 4c 20 3d 20 34 3b 20 0d 0a 09 69 6e 74 20 41 L = 4; ...int A
1110: 20 3d 20 32 31 3b 20 0d 0a 09 69 6e 74 20 42 20 = 21; ...int B
1120: 3d 20 33 30 3b 20 0d 0a 09 69 6e 74 20 5f 20 3d = 30; ...int _ =
1130: 20 32 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 2; ..END..CASE(
1140: 32 29 0d 0a 09 73 74 72 69 6e 67 20 69 6e 74 65 2)...string inte
1150: 6e 73 69 74 79 5f 5b 5d 20 3d 20 7b 22 31 31 31 nsity_[] = {"111
1160: 22 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c "};... vector <
1170: 73 74 72 69 6e 67 3e 20 69 6e 74 65 6e 73 69 74 string> intensit
1180: 79 28 69 6e 74 65 6e 73 69 74 79 5f 2c 20 69 6e y(intensity_, in
1190: 74 65 6e 73 69 74 79 5f 2b 73 69 7a 65 6f 66 28 tensity_+sizeof(
11a0: 69 6e 74 65 6e 73 69 74 79 5f 29 2f 73 69 7a 65 intensity_)/size
11b0: 6f 66 28 2a 69 6e 74 65 6e 73 69 74 79 5f 29 29 of(*intensity_))
11c0: 3b 20 0d 0a 09 69 6e 74 20 4d 20 3d 20 32 3b 20 ; ...int M = 2;
11d0: 0d 0a 09 69 6e 74 20 4c 20 3d 20 32 3b 20 0d 0a ...int L = 2; ..
11e0: 09 69 6e 74 20 41 20 3d 20 32 3b 20 0d 0a 09 69 .int A = 2; ...i
11f0: 6e 74 20 42 20 3d 20 33 3b 20 0d 0a 09 69 6e 74 nt B = 3; ...int
1200: 20 5f 20 3d 20 30 3b 20 0d 0a 45 4e 44 0d 0a 43 _ = 0; ..END..C
1210: 41 53 45 28 33 29 0d 0a 09 73 74 72 69 6e 67 20 ASE(3)...string
1220: 69 6e 74 65 6e 73 69 74 79 5f 5b 5d 20 3d 20 7b intensity_[] = {
1230: 22 35 39 30 35 39 22 2c 20 22 31 31 30 22 2c 20 "59059", "110",
1240: 22 31 31 33 32 32 33 30 32 33 31 22 7d 3b 0d 0a "1132230231"};..
1250: 09 20 20 76 65 63 74 6f 72 20 3c 73 74 72 69 6e . vector <strin
1260: 67 3e 20 69 6e 74 65 6e 73 69 74 79 28 69 6e 74 g> intensity(int
1270: 65 6e 73 69 74 79 5f 2c 20 69 6e 74 65 6e 73 69 ensity_, intensi
1280: 74 79 5f 2b 73 69 7a 65 6f 66 28 69 6e 74 65 6e ty_+sizeof(inten
1290: 73 69 74 79 5f 29 2f 73 69 7a 65 6f 66 28 2a 69 sity_)/sizeof(*i
12a0: 6e 74 65 6e 73 69 74 79 5f 29 29 3b 20 0d 0a 09 ntensity_)); ...
12b0: 69 6e 74 20 4d 20 3d 20 31 3b 20 0d 0a 09 69 6e int M = 1; ...in
12c0: 74 20 4c 20 3d 20 35 3b 20 0d 0a 09 69 6e 74 20 t L = 5; ...int
12d0: 41 20 3d 20 31 30 3b 20 0d 0a 09 69 6e 74 20 42 A = 10; ...int B
12e0: 20 3d 20 32 30 3b 20 0d 0a 09 69 6e 74 20 5f 20 = 20; ...int _
12f0: 3d 20 36 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 = 6; ..END..CASE
1300: 28 34 29 0d 0a 09 73 74 72 69 6e 67 20 69 6e 74 (4)...string int
1310: 65 6e 73 69 74 79 5f 5b 5d 20 3d 20 7b 22 31 31 ensity_[] = {"11
1320: 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 1111111111111111
1330: 31 31 31 31 31 31 22 2c 20 22 31 31 31 31 31 31 111111", "111111
1340: 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 1111111111111111
1350: 31 31 31 31 31 22 2c 20 22 31 31 31 31 31 31 31 11111", "1111111
1360: 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 1111111111111111
1370: 22 7d 3b 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c "};... vector <
1380: 73 74 72 69 6e 67 3e 20 69 6e 74 65 6e 73 69 74 string> intensit
1390: 79 28 69 6e 74 65 6e 73 69 74 79 5f 2c 20 69 6e y(intensity_, in
13a0: 74 65 6e 73 69 74 79 5f 2b 73 69 7a 65 6f 66 28 tensity_+sizeof(
13b0: 69 6e 74 65 6e 73 69 74 79 5f 29 2f 73 69 7a 65 intensity_)/size
13c0: 6f 66 28 2a 69 6e 74 65 6e 73 69 74 79 5f 29 29 of(*intensity_))
13d0: 3b 20 0d 0a 09 69 6e 74 20 4d 20 3d 20 31 32 3b ; ...int M = 12;
13e0: 20 0d 0a 09 69 6e 74 20 4c 20 3d 20 38 3b 20 0d ...int L = 8; .
13f0: 0a 09 69 6e 74 20 41 20 3d 20 34 3b 20 0d 0a 09 ..int A = 4; ...
1400: 69 6e 74 20 42 20 3d 20 32 37 30 30 3b 20 0d 0a int B = 2700; ..
1410: 09 69 6e 74 20 5f 20 3d 20 34 31 38 36 32 39 39 .int _ = 4186299
1420: 34 38 3b 20 0d 0a 45 4e 44 0d 0a 2f 2a 0d 0a 43 48; ..END../*..C
1430: 41 53 45 28 35 29 0d 0a 09 73 74 72 69 6e 67 20 ASE(5)...string
1440: 69 6e 74 65 6e 73 69 74 79 5f 5b 5d 20 3d 20 3b intensity_[] = ;
1450: 0d 0a 09 20 20 76 65 63 74 6f 72 20 3c 73 74 72 ... vector <str
1460: 69 6e 67 3e 20 69 6e 74 65 6e 73 69 74 79 28 69 ing> intensity(i
1470: 6e 74 65 6e 73 69 74 79 5f 2c 20 69 6e 74 65 6e ntensity_, inten
1480: 73 69 74 79 5f 2b 73 69 7a 65 6f 66 28 69 6e 74 sity_+sizeof(int
1490: 65 6e 73 69 74 79 5f 29 2f 73 69 7a 65 6f 66 28 ensity_)/sizeof(
14a0: 2a 69 6e 74 65 6e 73 69 74 79 5f 29 29 3b 20 0d *intensity_)); .
14b0: 0a 09 69 6e 74 20 4d 20 3d 20 3b 20 0d 0a 09 69 ..int M = ; ...i
14c0: 6e 74 20 4c 20 3d 20 3b 20 0d 0a 09 69 6e 74 20 nt L = ; ...int
14d0: 41 20 3d 20 3b 20 0d 0a 09 69 6e 74 20 42 20 3d A = ; ...int B =
14e0: 20 3b 20 0d 0a 09 69 6e 74 20 5f 20 3d 20 3b 20 ; ...int _ = ;
14f0: 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 36 29 0d 0a ..END..CASE(6)..
1500: 09 73 74 72 69 6e 67 20 69 6e 74 65 6e 73 69 74 .string intensit
1510: 79 5f 5b 5d 20 3d 20 3b 0d 0a 09 20 20 76 65 63 y_[] = ;... vec
1520: 74 6f 72 20 3c 73 74 72 69 6e 67 3e 20 69 6e 74 tor <string> int
1530: 65 6e 73 69 74 79 28 69 6e 74 65 6e 73 69 74 79 ensity(intensity
1540: 5f 2c 20 69 6e 74 65 6e 73 69 74 79 5f 2b 73 69 _, intensity_+si
1550: 7a 65 6f 66 28 69 6e 74 65 6e 73 69 74 79 5f 29 zeof(intensity_)
1560: 2f 73 69 7a 65 6f 66 28 2a 69 6e 74 65 6e 73 69 /sizeof(*intensi
1570: 74 79 5f 29 29 3b 20 0d 0a 09 69 6e 74 20 4d 20 ty_)); ...int M
1580: 3d 20 3b 20 0d 0a 09 69 6e 74 20 4c 20 3d 20 3b = ; ...int L = ;
1590: 20 0d 0a 09 69 6e 74 20 41 20 3d 20 3b 20 0d 0a ...int A = ; ..
15a0: 09 69 6e 74 20 42 20 3d 20 3b 20 0d 0a 09 69 6e .int B = ; ...in
15b0: 74 20 5f 20 3d 20 3b 20 0d 0a 45 4e 44 0d 0a 2a t _ = ; ..END..*
15c0: 2f 0d 0a 7d 0d 0a 2f 2f 20 45 4e 44 20 43 55 54 /..}..// END CUT
15d0: 20 48 45 52 45 0d 0a HERE..