#include <bits/stdc++.h>

#define forn(i, n) for (int i = 0; i < int(n); ++i)
#define fore(i, l, r) for (int i = int(l); i < int(r); ++i)

using namespace std;

const int N = int(2e7) + 15;

int n, p, b0, A, B;

bool read(){
	if (!(cin >> n >> p >> b0 >> A >> B))
		return false;
	return true;
}

int add(int a, int b){
	a += b;
	if (a >= p) a -= p;
	if (a < 0) a += p;
	return a;
}

int mul(int a, int b){
	return a * 1ll * b % p;
}

int binpow(int a, int b){
	int res = 1;
	while (b){
		if (b & 1)
			res = mul(res, a);
		a = mul(a, a);
		b >>= 1;
	}
	return res;
}

vector<int> a;

int ans, cnt;

void brute(int i, int bal, int cur){
	if (i == 2 * n){
		if (bal == 0){
			cnt = add(cnt, 1);
			ans = add(ans, cur);
		}
		return;
	}
	brute(i + 1, bal + 1, add(cur, -a[i]));
	if (bal > 0) brute(i + 1, bal - 1, add(cur, a[i]));
}

vector<int> fact, rfact;

void init(){
	fact.resize(2 * n + 1);
	fact[0] = 1;
	fore(i, 1, 2 * n + 1)
		fact[i] = mul(fact[i - 1], i);
	rfact.resize(2 * n + 1);
	rfact[2 * n] = binpow(fact[2 * n], p - 2);
	for (int i = 2 * n - 1; i >= 0; --i)
		rfact[i] = mul(rfact[i + 1], i + 1);
}

int cnk(int n, int k){
	if (k < 0 || n < 0 || k > n) return 0;
	return mul(fact[n], mul(rfact[k], rfact[n - k]));
}

int cat(int n){
	if (n < 0) return 0;
	return mul(fact[2 * n], mul(rfact[n], rfact[n + 1]));
}

void solve(){
	init();
	int b = b0;
	a.assign(2 * n + 1, 0);
	fore(i, 1, 2 * n + 1){
		b = add(mul(A, b), B);
		a[i] = add(a[i - 1], add(b, 1));
	}
	a.erase(a.begin());
	
	vector<int> na;
	//int ans = 0;
	int sum = 0;
	forn(i, n){
		sum = add(sum, mul(cat(i), cat(n - 1 - i)));
		na.push_back(sum);
		na.push_back(sum);
	}
	na.pop_back();
	reverse(na.begin(), na.end());
	na.resize(n);
	//for (int x : na) cerr << x << " ";
	//cerr << endl;
	
	forn(i, n){
		int cl = add(cat(n), -na[i]);
		na[i] = add(cl, -na[i]);
	}
	
	//for (int &x : na) x = add(0, -x);
	forn(i, n) na.push_back(add(0, -na[n - 1 - i]));
	
	//for (int x : na) cerr << x << " ";
	//cerr << endl;
	
	//	int cl = add(cat(n), -op);
	//	ans = add(ans, )
	//}
	
	/*for (int i = 0; i < n; i += 2){
		int sum = 0;//cat(n);
		for (int len = 0; i + 2 + 2 * len <= 2 * n; ++len){
			//sum = add(sum, -mul(cat(i / 2), mul(cat((2 * n - 2 * len - 2 - i) / 2), cat(len))));
			//sum = add(sum, -mul(cat(len), cat(n - len - 1)));
			sum = add(sum, mul(cat(i / 2), mul(cat((2 * n - 2 * len - 2 - i) / 2), cat(len))));
		}
		cerr << sum << " ";
	}
	cerr << endl;*/
	/*vector<vector<int>> dp(2 * n + 1, vector<int>(2 * n + 1, 0));
	dp[0][0] = 1;
	forn(i, 2 * n){
		forn(j, i + 1){
			dp[i + 1][j + 1] = add(dp[i + 1][j + 1], dp[i][j]);
			if (j > 0) dp[i + 1][j - 1] = add(dp[i + 1][j - 1], dp[i][j]);
		}
	}
	forn(i, n * 2){
		//cerr << add(cat(n), -dp[2 * n - i][i]) << " ";
		cerr << dp[2 * n - i][i] << " ";
	}
	cerr << endl;*/
	
	int ans = 0;
	forn(i, 2 * n) ans = add(ans, mul(a[i], na[i]));
	ans = mul(ans, binpow(cat(n), p - 2));
	cout << ans << "\n";
}

int main(){
	cin.tie(0);
	ios::sync_with_stdio(false);
	int t;
	cin >> t;
	while (t--){
		read();
		solve();
	}
}