//
#include <bits/stdc++.h>
using namespace std;

#define fi first
#define se second
#define all(x) x.begin(), x.end()
#define lch (o << 1)
#define rch (o << 1 | 1)

typedef double db;
typedef long long ll;
typedef unsigned int ui;
typedef pair<int, int> pint;
typedef tuple<int, int, int> tint;

const int MOD = 1e9 + 7;
const int INF = 0x3f3f3f3f;
const ll INF_LL = 0x3f3f3f3f3f3f3f3f;

void ExGCD(ll a, ll b, ll &x, ll &y, ll &g) {
	if (!b)
		x = 1, y = 0, g = a;
	else
		ExGCD(b, a % b, y, x, g), y -= x * (a / b);
}

ll divCeil(ll upp, ll low) {
	if (upp < 0)
		return upp / low;
	else
		return (upp + low - 1) / low;
}

ll divFloor(ll upp, ll low) {
	if (upp > 0)
		return upp / low;
	else
		return (upp - low + 1) / low;
}

void solve() {
	ll a, b, c, s;
	cin >> a >> b >> c >> s;

	if (s / a < s / c) swap(a, c);
	if (s / b < s / c) swap(b, c);

	ll x, y, g;
	ExGCD(a, b, x, y, g);

	ll sum = 0;
	ll zUpp = s / c;
	for (ll z = 0; z <= zUpp; z++) {
		ll p = s - z * c;
		if (p == 0) sum++;
		if (p <= 0) continue;
		if (p % g) continue;

		ll L = divCeil(-x * (p / g), b / g);
		ll R = divFloor(y * (p / g), a / g);
		if (L <= R) sum += R - L + 1;
	}
	cout << sum % MOD << endl;
}

int main() {
	ios::sync_with_stdio(0);

	int t;
	cin >> t;
	while (t--) solve();

	return 0;
}