Hex Artifact Content
Not logged in

Artifact d2930c13625ebb9212dca232b7282a48d8458bb6:


0000: 0a 2f 2f 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d  .//-------------
0010: 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d  ----------------
0020: 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d  ----------------
0030: 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d  ----------------
0040: 0a 2f 2f 20 28 4d 69 6e 69 6d 75 6d 29 20 53 65  .// (Minimum) Se
0050: 74 20 43 6f 76 65 72 0a 2f 2f 20 20 20 4e 50 2d  t Cover.//   NP-
0060: 63 6f 6d 70 6c 65 74 65 0a 2f 2f 0a 2f 2f 20 41  complete.//.// A
0070: 73 73 75 6d 70 74 69 6f 6e 3a 20 22 6d 61 73 6b  ssumption: "mask
0080: 22 20 6d 75 73 74 20 62 65 20 64 6f 77 6e 77 61  " must be downwa
0090: 72 64 20 63 6c 6f 73 65 72 64 2c 0a 2f 2f 20 20  rd closerd,.//  
00a0: 20 69 2e 65 2e 2c 20 58 20 73 75 62 73 65 74 20   i.e., X subset 
00b0: 59 20 26 20 59 20 69 6e 20 6d 61 73 6b 20 3d 3e  Y & Y in mask =>
00c0: 20 58 20 69 6e 20 6d 61 73 6b 0a 2f 2f 0a 2f 2f   X in mask.//.//
00d0: 20 56 65 72 69 66 69 65 64 20 62 79 0a 2f 2f 20   Verified by.// 
00e0: 20 20 2d 20 47 6f 6f 67 6c 65 20 43 6f 64 65 20    - Google Code 
00f0: 4a 61 6d 20 30 39 20 52 6f 75 6e 64 32 20 43 0a  Jam 09 Round2 C.
0100: 2f 2f 20 20 20 20 20 28 54 68 65 20 22 4f 6e 6c  //     (The "Onl
0110: 79 22 20 70 61 72 74 20 6e 65 76 65 72 20 74 65  y" part never te
0120: 73 74 65 64 2c 20 69 6e 20 74 68 65 20 74 65 73  sted, in the tes
0130: 74 20 63 61 73 65 73 20 6f 66 0a 2f 2f 20 20 20  t cases of.//   
0140: 20 20 20 74 68 65 20 70 72 6f 62 6c 65 6d 2c 20     the problem, 
0150: 74 68 61 74 20 6f 70 74 69 6d 69 7a 61 74 69 6f  that optimizatio
0160: 6e 20 68 61 64 20 6e 65 76 65 72 20 75 73 65 64  n had never used
0170: 2e 0a 2f 2f 20 20 20 20 20 20 42 65 20 63 61 72  ..//      Be car
0180: 65 66 75 6c 20 74 6f 20 75 73 65 20 74 68 61 74  eful to use that
0190: 3b 20 74 6f 20 62 65 20 6f 6e 20 74 68 65 20 73  ; to be on the s
01a0: 61 66 65 72 20 73 69 64 65 2c 0a 2f 2f 20 20 20  afer side,.//   
01b0: 20 20 20 63 6f 6e 73 69 64 65 72 20 75 73 69 6e     consider usin
01c0: 67 20 74 68 65 20 3c 74 72 75 65 2c 66 61 6c 73  g the <true,fals
01d0: 65 3e 20 6d 6f 64 65 2e 29 0a 2f 2f 2d 2d 2d 2d  e> mode.).//----
01e0: 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d  ----------------
01f0: 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d  ----------------
0200: 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d  ----------------
0210: 2d 2d 2d 2d 2d 2d 2d 2d 2d 0a 0a 69 6e 74 20 6d  ---------..int m
0220: 69 6e 43 6f 76 65 72 5f 65 78 68 61 75 73 74 69  inCover_exhausti
0230: 76 65 28 20 69 6e 74 20 67 6f 61 6c 2c 20 63 6f  ve( int goal, co
0240: 6e 73 74 20 73 65 74 3c 69 6e 74 3e 26 20 6d 61  nst set<int>& ma
0250: 73 6b 20 29 0a 7b 0a 09 74 79 70 65 64 65 66 20  sk ).{..typedef 
0260: 69 6e 74 20 53 54 41 54 45 3b 0a 0a 09 76 65 63  int STATE;...vec
0270: 74 6f 72 3c 53 54 41 54 45 3e 20 51 28 20 31 2c  tor<STATE> Q( 1,
0280: 20 30 20 29 3b 0a 09 73 65 74 3c 53 54 41 54 45   0 );..set<STATE
0290: 3e 20 76 69 73 69 74 65 64 3b 0a 09 66 6f 72 28  > visited;..for(
02a0: 69 6e 74 20 73 74 65 70 3d 30 3b 20 73 74 65 70  int step=0; step
02b0: 3c 3d 6d 61 73 6b 2e 73 69 7a 65 28 29 3b 20 2b  <=mask.size(); +
02c0: 2b 73 74 65 70 29 0a 09 7b 0a 09 09 76 65 63 74  +step)..{...vect
02d0: 6f 72 3c 53 54 41 54 45 3e 20 51 32 3b 0a 09 09  or<STATE> Q2;...
02e0: 66 6f 72 28 69 6e 74 20 69 3d 30 3b 20 69 3c 51  for(int i=0; i<Q
02f0: 2e 73 69 7a 65 28 29 3b 20 2b 2b 69 29 0a 09 09  .size(); ++i)...
0300: 7b 0a 09 09 09 53 54 41 54 45 20 73 20 3d 20 51  {....STATE s = Q
0310: 5b 69 5d 3b 0a 09 09 09 69 66 28 20 73 20 3d 3d  [i];....if( s ==
0320: 20 67 6f 61 6c 20 29 0a 09 09 09 09 72 65 74 75   goal ).....retu
0330: 72 6e 20 73 74 65 70 3b 0a 09 09 09 66 6f 72 28  rn step;....for(
0340: 73 65 74 3c 69 6e 74 3e 3a 3a 63 6f 6e 73 74 5f  set<int>::const_
0350: 69 74 65 72 61 74 6f 72 20 69 74 3d 6d 61 73 6b  iterator it=mask
0360: 2e 62 65 67 69 6e 28 29 3b 20 69 74 21 3d 6d 61  .begin(); it!=ma
0370: 73 6b 2e 65 6e 64 28 29 3b 20 2b 2b 69 74 29 0a  sk.end(); ++it).
0380: 09 09 09 09 69 66 28 20 21 76 69 73 69 74 65 64  ....if( !visited
0390: 2e 63 6f 75 6e 74 28 20 73 7c 2a 69 74 20 29 20  .count( s|*it ) 
03a0: 29 0a 09 09 09 09 7b 0a 09 09 09 09 09 76 69 73  ).....{......vis
03b0: 69 74 65 64 2e 69 6e 73 65 72 74 28 20 73 7c 2a  ited.insert( s|*
03c0: 69 74 20 29 3b 0a 09 09 09 09 09 51 32 2e 70 75  it );......Q2.pu
03d0: 73 68 5f 62 61 63 6b 28 20 73 7c 2a 69 74 20 29  sh_back( s|*it )
03e0: 3b 0a 09 09 09 09 7d 0a 09 09 7d 0a 09 09 51 2e  ;.....}...}...Q.
03f0: 73 77 61 70 28 51 32 29 3b 0a 09 7d 0a 09 72 65  swap(Q2);..}..re
0400: 74 75 72 6e 20 2d 31 3b 0a 7d 0a 0a 74 65 6d 70  turn -1;.}..temp
0410: 6c 61 74 65 3c 62 6f 6f 6c 20 53 75 62 73 2c 20  late<bool Subs, 
0420: 62 6f 6f 6c 20 4f 6e 6c 79 3e 0a 69 6e 74 20 6d  bool Only>.int m
0430: 69 6e 43 6f 76 65 72 28 20 69 6e 74 20 67 6f 61  inCover( int goa
0440: 6c 2c 20 73 65 74 3c 69 6e 74 3e 20 6d 61 73 6b  l, set<int> mask
0450: 20 29 0a 7b 0a 09 69 66 28 20 53 75 62 73 20 29   ).{..if( Subs )
0460: 0a 09 09 66 6f 72 28 73 65 74 3c 69 6e 74 3e 3a  ...for(set<int>:
0470: 3a 69 74 65 72 61 74 6f 72 20 69 74 3d 6d 61 73  :iterator it=mas
0480: 6b 2e 62 65 67 69 6e 28 29 3b 20 69 74 21 3d 6d  k.begin(); it!=m
0490: 61 73 6b 2e 65 6e 64 28 29 3b 20 29 0a 09 09 7b  ask.end(); )...{
04a0: 0a 09 09 09 62 6f 6f 6c 20 6e 65 65 64 65 64 20  ....bool needed 
04b0: 3d 20 74 72 75 65 3b 0a 09 09 09 66 6f 72 28 69  = true;....for(i
04c0: 6e 74 20 6d 3d 31 3b 20 6d 3c 3d 67 6f 61 6c 3b  nt m=1; m<=goal;
04d0: 20 6d 3c 3c 3d 31 29 0a 09 09 09 09 69 66 28 20   m<<=1).....if( 
04e0: 28 67 6f 61 6c 26 6d 29 20 26 26 20 21 28 2a 69  (goal&m) && !(*i
04f0: 74 26 6d 29 20 26 26 20 6d 61 73 6b 2e 63 6f 75  t&m) && mask.cou
0500: 6e 74 28 2a 69 74 7c 6d 29 29 0a 09 09 09 09 09  nt(*it|m))......
0510: 7b 6e 65 65 64 65 64 20 3d 20 66 61 6c 73 65 3b  {needed = false;
0520: 20 62 72 65 61 6b 3b 7d 0a 0a 09 09 09 69 66 28   break;}.....if(
0530: 20 6e 65 65 64 65 64 20 29 0a 09 09 09 09 2b 2b   needed ).....++
0540: 69 74 3b 0a 09 09 09 65 6c 73 65 0a 09 09 09 09  it;....else.....
0550: 6d 61 73 6b 2e 65 72 61 73 65 28 69 74 2b 2b 29  mask.erase(it++)
0560: 3b 0a 09 09 7d 0a 0a 09 69 66 28 20 4f 6e 6c 79  ;...}...if( Only
0570: 20 29 0a 09 09 66 6f 72 28 73 65 74 3c 69 6e 74   )...for(set<int
0580: 3e 3a 3a 69 74 65 72 61 74 6f 72 20 69 74 3d 6d  >::iterator it=m
0590: 61 73 6b 2e 62 65 67 69 6e 28 29 3b 20 69 74 21  ask.begin(); it!
05a0: 3d 6d 61 73 6b 2e 65 6e 64 28 29 3b 20 29 0a 09  =mask.end(); )..
05b0: 09 7b 0a 09 09 09 69 6e 74 20 6f 6e 6c 79 42 79  .{....int onlyBy
05c0: 4d 65 20 3d 20 67 6f 61 6c 20 26 20 2a 69 74 3b  Me = goal & *it;
05d0: 0a 09 09 09 66 6f 72 28 73 65 74 3c 69 6e 74 3e  ....for(set<int>
05e0: 3a 3a 69 74 65 72 61 74 6f 72 20 6a 74 3d 6d 61  ::iterator jt=ma
05f0: 73 6b 2e 62 65 67 69 6e 28 29 3b 20 6a 74 21 3d  sk.begin(); jt!=
0600: 6d 61 73 6b 2e 65 6e 64 28 29 3b 20 2b 2b 6a 74  mask.end(); ++jt
0610: 29 0a 09 09 09 09 6f 6e 6c 79 42 79 4d 65 20 26  ).....onlyByMe &
0620: 3d 20 7e 2a 6a 74 3b 0a 09 0a 09 09 09 69 66 28  = ~*jt;......if(
0630: 20 21 6f 6e 6c 79 42 79 4d 65 20 29 0a 09 09 09   !onlyByMe )....
0640: 09 2b 2b 69 74 3b 0a 09 09 09 65 6c 73 65 0a 09  .++it;....else..
0650: 09 09 09 6d 61 73 6b 2e 65 72 61 73 65 28 69 74  ...mask.erase(it
0660: 2b 2b 29 2c 20 67 6f 61 6c 20 26 3d 20 7e 6f 6e  ++), goal &= ~on
0670: 6c 79 42 79 4d 65 3b 0a 09 09 7d 0a 0a 09 72 65  lyByMe;...}...re
0680: 74 75 72 6e 20 6d 69 6e 43 6f 76 65 72 5f 65 78  turn minCover_ex
0690: 68 61 75 73 74 69 76 65 28 67 6f 61 6c 2c 20 6d  haustive(goal, m
06a0: 61 73 6b 29 3b 0a 7d 0a                          ask);.}.