#include <bits/stdc++.h>
static char buf[1000000],*paa=buf,*pd=buf;
static char buf2[1000000],*pp=buf2;
#define getchar() paa==pd&&(pd=(paa=buf)+fread(buf,1,1000000,stdin),paa==pd)?EOF:*paa++
inline void pc(char ch){
	if(pp-buf2==1000000) fwrite(buf2,1,1000000,stdout),pp=buf2;
	*pp++=ch;
}
inline void pcc(){
	fwrite(buf2,1,pp-buf2,stdout);
	pp=buf2;
}
inline int read(void) {
	int w=1;
	register int x(0);register char c(getchar());
	while(c<'0'||c>'9') {if(c=='-') w=-1;c=getchar();}
	while(c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return w*x;
}
void write(int x) {
	if(x<0) pc('-'),x=-x;
	static int sta[20];
	int top=0;
	do {
		sta[top++]=x%10,x/=10;
	} while(x);
	while(top) pc(sta[--top]+48);
}
void we(int x) {
	write(x);
	pc('\n');
}
#define pai 3.141592653589793238462643383279502884197169399375105820
#define md(x) ((x)>=MOD?(x)-MOD:(x))
int rev[4000009],p[1000009],iv[1000009];
const int g=3,MOD=998244353,MAX=1000005;
struct complex{
    double x,y;
    complex(double x=0,double y=0):x(x),y(y){}
    inline complex operator+(complex b) {
		return complex(x+b.x,y+b.y);
    }
    inline complex operator-(complex b) {
        return complex(x-b.x,y-b.y);
    }
    inline complex operator*(complex b) {
        complex res;
		res.x=x*b.x-y*b.y;
		res.y=x*b.y+y*b.x;
		return res;
    }
};
inline int pw(int a,int b) {
	int as=1;
	while(b) {
		if(b&1) as=1ll*as*a%MOD;
		a=1ll*a*a%MOD;
		b>>=1;
	}
	return as;
}
inline int ni(int a) {
	return pw(a,MOD-2);
}
void binom(void) {
	p[0]=1;
	for(int i=1;i<=MAX;i++) p[i]=1ll*p[i-1]*i%MOD;
	iv[MAX]=ni(p[MAX]);
	for(int i=MAX-1;i>=0;i--) iv[i]=1ll*iv[i+1]*(i+1)%MOD;
}
struct poly{
	std::vector<int> f;
	inline int sz(void) {
		return f.size()-1;
	}
	inline void input(int sz) {
		for(int i=0;i<=sz;i++) {
			f.push_back(read());
		}
	}
	inline void output() {
		for(int i=0;i<f.size();i++) write(f[i]),pc(' ');
		pcc();
	}
	inline void FFT(std::vector<complex>& a,int t,int lim) {
		for(int i=1;i<=lim;i++) if(i<rev[i]) std::swap(a[i],a[rev[i]]);
		complex rt,w,x,y;
		for(int mid=1;mid<lim;mid<<=1) {
			int r=mid<<1;
			rt=complex(cos(pai/mid),t*sin(pai/mid));
			for(int j=0;j<lim;j+=r) {
				w=complex(1,0);
				for(int k=0;k<mid;k++) {
					x=a[j|k];
					y=w*a[j|k|mid];
					a[j|k]=x+y;
					a[j|k|mid]=x-y;
					w=w*rt;
				}
			}
		}
		if(t==1) return;
		for(int i=0;i<=lim;i++) a[i].x/=lim,a[i].y/=lim; 
	}
	inline void NTT(std::vector<int>& a,int t,int lim,int g) {
		for(int i=1;i<=lim;i++) if(i<rev[i]) std::swap(a[i],a[rev[i]]);
		int rt,w,x,y;
		for(int mid=1;mid<lim;mid<<=1) {
			int r=mid<<1;
			rt=pw(g,(MOD-1)/r);
			if(t==-1) rt=ni(rt);
			for(int j=0;j<lim;j+=r) {
				w=1;
				for(int k=0;k<mid;k++) {
					x=a[j|k];
					y=1ll*w*a[j|k|mid]%MOD;
					a[j|k]=((x+y>=MOD)?(x+y-MOD):(x+y));
					a[j|k|mid]=((x-y<0)?(x-y+MOD):(x-y));
					w=1ll*w*rt%MOD;
				}
			}
		}
		if(t==1) return;
		int qq=ni(lim);
		for(int i=0;i<=lim;i++) a[i]=1ll*a[i]*qq%MOD;
	}
/*	inline poly operator*(poly a) {
		int n=std::max(f.size()-1,a.f.size()-1),lim=1,x=-1;
		while(lim<=(n<<1)) {
			lim<<=1;
			x++;
		}
		for(int i=1;i<=lim;i++) rev[i]=(rev[i>>1]>>1)|((i&1)<<x);
		std::vector<complex> tp;tp.resize(lim+1);
		for(int i=0;i<f.size();i++) tp[i].x=f[i];
		for(int i=0;i<a.f.size();i++) tp[i].y=a.f[i];	
		FFT(tp,1,lim);
		for(int i=0;i<=lim;i++) tp[i]=tp[i]*tp[i];
		FFT(tp,-1,lim);
		poly c;c.f.resize(f.size()-1+a.f.size());
		for(int i=0;i<f.size()-1+a.f.size();i++) {
			c.f[i]=tp[i].y/2+0.5;
			c.f[i]%=MOD;
		}
		return c;
	}*/
	inline poly operator+(poly a) {
		int n=std::max(f.size()-1,a.f.size()-1);
		a.f.resize(n+1);
		for(int i=0;i<f.size();i++) {
			a.f[i]=(a.f[i]+f[i])>=MOD?(a.f[i]+f[i])-MOD:(a.f[i]+f[i]); 
		}
		return a;
	}
	inline poly operator+(int x) {
		poly a;a.f=f;a.f[0]=(a.f[0]+x>=MOD)?(a.f[0]+x-MOD):(a.f[0]+x); 
		return a;
	}
	inline poly operator-(poly a) {
		int n=std::max(f.size()-1,a.f.size()-1);
		a.f.resize(n+1);
		for(int i=0;i<f.size();i++) {
			a.f[i]=(f[i]-a.f[i])<0?(f[i]-a.f[i])+MOD:(f[i]-a.f[i]); 
		}
		for(int i=f.size();i<=n;i++) {
			a.f[i]=(a.f[i]==0?0:(MOD-a.f[i]));
		}
		return a;
	}
	inline poly operator-(int x) {
		poly a;a.f=f;a.f[0]=(a.f[0]-x<0)?(a.f[0]-x+MOD):(a.f[0]-x); 
		return a;
	}
	inline poly operator*(int a) {
		poly tmp;tmp.f=f;
		for(int i=0;i<f.size();i++) tmp.f[i]=1ll*tmp.f[i]*a%MOD;
		return tmp;
	} 
	inline poly operator/(int a) {
		poly tmp;tmp.f=f;a=ni(a);
		for(int i=0;i<f.size();i++) tmp.f[i]=1ll*tmp.f[i]*a%MOD;
		return tmp;
	} 
	inline poly operator*(poly a) {
		if(a.f.size()<=16||f.size()<=16) {
			poly tp;
			tp.f.resize(f.size()-1+a.f.size());
			for(int i=0;i<f.size();i++) {
				for(int j=0;j<a.f.size();j++) {
					int as=1ll*f[i]*a.f[j]%MOD;
					tp.f[i+j]=(tp.f[i+j]+as)>=MOD?(tp.f[i+j]+as)-MOD:(tp.f[i+j]+as);
				}
			}
			return tp;
		}
		int n=std::max(f.size()-1,a.f.size()-1),lim=1,x=-1,ss=f.size()-1+a.f.size();
		while(lim<=(n<<1)) {
			lim<<=1;
			x++;
		}
		for(int i=1;i<=lim;i++) rev[i]=(rev[i>>1]>>1)|((i&1)<<x);
		poly tmp;tmp.f=f;
		tmp.f.resize(lim+1);
		a.f.resize(lim+1);
		NTT(tmp.f,1,lim,g);
		NTT(a.f,1,lim,g);
		for(int i=0;i<=lim;i++) tmp.f[i]=1ll*tmp.f[i]*a.f[i]%MOD;
		NTT(tmp.f,-1,lim,g);
		tmp.f.resize(ss);
		return tmp;		
	} 
	inline poly dao(void) {
		poly tmp;
		for(int i=1;i<f.size();i++) {
			tmp.f.push_back(1ll*f[i]*i%MOD);
		}
		return tmp;
	} 
	inline poly ivdao(void) {
		poly tmp;
		tmp.f.push_back(0);
		for(int i=0;i<f.size();i++) {
			tmp.f.push_back(1ll*f[i]*iv[i+1]%MOD*p[i]%MOD);
		}
		return tmp;
	} 
	inline poly sqr(void) {
		if(f.size()<=16) {
			poly tp;
			tp.f.resize(f.size()-1+f.size());
			for(int i=0;i<f.size();i++) {
				for(int j=0;j<f.size();j++) {
					int as=1ll*f[i]*f[j]%MOD;
					tp.f[i+j]=(tp.f[i+j]+as)>=MOD?(tp.f[i+j]+as)-MOD:(tp.f[i+j]+as);
				}
			}
			return tp;
		}
		int n=f.size()-1,lim=1,x=-1,ss=f.size()-1+f.size();
		while(lim<=(n<<1)) {
			lim<<=1;
			x++;
		}
		for(int i=1;i<=lim;i++) rev[i]=(rev[i>>1]>>1)|((i&1)<<x);
		poly tmp;tmp.f=f;
		tmp.f.resize(lim+1);
		NTT(tmp.f,1,lim,g);
		for(int i=0;i<=lim;i++) tmp.f[i]=1ll*tmp.f[i]*tmp.f[i]%MOD;
		NTT(tmp.f,-1,lim,g);
		tmp.f.resize(ss);
		return tmp;		
	} 
	inline poly inv(void) {
		poly tmp,tmp3;tmp3.f.push_back(f[0]);
		tmp.f.push_back(ni(f[0]));
		int s=f.size(),qq=0,a[30];
		a[++qq]=s;
		while(s>1) a[++qq]=(s+1)/2,s=a[qq];
		for(int i=qq-1;i>=1;i--) {
			for(int j=a[i+1];j<a[i];j++) tmp3.f.push_back(f[j]);
			tmp=tmp*2-tmp3*tmp.sqr();
			tmp.f.resize(a[i]);
		}
		return tmp;
	}
	inline poly ln(void) {
		poly qq=(this->dao()*this->inv()).ivdao();
		qq.f.resize(f.size());
		return qq;
	}
	inline poly exp(void) {
		poly tmp,tmp3;tmp3.f.push_back(f[0]);
		tmp.f.push_back(1);
		int s=f.size(),a[30],qq=0;
		a[++qq]=s;
		while(s>1) a[++qq]=(s+1)/2,s=a[qq];
		for(int i=qq-1;i>=1;i--) {
			for(int j=a[i+1];j<a[i];j++) tmp3.f.push_back(f[j]);
			poly tmp4=tmp;
			tmp.f.resize(a[i]);
			tmp=tmp4*(tmp3-tmp.ln()+1);
			tmp.f.resize(a[i]);
		}
		return tmp;
	}
	inline poly sqrt(void) {
		poly tmp,tmp3;tmp3.f.push_back(f[0]);
		tmp.f.push_back(1);
		int s=f.size(),a[30],qq=0;
		a[++qq]=s;
		while(s>1) a[++qq]=(s+1)/2,s=a[qq];
		for(int i=qq-1;i>=1;i--) {
			for(int j=a[i+1];j<a[i];j++) tmp3.f.push_back(f[j]);
			poly tmp5=tmp3+tmp.sqr();
			tmp.f.resize(a[i]);
			tmp=tmp5*(tmp*2).inv();
			tmp.f.resize(a[i]);
		}
		return tmp;		
	}
	inline poly ksm(int x) {
		if(x<64) {
			poly a,b;
			a.f=f;
			b.f.push_back(1); 
			while(x) {
				if(x&1) b=b*a;
				b.f.resize(f.size());
				x>>=1;
				a=a.sqr();
				a.f.resize(f.size());
			}
			return b;
		}
		poly tmp;tmp.f=f;
		tmp=tmp.ln();
		tmp=tmp*x;
		tmp=tmp.exp();
		return tmp;
	}
};
poly sv(int l,int r,int m) {
	if(l==r) {
		if(l==m) {
			poly tp;tp.f.push_back(2);
			return tp;
		}
		poly tp;tp.f.resize(l+1);
		tp.f[l]=tp.f[0]=1;
		return tp;
	}
	poly tp=sv(l,(l+r)>>1,m)*sv((l+r+2)>>1,r,m);
	if(tp.f.size()>m) {
		for(int i=m;i<tp.f.size();i++) {
			tp.f[i-m]=md(tp.f[i-m]+tp.f[i]);
		}
		tp.f.resize(m);
	}
	return tp;
}
inline int sv(long long a,long long m) {
	a/=m;
	poly tp;
	tp.f.push_back(2);
	for(int i=1;i<m;i++) tp.f.push_back(0);
	for(int i=1;i<m;i++) {
		poly as=tp;
		for(int j=0;j<m;j++) {
			as.f[(i+j>=m)?(i+j-m):(i+j)]=md(as.f[(i+j>=m)?(i+j-m):(i+j)]+tp.f[j]);
		}
		tp=as;
	}
	poly as;
	as.f.push_back(1);
	printf("ok\n");
	while(a) {
		printf("%lld\n",a);
		if(a&1) {
			as=as*tp;
			if(as.f.size()>m) {
				for(int i=m;i<as.f.size();i++) {
					as.f[i-m]=md(as.f[i-m]+as.f[i]);
				}
				as.f.resize(m);
			}
		}
		tp=tp.sqr();
		if(tp.f.size()>m) {
			for(int i=m;i<tp.f.size();i++) {
				tp.f[i-m]=md(tp.f[i-m]+tp.f[i]);
			}
			tp.f.resize(m);
		}
		a>>=1;
	}
	return as.f[0];
}
signed main(void) {
	int T,m,a,b,c,d;
	scanf("%d",&T);
	binom();
	while(T--) {
		scanf("%d %d %d %d %d",&m,&a,&b,&c,&d);
		if(a==0||b==0) a=0;
		else a=std::__gcd(a,b);
		if(c==0||d==0) c=0;
		else c=std::__gcd(c,d);
		long long a1=1,a2=1;
		for(int i=1;i<=a;i++) a1=a1*m;
		for(int i=1;i<=c;i++) a2=a2*m;
		if(a==0) a1=0;
		printf("%d\n",sv(a2-a1,m));
	}
}