#include #define rep(i, a) for (int i = (int)0; i < (int)a; ++i) #define rrep(i, a) for (int i = (int)a - 1; i >= 0; --i) #define REP(i, a, b) for (int i = (int)a; i < (int)b; ++i) #define RREP(i, a, b) for (int i = (int)a - 1; i >= b; --i) #define pb push_back #define eb emplace_back #define all(x) x.begin(), x.end() #define rall(x) x.rbegin(), x.rend() #define popcount __builtin_popcount using ll = long long; constexpr ll mod = 1e9 + 7; constexpr ll INF = 1LL << 60; template inline bool chmin(T &a, T b) { if (a > b) { a = b; return true; } return false; } template inline bool chmax(T &a, T b) { if (a < b) { a = b; return true; } return false; } ll gcd(ll n, ll m) { ll tmp; while (m != 0) { tmp = n % m; n = m; m = tmp; } return n; } ll lcm(ll n, ll m) { return abs(n) / gcd(n, m) * abs(m); //gl=xy } using namespace std; ll x,y; bool check(ll a,ll b,ll c){ return x-a*b==c&&y-a*c==b; } using P=pair; void solve(){ int t; cin>>t; while(t--){ cin>>x>>y; set>>st; ll ans=0; for(ll a=2;a*a<=min(x,y)+1;++a){ if((a*x-y)%(a*a-1))continue; ll b=(a*x-y)/(a*a-1),c=x-a*b; if(b>0&&c>0){ if(check(a,b,c)&&!st.count(pair(a,P(b,c)))){ ++ans; st.emplace(a,P{b,c}); } if(check(b,a,c)&&!st.count(pair(b,P(a,c)))){ ++ans; st.emplace(b,P{a,c}); } if(check(a,c,b)&&!st.count(pair(a,P(c,b)))){ ++ans; st.emplace(a,P{c,b}); } if(check(b,c,a)&&!st.count(pair(b,P(c,a)))){ ++ans; st.emplace(b,P{c,a}); } if(check(c,a,b)&&!st.count(pair(c,P(a,b)))){ ++ans; st.emplace(c,P{a,b}); } if(check(c,b,a)&&!st.count(pair(c,P(b,a)))){ ++ans; st.emplace(c,P{b,a}); } } } // for(auto x:st){ // cout<