結果
問題 | No.800 四平方定理 |
ユーザー |
|
提出日時 | 2019-03-22 21:35:45 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 133 ms / 2,000 ms |
コード長 | 874 bytes |
コンパイル時間 | 1,515 ms |
コンパイル使用メモリ | 167,572 KB |
実行使用メモリ | 64,384 KB |
最終ジャッジ日時 | 2024-09-19 02:55:17 |
合計ジャッジ時間 | 4,218 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define rep(i,j,n) for(int i=j;i<n;++i) #define all(i) i.begin(),i.end() #define rall(i) i.rbegin(), i.rend() #define INF 1e9 #define LINF 1e18 const int mod = 1e9 + 7; typedef long long i64; typedef pair<int, int> pi; template <class T> using vt = vector<T>; template <class T> using vvt = vector<vector<T>>; i64 gcd(i64 n, i64 m) {return (m == 0? n : gcd(m, n % m));} i64 lcd(i64 n, i64 m) {return (n / gcd(n, m) * m);} int dx[] = {1, 0, -1, 0}; int dy[] = {0, 1, 0, -1}; i64 cnt[10000000]; int main() { cin.tie(0); ios::sync_with_stdio(false); int n, d; cin >> n >> d; rep(x, 1, n + 1) { rep(y, 1, n + 1) { cnt[x * x + y * y]++; } } i64 ans = 0; rep(w, 1, n + 1) { rep(z, 1, n + 1) { if(w * w - z * z + d >= 0) ans += cnt[w * w - z * z + d]; } } cout << ans << endl; }