#include <stdio.h>
#include <malloc.h>
#include <omp.h>

long long factorial(long n)
{
   long long i,out;
   out = 1;
   for (i=1; i<n+1; i++) out *= i;
   return(out);
}

int main(int argc, char **argv)
{
  int i,j,threads;
  long long *x;
  long long n=12;

  /* Set number of threads equal to argv[1] if present */
  if (argc > 1)
    {
      threads = atoi(argv[1]);
      if (omp_get_dynamic())
        {
          omp_set_dynamic(0);
          printf("called omp_set_dynamic(0)\n");
        }
      omp_set_num_threads(threads);
    }
  printf("%d threads\n",omp_get_max_threads());

  x = (long long *) malloc(n * sizeof(long));
  for (i=0;i<n;i++) x[i]=factorial(i);
  j=0;
  /* Is the output the same if the following line is commented out? */
  #pragma omp parallel for firstprivate(x,j)
  for (i=1; i<n; i++)
    {
      j += i;
      x[i] = j*x[i-1];
    }
  for (i=0; i<n; i++)
    printf("factorial(%2d)=%14lld   x[%2d]=%14lld\n",i,factorial(i),i,x[i]);
  return 0;
}