#include <omp.h>
#include <iostream>
using namespace std;

#define SIZE 100000000
int A[SIZE];

int main()
{
	int sum = 0;

        #pragma omp parallel
        {
                int total_threads = omp_get_num_threads();
                int myid = omp_get_thread_num();
                if (myid == 0) cout << total_threads << " threads\n";
        }

        #pragma omp parallel for
	for (int i=0; i<SIZE; i++)
		A[i] = 1;
	
	#pragma omp parallel for reduction(+: sum)
	for (int i=0; i<SIZE; i++)
		sum = sum + A[i];

	cout << "sum = " << sum << "\n";
        return 0;
}