N, M, P, Q = map(int, input().split())
XF = [list(map(int, input().split())) for _ in range(Q)]
for x, f in XF:
    if x%P==0:
        if f==0:
            print(M)
        else:
            print(0)
    else:
        invx = pow(x, P-2, P)
        y = f * invx % P
        cnt = (M-y)//P - (-y)//P
        print(cnt)