#include using namespace std; using ll=long long; #define int ll #define rng(i,a,b) for(int i=int(a);i=int(a);i--) #define per(i,b) gnr(i,0,b) #define pb push_back #define eb emplace_back #define a first #define b second #define bg begin() #define ed end() #define all(x) x.bg,x.ed #define si(x) int(x.size()) #ifdef LOCAL #define dmp(x) cerr<<__LINE__<<" "<<#x<<" "< void chmax(t&a,u b){if(a void chmin(t&a,u b){if(b using vc=vector; template using vvc=vc>; using pi=pair; using vi=vc; template ostream& operator<<(ostream& os,const pair& p){ return os<<"{"< ostream& operator<<(ostream& os,const vc& v){ os<<"{"; for(auto e:v)os< void dmpr(ostream&os,const T&t,const Args&... args){ os< ostream& operator<<(ostream&os,const array&a){ return os<(all(a)); } template void print_tuple(ostream&,const T&){ } template void print_tuple(ostream&os,const T&t){ if(i)os<<","; os<(t); print_tuple(os,t); } template ostream& operator<<(ostream&os,const tuple&t){ os<<"{"; print_tuple<0,tuple,Args...>(os,t); return os<<"}"; } template void print(t x,int suc=1){ cout<>i; return i; } vi readvi(int n,int off=0){ vi v(n); rep(i,n)v[i]=read()+off; return v; } template void print(const vector&v,int suc=1){ rep(i,v.size()) print(v[i],i==int(v.size())-1?suc:2); } string readString(){ string s; cin>>s; return s; } template T sq(const T& t){ return t*t; } //#define CAPITAL void yes(bool ex=true){ #ifdef CAPITAL cout<<"YES"<<"\n"; #else cout<<"Yes"<<"\n"; #endif if(ex)exit(0); } void no(bool ex=true){ #ifdef CAPITAL cout<<"NO"<<"\n"; #else cout<<"No"<<"\n"; #endif if(ex)exit(0); } void possible(bool ex=true){ #ifdef CAPITAL cout<<"POSSIBLE"<<"\n"; #else cout<<"Possible"<<"\n"; #endif if(ex)exit(0); } void impossible(bool ex=true){ #ifdef CAPITAL cout<<"IMPOSSIBLE"<<"\n"; #else cout<<"Impossible"<<"\n"; #endif if(ex)exit(0); } constexpr ll ten(int n){ return n==0?1:ten(n-1)*10; } const ll infLL=LLONG_MAX/3; #ifdef int const int inf=infLL; #else const int inf=INT_MAX/2-100; #endif int topbit(signed t){ return t==0?-1:31-__builtin_clz(t); } int topbit(ll t){ return t==0?-1:63-__builtin_clzll(t); } int botbit(signed a){ return a==0?32:__builtin_ctz(a); } int botbit(ll a){ return a==0?64:__builtin_ctzll(a); } int popcount(signed t){ return __builtin_popcount(t); } int popcount(ll t){ return __builtin_popcountll(t); } bool ispow2(int i){ return i&&(i&-i)==i; } ll mask(int i){ return (ll(1)< void mkuni(vc&v){ sort(all(v)); v.erase(unique(all(v)),v.ed); } ll rand_int(ll l, ll r) { //[l, r] #ifdef LOCAL static mt19937_64 gen; #else static mt19937_64 gen(chrono::steady_clock::now().time_since_epoch().count()); #endif return uniform_int_distribution(l, r)(gen); } template void myshuffle(vc&a){ rep(i,si(a))swap(a[i],a[rand_int(0,i)]); } template int lwb(const vc&v,const t&a){ return lower_bound(all(v),a)-v.bg; } int gcd(int a,int b){return b?gcd(b,a%b):a;} template t pow_mod(t x,t n,t m){ t r=1; while(n){ if(n&1)r=(r*x)%m; x=(x*x)%m; n>>=1; } return r; } bool is_prime(ll n) { if (n <= 1) return false; if (n == 2) return true; if (n % 2 == 0) return false; ll d = n - 1; while (d % 2 == 0) d /= 2; for (ll a : {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37}) { if (n <= a) break; ll t = d; ll y = pow_mod<__int128_t>(a, t, n); // over while (t != n - 1 && y != 1 && y != n - 1) { y = __int128_t(y) * y % n; // flow t <<= 1; } if (y != n - 1 && t % 2 == 0) { return false; } } return true; } ll pollard_single(ll n) { auto f = [&](ll x) { return (__int128_t(x) * x + 1) % n; }; if (is_prime(n)) return n; if (n % 2 == 0) return 2; ll st = 0; while (true) { st++; ll x = st, y = f(x); while (true) { ll p = gcd((y - x + n), n); if (p == 0 || p == n) break; if (p != 1) return p; x = f(x); y = f(f(y)); } } } vc pollard(ll n) { if (n == 1) return {}; ll x = pollard_single(n); if (x == n) return {x}; vc le = pollard(x); vc ri = pollard(n / x); le.insert(le.end(), ri.begin(), ri.end()); return le; } vc to_freq(vi a){ sort(all(a)); vc res; for(auto v:a){ if(res.empty()||res.back().a!=v)res.eb(v,0); res.back().b++; } return res; } vi enum_divs(int n){ auto f=to_freq(pollard(n)); dmp(f); vi res; auto rec=[&](auto self,int i,int d)->void{ if(i==si(f))res.pb(d); else{ rep(_,f[i].b+1){ self(self,i+1,d); d*=f[i].a; } } }; rec(rec,0,1); sort(all(res)); return res; } signed main(){ cin.tie(0); ios::sync_with_stdio(0); cout<>n; vi ds=enum_divs(n*n); while(ds.back()>n)ds.pop_back(); n=n*n; int ans=0; rep(i,si(ds)){ rng(j,i,si(ds)){ int x=ds[i],y=ds[j]; if(n%(x*y)==0){ int rem=n/(x*y); int p=(x+y); int lw=j,up=si(ds); while(up-lw>1){ const int mid=(lw+up)/2; if(ds[mid]*(ds[mid]+p)<=rem)lw=mid; else up=mid; } if(ds[lw]*(ds[lw]+p)==rem){ ans++; } } } } print(ans); }