# include <bits/stdc++.h>

using namespace std;

long n, m, p;

inline long pow(long a, long n) {
  long r(1);
  while (n) {
    r = r * (n & 1 ? a : 1) % p;
    a = a * a % p, n >>= 1;
  }
  return r;
}
inline long inv(long x) {
  x = (x % p + p) % p;
  if (!x) return 0;
  long res(1);
  while (x != 1)
    res = p - res * (p / x) % p, x = p % x;
  return res;
}

inline long s(long const n, long const c) {
  if (c == 1) return n * (n + 1) / 2 % p;
  long const v0((c - pow(c, n + 1)) * inv(1 - c) % p);
  long const v1(n * pow(c, n + 1) % p);
  // cerr << __func__ << ' ' << n << ' ' << c << ' ';
  // cerr << (v0 - v1) * inv(1 - c) % p << endl;
  return (v0 - v1) * inv(1 - c) % p;
}
inline void tr(long& a, long& b, long l, long r, long c) {
  // cerr << __func__ << ' ' << a << ' ' << b << ' ';
  // cerr << l << ' ' << r << ' ' << c << endl;
  b = (b + a * c % p * pow(c + 1, l - 1) % p * s(r - 1, inv(c + 1))) % p;
  b = (b - a * c % p * pow(c + 1, l - 1) % p * s(l - 1, inv(c + 1))) % p;
  a = a * pow(inv(c + 1), r - l) % p;
  // cerr << "--->" << ' ' << a << ' ' << b << endl;
}

int main() {
  long t;
  cin >> t;
  for (long i(1); i <= t; ++i) {
    cin >> n >> m >> p;
    if (n == 1) {
      cout << 0 << endl;
      continue;
    }
    long a(1), b(0);
    if (m < n - 1) {
      tr(a, b, 1, n - m, 0);
      tr(a, b, n - m, n, 1);
    } else {
      long const q((m - (n - 2)) / (n - 1) % p);
      long const r((m - (n - 2)) % (n - 1) % p);
      long const p(n - r);
      tr(a, b, 1, 2, q + 0);
      tr(a, b, 2, p, q + 1);
      tr(a, b, p, n, q + 2);
    }
    b = (b + (n - 1) * a) % p;
    // cerr << a << ' ' << b << endl;
    cout << (b + p) * inv(a) % p << endl;
  }
}