import std;

void main () {
    long N, M, P, Q;
    {
        auto buf = readln.split.to!(long[]);
        N = buf[0], M = buf[1], P = buf[2], Q = buf[3];
    }
    solve(N, M, P, Q);
}

void solve (long N, long M, long P, long Q) {
    foreach (_; 0..Q) {
        auto buf = readln.split.to!(long[]);
        long x = buf[0], f = buf[1];

        if (x % P != 0) {
        // x^{-1}が存在
            long inv_x = modinv(x, P);
            long modular = (f * inv_x) % P;
            // y == modular(MOD P)

            if (M < modular) {
                writeln(0);
            } else if (modular != 0) {
                writeln(1 + (M - modular)/P);
            } else {
                writeln((M - modular)/P);
            }
        } else {
        // x^{-1}が存在しない(任意のyについて親交度がf_iである)
            writeln(M);
        }
    }
}

long modpow (long x, long a, const long MOD) {
    if (a == 1) {
        return x;
    }
    if (a % 2 == 1) {
        return (x * modpow(x, a-1, MOD)) % MOD;
    }
    long t = modpow(x, a/2, MOD);
    return (t * t) % MOD;
}

long modinv (long x, const long MOD) {
    return modpow(x, MOD - 2, MOD);
}