#include<bits/stdc++.h>
using namespace std;
using ll=long long;

vector<vector<ll>> mul(vector<vector<ll>> v1,vector<vector<ll>> v2,ll mod){
  ll s=v1.size();
  vector<vector<ll>> res(s,vector<ll>(s,0));
  for(ll i=0;i<s;i++){
    for(ll j=0;j<s;j++){
      for(ll k=0;k<s;k++){
        res[i][j]+=v1[i][k]*v2[k][j];
        res[i][j]+=mod;
        res[i][j]%=mod;
      }
    }
  }
  return res;
}

int main(){
  ll n,b,q;
  cin>>n>>b>>q;
  vector<vector<vector<ll>>> a(n,vector<vector<ll>>(2,vector<ll>(2)));
  for(ll i=0;i<n;i++){
    for(ll j=0;j<2;j++){
      for(ll k=0;k<2;k++){
        cin>>a[i][j][k];
      }
    }
  }
  vector<vector<vector<ll>>> sum(n,vector<vector<ll>>(2,vector<ll>(2)));
  sum[0][0][0]=a[0][0][0];
  sum[0][0][1]=a[0][0][1];
  sum[0][1][0]=a[0][1][0];
  sum[0][1][1]=a[0][1][1];
  for(ll i=1;i<n;i++){
    sum[i]=mul(a[i],sum[i-1],b);
  }
  while(q--){
    ll l,r,x,y;
    cin>>l>>r>>x>>y;
    cout<<(sum[max((ll)0,min(n-1,r-l))][0][0]%b*x+sum[max((ll)0,min(n-1,r-l))][0][1]*y%b)%b<<" ";
    cout<<(sum[max((ll)0,min(n-1,r-l))][1][0]%b*x+sum[max((ll)0,min(n-1,r-l))][1][1]*y%b)%b<<endl;
  }
}