#pragma GCC optimize("Ofast") // #pragma GCC optimize("unroll-loops") // #pragma GCC push_options // #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native") #include #include #include using namespace::std; __attribute__((constructor))void init(){cin.tie(0);ios::sync_with_stdio(false);cout< #include #include #include // #include // #include // namespace mp = boost::multiprecision; // typedef mp::number> cdouble; // typedef mp::cpp_int cint; templateusing pbds=__gnu_pbds::tree,__gnu_pbds::rb_tree_tag,__gnu_pbds::tree_order_statistics_node_update>; templateusing pbds_map=__gnu_pbds::tree,__gnu_pbds::rb_tree_tag,__gnu_pbds::tree_order_statistics_node_update>; templateusing hash_map=__gnu_pbds::gp_hash_table; templateusing pqueue =__gnu_pbds::priority_queue,__gnu_pbds::rc_binomial_heap_tag>; typedef long long lint; #define INF (1LL<<60) #define IINF (1<<30) #define EPS (1e-10) #define MOD 1000000007LL //#define MOD 998244353LL typedef vector vec; typedef vector> mat; typedef vector>> mat3; typedef vector svec; typedef vector> smat; templateinline void numout(T t){bool f=0;for(auto i:t){cout<<(f?" ":"")<inline void numout2(T t){for(auto i:t)numout(i);} templateinline void output(T t){bool f=0;for(auto i:t){cout<<(f?" ":"")<inline void output2(T t){for(auto i:t)output(i);} templateinline void _output(T t){bool f=0;for(lint i=0;iinline void _output2(T t){for(lint i=0;i=0;--i) #define rrepi(i,a,b) for(lint i=lint(b)-1;i>=lint(a);--i) #define irep(i) for(lint i=0;;++i) #define all(n) begin(n),end(n) #define dist(a,b,c,d) sqrt(pow(a-c,2)+pow(b-d,2)) inline lint gcd(lint A,lint B){return B?gcd(B,A%B):A;} inline lint lcm(lint A,lint B){return A/gcd(A,B)*B;} // inline cint cgcd(cint A,cint B){return B?cgcd(B,A%B):A;} // inline cint clcm(cint A,cint B){return A/cgcd(A,B)*B;} //inline bool chmin(auto& s,const auto& t){bool res=s>t;s=min(s,t);return res;} //inline bool chmax(auto& s,const auto& t){bool res=s dx={1,0,-1,0,1,1,-1,-1}; const vector dy={0,1,0,-1,1,-1,1,-1}; #define SUM(v) accumulate(all(v),0LL) //auto call=[&](auto f,auto... args){return f(f,args...);}; lint bs(lint mn,lint mx,function func){ mn--; mx++; while(mx-mn>1){ lint mid=(mn+mx)/2; if(!func(mid))mx=mid; else mn=mid; } return mn; } int main(){ lint n; cin>>n; vec v; for(lint i=1;i<=n;i++){ if(n*n%i)continue; v.push_back(i); } lint ans=0; rep(i,v.size())rep(j,i+1){ if((n*n)%(v[i]*v[j]))continue; lint k=(n*n)/(v[i]*v[j]); lint t=v[i]+v[j]; lint d2=t*t+k*4; lint d=bs(0,d2,[&](lint x){return x*x<=d2;}); if(d*d!=d2)continue; if((-t+d)%2)continue; if((-t+d)/2