import numpy as np
def main():
	n,b,q=map(int,input().split())
	arr = []
	inv_arr = []
	for _ in range(n):
		a11,a12=map(int,input().split())
		a21,a22=map(int,input().split())
		a = np.matrix([[a11,a12],[a21,a22]])
		arr.append(a)
		inv_a = np.matrix([[a22,-a12],[-a21,a11]])
		inv_arr.append(inv_a)
	
	mod_b = lambda x: x % b

	acc = [np.eye(2, dtype=int)]
	inv_acc = [np.eye(2, dtype=int)]
	for i in range(n):
		a = arr[i] * acc[-1]
		a = mod_b(a)
		acc.append(a)
		
		inv_a = inv_acc[-1] * inv_arr[i]
		inv_a = mod_b(inv_a)
		inv_acc.append(inv_a)

	for _ in range(q):
		l,r,x,y=map(int,input().split())
		a1 = acc[r]
		inv_a2 = inv_acc[l]
		a = a1 * inv_a2
		a = mod_b(a)
		zw = a * np.matrix([[x],[y]])
		z = zw[0,0]
		w = zw[1,0]
		
		z = z%b
		w = w%b
		print(z,w)

if __name__ == "__main__":
	main()