Pagini recente » Cod sursa (job #3230154) | Borderou de evaluare (job #2912962) | Cod sursa (job #3279391) | Cod sursa (job #2566000) | Cod sursa (job #1988453)
#include<fstream>
#include<deque>
using namespace std;
class InParser {
private:
FILE *fin;
char *buff;
int sp;
char read_ch() {
++sp;
if (sp == 4096) {
sp = 0;
fread(buff, 1, 4096, fin);
}
return buff[sp];
}
public:
InParser(const char* nume) {
fin = fopen(nume, "r");
buff = new char[4096]();
sp = 4095;
}
InParser& operator >> (int &n) {
char c;
while (!isdigit(c = read_ch()) && c != '-');
int sgn = 1;
if (c == '-') {
n = 0;
sgn = -1;
} else {
n = c - '0';
}
while (isdigit(c = read_ch())) {
n = 10 * n + c - '0';
}
n *= sgn;
return *this;
}
InParser& operator >> (long long &n) {
char c;
n = 0;
while (!isdigit(c = read_ch()) && c != '-');
long long sgn = 1;
if (c == '-') {
n = 0;
sgn = -1;
} else {
n = c - '0';
}
while (isdigit(c = read_ch())) {
n = 10 * n + c - '0';
}
n *= sgn;
return *this;
}
} in( "struti.in" );
ofstream out("struti.out");
int n,m,p,h,i,j,minim,l,a,b,stv,drv,stw,drw;
int mat[1001][1001],mil[1001][1001],mic[1001][1001],mal[1001][1001],mac[1001][1001];
int w[10001],v[10001];
int main(){
in >> n >> m >> p;
for( i = 1; i <= n; i ++ ){
for( j = 1; j <= m; j ++ ){
in >> mat[i][j];
}
}
for( h = 1; h <= p; h ++ ){
in >> a >> b;
minim = 8100;
l = 0;
for( i = 1; i <= n; i ++ ){
stw = 1; drw = 0;
stv = 1; drv = 0;
for( j = 1; j <= m; j ++ ){
while( drw >= stw && mat[i][j] < mat[i][w[drw]] ){
drw --;
}
drw++;
w[drw] = j;
if(w[stw] <= j - b ){
stw++;
}
mil[i][j] = mat[i][w[stw]];
while( drv >= stv && mat[i][j] > mat[i][v[drv]] ){
drv--;
}
drv++;
v[drv] = j;
if(v[stv] <= j - b ){
stv++;
}
mal[i][j] = mat[i][v[stv]];
}
}
for( j = 1; j <= m; j ++ ){
stw = 1; drw = 0;
stv = 1; drv = 0;
for( i = 1; i <= n; i ++ ){
while( drw >= stw && mil[i][j] < mil[w[drw]][j] ){
drw--;
}
drw++;
w[drw] = i;
if(w[stw] <= i - a ){
stw++;
}
mic[i][j] = mil[w[stw]][j];
while( drv >= stv && mal[i][j] > mal[v[drv]][j] ){
drv--;
}
drv++;
v[drv] = i;
if(v[stv] <= i - a ){
stv++;
}
mac[i][j] = mal[v[stv]][j];
}
}
for( i = a; i <= n; i ++ ){
for( j = b; j <= m; j ++ ){
if( mac[i][j] - mic[i][j] < minim ){
minim = mac[i][j] - mic[i][j];
l = 0;
}
if( mac[i][j] - mic[i][j] == minim ){
l++;
}
}
}
if( a != b ){
for( i = 1; i <= n; i ++ ){
stw = 1; drw = 0;
stv = 1; drv = 0;
for( j = 1; j <= m; j ++ ){
while( drw >= stw && mat[i][j] < mat[i][w[drw]] ){
drw --;
}
drw++;
w[drw] = j;
if(w[stw] <= j - a ){
stw++;
}
mil[i][j] = mat[i][w[stw]];
while( drv >= stv && mat[i][j] > mat[i][v[drv]] ){
drv--;
}
drv++;
v[drv] = j;
if(v[stv] <= j - a ){
stv++;
}
mal[i][j] = mat[i][v[stv]];
}
}
for( j = 1; j <= m; j ++ ){
stw = 1; drw = 0;
stv = 1; drv = 0;
for( i = 1; i <= n; i ++ ){
while( drw >= stw && mil[i][j] < mil[w[drw]][j] ){
drw--;
}
drw++;
w[drw] = i;
if(w[stw] <= i - b ){
stw++;
}
mic[i][j] = mil[w[stw]][j];
while( drv >= stv && mal[i][j] > mal[v[drv]][j] ){
drv--;
}
drv++;
v[drv] = i;
if(v[stv] <= i - b ){
stv++;
}
mac[i][j] = mal[v[stv]][j];
}
}
for( i = b; i <= n; i ++ ){
for( j = a; j <= m; j ++ ){
if( mac[i][j] - mic[i][j] < minim ){
minim = mac[i][j] - mic[i][j];
l = 0;
}
if( mac[i][j] - mic[i][j] == minim ){
l++;
}
}
}
}
out<<minim <<" "<<l<<"\n";
}
return 0;
}