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

int main() {
	const int MAX_I = 5;
	const int MAX_J = 10;
	const int MAX_THREADS = 4;
	int nbr_processed[MAX_THREADS];
	
	for (int i=0; i<MAX_THREADS; ++i) nbr_processed[i] = 0;
	
	#pragma omp parallel for 
	for (int i=0; i<MAX_I; ++i) {
		for (int j=0; j<MAX_J; ++j) {
			#pragma omp critical
			{
				cout << i << " " << j << " " << omp_get_thread_num() << endl;
				++nbr_processed[omp_get_thread_num() ];
			}
		}
	}
	
	for (int i=0; i<MAX_THREADS; ++i) {
		cout << "thread " << i << " processed " <<  nbr_processed[i] << " elements" << endl;
	}
	
	return EXIT_SUCCESS;
}