結果
| 問題 |
No.800 四平方定理
|
| コンテスト | |
| ユーザー |
Bwambocos
|
| 提出日時 | 2019-03-17 22:54:40 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 821 bytes |
| コンパイル時間 | 1,591 ms |
| コンパイル使用メモリ | 168,040 KB |
| 実行使用メモリ | 13,756 KB |
| 最終ジャッジ日時 | 2024-07-08 01:26:01 |
| 合計ジャッジ時間 | 6,506 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | -- * 3 |
| other | AC * 1 WA * 9 TLE * 1 -- * 19 |
ソースコード
#include "bits/stdc++.h"
#define in std::cin
#define out std::cout
#define rep(i,N) for(LL i=0;i<N;++i)
typedef long long int LL;
int main()
{
LL N, D;
in >> N >> D;
LL ans = 0;
for (LL x = 1; x <= N; ++x)
{
for (LL y = 1; y <= N; ++y)
{
LL C = x * x + y * y;
LL M = D - C;
if (M == 0)
{
if (x == y) ans += (N - 1) * 3 + 1;
else ans += (N - 2) * 6 + 2 * 3;
}
else
{
for (LL d = 1; d*d <= std::abs(M); ++d)
{
if (M % d != 0) continue;
LL a = d, b = M / d;
if (a < b) std::swap(a, b);
if ((a + b) % 2 == 0)
{
LL z = (a + b) / 2;
LL w = a - z;
if (b > 0 && z < w) std::swap(z, w);
if (b < 0 && z > w) std::swap(z, w);
if (1 <= z && z <= N && 1 <= w && w <= N) ++ans;
}
}
}
}
}
out << ans << std::endl;
}
Bwambocos