#include using namespace std; typedef signed long long ll; #undef _P #define _P(...) (void)printf(__VA_ARGS__) #define FOR(x,to) for(x=0;x<(to);x++) #define FORR(x,arr) for(auto& x:arr) #define ITR(x,c) for(__typeof(c.begin()) x=c.begin();x!=c.end();x++) #define ALL(a) (a.begin()),(a.end()) #define ZERO(a) memset(a,0,sizeof(a)) #define MINUS(a) memset(a,0xff,sizeof(a)) //------------------------------------------------------- //A051584 ll N; set H[1010100]; ll issq(ll V) { ll a=sqrt(V); if(a*a==V) return a; if((a-1)*(a-1)==V) return a-1; if((a+1)*(a+1)==V) return a+1; return -1; } void solve() { int i,j,k,l,r,x,y; string s; cin>>N; ll M=N*N; vector A; for(i=1;1LL*i*i<=M;i++) if(M%i==0) { A.push_back(i); if(1LL*i*i!=M) A.push_back(M/i); } sort(ALL(A)); int num=0; FOR(x,A.size()) { if(A[x]>N) break; for(y=x;yN) break; if(M%(A[x]*A[y])) continue; ll c=-M/(A[x]*A[y]); ll b=A[x]+A[y]; ll d=b*b-4*c; if(d<0) continue; ll v=issq(d); if(v<0) continue; v-=b; if(v<0) continue; if(v%2) continue; if(A[y]<=v/2 && v/2