Cod sursa(job #1094658)

Utilizator andrettiAndretti Naiden andretti Data 29 ianuarie 2014 18:07:47
Problema Struti Scor 100
Compilator cpp Status done
Runda Teme Pregatire ACM Unibuc 2013 Marime 2.69 kb
#include<stdio.h>
#include<algorithm>
#include<deque>
#define pb push_back
#define popb pop_back
#define popf pop_front
#define mp make_pair
#define x first
#define y second
#define maxn 1005
#define maxb 8192
#define inf 0x3f3f3f3f
using namespace std;

int n,m,p;
int a[maxn][maxn],mx[maxn][maxn],mi[maxn][maxn];
int sol=inf,nr;
deque < pair<int,int> > qmin,qmax;
char buff[maxb];
int pos=0;

int get_int()
{
    int num=0;
    while(buff[pos]<'0' || buff[pos]>'9')
     if(++pos==maxb) fread(buff,1,maxb,stdin),pos=0;

    while(buff[pos]>='0' && buff[pos]<='9')
    {
         num=num*10+buff[pos]-'0';
         if(++pos==maxb) fread(buff,1,maxb,stdin),pos=0;
    }
    return num;
}

void read()
{
    n=get_int(); m=get_int(); p=get_int();
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++)
      a[i][j]=get_int();
}

void det(int DX,int DY)
{
    for(int i=1;i<=n;i++)
    {
        qmin.clear();
        qmax.clear();
        for(int j=1;j<=m;j++)
        {
            while(!qmin.empty() && a[qmin.back().x][qmin.back().y]>=a[i][j]) qmin.popb();
            qmin.pb(mp(i,j));
            while(!qmax.empty() && a[qmax.back().x][qmax.back().y]<=a[i][j]) qmax.popb();
            qmax.pb(mp(i,j));

            if(j>=DY)
            {
                if(qmin.front().y==j-DY) qmin.popf();
                if(qmax.front().y==j-DY) qmax.popf();
                mi[i][j]=a[qmin.front().x][qmin.front().y];
                mx[i][j]=a[qmax.front().x][qmax.front().y];
            }
        }
    }

   int maxx,minn;
   for(int j=DY;j<=m;j++)
   {
       qmin.clear();
       qmax.clear();
       for(int i=1;i<=n;i++)
       {
           while(!qmin.empty() && mi[qmin.back().x][qmin.back().y]>=mi[i][j]) qmin.popb();
           qmin.pb(mp(i,j));
           while(!qmax.empty() && mx[qmax.back().x][qmax.back().y]<=mx[i][j]) qmax.popb();
           qmax.pb(mp(i,j));

           if(i>=DX)
           {
               if(qmin.front().x==i-DX) qmin.popf();
               if(qmax.front().x==i-DX) qmax.popf();
               maxx=mx[qmax.front().x][qmax.front().y];
               minn=mi[qmin.front().x][qmin.front().y];

               if(sol>maxx-minn) sol=maxx-minn,nr=1;
               else
                if(sol==maxx-minn) nr++;
           }
       }
   }
}

void solve()
{
    int X,Y;
    for(int i=1;i<=p;i++)
    {
        X=get_int(); Y=get_int();
        sol=inf; nr=0;;
        det(X,Y); if(X!=Y) det(Y,X);
        printf("%d %d\n",sol,nr);
    }
}

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

    read();
    solve();

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