#define W 32 
#define H 32

struct Box
{
  unsigned int r[32] ;
};

#define littleW 3
#define bigW    23

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

#define littleSection (1<<littleW)-1
// the bit string defining a little ball, in horizontal section
#define bigSection (1<<bigW)-1
// the bit string defining a big ball, in horizontal 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;
  cout << "littleSection = " << littleSection << endl;
  for ( int i = 0 ; i < littleOptions ; i ++ ) {
    for ( int h=0; h < H ; h ++ ) {
      little[i].r[h]=0;
    }
    int ix = i % (W-littleW+1) ; 
    int iy = i / (W-littleW+1) ;
    for ( int row = iy ; row < iy +littleW ; row ++ ) {
      little[i].r[row] = littleSection << ix ;
    }
  }
  for ( int i = 0 ; i < bigOptions ; i ++ ) {
    for ( int h=0; h < H ; h ++ ) {
      big[i].r[h]=0;
    }
    int ix = i % (W-bigW+1) ; 
    int iy = i / (W-bigW+1) ;
    for ( int row = iy ; row < iy +bigW ; row ++ ) {
      big[i].r[row] = bigSection << ix ;
    }
  }
  if(0){
    for ( int j = 0 ; j < bigOptions/100 ; j ++ ) {
      cout << j << " th state of big dude" << endl;
      boxout(big[j]);
    }
    for ( int i = 0 ; i < littleOptions/100 ; 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 ;       // counts number of microstates for this j state
    for ( int i = 0 ; i < littleOptions ; i ++ ) {
      Box bi = little[i];
      int ijcompatible=1;
      for ( int h=0; ijcompatible && (h<H) ; h++ ) {
	if ( bi.r[h] & bj.r[h] ) {
	  // collision
	  ijcompatible=0;
	} else {
	  // keep going
	}
      }
      if ( ijcompatible ) 
	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 ;       // counts number of microstates for this j state
    for ( int i = 0 ; i < littleOptions ; i ++ ) {
      Box state = bj ;
      Box bi = little[i];
      int ijcompatible=1;
      for ( int h=0; ijcompatible && (h<H) ; h++ ) {
	if ( bj.r[h] & bi.r[h] ) {
	  // collision
	  ijcompatible=0;
	} else {
	  state.r[h] |= bi.r[h] ;
	  //	  state.r[h] = bj.r[h] | state.r[h] ;
	}
      }
      if ( ijcompatible ) {
	for ( int k = 0 ; k < littleOptions ; k ++ ) {
	  Box bk = little[k];
	  int ijkcompatible=1;
	  for ( int h=0; ijkcompatible && (h<H) ; h++ ) {
	    if ( state.r[h] & bk.r[h] ) {
	      // collision
	      ijkcompatible=0;
	    } else {
	      // keep going    stateBeforek.r[h] |= bk.r[h] ;
	    }
	  }
	  if ( ijkcompatible ) {
	    ok ++ ;
	  }
	}
      }
    }
    cout << " " << ok ;
    if ( !((j+1)%(W-bigW+1)) ) cout << endl ;
  }

  return(0);
}

void boxout( Box bj ) {
  for ( int h = 0 ; h < H ; h ++ ) {
    bitsout ( bj.r[h] ) ;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 ;
  }
}