Cod sursa(job #2531351)

Utilizator alex2209alexPavel Alexandru alex2209alex Data 26 ianuarie 2020 10:21:03
Problema Struti Scor 100
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 3.68 kb
#include <bits/stdc++.h>

using namespace std;
//ifstream f("struti.in");
//ofstream g("struti.out");
FILE *input = fopen("struti.in","r");
FILE *output = fopen("struti.out","w");
//------------------------------------
///Globale
int raspuns2;
short n,m,v[1001][1001],raspuns,nr,v2[1001][1001],v3[1001][1001],in[1001],sf[1001],in2[1001],sf2[1001];
deque<pair<short,short>>dq2,dq3;
pair<int,int>mi[1001][1001],ma[1001][1001];
//------------------------------------
///Functii
void citire();
//------------------------------------
int main()
{
    ios::sync_with_stdio(false);
	citire();
	return 0;
}
//------------------------------------
void afisare()
{
    fprintf(output,"%hu %d\n",raspuns,raspuns2);
}
//------------------------------------
void calculare_raspuns(short c)
{
	for(short i = 1; i <= nr; ++i)
	{
		dq2.clear();
		dq3.clear();
		for(short j = 1; j <= c; ++j)
		{
			while(!dq2.empty() && v2[i][j] < dq2.back().first)
				dq2.pop_back();
			dq2.push_back({v2[i][j],j});
			while(!dq3.empty() && v3[i][j] > dq3.back().first)
				dq3.pop_back();
			dq3.push_back({v3[i][j],j});
		}
		short r = dq3.front().first - dq2.front().first;
		if(r < raspuns)
		{
			raspuns = r;
			raspuns2 = 1;
		}
		else if(r == raspuns)
			raspuns2++;
		for(short j = c + 1; j <= m; ++j)
		{
			while(!dq2.empty() && v2[i][j] < dq2.back().first)
				dq2.pop_back();
			dq2.push_back({v2[i][j],j});
			while(!dq3.empty() && v3[i][j] > dq3.back().first)
				dq3.pop_back();
			dq3.push_back({v3[i][j],j});
			if(dq2.front().second == j - c)
				dq2.pop_front();
			if(dq3.front().second == j - c)
				dq3.pop_front();
			short r = dq3.front().first - dq2.front().first;
			if(r < raspuns)
			{
				raspuns = r;
				raspuns2 = 1;
			}
			else if(r == raspuns)
				raspuns2++;
		}
	}
}
//------------------------------------
void calculare_minim_maxim_pe_linii(short l)
{
	nr = 0;
	for(short i = 1; i <= m; ++i)
	{
		in[i] = 1;
        sf[i] = 0;
        in2[i] = 1;
        sf2[i] = 0;
	}
	for(short j = 1; j <= m; ++j)
		for(short i = 1; i <= l; ++i)
		{
			while(sf[j] >= in[j] && v[i][j] <= mi[sf[j]][j].first)
				sf[j]--;
            sf[j]++;
            mi[sf[j]][j] = {v[i][j],i};
			while(sf2[j] >= in2[j] && v[i][j] >= ma[sf2[j]][j].first)
				sf2[j]--;
            sf2[j]++;
            ma[sf2[j]][j] = {v[i][j],i};
		}
	nr++;
	for(short j = 1; j <= m; ++j)
		v2[nr][j] = mi[1][j].first;
    for(short j = 1; j <= m; ++j)
		v3[nr][j] = ma[1][j].first;
	for(short i = l + 1; i <= n; ++i)
	{
		nr++;
		for(short j = 1; j <= m; ++j)
		{
			while(sf[j] >= in[j] && v[i][j] <= mi[sf[j]][j].first)
				sf[j]--;
            sf[j]++;
            mi[sf[j]][j] = {v[i][j],i};
			if(mi[in[j]][j].second == i - l)
				in[j]++;
            while(sf2[j] >= in2[j] && v[i][j] >= ma[sf2[j]][j].first)
				sf2[j]--;
            sf2[j]++;
            ma[sf2[j]][j] = {v[i][j],i};
			if(ma[in2[j]][j].second == i - l)
				in2[j]++;
			v2[nr][j] = mi[in[j]][j].first;
			v3[nr][j] = ma[in2[j]][j].first;
		}
	}
}
//------------------------------------
void rezolvare(short l,short c)
{
	if(l <= n && c <= m)
	{
		calculare_minim_maxim_pe_linii(l);
		calculare_raspuns(c);
	}
	if(l == c)
		return;
	swap(l,c);
	if(l <= n && c <= m)
	{
		calculare_minim_maxim_pe_linii(l);
		calculare_raspuns(c);
	}
}
//------------------------------------
void citire()
{
	short p;
	fscanf(input,"%hu%hu%hu",&n,&m,&p);
	for(short i = 1; i <= n; ++i)
		for(short j = 1; j <= m; ++j)
			fscanf(input,"%hu",&v[i][j]);
	while(p--)
	{
		raspuns = 8001;
		raspuns2 = 0;
		short l,c;
		fscanf(input,"%hu%hu",&l,&c);
		rezolvare(l,c);
		afisare();
	}
}