def floor_sum(n, m, a, b): """ return \\sum_{i=0}^{n-1} ((a*i+b)//m) """ ret = 0 while True: if a >= m: ret += (n - 1) * n // 2 * (a // m) a %= m if b >= m: ret += n * (b // m) b %= m y_max = (a * n + b) // m if y_max == 0: return ret x_max = y_max * m - b ret += (n - (x_max + a - 1) // a) * y_max n, m, a, b = y_max, a, m, -x_max % a MOD = 998244353 def solve(): n, m, l, r = map(int, input().split()) ans = floor_sum(r + 1, n - 1, 1, 0) ans -= floor_sum(l, n - 1, 1, 0) l, r = m - r, m - l ans += floor_sum(r + 1, n - 1, 1, 0) ans -= floor_sum(l, n - 1, 1, 0) ans += r - l + 1 print(ans % MOD) for _ in range(int(input())): solve()