結果

問題 No.800 四平方定理
ユーザー HaraTakashiHaraTakashi
提出日時 2022-01-29 11:04:09
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 80 ms / 2,000 ms
コード長 1,428 bytes
コンパイル時間 2,111 ms
コンパイル使用メモリ 201,128 KB
実行使用メモリ 65,628 KB
最終ジャッジ日時 2023-08-30 08:26:34
合計ジャッジ時間 4,298 ms
ジャッジサーバーID
(参考情報)
judge15 / judge12
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,376 KB
testcase_01 AC 3 ms
4,376 KB
testcase_02 AC 3 ms
4,376 KB
testcase_03 AC 3 ms
4,412 KB
testcase_04 AC 2 ms
4,376 KB
testcase_05 AC 2 ms
4,376 KB
testcase_06 AC 3 ms
4,444 KB
testcase_07 AC 3 ms
4,424 KB
testcase_08 AC 3 ms
4,380 KB
testcase_09 AC 3 ms
4,440 KB
testcase_10 AC 30 ms
34,228 KB
testcase_11 AC 34 ms
38,012 KB
testcase_12 AC 35 ms
37,096 KB
testcase_13 AC 30 ms
35,424 KB
testcase_14 AC 33 ms
37,860 KB
testcase_15 AC 37 ms
37,604 KB
testcase_16 AC 33 ms
38,328 KB
testcase_17 AC 33 ms
38,328 KB
testcase_18 AC 37 ms
37,896 KB
testcase_19 AC 36 ms
38,112 KB
testcase_20 AC 1 ms
4,380 KB
testcase_21 AC 2 ms
4,376 KB
testcase_22 AC 38 ms
38,316 KB
testcase_23 AC 80 ms
65,444 KB
testcase_24 AC 71 ms
65,628 KB
testcase_25 AC 78 ms
65,596 KB
testcase_26 AC 1 ms
4,376 KB
testcase_27 AC 1 ms
4,380 KB
testcase_28 AC 71 ms
65,372 KB
testcase_29 AC 77 ms
65,428 KB
testcase_30 AC 72 ms
65,212 KB
testcase_31 AC 67 ms
61,008 KB
testcase_32 AC 75 ms
65,620 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using P = pair<int, int>;
using T = tuple<int, int, int>;
#define al(a) a.begin(), a.end()
#define ral(a) a.rbegin(), a.rend()
#define sz(a) (int)a.size()
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define rrep(i, n) for (int i = 1; i <= (n); ++i)
#define drep(i, n) for (int i = (n)-1; i >= 0; --i)
#define db(a, b) cout << #a << ": " << a << " " << #b << ": " << b << endl;

int main() {
  ll n, d;
  cin >> n >> d;
  const int MX = n * n + n * n;
  vector<int> cnt1(MX + 1), cnt2(MX + 1);
  for (int x = 1; x <= n; ++x) {
    for (int y = 1; y <= n; ++y) {
      int s = x * x + y * y;
      cnt1[s]++;
    }
  }

  ll ans = 0;
  for (int w = 1; w <= n; ++w) {
    for (int z = 1; z <= n; ++z) {
      int t = w * w - z * z + d;
      if (t < 0) break;
      if (t > MX) continue;
      cnt2[t]++;
    }
  }

  rrep(i, MX) ans += ll(cnt1[i]) * cnt2[i];
  cout << ans << endl;
  return 0;
}
// int main() {
//   ll n, d;
//   cin >> n >> d;
//   unordered_map<int, int> mp;
//   for (int x = 1; x <= n; ++x) {
//     for (int y = 1; y <= n; ++y) {
//       int s = x * x + y * y;
//       mp[s]++;
//     }
//   }

//   ll ans = 0;
//   for (int z = 1; z <= n; ++z) {
//     for (int w = 1; w <= n; ++w) {
//       int t = w * w - z * z + d;
//       if (t < 0) break;
//       ans += mp[t];
//     }
//   }
//   cout << ans << endl;
//   return 0;
// }
0