def floor_sum(n, m, a, b):
	ans = 0
	while True:
		if a >= m:
			ans += (n * (n - 1) >> 1) * (a // m)
			a %= m
		if b >= m:
			ans += n * (b // m)
			b %= m
		y = a * n + b
		if y < m:
			return ans
		n, b, m, a = y // m, y % m, a, m
for _ in range(int(input())):
	n, d, m, s = map(int, input().split())
	p, dm = 1 << s, d * m
	if p != dm:
		n = min(n, d * p // abs(dm - p))
		n -= abs(floor_sum(n + 1, p, m, 0) - floor_sum(n + 1, d, 1, 0))
	print(n)