結果
| 問題 | 
                            No.800 四平方定理
                             | 
                    
| コンテスト | |
| ユーザー | 
                             | 
                    
| 提出日時 | 2020-01-14 05:34:19 | 
| 言語 | C++14  (gcc 13.3.0 + boost 1.87.0)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 1,242 ms / 2,000 ms | 
| コード長 | 1,330 bytes | 
| コンパイル時間 | 1,058 ms | 
| コンパイル使用メモリ | 120,484 KB | 
| 実行使用メモリ | 36,304 KB | 
| 最終ジャッジ日時 | 2024-12-24 00:51:47 | 
| 合計ジャッジ時間 | 17,284 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge4 / judge5 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 3 | 
| other | AC * 30 | 
ソースコード
#include <cstdio>
#include <iostream>
#include <string>
#include <sstream>
#include <stack>
#include <algorithm>
#include <cmath>
#include <queue>
#include <map>
#include <set>
#include <cstdlib>
#include <bitset>
#include <tuple>
#include <assert.h>
#include <deque>
#include <bitset>
#include <iomanip>
#include <limits>
#include <chrono>
#include <random>
#include <array>
#include <unordered_map>
#include <functional>
#include <complex>
template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; }
template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; }
const long long MAX = 5100000;
const long long INF = 1LL << 60;
//const long long MOD = 1000000007LL;
const long long mod = 998244353LL;
using namespace std;
typedef unsigned long long ull;
typedef long long ll;
int main()
{
	/*
	cin.tie(nullptr);
	ios::sync_with_stdio(false);
	*/
	ll N, D; scanf("%lld %lld", &N, &D);
	vector<ll> t;
	for (ll i = 1; i <= N; i++) {
		for (ll j = 1; j <= N; j++) {
			t.push_back(i*i + j * j);
		}
	}
	sort(t.begin(), t.end());
	ll res = 0;
	for (ll i = 1; i <= N; i++) {
		for (ll j = 1; j <= N; j++) {
			ll tmp = j * j - i * i + D;
			res += upper_bound(t.begin(), t.end(), tmp) - lower_bound(t.begin(), t.end(), tmp);
		}
	}
	cout << res << endl;
	return 0;
}