import sys
input = sys.stdin.readline

N,mod,Q=map(int,input().split())

def prod(A,B,k,l,m):# A:k*l,B:l*m
    C=[[None for i in range(m)] for j in range(k)]

    for i in range(k):
        for j in range(m):
            ANS=0
            for pl in range(l):
                ANS=(ANS+A[i][pl]*B[pl][j])%mod

            C[i][j]=ANS

    return C

A=[]

for i in range(N):
    x,y=map(int,input().split())
    z,w=map(int,input().split())

    A.append([[x,y],[z,w]])


S=[[[1,0],[0,1]]]

for i in range(N):
    a=A[i]

    S.append(prod(a,S[-1],2,2,2))

for tests in range(Q):
    l,r,x,y=map(int,input().split())
    if l==r:
        print(x,y)
        continue

    k=S[l]
    k0,k1=k[0]
    k2,k3=k[1]

    L= prod(S[r],[[k3,-k1],[-k2,k0]],2,2,2)

    #print(S[r],k,[[k3,-k1],[-k2,k0]],L)

    k0,k1=L[0]
    k2,k3=L[1]

    print(x*k0+y*k1,x*k2+y*k3)