#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long ll; typedef unsigned int ui; const ll mod = 1000000007; typedef long double ld; const ld INF = (ld)10000000000000; typedef pair P; #define stop char nyaa;cin>>nyaa; #define rep(i,n) for(int i=0;i=0;i--) #define Rep(i,sta,n) for(int i=sta;i=1;i--) #define Rep1(i,sta,n) for(int i=sta;i<=n;i++) typedef complex Point; const ld eps = 1e-8; const ld pi = acos(-1.0); typedef pair LDP; typedef pair LP; int cnt[7000000]; int main() { int n, d; cin >> n >> d; int ans = 0; rep1(i, n) { rep1(j, n) { if (i*i - j * j + d >= 0) { cnt[i*i - j * j + d]++; } } } rep1(i, n) { rep1(j, n) { ans += cnt[i*i + j * j]; } } cout << ans << endl; //stop return 0; }