#include <bits/stdc++.h>

constexpr int mod = 1e9 + 7;
int add(int a, int b) {
    return a + b - mod * (a + b >= mod);
}
int add(int a, auto... args) {
    return add(a, add(args...));
}
void add_to(int& a, int b) {
    a = add(a, b);
}
int mul(int a, int b) {
    return a * 1ULL * b % mod;
}
int mul(int a, auto... args) {
    return mul(a, mul(args...));
}
int power(int b, int e) {
    int r = 1;
    for (; e > 0; e >>= 1) {
        if (e & 1)
            r = mul(r, b);
        b = mul(b, b);
    }
    return r;
}

struct ST {
    struct Data {
        int len;
        int cnt0, cnt1;
        int fuck, fuck_l, fuck_r;

        Data(int val = -1) {
            if (val == -1) {
                len = 0;
                return;
            }
            len = 1;
            cnt0 = 2 - val;
            cnt1 = val;
            fuck = 0;
            fuck_l = 0;
            fuck_r = 0;
        }
    };

    Data merge(Data a, Data b) {
        if (a.len == 0 || b.len == 0) {
            return a.len == 0 ? b : a;
        }

        Data d(-1);

        d.len = a.len + b.len;
        d.cnt0 = mul(a.cnt0, b.cnt0);
        d.cnt1 = add(mul(a.cnt1, b.cnt1), mul(a.cnt0, b.cnt1), mul(a.cnt1, b.cnt0));
        d.fuck = add(
            mul(a.fuck, add(b.cnt0, b.cnt1)),
            mul(b.fuck, add(a.cnt0, a.cnt1)),
            mul(a.cnt1, b.fuck_l),
            mul(b.cnt1, a.fuck_r));

        d.fuck_l = add(mul(a.cnt0, a.len, b.cnt1), mul(a.cnt0, b.fuck_l), mul(a.fuck_l, add(b.cnt0, b.cnt1)));
        d.fuck_r = add(mul(b.cnt0, b.len, a.cnt1), mul(b.cnt0, a.fuck_r), mul(b.fuck_r, add(a.cnt0, a.cnt1)));

        return d;
    }

    std::vector<Data> data;
    int size;

    ST(int n) {
        size = std::bit_ceil((size_t)n);
        data.assign(size * 2, Data(-1));
        for (int i = 0; i < n; i++) {
            data[size + i] = Data(0);
        }
        // size = 1 << std::__lg(std::max)
    }

    void update(int i, int val) {
        i += size;
        data[i] = Data(val);
        for (i >>= 1; i > 0; i >>= 1) {
            data[i] = merge(data[2 * i], data[2 * i + 1]);
        }
    }

    int get() {
        return data[1].fuck;
    }
};

int32_t main() {
    std::ios_base::sync_with_stdio(false);
    std::cin.tie(nullptr);

    int n;
    std::cin >> n;
    std::vector<std::array<int, 2>> input(n);
    for (int i = 0; i < 2; i++) {
        for (int j = 0; j < n; j++) {
            std::cin >> input[j][i];
        }
    }
    std::vector<int> vec(2 * n);
    std::iota(vec.begin(), vec.end(), 0);
    std::sort(vec.rbegin(), vec.rend(), [&](int a, int b) {
        return input[a / 2][a % 2] < input[b / 2][b % 2];
    });

    ST st(n);
    std::vector<int> cnt(n);
    int ans = 0;
    for (int i = 1; i < vec.size(); i++) {
        int a = vec[i - 1], b = vec[i];
        st.update(a / 2, ++cnt[a / 2]);

        int dt = input[a / 2][a % 2] - input[b / 2][b % 2];
        int ct = st.get();
        ans = add(ans, mul(dt, ct));
    }

    std::cout << ans << "\n";

    return 0;
}