結果

問題 No.800 四平方定理
コンテスト
ユーザー Radperia
提出日時 2019-03-17 23:58:41
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 123 ms / 2,000 ms
コード長 1,020 bytes
コンパイル時間 1,450 ms
コンパイル使用メモリ 166,224 KB
実行使用メモリ 64,512 KB
最終ジャッジ日時 2024-07-08 07:37:22
合計ジャッジ時間 3,745 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define FOR(i, a, b) for (ll i = (a); i < (b); ++i)
#define REP(i, n) for (ll i = 0, i##_len = (n); i < i##_len; ++i)
#define REPR(i, n) for (ll i = n; i >= 0; i--)
#define INF 1145141919
int dy[] = {0, 0, 1, -1};
int dx[] = {1, -1, 0, 0};
int ny, nx;
typedef long long int ll;
using namespace std;

ll gcd(ll a, ll b) { return b ? gcd(b, a % b) : a; }

ll lcm(ll m, ll n)
{
    if ((0 == m) || (0 == n))
        return 0;
    return ((m / gcd(m, n)) * n);
}

ll llpow(ll x, ll y)
{
    ll ans = 1;
    REP(i, y)
    ans *= x;
    return ans;
}

int ctoi(char c)
{
    if (c >= '0' && c <= '9')
    {
        return c - '0';
    }
    return 0;
}
ll ans = 0;
ll cnt[100000000];
////////////////////////////////////////////////////////////////////
int main()
{
    ll N, D; cin >> N >> D;

    FOR(x,1,N+1){
        FOR(y,1,N+1) cnt[x*x+y*y]++;
    }
    FOR(z,1,N+1) FOR(w,1,N+1){
        ll res = D + w*w - z*z;
        if(res>0) ans+=cnt[res];
    }

    cout << ans << endl;
    return 0;

}
0