Вычисление коэффициентов, которые позволяют отобрать из 20000 матриц всего лишь от 800 до 1600 матриц

#include <stdio.h>                
#include <locale.h>
#include <ctype.h>
#include <string.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
int main()
{
    //srand(time(NULL)); 
    int l=20000,k=34; 
    int c=2000;                             
unsigned int timer;
int p=1,n,i1,i2,i3,q;
int flag,flag1,flag2=0; 
double dd,D0,D; 
double a;
for (n=2; n<=100; n++){
int H[k][n];
for (a=0.90; a>=0.5; a=a-0.00015){
p=1;
D0=a*sqrt(k*n);
    int ***M=(int ***)malloc(c * sizeof(*M));
    for (i2=0; i2<c; i2++)
    {
        M[i2]=(int **)malloc(k * sizeof(*M[i2]));
        for (i3=0; i3<k; i3++)
        {
        M[i2][i3]=(int *)malloc(n * sizeof(*M[i2][i3]));
        }
    }
    
      for (i2=0; i2<k; i2++)   
      {
        for (i3=0; i3<n; i3++)
        {
         M[0][i2][i3]==rand()%2;         
        }
      }
 
for (i1=1; i1<l; i1++)       
{
   for (i2=0; i2<k; i2++)  
      {
        for (i3=0; i3<n; i3++)
        {
        H[i2][i3]=rand()%2;    
    }
  }
    
flag=999;
flag1=0;
 
  for (q=0; q<p; q++)
  { 
      dd=0; 
   for (i2=0; i2<k; i2++)  
      { 
        for (i3=0; i3<n; i3++)
       {
        dd=dd+pow((M[q][i2][i3]-H[i2][i3]),2);
       }
       D=sqrt(dd);
        if (D>=D0) break;           
      }
                    
  D=sqrt(dd);  
  if (D<D0)  
    { 
    flag=0;
    break;
    }      
  }  
      
  if ((flag==999))       
  {
    for (i2=0; i2<k; i2++)  
      {
        for (i3=0; i3<n; i3++)
       {
      M[p][i2][i3]=H[i2][i3];
       }
            
      }
            
  p++;                   
  } 
if (p>1600) break; 
}
 
 
for (i2=0; i2<c; i2++)
{
for (i3=0; i3<k; i3++)
{
 free(M[i2][i3]);
}  
free(M[i2]);
}
free(M);
 
if ((p>=800)&&(p<=1600)) {flag2=999;break;} }
if(flag2==999){
printf("n = %d\n",n);
printf("a = %.5f\n",a);   
printf("number of matrices = %d\n",p);
printf("\n");} flag2=0;}
 
 
timer=(clock() / CLOCKS_PER_SEC); 
printf("time = %d\n",timer);
return 0;
}