# include <cstdlib>
# include <iostream>
# include <omp.h>
using namespace std;
int main ( int argc, char *argv[] )
{
  int i, j, total, n = 1000, total = 0;
  bool prime;

# pragma omp parallel \
  shared ( n ) \
  private ( i, prime, j ) \
  reduction ( + : total )

# pragma omp for
  for ( i = 2; i <= n; i++ )
  {
    prime = true;
    for ( j = 2; j < i; j++ )
    {
      if ( i % j == 0 )
      {
        prime = false;
        break;
      }
    }
    if ( prime ) { total = total + i; }
  }
  cout << "PRIME_SUM(2:" << n << ") = " << total << "\n";
  return 0;
}