#include <stdio.h>
#include <assert.h>
#include <math.h>
typedef unsigned long int L;
typedef struct{L q; L r;} res;
res divq(res, L);
res mulq(L, L);
res addq(res, res);
#define sz 4000
#define T 20000
L e[sz];

void div(int r, int j){int k; res d = {0, r};
  for(k=0; k<sz; ++k) {d = divq((res){d.r, e[k]}, j); e[k] = d.q;}}

int main(){
float dt = (T*log(T) - T + log(2.*3.14159265*T)/2)/log(10),
      da = 64*log(2)*sz/log(10);
  printf("Enough terms for %5.1f digits.\n"
         "Enough array for %5.1f digits.\n", dt, da);
  {int j=sz; while(j--) e[j]=0;}
  {int j=T; while(j--) div(1, j+2);}
  if(0) {int j; for(j=0; j<sz; ++j) printf("%016lX", e[j]);} // Print hex
  if(1) {printf("2."); int j=(dt<da?dt:da)/13; res p = {0, 0};
    while(j--) {int k=sz; while(k--) {
       p = addq(mulq(e[k], 10000000000000), (res){0, p.q}); e[k] = p.r;}
     printf("%013ld", p.q);} printf("\n");} // Print decimal.
return 0;} 

// a68g -p "long long exp(1)" -prec 7000