#include<cstdio>
#include<cassert>
#include<cstring>
#include<map>
#include<set>
#include<time.h>
#include<algorithm>
#include<stack>
#include<queue>
#include<utility>
#include<cmath>
#include<iostream>
#include<string>
#include<vector>
using namespace std;
long long M;
long pow(long long a, long long n, long long M) {
long long ret = 1, unit = a % M;
for(long long y = 0; (1ll<<y) <= n; ++y) {
if(n & (1ll << y))
ret = (ret * unit) % M;
unit = (unit * unit) % M;
}
return ret;
}
long long mod_inverse(long long a, long long M) {
return pow(a, M-2, M);
}
vector<long long> eliminate(vector<vector<long long> >& A, vector<long long>& b) {
long long n = A.size();
for(long long y = 0; y < n; ++y) {
for(long long y2 = y; y2 < n; ++y2) {
if(A[y2][y]) {
if(y != y2) {
swap(A[y], A[y2]);
swap(b[y], b[y2]);
}
break;
}
}
if(A[y][y] == 0) return vector<long long>();
long long div = mod_inverse(A[y][y], M);
for(long long x = y; x < n; ++x)
A[y][x] = (A[y][x] * div) % M;
b[y] = (b[y] * div) % M;
for(long long y2 = y+1; y2 < n; ++y2) {
long long mult = A[y2][y];
for(long long x = y; x < n; ++x)
A[y2][x] = (A[y2][x] - (A[y][x] * mult) % M + M) % M;
b[y2] = (b[y2] - (b[y] * mult) % M + M) % M;
}
}
for(long long y = n-1; y >= 0; --y) {
for(long long y2 = y-1; y2 >= 0; --y2) {
long long mult = A[y2][y];
for(long long x = y; x < n; ++x)
A[y2][x] = (A[y2][x] - (A[y][x] * mult) % M + M) % M;
b[y2] = (b[y2] - (b[y] * mult) % M + M) % M;
}
}
return b;
}
int main() {
long long t;
cin >> M >> t;
for(long long cc = 0; cc < t; ++cc) {
vector<vector<long long> > A;
{
long long a, b, c, d;
cin >> a >> b >> c >> d;
A.push_back(vector<long long>{a, -b, -c, -d});
A.push_back(vector<long long>{b, a, -d, c});
A.push_back(vector<long long>{c, d, a, -b});
A.push_back(vector<long long>{d, -c, b, a});
}
vector<long long> b = { 1, 0, 0, 0 };
// for(long long i = 0; i < 4; ++i) {
// for(long long j = 0; j < 4; ++j) {
// printf("%d ", A[i][j]);
// }
// printf("\n");
// }
vector<long long> sol = eliminate(A, b);
if(sol.empty())
printf("0 0 0 0\n");
else
printf("%d %d %d %d\n", sol[0], sol[1], sol[2], sol[3]);
}
}
