#include <bits/stdc++.h>
#define fst first
#define snd second
#define fore(i,a,b) for(int i=a,ThxDem=b;i<ThxDem;++i)
#define pb push_back
#define ALL(s) s.begin(),s.end()
#define FIN ios_base::sync_with_stdio(0);cin.tie(0);cout.tie(0)
#define SZ(s) int(s.size())
using namespace std;
typedef long long ll;
typedef pair<int,int> ii;

// MAXN must be power of 2 !!
// MOD-1 needs to be a multiple of MAXN !!
// big mod and primitive root for NTT:
typedef ll tf;
typedef vector<tf> poly;
const tf MOD=2305843009255636993,RT=5,MAXN=1<<21;
// FFT
struct CD {
	double r,i;
	CD(double r=0, double i=0):r(r),i(i){}
	double real()const{return r;}
	void operator/=(const int c){r/=c, i/=c;}
};
CD operator*(const CD& a, const CD& b){
	return CD(a.r*b.r-a.i*b.i,a.r*b.i+a.i*b.r);}
CD operator+(const CD& a, const CD& b){return CD(a.r+b.r,a.i+b.i);}
CD operator-(const CD& a, const CD& b){return CD(a.r-b.r,a.i-b.i);}
const double pi=acos(-1.0);
// NTT
/*
struct CD {
	tf x;
	CD(tf x):x(x){}
	CD(){}
};
CD operator*(const CD& a, const CD& b){return CD(mulmod(a.x,b.x));}
CD operator+(const CD& a, const CD& b){return CD(addmod(a.x,b.x));}
CD operator-(const CD& a, const CD& b){return CD(submod(a.x,b.x));}
vector<tf> rts(MAXN+9,-1);
CD root(int n, bool inv){
	tf r=rts[n]<0?rts[n]=pm(RT,(MOD-1)/n):rts[n];
	return CD(inv?pm(r,MOD-2):r);
}
*/
CD cp1[MAXN+9],cp2[MAXN+9];
int R[MAXN+9];
void dft(CD* a, int n, bool inv){
	fore(i,0,n)if(R[i]<i)swap(a[R[i]],a[i]);
	for(int m=2;m<=n;m*=2){
		double z=2*pi/m*(inv?-1:1); // FFT
		CD wi=CD(cos(z),sin(z)); // FFT
		// CD wi=root(m,inv); // NTT
		for(int j=0;j<n;j+=m){
			CD w(1);
			for(int k=j,k2=j+m/2;k2<j+m;k++,k2++){
				CD u=a[k];CD v=a[k2]*w;a[k]=u+v;a[k2]=u-v;w=w*wi;
			}
		}
	}
	if(inv)fore(i,0,n)a[i]/=n; // FFT
	//if(inv){ // NTT
	//	CD z(pm(n,MOD-2)); // pm: modular exponentiation
	//	fore(i,0,n)a[i]=a[i]*z;
	//}
}
poly multiply(poly& p1, poly& p2){
	int n=p1.size()+p2.size()+1;
	int m=1,cnt=0;
	while(m<=n)m+=m,cnt++;
	fore(i,0,m){R[i]=0;fore(j,0,cnt)R[i]=(R[i]<<1)|((i>>j)&1);}
	fore(i,0,m)cp1[i]=0,cp2[i]=0;
	fore(i,0,p1.size())cp1[i]=p1[i];
	fore(i,0,p2.size())cp2[i]=p2[i];
	dft(cp1,m,false);dft(cp2,m,false);
	fore(i,0,m)cp1[i]=cp1[i]*cp2[i];
	dft(cp1,m,true);
	poly res;
	n-=2;
	fore(i,0,n)res.pb((tf)floor(cp1[i].real()+0.5)); // FFT
	//fore(i,0,n)res.pb(cp1[i].x); // NTT
	return res;
}

poly add(poly p, poly q){
	int n=SZ(p),m=SZ(q);
	poly ans(max(n,m));
	fore(i,0,n) ans[i]+=p[i];
	fore(i,0,m) ans[i]+=q[i];
	return ans;
}

poly sub(poly p, poly q){
	int n=SZ(p),m=SZ(q);
	poly ans(max(n,m));
	fore(i,0,n) ans[i]+=p[i];
	fore(i,0,m) ans[i]-=q[i];
	return ans;
}

int main(){FIN;
	int n1,m1; cin>>n1>>m1;
	vector<vector<int>> a(n1,vector<int>(m1));
	fore(i,0,n1) fore(j,0,m1) cin>>a[i][j];
	
	int n2,m2; cin>>n2>>m2;
	vector<vector<int>> b(n2,vector<int>(m2));
	fore(i,0,n2) fore(j,0,m2) cin>>b[i][j];
	
	
	poly p;
	fore(i,0,n1){
		fore(j,0,m1) p.pb(a[i][j]);
		if(i!=n1-1) fore(j,0,m2-m1) p.pb(0);
	}
	
	poly q;
	fore(i,0,n2) fore(j,0,m2) q.pb(b[i][j]);
	
	// (p-q)^2 = p^2 - 2*p*q + q^2
	
	ll sum2=0;
	for(auto &x:p) sum2+=x*x;
	auto sq=q;
	for(auto &x:sq) x=x*x;
	fore(i,1,SZ(sq)) sq[i]+=sq[i-1];
	
	poly qr=p;
	reverse(ALL(qr));
	
	poly now=multiply(q,qr);
	poly me;
	fore(i,0,SZ(q)-SZ(p)+1){
	  ll val=sq[i+SZ(p)-1]; if(i) val-=sq[i-1];
	  val-=2ll*now[SZ(p)+i-1];
	  val+=sum2;
	  me.pb(val);
	}
	
	//(p0*q^2)
	poly p0=p;
	for(auto &x:p0) x=x==0;
	for(auto &x:q) x=x*x;
	reverse(ALL(p0));
	poly now2=multiply(p0,q);
	poly he;
	fore(i,0,SZ(q)-SZ(p0)+1) he.pb(now2[SZ(p)+i-1]);

	// (p-q)^2 - (p0*q^2)
	poly ans=sub(me,he);
	
	vector<ii> res;
	
	fore(i,0,n2) fore(j,0,m2) if(i+n1<=n2 && j+m1<=m2 && ans[m2*i+j]==0){
		res.pb({i,j});
	}
	
	cout<<SZ(res)<<"\n";
	for(auto x:res)cout<<x.fst+1<<" "<<x.snd+1<<"\n";
}