#include #include #include using namespace std; long long N, D, sum; vectorE1, E2; int main() { cin >> N >> D; for (int i = 1; i <= N; i++) { for (int j = 1; j <= N; j++) { E1.push_back(i*i + j * j); E2.push_back(i*i - j * j); } } sort(E1.begin(), E1.end()); sort(E2.begin(), E2.end()); for (int i = 0; i < E1.size(); i++) { int target = D - E1[i]; int pos1 = lower_bound(E2.begin(), E2.end(), target) - E2.begin(); int pos2 = lower_bound(E2.begin(), E2.end(), target + 1) - E2.begin(); sum += 1LL * (pos2 - pos1); } cout << sum << endl; return 0; }