cp's OEIS Frontend

C

/* To compile: gcc -Wall -O2 A227644.c -o A227644 -lgmp */
#include 
#include 
#include 
int main()
{
   int bsz=256, a=0;
   mpz_t *f, t;
   f = malloc(sizeof(mpz_t) * bsz);
   mpz_init(t); mpz_init(f[0]); mpz_set_ui(f[0], 1);
   while (1)
   {
       a += 1;
       if (a == bsz)
       {
           bsz *= 2;
           f = (mpz_t *) realloc(f, sizeof(mpz_t) * bsz);
       }
       mpz_init(f[a]);
       mpz_mul_ui(f[a], f[a-1], a);
       for (int i=1; i<=a; i++)
       {
           mpz_add(t, f[a], f[i]);
           if (mpz_perfect_power_p(t))
           {
               gmp_printf("%Zd, ", t);
               fflush(stdout);
           }
       }
   }
   return 0;
}