#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define rep(i, n) for (int i = 0; i < (int)(n); ++i) //#define cerr if(false) cerr #ifdef DEBUG #define show(...) cerr << #__VA_ARGS__ << " = ", debug(__VA_ARGS__); #else #define show(...) 42 #endif using namespace std; using ll = long long; using pii = pair; template ostream& operator<<(ostream& os, pair a) { os << '(' << a.first << ',' << a.second << ')'; return os; } template ostream& operator<<(ostream& os, vector v) { for (auto x : v) os << x << ' '; return os; } void debug() { cerr << '\n'; } template void debug(H a, T... b) { cerr << a; if (sizeof...(b)) cerr << ", "; debug(b...); } template vector divisor(T n) { vector res; for (T i = 1; i * i <= n; i++) { if (n % i == 0) { res.push_back(i); if (i * i != n) { res.push_back(n / i); } } } sort(res.begin(), res.end()); return res; } template map prime_factor(T n) { map res; for (T i = 2; i * i <= n; i++) { while (n % i == 0) { res[i]++; n /= i; } } if (n != 1) { res[n]++; } return res; } template T gcd(T a, T b) { while (b) { a %= b; swap(a, b); } return a; } template T lcm(T a, T b) { return a / gcd(a, b) * b; } int main(){ ll n; cin >> n; auto d = divisor(n * n); ll ans = 0; set>st; for(auto x : d){ for(auto y : d){ if((n * n) % (x * y))continue; //z(z + x + y) = n * n / x / y ll b = x + y; ll c = - n * n / x / y; if(b*b-4*c<0)continue; ll r = (ll) sqrt(b*b-4*c); if(r * r < (b*b-4*c))r++; if(r * r > (b*b-4*c))r--; if(r * r != (b*b-4*c))continue; if(abs(-b+r)%2)continue; auto v = vector{x,y,(-b+r)/2}; sort(v.begin(), v.end()); st.insert(v); } } cout << st.size() << endl; }