n,m,p,q = map(int,input().split())
xy = [[0]*m for i in range (n)]
for i in range (n):
    for j in range (m):
        xy[i][j] = ((i+1)*(j+1))%p
for i in range (q):
    x,f = map(int,input().split())
    print(xy[x].count(f))