Artifact cbc7da6c8638afb4040f558e44a184a6b83483eb:
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 50 72 6f 62 P;....class Prob
01a0: 6c 65 6d 53 65 74 73 20 7b 20 70 75 62 6c 69 63 lemSets { public
01b0: 3a 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 20 6d 61 :...long long ma
01c0: 78 53 65 74 73 28 6c 6f 6e 67 20 6c 6f 6e 67 20 xSets(long long
01d0: 45 2c 20 6c 6f 6e 67 20 6c 6f 6e 67 20 45 4d 2c E, long long EM,
01e0: 20 6c 6f 6e 67 20 6c 6f 6e 67 20 4d 2c 20 6c 6f long long M, lo
01f0: 6e 67 20 6c 6f 6e 67 20 4d 48 2c 20 6c 6f 6e 67 ng long MH, long
0200: 20 6c 6f 6e 67 20 48 29 0d 0a 09 7b 0d 0a 09 09 long H)...{....
0210: 4c 4c 20 4c 3d 30 2c 20 52 3d 30 78 37 66 66 66 LL L=0, R=0x7fff
0220: 66 66 66 66 66 66 66 66 66 66 66 66 4c 4c 3b 20 ffffffffffffLL;
0230: 2f 2f 20 63 61 6e 20 6d 61 6b 65 20 5b 4c 2c 52 // can make [L,R
0240: 29 0d 0a 09 09 77 68 69 6c 65 28 52 2d 4c 3e 31 )....while(R-L>1
0250: 29 20 7b 0d 0a 09 09 09 4c 4c 20 43 20 3d 20 4c ) {.....LL C = L
0260: 2b 28 52 2d 4c 29 2f 32 3b 0d 0a 09 09 09 28 63 +(R-L)/2;.....(c
0270: 61 6e 28 45 2c 45 4d 2c 4d 2c 4d 48 2c 48 2c 43 an(E,EM,M,MH,H,C
0280: 29 20 3f 20 4c 20 3a 20 52 29 20 3d 20 43 3b 0d ) ? L : R) = C;.
0290: 0a 09 09 7d 0d 0a 09 09 72 65 74 75 72 6e 20 4c ...}....return L
02a0: 3b 0d 0a 09 7d 0d 0a 0d 0a 09 62 6f 6f 6c 20 63 ;...}.....bool c
02b0: 61 6e 28 6c 6f 6e 67 20 6c 6f 6e 67 20 45 2c 20 an(long long E,
02c0: 6c 6f 6e 67 20 6c 6f 6e 67 20 45 4d 2c 20 6c 6f long long EM, lo
02d0: 6e 67 20 6c 6f 6e 67 20 4d 2c 20 6c 6f 6e 67 20 ng long M, long
02e0: 6c 6f 6e 67 20 4d 48 2c 20 6c 6f 6e 67 20 6c 6f long MH, long lo
02f0: 6e 67 20 48 2c 20 4c 4c 20 6e 50 72 6f 62 73 29 ng H, LL nProbs)
0300: 20 7b 0d 0a 09 09 4c 4c 20 65 4e 65 65 64 20 3d {....LL eNeed =
0310: 20 6d 61 78 28 30 4c 4c 2c 20 6e 50 72 6f 62 73 max(0LL, nProbs
0320: 2d 45 29 3b 0d 0a 09 09 4c 4c 20 6d 4e 65 65 64 -E);....LL mNeed
0330: 20 3d 20 6d 61 78 28 30 4c 4c 2c 20 6e 50 72 6f = max(0LL, nPro
0340: 62 73 2d 4d 29 3b 0d 0a 09 09 4c 4c 20 68 4e 65 bs-M);....LL hNe
0350: 65 64 20 3d 20 6d 61 78 28 30 4c 4c 2c 20 6e 50 ed = max(0LL, nP
0360: 72 6f 62 73 2d 48 29 3b 0d 0a 09 09 69 66 28 45 robs-H);....if(E
0370: 4d 20 3c 20 65 4e 65 65 64 29 20 72 65 74 75 72 M < eNeed) retur
0380: 6e 20 66 61 6c 73 65 3b 0d 0a 09 09 45 4d 20 2d n false;....EM -
0390: 3d 20 65 4e 65 65 64 3b 0d 0a 09 09 69 66 28 4d = eNeed;....if(M
03a0: 48 20 3c 20 68 4e 65 65 64 29 20 72 65 74 75 72 H < hNeed) retur
03b0: 6e 20 66 61 6c 73 65 3b 0d 0a 09 09 4d 48 20 2d n false;....MH -
03c0: 3d 20 68 4e 65 65 64 3b 0d 0a 09 09 72 65 74 75 = hNeed;....retu
03d0: 72 6e 20 6d 4e 65 65 64 20 3c 3d 20 45 4d 2b 4d rn mNeed <= EM+M
03e0: 48 3b 0d 0a 09 7d 0d 0a 7d 3b 0d 0a 0d 0a 2f 2f H;...}..};....//
03f0: 20 42 45 47 49 4e 20 43 55 54 20 48 45 52 45 0d BEGIN CUT HERE.
0400: 0a 23 69 6e 63 6c 75 64 65 20 3c 63 74 69 6d 65 .#include <ctime
0410: 3e 0d 0a 64 6f 75 62 6c 65 20 73 74 61 72 74 5f >..double start_
0420: 74 69 6d 65 3b 20 73 74 72 69 6e 67 20 74 69 6d time; string tim
0430: 65 72 28 29 0d 0a 20 7b 20 6f 73 74 72 69 6e 67 er().. { ostring
0440: 73 74 72 65 61 6d 20 6f 73 3b 20 6f 73 20 3c 3c stream os; os <<
0450: 20 22 20 28 22 20 3c 3c 20 69 6e 74 28 28 63 6c " (" << int((cl
0460: 6f 63 6b 28 29 2d 73 74 61 72 74 5f 74 69 6d 65 ock()-start_time
0470: 29 2f 43 4c 4f 43 4b 53 5f 50 45 52 5f 53 45 43 )/CLOCKS_PER_SEC
0480: 2a 31 30 30 30 29 20 3c 3c 20 22 20 6d 73 65 63 *1000) << " msec
0490: 29 22 3b 20 72 65 74 75 72 6e 20 6f 73 2e 73 74 )"; return os.st
04a0: 72 28 29 3b 20 7d 0d 0a 74 65 6d 70 6c 61 74 65 r(); }..template
04b0: 3c 74 79 70 65 6e 61 6d 65 20 54 3e 20 6f 73 74 <typename T> ost
04c0: 72 65 61 6d 26 20 6f 70 65 72 61 74 6f 72 3c 3c ream& operator<<
04d0: 28 6f 73 74 72 65 61 6d 26 20 6f 73 2c 20 63 6f (ostream& os, co
04e0: 6e 73 74 20 76 65 63 74 6f 72 3c 54 3e 26 20 76 nst vector<T>& v
04f0: 29 0d 0a 20 7b 20 6f 73 20 3c 3c 20 22 7b 20 22 ).. { os << "{ "
0500: 3b 0d 0a 20 20 20 66 6f 72 28 74 79 70 65 6e 61 ;.. for(typena
0510: 6d 65 20 76 65 63 74 6f 72 3c 54 3e 3a 3a 63 6f me vector<T>::co
0520: 6e 73 74 5f 69 74 65 72 61 74 6f 72 20 69 74 3d nst_iterator it=
0530: 76 2e 62 65 67 69 6e 28 29 3b 20 69 74 21 3d 76 v.begin(); it!=v
0540: 2e 65 6e 64 28 29 3b 20 2b 2b 69 74 29 0d 0a 20 .end(); ++it)..
0550: 20 20 6f 73 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 os << '\"' <<
0560: 2a 69 74 20 3c 3c 20 27 5c 22 27 20 3c 3c 20 28 *it << '\"' << (
0570: 69 74 2b 31 3d 3d 76 2e 65 6e 64 28 29 20 3f 20 it+1==v.end() ?
0580: 22 22 20 3a 20 22 2c 20 22 29 3b 20 6f 73 20 3c "" : ", "); os <
0590: 3c 20 22 20 7d 22 3b 20 72 65 74 75 72 6e 20 6f < " }"; return o
05a0: 73 3b 20 7d 0d 0a 76 6f 69 64 20 76 65 72 69 66 s; }..void verif
05b0: 79 5f 63 61 73 65 28 63 6f 6e 73 74 20 6c 6f 6e y_case(const lon
05c0: 67 20 6c 6f 6e 67 26 20 45 78 70 65 63 74 65 64 g long& Expected
05d0: 2c 20 63 6f 6e 73 74 20 6c 6f 6e 67 20 6c 6f 6e , const long lon
05e0: 67 26 20 52 65 63 65 69 76 65 64 29 20 7b 0d 0a g& Received) {..
05f0: 20 62 6f 6f 6c 20 6f 6b 20 3d 20 28 45 78 70 65 bool ok = (Expe
0600: 63 74 65 64 20 3d 3d 20 52 65 63 65 69 76 65 64 cted == Received
0610: 29 3b 0d 0a 20 69 66 28 6f 6b 29 20 63 65 72 72 );.. if(ok) cerr
0620: 20 3c 3c 20 22 50 41 53 53 45 44 22 20 3c 3c 20 << "PASSED" <<
0630: 74 69 6d 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b timer() << endl;
0640: 20 20 65 6c 73 65 20 7b 20 63 65 72 72 20 3c 3c else { cerr <<
0650: 20 22 46 41 49 4c 45 44 22 20 3c 3c 20 74 69 6d "FAILED" << tim
0660: 65 72 28 29 20 3c 3c 20 65 6e 64 6c 3b 0d 0a 20 er() << endl;..
0670: 63 65 72 72 20 3c 3c 20 22 5c 74 6f 3a 20 5c 22 cerr << "\to: \"
0680: 22 20 3c 3c 20 45 78 70 65 63 74 65 64 20 3c 3c " << Expected <<
0690: 20 27 5c 22 27 20 3c 3c 20 65 6e 64 6c 20 3c 3c '\"' << endl <<
06a0: 20 22 5c 74 78 3a 20 5c 22 22 20 3c 3c 20 52 65 "\tx: \"" << Re
06b0: 63 65 69 76 65 64 20 3c 3c 20 27 5c 22 27 20 3c ceived << '\"' <
06c0: 3c 20 65 6e 64 6c 3b 20 7d 20 7d 0d 0a 23 64 65 < endl; } }..#de
06d0: 66 69 6e 65 20 43 41 53 45 28 4e 29 20 7b 63 65 fine CASE(N) {ce
06e0: 72 72 20 3c 3c 20 22 54 65 73 74 20 43 61 73 65 rr << "Test Case
06f0: 20 23 22 20 3c 3c 20 4e 20 3c 3c 20 22 2e 2e 2e #" << N << "...
0700: 22 20 3c 3c 20 66 6c 75 73 68 3b 20 73 74 61 72 " << flush; star
0710: 74 5f 74 69 6d 65 3d 63 6c 6f 63 6b 28 29 3b 0d t_time=clock();.
0720: 0a 23 64 65 66 69 6e 65 20 45 4e 44 09 20 76 65 .#define END. ve
0730: 72 69 66 79 5f 63 61 73 65 28 5f 2c 20 50 72 6f rify_case(_, Pro
0740: 62 6c 65 6d 53 65 74 73 28 29 2e 6d 61 78 53 65 blemSets().maxSe
0750: 74 73 28 45 2c 20 45 4d 2c 20 4d 2c 20 4d 48 2c ts(E, EM, M, MH,
0760: 20 48 29 29 3b 7d 0d 0a 69 6e 74 20 6d 61 69 6e H));}..int main
0770: 28 29 7b 0d 0a 0d 0a 43 41 53 45 28 30 29 0d 0a (){....CASE(0)..
0780: 09 6c 6f 6e 67 20 6c 6f 6e 67 20 45 20 3d 20 32 .long long E = 2
0790: 4c 4c 3b 20 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 LL; ...long long
07a0: 20 45 4d 20 3d 20 32 4c 4c 3b 20 0d 0a 09 6c 6f EM = 2LL; ...lo
07b0: 6e 67 20 6c 6f 6e 67 20 4d 20 3d 20 31 4c 4c 3b ng long M = 1LL;
07c0: 20 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 20 4d 48 ...long long MH
07d0: 20 3d 20 32 4c 4c 3b 20 0d 0a 09 6c 6f 6e 67 20 = 2LL; ...long
07e0: 6c 6f 6e 67 20 48 20 3d 20 32 4c 4c 3b 20 0d 0a long H = 2LL; ..
07f0: 09 6c 6f 6e 67 20 6c 6f 6e 67 20 5f 20 3d 20 33 .long long _ = 3
0800: 4c 4c 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 LL; ..END..CASE(
0810: 31 29 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 20 45 1)...long long E
0820: 20 3d 20 31 30 30 4c 4c 3b 20 0d 0a 09 6c 6f 6e = 100LL; ...lon
0830: 67 20 6c 6f 6e 67 20 45 4d 20 3d 20 31 30 30 4c g long EM = 100L
0840: 4c 3b 20 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 20 L; ...long long
0850: 4d 20 3d 20 31 30 30 4c 4c 3b 20 0d 0a 09 6c 6f M = 100LL; ...lo
0860: 6e 67 20 6c 6f 6e 67 20 4d 48 20 3d 20 30 4c 4c ng long MH = 0LL
0870: 3b 20 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 20 48 ; ...long long H
0880: 20 3d 20 30 4c 4c 3b 20 0d 0a 09 6c 6f 6e 67 20 = 0LL; ...long
0890: 6c 6f 6e 67 20 5f 20 3d 20 30 4c 4c 3b 20 0d 0a long _ = 0LL; ..
08a0: 45 4e 44 0d 0a 43 41 53 45 28 32 29 0d 0a 09 6c END..CASE(2)...l
08b0: 6f 6e 67 20 6c 6f 6e 67 20 45 20 3d 20 36 35 37 ong long E = 657
08c0: 4c 4c 3b 20 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 LL; ...long long
08d0: 20 45 4d 20 3d 20 36 35 37 4c 4c 3b 20 0d 0a 09 EM = 657LL; ...
08e0: 6c 6f 6e 67 20 6c 6f 6e 67 20 4d 20 3d 20 36 35 long long M = 65
08f0: 37 4c 4c 3b 20 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 7LL; ...long lon
0900: 67 20 4d 48 20 3d 20 36 35 37 4c 4c 3b 20 0d 0a g MH = 657LL; ..
0910: 09 6c 6f 6e 67 20 6c 6f 6e 67 20 48 20 3d 20 36 .long long H = 6
0920: 35 37 4c 4c 3b 20 0d 0a 09 6c 6f 6e 67 20 6c 6f 57LL; ...long lo
0930: 6e 67 20 5f 20 3d 20 31 30 39 35 4c 4c 3b 20 0d ng _ = 1095LL; .
0940: 0a 45 4e 44 0d 0a 43 41 53 45 28 33 29 0d 0a 09 .END..CASE(3)...
0950: 6c 6f 6e 67 20 6c 6f 6e 67 20 45 20 3d 20 31 4c long long E = 1L
0960: 4c 3b 20 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 20 L; ...long long
0970: 45 4d 20 3d 20 32 4c 4c 3b 20 0d 0a 09 6c 6f 6e EM = 2LL; ...lon
0980: 67 20 6c 6f 6e 67 20 4d 20 3d 20 33 4c 4c 3b 20 g long M = 3LL;
0990: 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 20 4d 48 20 ...long long MH
09a0: 3d 20 34 4c 4c 3b 20 0d 0a 09 6c 6f 6e 67 20 6c = 4LL; ...long l
09b0: 6f 6e 67 20 48 20 3d 20 35 4c 4c 3b 20 0d 0a 09 ong H = 5LL; ...
09c0: 6c 6f 6e 67 20 6c 6f 6e 67 20 5f 20 3d 20 33 4c long long _ = 3L
09d0: 4c 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 28 34 L; ..END..CASE(4
09e0: 29 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 20 45 20 )...long long E
09f0: 3d 20 31 30 30 30 30 30 30 30 30 30 30 30 30 30 = 10000000000000
0a00: 30 30 30 30 30 4c 4c 3b 20 0d 0a 09 6c 6f 6e 67 00000LL; ...long
0a10: 20 6c 6f 6e 67 20 45 4d 20 3d 20 31 30 30 30 30 long EM = 10000
0a20: 30 30 30 30 30 30 30 30 30 30 30 30 30 30 4c 4c 00000000000000LL
0a30: 3b 20 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 20 4d ; ...long long M
0a40: 20 3d 20 31 30 30 30 30 30 30 30 30 30 30 30 30 = 1000000000000
0a50: 30 30 30 30 30 30 4c 4c 3b 20 0d 0a 09 6c 6f 6e 000000LL; ...lon
0a60: 67 20 6c 6f 6e 67 20 4d 48 20 3d 20 31 30 30 30 g long MH = 1000
0a70: 30 30 30 30 30 30 30 30 30 30 30 30 30 30 30 4c 000000000000000L
0a80: 4c 3b 20 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 20 L; ...long long
0a90: 48 20 3d 20 31 30 30 30 30 30 30 30 30 30 30 30 H = 100000000000
0aa0: 30 30 30 30 30 30 30 4c 4c 3b 20 0d 0a 09 6c 6f 0000000LL; ...lo
0ab0: 6e 67 20 6c 6f 6e 67 20 5f 20 3d 20 31 36 36 36 ng long _ = 1666
0ac0: 36 36 36 36 36 36 36 36 36 36 36 36 36 36 36 4c 666666666666666L
0ad0: 4c 3b 20 0d 0a 45 4e 44 0d 0a 2f 2a 0d 0a 43 41 L; ..END../*..CA
0ae0: 53 45 28 35 29 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e SE(5)...long lon
0af0: 67 20 45 20 3d 20 4c 4c 3b 20 0d 0a 09 6c 6f 6e g E = LL; ...lon
0b00: 67 20 6c 6f 6e 67 20 45 4d 20 3d 20 4c 4c 3b 20 g long EM = LL;
0b10: 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 20 4d 20 3d ...long long M =
0b20: 20 4c 4c 3b 20 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e LL; ...long lon
0b30: 67 20 4d 48 20 3d 20 4c 4c 3b 20 0d 0a 09 6c 6f g MH = LL; ...lo
0b40: 6e 67 20 6c 6f 6e 67 20 48 20 3d 20 4c 4c 3b 20 ng long H = LL;
0b50: 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 20 5f 20 3d ...long long _ =
0b60: 20 4c 4c 3b 20 0d 0a 45 4e 44 0d 0a 43 41 53 45 LL; ..END..CASE
0b70: 28 36 29 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 20 (6)...long long
0b80: 45 20 3d 20 4c 4c 3b 20 0d 0a 09 6c 6f 6e 67 20 E = LL; ...long
0b90: 6c 6f 6e 67 20 45 4d 20 3d 20 4c 4c 3b 20 0d 0a long EM = LL; ..
0ba0: 09 6c 6f 6e 67 20 6c 6f 6e 67 20 4d 20 3d 20 4c .long long M = L
0bb0: 4c 3b 20 0d 0a 09 6c 6f 6e 67 20 6c 6f 6e 67 20 L; ...long long
0bc0: 4d 48 20 3d 20 4c 4c 3b 20 0d 0a 09 6c 6f 6e 67 MH = LL; ...long
0bd0: 20 6c 6f 6e 67 20 48 20 3d 20 4c 4c 3b 20 0d 0a long H = LL; ..
0be0: 09 6c 6f 6e 67 20 6c 6f 6e 67 20 5f 20 3d 20 4c .long long _ = L
0bf0: 4c 3b 20 0d 0a 45 4e 44 0d 0a 2a 2f 0d 0a 7d 0d L; ..END..*/..}.
0c00: 0a 2f 2f 20 45 4e 44 20 43 55 54 20 48 45 52 45 .// END CUT HERE
0c10: 0d 0a ..