結果
| 問題 |
No.800 四平方定理
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2021-08-19 23:25:49 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 106 ms / 2,000 ms |
| コード長 | 800 bytes |
| コンパイル時間 | 1,384 ms |
| コンパイル使用メモリ | 169,648 KB |
| 実行使用メモリ | 81,396 KB |
| 最終ジャッジ日時 | 2024-10-13 01:19:57 |
| 合計ジャッジ時間 | 4,075 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 30 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
#define all(A) A.begin(),A.end()
using vll = vector<ll>;
#define rep(i, n) for (long long i = 0; i < (long long)(n); i++)
using Graph = vector<vector<ll>>;
ll gcd(ll(a), ll(b)) {
ll c = a;
while (a % b != 0) {
c = a % b;
a = b;
b = c;
}
return b;
}
int main() {
ll N,D;
cin>>N>>D;
//map<ll,ll> XY,ZW;
vll XY(2000*2000+D+5,0);
vll ZW(2000*2000+D+5,0);
for(ll X=1;X<=N;X++){
for(ll Y=1;Y<=N;Y++){
if(X*X+Y*Y<=N*N+D){
XY[X*X+Y*Y]++;
}
if(X*X+D-Y*Y>=2){
ZW[X*X+D-Y*Y]++;
}
}
}
ll an=0;
rep(i,2000*2000+D+2){
an+=XY[i]*ZW[i];
}
cout<<an<<endl;
}