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

#define SIZE 100000000
int A[SIZE];

int main()
{
	int sum = 0, total;
        int lsum[100];
        for (int i = 0; i<100; i++) lsum[i] = 0;

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

        #pragma omp parallel for 
	for (int i=0; i<SIZE; i++)
		A[i] = 1;
	
	#pragma omp parallel
	{
                int myid = omp_get_thread_num();
                #pragma omp for
		for (int i=0; i<SIZE; i++) 
			lsum[myid] = lsum[myid] + A[i];
        }

        for (int i = 0; i<total; i++) sum += lsum[i];

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