Calculate factorials using gmp

fact.c

/*
 * $ cc -O3 -Wall -pedantic -Wextra fact.c -lgmp
 */

#include <stdio.h>
#include <gmp.h>

static void
factorial(long n, mpz_t r)
{ 
	mpz_init_set_si(r, 1);
	for (; n > 1; n--) {
		mpz_mul_si(r, r, n);
	}
}

int
main (void)
{
	int i;
	mpz_t r;

	for (i=0; i<14888; i++) {
		factorial(i, r);
		gmp_printf ("%Zd\n", r);
	}

	return 0;
}