#include<bits/stdc++.h>
using namespace std;

using i64 = long long;
using u64 = unsigned long long;

struct debug {
#define contPrint { *this << "["; \
        int f = 0; for(auto it : x) { *this << (f?", ":""); *this << it; f = 1;} \
        *this << "]"; return *this;}
 
    ~debug(){cerr << endl;}
    template<class c> debug& operator<<(c x) {cerr << x; return *this;}
    template<class c, class d>
    debug& operator<<(pair<c, d> x) {*this << "(" << x.first << ", " << x.second << ")"; 
        return *this;}
    template<class c> debug& operator<<(vector<c> x) contPrint;
#undef contPrint
};

#define dbg(x) "[" << #x << ": " << x << "]  "
#define Wa() cerr << "[LINE: " << __LINE__ << "] -> "; debug() << 
#define FASTIO ios_base::sync_with_stdio(false); cin.tie(NULL);

using ld = long double;
using cd = complex<ld>;
const ld PI = acos(-1.0);


void fft(vector<cd>& a, bool invert)
{
    int n = (int)a.size();

    for(int i = 1, j = 0; i < n; i++) {
        int bit = n >> 1;
        for(; j & bit; bit >>= 1) j ^= bit;
        j ^= bit;
        if(i < j) swap(a[i], a[j]);
    }

    for(int len = 2; len <= n; len <<= 1) {
        ld ang = 2 * PI / len * (invert ? -1 : 1);
        cd wlen(cosl(ang), sinl(ang));
        for(int i = 0; i < n; i += len) {
            cd w(1);
            for(int j = 0; j < len / 2; j++) {
                cd u = a[i + j], v = a[i + j + len/2] * w;
                a[i + j] = u + v;
                a[i + j + len/2] = u - v;
                w *= wlen;
            }
        }
    }
    if(invert) {
        for(cd& x : a) {
            x /= n;
        }
    }
}

vector<i64> multiply(vector<i64> const& a, vector<i64> const& b) {
    vector<cd> fa(a.begin(), a.end()), fb(b.begin(), b.end());
    int n = 1;
    while(n < (int)a.size() + b.size()) {
        n <<= 1;
    }
    fa.resize(n);
    fb.resize(n);

    fft(fa, false);
    fft(fb, false);
    for(int i = 0; i < n; i++) fa[i] *= fb[i];
    fft(fa, true);

    vector<i64> res(n);
    for(int i = 0; i < n; i++) res[i] = round(fa[i].real());
    return res;
}

int main() 
{
    int n1, m1;
    scanf("%d %d", &n1, &m1);
    int motif[n1][m1];
    for(int i = 0; i < n1; i++) {
        for(int j = 0; j < m1; j++) {
            scanf("%d", &motif[i][j]);
        }
    }

    int n2, m2;
    scanf("%d %d", &n2, &m2);
    int mosaic[n2][m2];
    for(int i = 0; i < n2; i++) {
        for(int j = 0; j < m2; j++) {
            scanf("%d", &mosaic[i][j]);
        }
    }

    if(n1 > n2 || m1 > m2) {
        puts("0");
        return 0;
    }

    vector<i64> motif_line((n1 - 1) * m2 + m1);
    for(int i = n1-1, start = 0; i >= 0; i--, start += m2) {
        for(int j = m1-1, k = 0; j >= 0; j--, k++) {
            motif_line[start + k] = motif[i][j];
        }
    }

    vector<i64> mosaic_line(n2 * m2);
    for(int i = 0, k = 0; i < n2; i++) {
        for(int j = 0; j < m2; j++) {
            mosaic_line[k++] = mosaic[i][j];
        }
    }

    vector<i64> motif_line_sq;
    vector<i64> motif_line_cb;
    for(int i : motif_line) {
        motif_line_sq.push_back(i * i);
        motif_line_cb.push_back(i * i * i);
    }

    vector<i64> mosaic_line_sq;
    vector<i64> mosaic_line_cb;
    for(int i : mosaic_line) {
        mosaic_line_sq.push_back(i * i);
        mosaic_line_cb.push_back(i * i * i);
    }

    auto res1 = multiply(mosaic_line_cb, motif_line);
    auto res2 = multiply(mosaic_line_sq, motif_line_sq);
    auto res3 = multiply(mosaic_line, motif_line_cb);

    vector<i64> res4(res1.size());

    for(int i = 0; i < (int)res1.size(); i++) {
        res4[i] = res1[i] - 2 * res2[i] + res3[i];
        // Wa() dbg(i) dbg(res1[i]) dbg(2*res2[i]) dbg(res3[i]);
    }

    vector<pair<int, int>> ans;
    for(int i = 0; i <= n2-n1; i++) {
        for(int j = 0; j <= m2-m1; j++) {
            // Wa() dbg(i * m2 + j + (n1 - 1) * m2 + m1 - 1) dbg(res4[i * m2 + j + (n1 - 1) * m2 + m1 - 1]);
            if(res4[i * m2 + j + (n1 - 1) * m2 + m1 - 1] == 0) {
                ans.push_back({i + 1, j + 1});
            }
        }
    }
    sort(ans.begin(), ans.end());

    printf("%d\n", (int)ans.size());
    for(auto it : ans) {
        printf("%d %d\n", it.first, it.second);
    }
}