1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64
| #include <iostream> using namespace std; const int N = 1000+10; int a,b,n; int s[N][N],q[N],hh=0,tt=-1,d; int anserA[N],anserB[N],anserC[N],anserD[N];
int rsMin[N][N],rsMax[N][N];
void get_max(int * p,int * f,int col){ hh=0,tt=-1,d=0; for(int i=0;i<col;i++){ if(i-n+1>q[hh]) hh++; while(hh<=tt&&p[q[tt]]<=p[i]) tt--; q[++tt] = i; if(i>=n-1) f[d++] = p[q[hh]]; } f[N-1] = d; }
void get_min(int * p,int * f,int col){ hh=0,tt=-1,d=0; for(int i=0;i<col;i++){ if(i-n+1>q[hh]) hh++; while(hh<=tt&&p[q[tt]]>=p[i]) tt--; q[++tt] = i; if(i>=n-1) f[d++] = p[q[hh]]; } f[N-1] = d; }
int main() { cin>>a>>b>>n;
for(int i=0;i<a;i++) for(int j=0;j<b;j++) cin>>s[i][j]; for(int i=0;i<a;i++){ get_max(s[i],rsMax[i],b); get_min(s[i],rsMin[i],b); }
int anser = 1e10;
for(int i=0;i<rsMax[0][N-1];i++){ for(int j=0;j<a;j++){ anserA[j] = rsMax[j][i]; anserB[j] = rsMin[j][i]; } get_max(anserA,anserC,a); get_min(anserB,anserD,a); for(int k=0;k<anserC[N-1];k++) { anser = min(anser,anserC[k]-anserD[k]); } } cout << anser; return 0; }
|