#include <cstdio>
#include <vector>
#include <ctime>
#include <cmath>
using namespace std;
int main ()
{
  unsigned N = 13107200, M = 15, jump = 625681;
  unsigned pos = 0;
  vector<double> a;
  double *pA, x=2.0;
  clock_t start, finish;
  a.resize (N);
  pA = &a[0];
  start = clock();
  for (unsigned i = 0; i < M; i++)
  {
    pos = 0;
    for (unsigned j = 0; j < N; j++)
    {
      pA[pos] = x;
      x += sqrt (x);
      pos = (pos + jump)%N;
    }
  }
  finish = clock ();
  printf ("%d, seconds: %d\n", finish-start, (finish - start) / CLOCKS_PER_SEC);
  return 0;
}