#include <stdio.h>
#include <stdlib.h>
#include <pthread.h>

#define N 5000
int S[N][N];
int A[N][N];
int B[N][N];

void *w (void *arg) {  
	int thid = (int)(long)arg;
	int i, j;
	
	//printf("td:%d\n", thid);
	for(i=thid*N/2; i < thid*N/2 + N/2; i++)	
		for(j=0; j<N; j++)
			S[i][j] = A[i][j] + B[i][j];
			
/* Caso N fosse um numero impar
	if(thid && N%2 )  
		for(j=0; j<N; j++)
			S[N-1][j] = A[N-1][j] + B[N-1][j];
*/			
	return NULL;
}

int main(int argc, char *argv[]){
  pthread_t p1, p2;
  int i, j;

  // incializaçao matrizes
  for(i=0; i < N; i++)	
		for(j=0; j<N; j++) {
			A[i][j] = rand()%100;
			B[i][j] = rand()%100;
		}
		
  pthread_create(&p1, NULL, w, (void *)0);       
  pthread_create(&p2, NULL, w, (void *)1);    
  pthread_join(p1, NULL); pthread_join(p2, NULL); 
/*
  for(i=0; i < N; i++) {
		for(j=0; j<N; j++) 
			printf("S[%d][%d]=%d == %d + %d\n", i, j, S[i][j], A[i][j], B[i][j]);
		printf("\n");
  }
*/
  return 0;
}