Cod sursa(job #1157323)

Utilizator andrettiAndretti Naiden andretti Data 28 martie 2014 13:44:12
Problema Plantatie Scor 50
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.42 kb
#include<stdio.h>
#include<algorithm>
#define maxn 505
#define maxlg 25
using namespace std;

int n,m;
int a[maxn][maxn];
int rmq[maxn][maxn][maxlg],lg[maxn],pw[maxlg];

void read()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
     for(int j=1;j<=n;j++)
      scanf("%d",&a[i][j]);
}

void RMQ()
{
    for(int i=1;i<=n;i++)
     for(int j=1;j<=n;j++)
      rmq[i][j][0]=a[i][j];

    for(int i=2;i<=n;i++) lg[i]=lg[i/2]+1;
    pw[0]=1; for(int i=1;i<=lg[n]+2;i++) pw[i]=pw[i-1]*2;

    int cnt;
    for(int k=1;pw[k]<=n;k++)
     for(int i=1;i+pw[k]-1<=n;i++)
      for(int j=1;j+pw[k]-1<=n;j++)
      {
          cnt=pw[k-1];
          rmq[i][j][k]=max(rmq[i][j][k-1],rmq[i][j+cnt][k-1]);
          rmq[i][j][k]=max(rmq[i][j][k],rmq[i+cnt][j][k-1]);
          rmq[i][j][k]=max(rmq[i][j][k],rmq[i+cnt][j+cnt][k-1]);
      }
}

int det(int x,int y,int k)
{
    int sol,cnt;
    int p=lg[k]; cnt=pw[p];
    sol=max(rmq[x][y][p],rmq[x][y+k-cnt][p]);
    sol=max(sol,rmq[x+k-cnt][y][p]);
    sol=max(sol,rmq[x+k-cnt][y+k-cnt][p]);
    return sol;
}

void solve()
{
    int x,y,z;
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d",&x,&y,&z);
        printf("%d\n",det(x,y,z));
    }
}

int main()
{
    freopen("plantatie.in","r",stdin);
    freopen("plantatie.out","w",stdout);

    read();
    RMQ();
    solve();

    fclose(stdin);
    fclose(stdout);
    return 0;
}