結果
問題 | No.800 四平方定理 |
ユーザー |
![]() |
提出日時 | 2022-01-29 11:04:09 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 83 ms / 2,000 ms |
コード長 | 1,428 bytes |
コンパイル時間 | 2,497 ms |
コンパイル使用メモリ | 194,076 KB |
最終ジャッジ日時 | 2025-01-27 17:46:24 |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
#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;// }