#include "stdio.h"




int ** tomjc( short *S, int n, int *F, int lenF, int k, int *lenoutput )
{
  int i,j,l,r;
  int lenS = n*k,nonempty=0,max=0,min=n*k,bigjump;
  int *allow,lenAllow,subset,*F2,*tempF;
  short *S2;
  int **output, **oldoutput, lenoldoutput;

  int temp;

  for ( i=0; i<lenS; i++ ) {
    if ( S[i]==1 ) {
      if ( i>max ) {
	max = i;
      }
      if ( i<min ) {
	min = i;
      }
      nonempty = 1;
    }
  }
  if ( !nonempty ) {
    output = (int**) malloc( 2* sizeof(int*) );
    *lenoutput = 1;
    output[0] = (int*) malloc( 2*n* sizeof(int) );
    for (i=0; i<2*n; i++) {
      output[0][i] = -1;
    }
    output[1] = NULL;
    return output;
  }

  bigjump = (int) ((max-min)/(k-1));

  tempF = (int*) malloc(bigjump* sizeof(int));
  allow = (int*) malloc(bigjump* sizeof(int));
  for ( i=1; i<=bigjump; i++ ) {
    tempF[i] = 0;
  }
  lenAllow = bigjump;
  for ( i=0; i<lenF; i++ ) {
    if ( F[i] <= bigjump ) {
      tempF[F[i]] = 1;
      lenAllow--;
    }
  }
  j=0;
  for ( i=1; i<=bigjump; i++ ) {
    if ( tempF[i] == 0 ) {
      allow[j] = i;
      j++;
    }
  }
  
  output = (int**) malloc( sizeof(int*) );
  *lenoutput = 0;
  output[0] = NULL;

  for ( i=0; i<lenAllow; i++ ) {
    r = allow[i];
    subset = 1;
    for ( j=0; j<k; j++ ) {
      if( S[min+j*r]==0 ) {
	subset = 0;
	break;
      }
    }
    if ( subset ) {
      S2 = (short*) malloc(lenS*sizeof(short));
      for ( j=0; j<lenS; j++ ) {
	S2[j] = S[j];
      }
      for ( j=0; j<k; j++ ) {
	S2[min+j*r] = 0;
      }

      F2 = (int*) malloc((lenF+1)*sizeof(int));
      for ( j=0; j<lenF; j++ ) {
	F2[j] = F[j];
      }
      F2[j] = r;
      
      oldoutput = tomjc( S2, n, F2, lenF+1, k, &lenoldoutput );

      temp = *lenoutput+lenoldoutput;
      
      /*
      free(S2);
      free(F2);
      */
      realloc( output, (*lenoutput+lenoldoutput+1)*sizeof(int*) );
      
      for ( j=0; j<lenoldoutput; j++ ) {
	for ( l=0; l<2*n; l++ ) {
	  if ( oldoutput[j][l]==-1 ) {
	    oldoutput[j][l] = min;
	    oldoutput[j][l+1] = r;
	    break;
	  }
	}
	output[*lenoutput+j] = oldoutput[j];
      }
      *lenoutput += lenoldoutput;
      output[*lenoutput] = NULL;
    }
  }

  /*
  free(allow);
  free(tempF);
  */
  
  return output;
    
      
}





int main( void )
{

  int n=3, k=2, l=0;
  short *S;
  int i,j,*F,lenS,lenF,**output;

  S = (short*) malloc( n*k*sizeof(short) );
  for ( i=0; i<n*k; i++ ) {
    S[i]=1;
  }
  lenS = n*k;
  
  lenF = 0;


  
  output = tomjc( S, n, F,lenF, k, &l );
  printf("%d\n",l);
  
  for ( i=0; i<l; i++ ) {
    for ( j=0; j<2*n; j++ ) {
      printf("%d,",output[i][j]);
    }
    printf("\n");
  }
    
  return 0;
}