#include <bits/stdc++.h>
using namespace std;
const int N=1e6+10;
const int Mod=998244353;
const int inv2=(Mod+1)/2;
typedef long long ll;

int qpow(int x,ll y) {
	int z=1;
	for(;y;y>>=1) {
		if(y&1) z=1ll*z*x%Mod;
		x=1ll*x*x%Mod;
	}
	return z;
}

int gcd(int x,int y) {
	if(!y) return x;
	return gcd(y,x%y);
}

int exgcd(int a,int b,int &x,int &y) {
	if(!b) {
		x=1,y=0;
		return a;
	}
	int z=exgcd(b,a%b,x,y);
	int t=y;
	y=x-a/b*y;
	x=t;
	return z;
}

int T;
ll n;
int a,b;

int solve(int a,int b,int p) {
	b=(p+b%p)%p;
	int g=gcd(a,p);
	if(b%g) return 0;
	int x,y;
	exgcd(a,p,x,y);
	x%=p/g;
	if(x<0) x+=p/g;
	x=(1ll*x*b/g)%(p/g);
	return x;
}

int main() {
	scanf("%d",&T);
	while(T--) {
		scanf("%lld%d%d",&n,&a,&b);
		int ans=0;
		for(int i=1;i<=min(n,(ll)b);i++) ans=(ans+solve(a%i,b%i,i))%Mod;
		if(n>b) {
			if(b%a==0) {
				int x=b/a;
				ans=(ans+(n-b)%Mod*x%Mod)%Mod;
			}
			else {
				for(int i=b+1;i<=b+a && i<=n;i++) {
					int z=solve(i,-b,a);
					if(b%a && !z) continue;
					ll f=(b+1ll*i*z)/a;
					ll c=(n-i)/a+1;
					ans=(ans+f%Mod*(c%Mod)%Mod+c%Mod*((c-1)%Mod)%Mod*inv2%Mod*z%Mod)%Mod;
				}
			}
		}
		printf("%d\n",ans);
	}
	return 0;
}