N,B,Q=map(int,input().split())
A1=[]
A2=[]
T=[]
T2=[]
for i in range(N):
  v=[]
  v2=[]
  for j in range(2):
    a,b=map(int,input().split())
    a%=B
    b%=B
    v.append([a,b])
    v2.append([0,0])
  v2[1][1]=v[0][0]
  v2[0][0]=v[1][1]
  v2[0][1]=-v[0][1]
  v2[1][0]=-v[1][0]
  A1.append(v)
  A2.append(v2)
  if i==0:
    T.append(v.copy())
    T2.append(v2.copy())
for w in range(1,N):
  r=[[0]*2 for i in range(2)]
  r2=[[0]*2 for i in range(2)]
  for i in range(2):
    for j in range(2):
      for k in range(2):
        r[i][j]+=A1[w][i][k]*T[w-1][k][j]
        r[i][j]%=B
        r2[i][j]+=T2[w-1][i][k]*A2[w][k][j]
        r2[i][j]%=B
  T.append(r.copy())
  T2.append(r2.copy())
for _ in range(Q):
  l,r,x,y=map(int,input().split())
  x%=B
  y%=B
  if r==0:
    print(x,y)
    continue
  p=[[0]*2 for i in range(2)]
  for i in range(2):
    for j in range(2):
      if l==0:
        p[i][j]=T[r-1][i][j]
      for k in range(2):
        if l>0:
          p[i][j]+=T[r-1][i][k]*T2[l-1][k][j]
          p[i][j]%=B
  z=p[0][0]*x+p[0][1]*y
  w=p[1][0]*x+p[1][1]*y
  z%=B
  w%=B
  print(z,w)