結果

問題 No.800 四平方定理
ユーザー siman
提出日時 2022-04-05 09:00:09
言語 C++17(clang)
(17.0.6 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 887 bytes
コンパイル時間 2,718 ms
コンパイル使用メモリ 138,984 KB
実行使用メモリ 105,036 KB
最終ジャッジ日時 2024-11-26 08:50:01
合計ジャッジ時間 10,066 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 22 RE * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cassert>
#include <cmath>
#include <algorithm>
#include <iostream>
#include <iomanip>
#include <climits>
#include <map>
#include <queue>
#include <set>
#include <cstring>
#include <vector>

using namespace std;
typedef long long ll;

const int MAX_N = 2000;
const int MAX_D = 1000000;

ll cnt1[2 * MAX_N * MAX_N + 10];
ll cnt2[MAX_N * MAX_N + MAX_D + 10];

int main() {
  int N, D;
  cin >> N >> D;
  memset(cnt1, 0, sizeof(cnt1));
  memset(cnt2, 0, sizeof(cnt2));

  for (int x = 1; x <= N; ++x) {
    for (int y = 1; y <= N; ++y) {
      int v = x * x + y * y;
      cnt1[v]++;
    }
  }

  for (int w = 1; w <= N; ++w) {
    for (int z = 1; z <= N; ++z) {
      int v = w * w - z * z + D;
      if (v <= 0) continue;
      cnt2[v]++;
    }
  }

  ll ans = 0;

  for (int v = 1; v <= 2 * N * N; ++v) {
    ans += cnt1[v] * cnt2[v];
  }

  cout << ans << endl;

  return 0;
}
0