#include #include using namespace std; #define FOR(i, begin, end) for(int i=(begin),i##_end_=(end);i #include using mint = atcoder::modint998244353; // [x^B] 1/(m!) * (x / (1 - x)^2 )^m (1 / (1 - x))^(W - m + 1) を m = 0, ..., W で列挙 // [x^B] 1/m! * x^m / (1 - x)^(W + 1 + m) // 1/m! [x^(B - m)] 1 / (1 - x)^(W + 1 + m) // 1/m! Comb(B + W, W + m) vector gen(int W, int B) { vector ret(W + 1); mint num = 1, den = 1; REP(_, W) { num *= B + W - _; den *= W - _; } ret.at(0) = num / den; FOR(m, 1, W + 1) { num *= B + W - W + 1 - m; den *= W + m; den *= m; ret.at(m) = num / den; } return ret; } int main() { cin.tie(nullptr); ios::sync_with_stdio(false); int H, W; cin >> H >> W; vector dp{1}; for (int i = 0; i < H; ++i) { int a, b; cin >> a >> b; auto f = gen(W, b); auto g = gen(W, a - 1); REP(i, g.size()) f.at(i) -= g.at(i); dp = atcoder::convolution(dp, f); dp.resize(W + 1); } mint ret = dp.back(); FOR(i, 1, W + 1) ret *= i; cout << ret.val() << '\n'; }