#include <bits/stdc++.h>
using namespace std;
#pragma GCC target "avx2"
#pragma GCC optimize "Ofast"
#include <immintrin.h>

typedef __m256i m256;
typedef unsigned uint;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
typedef vector<int> vi;
typedef vector<pii> vpii;
typedef complex<double> cpx;
template <typename T> using minPq = priority_queue<T, vector<T>, greater<T>>;
#define ms(a, x) memset(a, x, sizeof(a))
#define pb push_back
#define fs first
#define sn second
#define ALL(v) begin(v), end(v)
#define SZ(v) ((int) (v).size())
#define lbv(v, x) (lower_bound(ALL(v), x) - (v).begin())
#define ubv(v, x) (upper_bound(ALL(v), x) - (v).begin())
template <typename T> inline void UNIQUE(vector<T> &v){sort(ALL(v)); v.resize(unique(ALL(v)) - v.begin());}
const int INF = 0x3f3f3f3f;
const ll LLINF = 0x3f3f3f3f3f3f3f3f;
const double PI = acos(-1);
#define FR(i, n) for(int i = 0; i < (n); i++)
#define FOR(i, a, b) for(int i = (a); i < (b); i++)
#define FORR(i, a, b) for(int i = (a); i >= (b); i--)
#define dbg(x) {cerr << #x << ' ' << x << endl;}
#define dbgArr(arr, n) {cerr << #arr; FR(_i, n) cerr << ' ' << (arr)[_i]; cerr << endl;}
template <typename T, typename U>
ostream& operator<<(ostream &os, pair<T, U> p){return os << "(" << p.fs << ", " << p.sn << ")";}

const int MOD = 998244353;
const int MOD_NINV = 998244351; // mod R
const int R_INV = 232013824; // mod M
const int ONE = (1ULL << 32) % MOD;

uint red(const uint &x){return min(x, x-MOD);}
uint redNeg(const uint &x){return min(x, x+MOD);}
uint redc(const ull &x){return red((x + (ull)(uint) ((ull)(uint) x*MOD_NINV) * MOD) >> 32);}
uint toMont(const uint &x){return ((ull) x << 32) % MOD;}
uint fromMont(const uint &x){return (ull) x * R_INV % MOD;}

#define vload(ptr) _mm256_loadu_si256((m256*) (ptr))
#define vstore(ptr, v) _mm256_storeu_si256((m256*) (ptr), v)
#define vrshift32(v) _mm256_bsrli_epi128(v, 4)

uint mulSum(int *arr, int *brr, int sz){
    m256 res = _mm256_setzero_si256();
    m256 mod_ninv = _mm256_set1_epi32(MOD_NINV);
    m256 mod = _mm256_set1_epi32(MOD);
    for(int i = 0; i < sz; i += 8){
        m256 a = vload(arr+i);
        m256 b = vload(brr+i);
        m256 v1 = _mm256_mul_epu32(a, b);
        m256 v2 = _mm256_mul_epu32(vrshift32(a), vrshift32(b));
        m256 p1 = _mm256_mul_epu32(_mm256_mul_epu32(v1, mod_ninv), mod);
        m256 p2 = _mm256_mul_epu32(_mm256_mul_epu32(v2, mod_ninv), mod);
        v1 = vrshift32(_mm256_add_epi64(v1, p1));
        v2 = _mm256_add_epi64(v2, p2);
        res = _mm256_add_epi32(res, _mm256_blend_epi32(v1, v2, 0b10101010));
        res = _mm256_min_epu32(res, _mm256_sub_epi32(res, mod));
        res = _mm256_min_epu32(res, _mm256_sub_epi32(res, mod));
    }

    int tmp[8];
    vstore(tmp, res);
    ull ans = 0;
    FR(i, 8) ans += tmp[i];
    return ans % MOD;
}

const int BSZ = 28;
const int MN = 505;
int n, m;
int xlrr[MN], xdrr[MN], xlPref[MN], xdPref[MN];
int ylrr[MN], ydrr[MN], ylPref[MN], ydPref[MN];
int ans[MN][MN];
int xBigArr[MN*MN], xSmallArr[BSZ][MN*BSZ], yArr[MN*MN*2];

void conv(int *arr, int *brr, int d, int sz){
    while((sz-d) & 7) sz++;
    
    brr[0] = arr[0];
    FR(i, sz) brr[i+1] = red(brr[i] + arr[i+1]);

    m256 mod = _mm256_set1_epi32(MOD);
    for(int i = sz-7; i > d; i -= 8){
        m256 v = _mm256_sub_epi32(vload(brr+i), vload(brr+i-d-1));
        v = _mm256_min_epu32(v, _mm256_add_epi32(v, mod));
        vstore(brr+i, v);
    }
}

void procBlk(int bl, int br){
    FR(i, xdrr[bl]+1) xSmallArr[0][i] = ONE;
    FOR(i, bl+1, br) conv(xSmallArr[i-1-bl], xSmallArr[i-bl], xdrr[i], xdPref[i+1]-xdPref[bl]);

    int ysz = xdPref[bl];
    FR(i, ysz+1) yArr[i] = xBigArr[ysz-i];

    FR(bp, m){
        conv(yArr, yArr, ydrr[bp], ysz += ydrr[bp]);

        FOR(ap, bl, br){
            int dif = xlPref[ap+1] - ylPref[bp+1] + xdPref[bl];
            int l = max(0, -dif);
            int r = min(ysz-dif, xdPref[ap+1]-xdPref[bl]) + 1;
            ans[ap][bp] = mulSum(xSmallArr[ap-bl] + l, yArr + dif + l, r-l);
        }
    }

    FR(i, ysz+1) yArr[i] = 0;
    FOR(i, bl, br) FR(j, xdPref[i+1]-xdPref[bl]+1) xSmallArr[i-bl][j] = 0;

    FOR(i, bl, br) conv(xBigArr, xBigArr, xdrr[i], xdPref[i+1]);
}

int main(){
    cin.sync_with_stdio(0); cin.tie(0);
    cin >> n >> m;
    FR(i, n){
        int l, r;
        cin >> l >> r;
        xlrr[i] = l;
        xdrr[i] = r-l;
        xlPref[i+1] = xlPref[i] + xlrr[i];
        xdPref[i+1] = xdPref[i] + xdrr[i];
    }
    FR(i, m){
        int l, r;
        cin >> l >> r;
        ylrr[i] = l;
        ydrr[i] = r-l;
        ylPref[i+1] = ylPref[i] + ylrr[i];
        ydPref[i+1] = ydPref[i] + ydrr[i];
    }

    xBigArr[0] = ONE;
    for(int bl = 0; bl < n; bl += BSZ) procBlk(bl, min(n, bl+BSZ));

    FR(i, n) FR(j, m) cout << fromMont(ans[i][j]) << " \n"[j == m-1];
}