#include <bits/stdc++.h>

using namespace std;

using ll = long long;
constexpr ll mod = 1e9 + 7;

constexpr ll N = 2e7 + 5;

int fact[N + 1], inv[N + 1];

ll C(ll n, ll k) {
    if (k < 0)return 0;
    return (ll) fact[n] * inv[k] % mod * inv[n - k] % mod;
}

ll obr(ll n) {
    return (ll) inv[n] * fact[n - 1];
}

ll T(ll n, ll k) {
    return C(n + k, n) * (n - k + 1) % mod * obr(n + 1);
}

ll modPow(ll x, ll n) {
    if (!n)return 1;
    ll u = modPow(x, n / 2);
    u = u * u % mod;
    if (n % 2)u = u * x % mod;
    return u;
}

void solve() {
    ll n, c;
    cin >> n >> c;
    vector<int> small(n + 1, 1), big(n + 1, 1);
    for (ll i = 1; i <= n; ++i)small[i] = (ll) (c - 1) * small[i - 1] % mod;
    for (ll i = 1; i <= n; ++i)big[i] = (ll) c * big[i - 1] % mod;
    ll res = 0;
    for (ll t = 1; t <= n; ++t) {
        ll ways = T(n - 1, n - t);
        (res += ways * big[t] % mod * small[n - t]) %= mod;
    }
    cout << res << "\n";
}

int main() {
    cin.tie(0)->sync_with_stdio(false);
    fact[0] = 1;
    for (ll i = 1; i <= N; ++i)fact[i] = (ll) fact[i - 1] * i % mod;
    inv[N] = modPow((ll) fact[N], mod - 2);
    for (ll i = N - 1; ~i; --i) {
        inv[i] = (ll) inv[i + 1] * (i + 1) % mod;
    }
    ll t;
    cin >> t;
    while (t--)solve();
}