#define _USE_MATH_DEFINES
#include <bits/stdc++.h>

using namespace std;

using i64 = long long;
#define forn(a, e) for (i64 a = 0; a < (i64)(e); a++)
#define forr(a, s, e) for (i64 a = s; a < (i64)(e); a++)
#define fore(e, a) for (auto& e : a)

#ifdef LOCAL
#define logv(a) {cerr << #a << " = "; fore(e, a) {cerr << e << " ";} cerr << "\n";}
#define logvp(a) {cerr << #a << " = "; fore(e, a) {cerr << "(" << e.first << ", " << e.second << ") ";} cerr << "\n";}
#define logvv(a) {cerr << #a << " = \n"; fore(r, a) { fore(e, r) {cerr << e << " ";} cerr << "\n";} }
#define logvf(a, field) {cerr << #a"."#field << " = \n"; fore(e, a) { cerr << e.field << " ";} cerr << "\n"; }
#define logvff(a, f1, f2) {cerr << #a".{"#f1 << ", "#f2 << "} = \n"; fore(e, a) { cerr << "(" << e.f1 <<", " << e.f2 << ") ";} cerr << "\n"; }
#define logs(a) cerr << #a << " = " << (a) << "\n";
#define logss(a, b) cerr << #a << " = " << (a) << ", " << #b << " = " << (b) << "\n";
#define logp(a) cerr << #a << " = " << "(" << a.first << ", " << a.second << ")" << "\n";
#define cond(pred, stmt) if (pred) { stmt }
#else
#define logv(a)
#define logvp(a)
#define logvv(a)
#define logvf(a, field)
#define logvff(a, f1, f2)
#define logs(a)
#define logss(a, b)
#define logp(a)
#define cond(pred, stmt)
#endif

using iip = pair<int, int>;
using llp = pair<i64, i64>;
using ivec = vector<int>;
using llvec = vector<i64>;
using svec = vector<string>;
template<typename T> using vec = vector<T>;

template<typename T, typename Dim>
auto make_vec(T value, Dim dim) { return vector<T>(dim, value); }
template<typename T, typename Dim1, typename... Dim>
auto make_vec(T value, Dim1 dim1, Dim... dims) { return make_vec(make_vec(value, dims...), dim1); }

template<typename T>
bool uax(T& v, const T& newv) { if (v < newv) { v = newv; return true; } else return false; }
template<typename T>
bool uin(T& v, const T& newv) { if (v > newv) { v = newv; return true; } else return false; }

template<typename T>
istream& operator>>(istream& is, vector<T>& c) { for (auto& e : c) is >> e; return is; }
template<typename T, size_t N>
istream& operator>>(istream& is, array<T, N>& c) { for (auto& e : c) is >> e; return is; }

template<typename ...T>
istream& read(T&... args) { return (cin >> ... >> args); }

static mt19937 rande(123123);
template<typename T>
T rand_int(T from, T to) { uniform_int_distribution<T> distr(from, to); return distr(rande); }

const i64 INF = 1e18;
const i64 M = 1e9 + 7;
const int MXB = 20;
using bin = bitset<MXB + 1>;


i64 bpow(i64 a, i64 p) {
    if (p == 1) {
        return a;
    }
    if (p == 0) {
        return 1;
    }

    i64 res = bpow(a, p / 2);
    res *= res;
    res %= M;

    if (p % 2) {
        res *= a;
        res %= M;
    }

    return res;
}

i64 inv(i64 a) {
    return bpow(a, M - 2);
}

int main() {
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);
    rande.seed(chrono::steady_clock::now().time_since_epoch().count());

    int n, c;
    while (read(n, c)) {
        ivec a(n);
        ivec b(n);
        read(a, b);

        auto dp = make_vec(0ll, n + 1, c + 1);
        dp[0][0] = 1;
        auto pp = make_vec(0ll, *max_element(b.begin(), b.end()) + 1, c + 1);
        auto ppi = pp;
        forr(i, 1, pp.size()) forn(j, c + 1) {
            pp[i][j] = bpow(i, j);
            ppi[i][j] = inv(pp[i][j]);
        }

        forr(i, 1, n + 1) {
            forr(j, a[i - 1], b[i - 1] + 1) {
                llvec cc(c + 2);
                forn(k, c + 1) {
                    cc[k + 1] = (cc[k] + dp[i - 1][k] * ppi[j][k] % M) % M;
                }
                forr(k, 0, c + 1) {
                    dp[i][k] += cc[k + 1] * pp[j][k] % M;
                    dp[i][k] %= M;
                }
            }
        }

        logvv(dp);
        cout << dp[n][c] << endl;
    }
}

Measured memory_bytes: 7176