#include <cstdio>
#include <cctype>
namespace std {struct istream {
#define M (1 << 22 | 3)char buf[M], *ch = buf - 1;inline istream() { fread(buf, 1, M, stdin); }inline istream& operator >> (int &x) {while (isspace(*++ch));for (x = *ch & 15; isdigit(*++ch); ) x = x * 10 + (*ch & 15);return *this;}inline istream& operator >> (long long &x) {while (isspace(*++ch));for (x = *ch & 15; isdigit(*++ch); ) x = x * 10 + (*ch & 15);return *this;}
#undef M} cin;struct ostream {
#define M (1 << 20 | 3)char buf[M], *ch = buf - 1;inline ostream& operator << (long long x) {if (!x) {*++ch = '0'; return *this;}static int S[20], *top; top = S;while (x) {*++top = x % 10 ^ 48; x /= 10;}for (; top != S; --top) *++ch = *top;return *this;}inline ostream& operator << (const char x) {*++ch = x; return *this;}inline ~ostream() { fwrite(buf, 1, ch - buf + 1, stdout); }
#undef M} cout;
}int Tim;
long long n, a, b, mod;inline void reduce(long long &x) { x += x >> 63 & mod; }
inline long long mul(long long x, long long y) {long long res = x * y - static_cast<long long> (static_cast<long double> (x) * y / mod + 0.5) * mod;return res + (res >> 63 & mod);
}struct Matrix {long long s00, s10, s11;Matrix() { }Matrix(long long __s00, long long __s10, long long __s11) : s00(__s00), s10(__s10), s11(__s11) { }
#define M(l, r) mul(s##l, rhs.s##r)inline void operator *= (const Matrix &rhs) {static long long __s00, __s10, __s11;__s00 = M(00, 00);reduce(__s10 = M(10, 00) + M(11, 10) - mod);__s11 = M(11, 11);s00 = __s00, s10 = __s10, s11 = __s11;}
#undef M
} ;long long calc(long long n) {a %= mod, b %= mod;Matrix base(a, 1, b), res(0, 1, b);for (; n; n >>= 1, base *= base) if (n & 1) res *= base;return res.s10;
}int main() {std::cin >> Tim;while (Tim --> 0) {std::cin >> n >> a >> b >> mod;std::cout << calc(n) << '\n';}return 0;
}