#pragma GCC optimize("O3")
#include<bits/stdc++.h> 
using namespace std;
using ll=long long;
using P=pair<ll,ll>;
template<class T> using V=vector<T>; 
#define fi first
#define se second
#define all(v) (v).begin(),(v).end()
const ll inf=(1e18);
constexpr ll mod=1000000007;
ll gcd(ll a,ll b) {return b ? gcd(b,a%b):a;}
ll lcm(ll c,ll d){return c/gcd(c,d)*d;}
struct __INIT{__INIT(){cin.tie(0);ios::sync_with_stdio(false);cout<<fixed<<setprecision(15);}} __init;
template<class T> bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; }
template<class T> bool chmin(T &a, const T &b) { if (a>b) { a=b; return 1; } return 0; }
int main(){
    ll n,d;
    cin>>n>>d;
    V<ll> a(n);
    for(ll i=0;i<n;i++){
        a[i]=(i+1)*(i+1);
    }
    V<ll> dat(n*n);
    int cnt=0;
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            dat[cnt++]=a[i]+a[j];
           // cout<<dat[cnt-1]<<endl;
        }
    }
    V<ll> dat2;
    cnt=0;
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            if(a[i]-a[j]+d>=2ll){
                dat2.emplace_back(a[i]-a[j]+d);
             //   cout<<dat2[cnt++]<<endl;
            }
        }
    }
    int m=dat2.size();
    sort(all(dat));
    ll ans=0;
   // cout<<m<<endl;
    for(int i=0;i<m;i++){
        ans+=upper_bound(all(dat),dat2[i])-lower_bound(all(dat),dat2[i]);
    }
    cout<<ans<<"\n";
}