結果

問題 No.800 四平方定理
ユーザー Thistle
提出日時 2019-03-28 22:51:12
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 970 ms / 2,000 ms
コード長 2,791 bytes
コンパイル時間 1,274 ms
コンパイル使用メモリ 123,636 KB
実行使用メモリ 66,004 KB
最終ジャッジ日時 2024-10-14 17:19:30
合計ジャッジ時間 16,033 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#define _USE_MATH_DEFINES
#include<iostream>
#include<string>
#include<queue>
#include<cmath>
#include<map>
#include<set>
#include<list>
#include<iomanip>
#include<vector>
#include<random>
#include<functional>
#include<algorithm>
#include<stack>
#include<cstdio>
#include<bitset>
#include<unordered_map>
#include<climits>
#include<fstream>
using namespace std;
///////////////////library zone!!!!!!!!!!!!!!!!!!!!!!!!!!!!
typedef long long ll;
typedef long double ld;
#define all(a) (a).begin(),(a).end()
#define EPS (1e-5)
const ll Mod = 1000000007;
const ll mod = 998244353;
struct H {
    ll x, y;
    bool operator<(const H &b) const {
        if (x != b.x) return x < b.x;
        return y < b.y;
    }
    bool operator>(const H &b) const {
        if (x != b.x) return x > b.x;
        return y > b.y;
    }
    bool operator==(const H &b) const { return x == b.x&&y == b.y; }
    bool operator!=(const H &b) const { return (*this) != b; }
};
struct P {
    ll pos, cost;
    bool operator<(const P &b) const { return cost < b.cost; }
    bool operator>(const P &b) const { return cost > b.cost; }
};
struct B {
    ll to, cost;
};
struct E {
    ll from, to, cost;
    bool operator<(const E &b) const { return cost < b.cost; }
    bool operator>(const E &b) const { return cost > b.cost; }
};
template<typename T, typename U>
void chmin(T &a, U b) {
    if (a > b) a = b;
}
template<typename T, typename U>
void chmax(T &a, U b) {
    if (a < b) a = b;
}
template<typename T>
T max_0(T a) {
    if (a < 0) return 0;
    return a;
}
template<typename T>
T min_0(T a) {
    if (a > 0) return 0;
    return a;
}
ll read() {
    ll u;
    ll k = scanf("%lld", &u);
    return u;
}
ll gcd(ll i, ll j) {
    if (i > j) swap(i, j);
    if (i == 0) return j;
    return gcd(j%i, i);
}
ll mod_pow(ll x, ll n, ll p) {
    ll res = 1;
    while (n > 0) {
        if (n & 1) res = res * x % p;
        x = x * x % p;
        n >>= 1;
    }
    return res;
}//x^n%p
vector<string>split(string s, char c = ' ') {
    vector<string>vec;
    string h = "";
    for (int i = 0; i < s.size();) {
        int j = i;
        while (j < s.size() && s[j] != ' ') {
            h += s[j++];
        }
        if (h.size() > 0) vec.push_back(h);
h = "";
i = j + 1;
    }
    return vec;
}
template<typename T>
void print_vec(vector<T>a) {
    for (int i = 0; i < a.size(); i++) {
        cout << (i == 0 ? "" : " ") << a[i];
    }
    cout << endl;
}
const ll Inf = 3023372036854775807;
const int inf = 1500000000;
#define int long long
//----------------------------------------------------
int n, d;
vector<int>a, b;
signed main() {
    cin >> n >> d;
    for (int i = 1; i <= n; i++) {
        for (int j = 1; j <= n; j++) {
            a.push_back(i*i + j * j);
            b.push_back(i*i - j * j);
        }
    }
    sort(all(a)); sort(all(b));
    //x*x+y*y+z*z=w*w+d
    //x*x+y*y=w*w-z*z+d
    int ans = 0;
    for (int i = 0; i < a.size(); i++) {
        ans += upper_bound(all(b), a[i] - d) - lower_bound(all(b), a[i] - d);
    }
    cout << ans << endl;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0