Алгоритм BBS(Блюма — Блюма — Шубa)

int x0=56;
 
int bbs()
{
 int M=111;
 static int x=x0;
 
 x=x*x % M;
 
 x0=x;
 return x;
}

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <iostream>
/////////////////
int x0 = 56;
int M = 2111;
int n = 4;
/////////////////
int *bits;
int k=0;
int c=0;
 
int bbs();
int * dec2bin(int);
int check_bits(int *);
void generate(int,int);
 
int main()
{
    using std::cout;
    using std::cin;
    using std::endl;
 
    setlocale(0,"");
 
     int numbers = 0;
 
     cout<<"Ââåäiòü êiëüêiñòü ÷èñåë: ";
     cin>>numbers;
     cout<<"Êiëüêiñòü áiò: ";
     cin>>n;
 
     generate(numbers,n);
    
 
 
     return 0;
}
 
int bbs()
{
    static int x=x0;
 
    x=x*x%M;
    c = x;
    x0 = x;
    
    return check_bits(dec2bin(x));
}
 
int * dec2bin(int z)
{   
    int i=1,j=1;
    int *t;
    t = (int*)calloc(100,sizeof(int));
    
    while(z>1)
    {
    *(t+i) = z%2;
    z=z/2; 
    i++;
    }
    
    k=i;
    bits = (int*)calloc(k,sizeof(int));
 
    int m = 1;
    *(bits+0)=z;
    for(j=i-1;j>0;j--,m++) {*(bits+m)=*(t+j);}
    
    free(t);
 
    return bits;
}
 
int check_bits(int *pbits)
{
    int bit_1 = 0;
    for(int i=0;i<k;i++) if(*(pbits+i) & 1) bit_1++;
 
    if(bit_1 % 2) return 0;
    else return 1;
}
 
void generate(int numbers,int n_bits)
{
    int *g = (int*)calloc(n_bits,sizeof(int));
    int rand = 0;
    int p;
    
    for(int j=0;j<numbers;j++)
    {
      p = n_bits-1;
      for(int i=0;i<n_bits;i++) 
      {
        *(g+i) = bbs();
        rand += *(g+i)*pow(2.0,p--);
      } 
 
      std::cout<<rand<<"; ";
      rand = 0;
    }
}