#include const double pi = 3.141592653589793238462643383279; using namespace std; //typedef //------------------------------------------ typedef vector VI; typedef vector VVI; typedef vector VS; typedef pair PII; typedef pair PLL; typedef pair TIII; typedef long long LL; typedef unsigned long long ULL; typedef vector VLL; typedef vector VVLL; //container util //------------------------------------------ #define ALL(a) (a).begin(),(a).end() #define RALL(a) (a).rbegin(), (a).rend() #define PB push_back #define MP make_pair #define SZ(a) int((a).size()) #define SQ(a) ((a)*(a)) #define EACH(i,c) for(typeof((c).begin()) i=(c).begin(); i!=(c).end(); ++i) #define EXIST(s,e) ((s).find(e)!=(s).end()) #define SORT(c) sort((c).begin(),(c).end()) //repetition //------------------------------------------ #define FOR(i,s,n) for(int i=s;i<(int)n;++i) #define REP(i,n) FOR(i,0,n) #define MOD 1000000007 #define rep(i, a, b) for(int i = a; i < (b); ++i) #define trav(a, x) for(auto& a : x) #define all(x) x.begin(), x.end() #define sz(x) (int)(x).size() typedef long long ll; typedef pair pii; typedef vector vi; #define chmin(x,y) x=min(x,y) #define chmax(x,y) x=max(x,y) //ここから編集 long long extGCD(long long a, long long b, long long &x, long long &y) { if (b == 0) { x = 1; y = 0; return a; } long long d = extGCD(b, a%b, y, x); y -= a/b * x; return d; } // 負の数にも対応した mod (a = -11 とかでも OK) inline long long mod(long long a, long long m) { return (a % m + m) % m; } // 逆元計算 (ここでは a と m が互いに素であることが必要) long long modinv(long long a, long long m) { long long x, y; extGCD(a, m, x, y); return mod(x, m); // 気持ち的には x % m だが、x が負かもしれないので } ll mul(ll a, ll b, ll c){ if((1LL<<60)/a < b) return (1LL<<60); if((1LL<<60)/(a*b) < c) return (1LL<<60); return a*b*c; } int main() { cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(5); ll N; cin >> N; vector elem; for(ll i=1; i <= N; i++){ if((N*N)%i == 0){ elem.push_back(i); } } ll ans = 0; int n = elem.size(); for(int i=0; iN*N) break; for(int k=j; k N*N) break; if(elem[i] * elem[j] * elem[k] * (elem[i]+elem[j]+elem[k]) == N*N) ans++; } } } cout << ans << endl; return 0; }