Cod sursa(job #514269)

Utilizator mihaipopa12Popa Mihai mihaipopa12 Data 18 decembrie 2010 11:44:05
Problema Plantatie Scor 10
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.43 kb
#include<stdio.h>

FILE*f=fopen("plantatie.in","r");
FILE*g=fopen("plantatie.out","w");

int N,Q,i,j,A[512][512],p2[12],M[512][512][11],E[15],R,L,e,k;

int max(int a,int b){
	if ( a > b )
		return a;
	return b;
}

int main () {
	
	fscanf(f,"%d %d",&N,&Q);
	
	for ( i = 1 ; i <= N ; ++i ){
		for ( j = 1 ; j <= N ; ++j ){
			fscanf(f,"%d",&A[i][j]);
			M[i][j][0] = A[i][j];
		}
		E[ i + 1 ] = E[ ( i + 1 ) / 2 ] + 1;
	}
	for ( i = 0 ; i <= N ; ++i ){
		for ( j = 0 ; j <= N ; ++j )
			for ( k = 1 ; k <= 10 ; ++k )
				M[i][j][k] = -1;
			
	}
	
	//M[i][j][k] = maximul din patratul cu coltul stanga-sus in i si j si latura 2 ^ k
	
	for ( i = 1, p2[0] = 1 ; p2[i-1] <= 512 ; p2[i] = 2 * p2[i-1] , ++i ) {};
	
	for ( k = 1 ; k <= E[N] ; ++k ){
		for ( i = 1 ; i <= N ; ++i ){
			for ( j = 1 ; j <= N ; ++j ){
				M[i][j][k] = M[i][j][k-1];
				if( i + p2[k-1] <= N )
					M[i][j][k] = max(M[i][j][k],M[i+p2[k]/2][j][k-1]);
				if( j+p2[k-1] <= N )
					M[i][j][k] = max(M[i][j][k],M[i][j+p2[k-1]][k-1]);
				if ( i + p2[k-1] <= N && j+p2[k-1] <= N )
					M[i][j][k] = max(M[i][j][k],M[i+p2[k-1]][j+p2[k-1]][k-1]);
			}
		}
	}
	
	while ( Q-- ) {
		fscanf(f,"%d %d %d",&i,&j,&L);
		e = E[L];
		R = 0;
		R = max(R,M[i][j][e]); 
		R = max(R,M[i+L-p2[e]][j][e]);
		R = max(R,M[i][j+L-p2[e]][e]);
		R = max(R,M[i+L-p2[e]][j+L-p2[e]][e]);
		fprintf(g,"%d\n",R);
	}
	
	
	fclose(f);
	fclose(g);
	
	return 0;
}