# include <cstdlib>
# include <iostream>
# include <omp.h>

using namespace std;

int main ( int argc, char *argv[] )
{
  int i, n = 10000000;
  double s, *x, *y;

  s = 123.456;

  x = new double[n];
  y = new double[n];

  for ( i = 0; i < n; i++ )
  {
    x[i] = ( double ) ( ( i + 1 ) % 17 );
    y[i] = ( double ) ( ( i + 1 ) % 31 );
  }

# pragma omp parallel

# pragma omp for
  for ( i = 0; i < n; i++ )
  {
    y[i] = y[i] + s * x[i];
  }

  delete [] x;
  delete [] y;

  return 0;
}