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

int main()
{
	int num_threads = 4;
	int sum = 0;
	
	#pragma omp parallel num_threads(num_threads)
	{
		int thread_id = omp_get_thread_num();
		int thread_sum = 0;
		
		// Perform some computation
		for (int i = thread_id; i <= 100; i += num_threads)
		{
			thread_sum += i;
		}
		
		// Use the "pv" primitive to update the shared sum variable
		#pragma omp critical
		{
			sum += thread_sum;
		}
	}
	
	std::cout << "Sum: " << sum << std::endl;
	
	return 0;
}