#define W 12 
struct Box
{
unsigned r0 :  W ;
unsigned r1 :  W ;
unsigned r2 :  W ;
unsigned r3 :  W ;
unsigned r4 :  W ;
unsigned r5 :  W ;
unsigned r6 :  W ;
unsigned r7 :  W ;
unsigned r8 :  W ;
unsigned r9 :  W ;
unsigned r10 :  W ;
unsigned r11 :  W ;
};

#define H 12

#define smallW 3
#define bigW   4

#define littleOptions ((W-smallW+1)*(H-smallW+1))
#define bigOptions    ((W-bigW+1)*(H-bigW+1))

#define littleSection 7
// the bit string defining a little ball, in section
#define bigSection 15
// the bit string defining a big ball, in section

#include <iostream>
using namespace std;

Box little[littleOptions] ;  
Box big[bigOptions] ;  

void bitsout( unsigned int n );
void boxout( Box bj ) ;

int main()
{
  cout << "bigOptions = " << bigOptions << endl;
  cout << "littleOptions = " << littleOptions << endl;
  for ( int i = 0 ; i < littleOptions ; i ++ ) {
    little[i].r0=0;    little[i].r1=0;    little[i].r2=0;    little[i].r3=0;
    little[i].r4=0;    little[i].r5=0;    little[i].r6=0;    little[i].r7=0;
    little[i].r8=0;    little[i].r9=0;    little[i].r10=0;    little[i].r11=0;
    int ix = i % (W-smallW+1) ; 
    int iy = i / (W-smallW+1) ;
    if ( iy==0)little[i].r2 =little[i].r1 = little[i].r0 = littleSection << ix ;
    if ( iy==1)little[i].r3 =little[i].r2 = little[i].r1 = littleSection << ix ;
    if ( iy==2)little[i].r4 =little[i].r3 = little[i].r2 = littleSection << ix ;
    if ( iy==3)little[i].r5 =little[i].r4 = little[i].r3 = littleSection << ix ;
    if ( iy==4)little[i].r6 =little[i].r5 = little[i].r4 = littleSection << ix ;
    if ( iy==5)little[i].r7 =little[i].r6 = little[i].r5 = littleSection << ix ;
    if ( iy==6)little[i].r8 =little[i].r7 = little[i].r6 = littleSection << ix ;
    if ( iy==7)little[i].r9 =little[i].r8 = little[i].r7 = littleSection << ix ;
    if ( iy==8)little[i].r10=little[i].r9 = little[i].r8 = littleSection << ix ;
    if ( iy==9)little[i].r11=little[i].r10=little[i].r9  = littleSection << ix ;
  }
  for ( int i = 0 ; i < bigOptions ; i ++ ) {
    big[i].r0=0;    big[i].r1=0;    big[i].r2=0;    big[i].r3=0;
    big[i].r4=0;    big[i].r5=0;    big[i].r6=0;    big[i].r7=0;
    int ix = i % (W-bigW+1) ; 
    int iy = i / (W-bigW+1) ;
    if ( iy==0 ) big[i].r3 =big[i].r2 =big[i].r1 = big[i].r0=bigSection << ix ;
    if ( iy==1 ) big[i].r4 =big[i].r3 =big[i].r2 = big[i].r1=bigSection << ix ;
    if ( iy==2 ) big[i].r5 =big[i].r4 =big[i].r3 = big[i].r2=bigSection << ix ;
    if ( iy==3 ) big[i].r6 =big[i].r5 =big[i].r4 = big[i].r3=bigSection << ix ;
    if ( iy==4 ) big[i].r7 =big[i].r6 =big[i].r5 = big[i].r4=bigSection << ix ;
    if ( iy==5 ) big[i].r8 =big[i].r7 =big[i].r6 = big[i].r5=bigSection << ix ;
    if ( iy==6 ) big[i].r9 =big[i].r8 =big[i].r7 = big[i].r6=bigSection << ix ;
    if ( iy==7 ) big[i].r10=big[i].r9 =big[i].r8 = big[i].r7=bigSection << ix ;
    if ( iy==8 ) big[i].r11=big[i].r10=big[i].r9 = big[i].r8=bigSection << ix ;
  }
  if(1){
    for ( int j = 0 ; j < bigOptions/10 ; j ++ ) {
      cout << j << " th state of big dude" << endl;
      boxout(big[j]);
    }
    for ( int i = 0 ; i < littleOptions/10 ; i ++ ) {
      cout << i << " th state of little dude " << endl;
      boxout( little[i] );
    }
  }
  cout << "Single small-particle states" <<endl;
  for ( int j = 0 ; j < bigOptions ; j ++ ) {
    Box bj = big[j];
    int ok = 0 ; 
    for ( int i = 0 ; i < littleOptions ; i ++ ) {
      Box bi = little[i];
      if ( (bi.r0 & bj.r0) ||  (bi.r1 & bj.r1) ||
	    (bi.r2 & bj.r2) ||	    (bi.r3 & bj.r3) ||
	    (bi.r4 & bj.r4) ||	    (bi.r5 & bj.r5) ||
	   (bi.r6 & bj.r6) ||  (bi.r7 & bj.r7)      ||
	   (bi.r8 & bj.r8) ||  (bi.r9 & bj.r9)      ||
	   (bi.r10& bj.r10) ||  (bi.r11 & bj.r11) ) {
	// collision
      } else {
	ok ++ ;
      }
    }
    cout << " " << ok ;
    if ( !((j+1)%(W-bigW+1)) ) cout << endl ;
  }
  cout << "Two small-particle states" <<endl;
  for ( int j = 0 ; j < bigOptions ; j ++ ) {
    Box bj = big[j];
    int ok = 0 ; 
    for ( int i = 0 ; i < littleOptions ; i ++ ) {
      Box bi = little[i];
      if ( (bi.r0 & bj.r0) ||  (bi.r1 & bj.r1) ||
	    (bi.r2 & bj.r2) ||	    (bi.r3 & bj.r3) ||
	    (bi.r4 & bj.r4) ||	    (bi.r5 & bj.r5) ||
	   (bi.r6 & bj.r6) ||  (bi.r7 & bj.r7)   ||
	   (bi.r8 & bj.r8) ||  (bi.r9 & bj.r9)      ||
	   (bi.r10& bj.r10) ||  (bi.r11 & bj.r11)  ) {
	// collision
      } else {
	for ( int ii = 0 ; ii < littleOptions ; ii ++ ) {
	  Box bk = little[ii];
	  if ( (bk.r0 & bj.r0) ||  (bk.r1 & bj.r1) ||
	       (bk.r2 & bj.r2) ||	    (bk.r3 & bj.r3) ||
	       (bk.r4 & bj.r4) ||	    (bk.r5 & bj.r5) ||
	       (bk.r6 & bj.r6) ||  (bk.r7 & bj.r7)   ||
	   (bk.r8 & bj.r8) ||  (bk.r9 & bj.r9)      ||
	   (bk.r10& bj.r10) ||  (bk.r11 & bj.r11) || 
	       (bi.r0 & bk.r0) ||  (bi.r1 & bk.r1) ||
	       (bi.r2 & bk.r2) ||	    (bi.r3 & bk.r3) ||
	       (bi.r4 & bk.r4) ||	    (bi.r5 & bk.r5) ||
	       (bi.r6 & bk.r6) ||  (bi.r7 & bk.r7)   ||
	       (bi.r8 & bk.r8) ||  (bi.r9 & bk.r9)      ||
	       (bi.r10& bk.r10) ||  (bi.r11 & bk.r11)  ) {
	    // collision
	  } else {
	    ok ++ ;
	  }
	}
      }
    }
    cout << " " << ok ;
    if ( !((j+1)%(W-bigW+1)) ) cout << endl ;
  }

  return(0);
}

void boxout( Box bj ) {
    bitsout ( bj.r0 ) ;cout<<endl;
    bitsout ( bj.r1 ) ;cout<<endl;
    bitsout ( bj.r2 ) ;cout<<endl;
    bitsout ( bj.r3 ) ;cout<<endl;
    bitsout ( bj.r4 ) ;cout<<endl;
    bitsout ( bj.r5 ) ;cout<<endl;
    bitsout ( bj.r6 ) ;cout<<endl;
    bitsout ( bj.r7 ) ;cout<<endl;
    bitsout ( bj.r8 ) ;cout<<endl;
    bitsout ( bj.r9 ) ;cout<<endl;
    bitsout ( bj.r10) ;cout<<endl;
    bitsout ( bj.r11) ;cout<<endl;
}

void bitsout( unsigned int m )
{
  int lastbit  ; 
  unsigned int two_to_power_i ;
  
  for ( int i = W-1 ; i >= 0 ; i -- ) {
    two_to_power_i = (1<<i) ;
    lastbit = ( two_to_power_i & m ) == 0 ? 0 : 1 ;
    cout << lastbit ;
  }
}