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

int g, per;
vector<int> st;
long long num(long long Y, int K, int H) {
    int x = Y % H;
    long long mx = Y / K;
    if(st[x] == -1 || st[x] > mx) return 0;
    return (mx - st[x]) / per + 1;
    
}

long long MOD = 1e9 + 7;

int main() {
    int T; cin >> T;
    while(T--) {
        int N, K, H;
        long long Y;
        cin >> N >> K >> H >> Y;
        if(N < K) swap(N, K);
        if(N < H) swap(N, H);

        g = __gcd(K, H), per = H / g;
        st.clear(), st.resize(H + 1, -1);
        for(int i = 0; i < per; i++) {
            st[((K * 1LL * i) % H + H) % H] = i;
        }
        long long ans = 0LL;
        for(int i = 0; i <= Y / N; i++) {
            ans += num(Y - N * 1LL * i, K, H);
        }
        ans %= MOD;
        cout << ans << "\n";
    }
}