#pragma GCC optimize("O3") #pragma GCC optimize(2) #pragma GCC optimize(3) #pragma GCC optimize("Ofast") #pragma GCC optimize("inline") #pragma GCC optimize("-fgcse") #pragma GCC optimize("-fgcse-lm") #pragma GCC optimize("-fipa-sra") #pragma GCC optimize("-ftree-pre") #pragma GCC optimize("-ftree-vrp") #pragma GCC optimize("-fpeephole2") #pragma GCC optimize("-ffast-math") #pragma GCC optimize("-fsched-spec") #pragma GCC optimize("unroll-loops") #pragma GCC optimize("-falign-jumps") #pragma GCC optimize("-falign-loops") #pragma GCC optimize("-falign-labels") #pragma GCC optimize("-fdevirtualize") #pragma GCC optimize("-fcaller-saves") #pragma GCC optimize("-fcrossjumping") #pragma GCC optimize("-fthread-jumps") #pragma GCC optimize("-funroll-loops") #pragma GCC optimize("-freorder-blocks") #pragma GCC optimize("-fschedule-insns") #pragma GCC optimize("inline-functions") #pragma GCC optimize("-ftree-tail-merge") #pragma GCC optimize("-fschedule-insns2") #pragma GCC optimize("-fstrict-aliasing") #pragma GCC optimize("-falign-functions") #pragma GCC optimize("-fcse-follow-jumps") #pragma GCC optimize("-fsched-interblock") #pragma GCC optimize("-fpartial-inlining") #pragma GCC optimize("no-stack-protector") #pragma GCC optimize("-freorder-functions") #pragma GCC optimize("-findirect-inlining") #pragma GCC optimize("-fhoist-adjacent-loads") #pragma GCC optimize("-frerun-cse-after-loop") #pragma GCC optimize("inline-small-functions") #pragma GCC optimize("-finline-small-functions") #pragma GCC optimize("-ftree-switch-conversion") #pragma GCC optimize("-foptimize-sibling-calls") #pragma GCC optimize("-fexpensive-optimizations") #pragma GCC optimize("inline-functions-called-once") #pragma GCC optimize("-fdelete-null-pointer-checks") #pragma GCC optimize("Ofast") #include using namespace std; vector v[333]; int main(void) { cin.tie(0); ios::sync_with_stdio(false); long long int A; cin >> A; for(int i=1;i<=4444;i++) { long long int val = i; val*=i; val*=i; val%=333; v[val].push_back(i); } long long int threshold = 4444; threshold*=4444; threshold*=4444; threshold*=4444; long long int res = 0; for(long long int x=1;x<=4444;x++) { for(long long int y=1;y<=4444;y++) { long long int val = (x*x*x*x) + (y*y*y*y); if(val > threshold) break; long long int temp = threshold - val; long long int z = pow(temp,0.25); while(1) { if(z*z*z*z > temp) { z-=1; break; } z++; } long long int val2 = (x*x*x)%333; val2 += ((y*y*y)%333); val2%=333; val2 = A - val2; val2%=333; if(val2 < 0) val2 += 333; int idx = upper_bound(v[val2].begin(),v[val2].end(),z) - v[val2].begin(); res += idx; } } cout << res << '\n'; return 0; }