#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
const ll mo=998244353;

struct node {
	ll x;
	ll y;
	ll id;
};
map<int,vector<node> > g;
ll ans[100005];

int main() {
	ios::sync_with_stdio(false);
	cin.tie(0), cout.tie(0);
	ll t; cin>>t;
	for(ll i=1;i<=t;i++) {
		ll x,y;
		cin>>x>>y;
		if(x<y) swap(x,y);
		g[(x+y)/(y+1)].push_back((node){x,y,i});
	}
	for(auto q:g) {
		int i=q.first;
		vector<vector<ll> > dp[2], sdp[2];
		dp[0].resize(2002);
		sdp[0].resize(2002);
		dp[1].resize(2002);
		sdp[1].resize(2002);
		ll lim=0;
		if(i==1) lim=2000;
		else lim=2000/(i-1);
		for(ll j=0;j<2;j++) {
			for(ll k=0;k<=2000;k++) {
				dp[j][k].resize(lim+2);
				sdp[j][k].resize(lim+2);
			}
		}
		dp[0][0][0]=dp[1][0][0]=1;
		sdp[0][0][0]=sdp[1][0][0]=1;
		for(ll x=0;x<=2000;x++) {
			for(ll y=0;y<=lim;y++) {
				if(x+y==0) continue;
				if(x>=i+1) dp[0][x][y]=(sdp[1][x-1][y]-sdp[1][x-i-1][y]+mo)%mo;
				else if(x>=1) dp[0][x][y]=sdp[1][x-1][y];
				if(y>=i+1) dp[1][x][y]=(sdp[0][x][y-1]-sdp[0][x][y-i-1]+mo)%mo;
				else if(y>=1) dp[1][x][y]=sdp[0][x][y-1];
				sdp[0][x][y]=((y-1>=0?sdp[0][x][y-1]:0)+dp[0][x][y])%mo;
				sdp[1][x][y]=((x-1>=0?sdp[1][x-1][y]:0)+dp[1][x][y])%mo;
			}
		}
		for(auto p:g[i]) ans[p.id]=(dp[0][p.x][p.y]+dp[1][p.x][p.y])%mo;
	}
	for(ll i=1;i<=t;i++) cout<<ans[i]<<'\n';
	return 0;
}