結果
問題 | No.800 四平方定理 |
ユーザー |
|
提出日時 | 2022-12-01 23:18:25 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 127 ms / 2,000 ms |
コード長 | 1,339 bytes |
コンパイル時間 | 2,513 ms |
コンパイル使用メモリ | 248,348 KB |
最終ジャッジ日時 | 2025-02-09 03:07:30 |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
#pragma region Macros#include <bits/extc++.h>using namespace std;using namespace __gnu_pbds;// using namespace __gnu_cxx;#define TO_STRING(var) # var#define pb emplace_back#define int ll#define endl '\n'using ll = long long;using ld = long double;const ld PI = acos(-1);const ld EPS = 1e-10;const int INF = 1 << 30;const ll INFL = 1LL << 62;// const int MOD = 998244353;const int MOD = 1000000007;int ceil(int x, int y) {return (x > 0 ? (x + y - 1) / y : x / y);}int POW(int x, int y) {if (y != (int)(y) or y < 0 or x != 0 && x != 1 && y > 64) {cout << "Error" << endl;return 0;}if (y == 0) return 1;if (y % 2 == 0) return POW(x * x, y / 2);return x * POW(x, y - 1);}__attribute__((constructor))void constructor() {ios::sync_with_stdio(false);cin.tie(nullptr);cout << fixed << setprecision(15);}#pragma endregionsigned main() {int N, D;cin >> N >> D;vector<int> mp1(2 * N * N + 1);for (int x = 1; x <= N; x++) {for (int y = 1; y <= N; y++) {mp1[x * x + y * y]++;}}vector<int> mp2(2 * N * N + 1 + D);for (int z = 1; z <= N; z++) {for (int w = 1; w <= N; w++) {if (w * w - z * z + D < 0) continue;mp2[w * w - z * z + D]++;}}int ans = 0;for (int i = 1; i <= 2 * N * N; i++) {ans += mp1[i] * mp2[i];}cout << ans << endl;}