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

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 BSZ = 28;
const int MN = 505;
const int MOD = 998244353;
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];

uint red(const uint &x){return min(x, x-MOD);}
uint redNeg(const uint &x){return min(x, x+MOD);}

void conv(int *arr, int *brr, int d, int sz){
    brr[0] = arr[0];
    FR(i, sz) brr[i+1] = red(brr[i] + arr[i+1]);
    FORR(i, sz, d+1) brr[i] = redNeg(brr[i] - brr[i-d-1]);
}

void procBlk(int bl, int br){
    FR(i, xdrr[bl]+1) xSmallArr[0][i] = 1;
    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){
            ll v = 0;
            int dif = xlPref[ap+1] - ylPref[bp+1] + xdPref[bl];
            FOR(i, max(0, -dif), min(ysz-dif, xdPref[ap+1]-xdPref[bl]) + 1){
                v += (ull) xSmallArr[ap-bl][i] * yArr[i+dif] % MOD;
            }
            ans[ap][bp] = v % MOD;
        }
    }

    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] = 1;
    for(int bl = 0; bl < n; bl += BSZ) procBlk(bl, min(n, bl+BSZ));

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