#include using namespace std; #define all(x) (x).begin(),(x).end() #define rep(i, n) for (int i = 0; i < (n); i++) #define chmin(x, y) (x) = min((x), (y)) #define chmax(x, y) (x) = max((x), (y)) #define endl "\n" typedef long long ll; typedef pair pii; typedef pair pll; template ostream &operator<<(ostream &os, const vector &vec) {os << "["; for (const auto &v : vec) {os << v << ","; } os << "]"; return os;} template ostream &operator<<(ostream &os, const pair &p) {os << "(" << p.first << ", " << p.second << ")"; return os;} void solve() { int N, D; cin >> N >> D; unordered_map A, B; for (int i = 1; i <= N; i++) { for (int j = 1; j <= N; j++) { A[i * i + j * j]++; B[i * i - j * j]++; } } ll ans = 0; for(auto b : B) { if (!A.count(b.first + D)) continue; ll c = A[b.first + D]; ans += b.second * c; } cout << ans << endl; } int main() { #ifdef LOCAL_ENV cin.exceptions(ios::failbit); #endif cin.tie(0); ios::sync_with_stdio(false); cout.setf(ios::fixed); cout.precision(16); solve(); }