結果

問題 No.800 四平方定理
ユーザー HaarHaar
提出日時 2019-03-17 21:39:55
言語 C++14
(gcc 11.2.0 + boost 1.78.0)
結果
AC  
実行時間 72 ms / 2,000 ms
コード長 1,739 bytes
コンパイル時間 1,352 ms
使用メモリ 38,280 KB
最終ジャッジ日時 2023-02-11 05:12:55
合計ジャッジ時間 3,865 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge13
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
testcase_00 AC 21 ms
38,240 KB
testcase_01 AC 21 ms
38,184 KB
testcase_02 AC 21 ms
38,280 KB
testcase_03 AC 20 ms
38,156 KB
testcase_04 AC 20 ms
38,156 KB
testcase_05 AC 20 ms
38,160 KB
testcase_06 AC 22 ms
38,280 KB
testcase_07 AC 21 ms
38,208 KB
testcase_08 AC 21 ms
38,192 KB
testcase_09 AC 20 ms
38,088 KB
testcase_10 AC 31 ms
38,188 KB
testcase_11 AC 34 ms
38,260 KB
testcase_12 AC 34 ms
38,156 KB
testcase_13 AC 33 ms
38,160 KB
testcase_14 AC 32 ms
38,276 KB
testcase_15 AC 36 ms
38,232 KB
testcase_16 AC 32 ms
38,088 KB
testcase_17 AC 36 ms
38,252 KB
testcase_18 AC 35 ms
38,160 KB
testcase_19 AC 36 ms
38,280 KB
testcase_20 AC 19 ms
38,156 KB
testcase_21 AC 21 ms
38,156 KB
testcase_22 AC 36 ms
38,264 KB
testcase_23 AC 66 ms
38,080 KB
testcase_24 AC 67 ms
38,152 KB
testcase_25 AC 64 ms
38,088 KB
testcase_26 AC 19 ms
38,184 KB
testcase_27 AC 21 ms
38,204 KB
testcase_28 AC 61 ms
38,276 KB
testcase_29 AC 72 ms
38,148 KB
testcase_30 AC 61 ms
38,200 KB
testcase_31 AC 58 ms
38,224 KB
testcase_32 AC 62 ms
38,192 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define FOR(v, a, b) for(int v = (a); v < (b); ++v)
#define FORE(v, a, b) for(int v = (a); v <= (b); ++v)
#define REP(v, n) FOR(v, 0, n)
#define REPE(v, n) FORE(v, 0, n)
#define REV(v, a, b) for(int v = (a); v >= (b); --v)
#define ALL(x) (x).begin(), (x).end()
#define ITR(it, c) for(auto it = (c).begin(); it != (c).end(); ++it)
#define RITR(it, c) for(auto it = (c).rbegin(); it != (c).rend(); ++it)
#define EXIST(c,x) ((c).find(x) != (c).end())
#define LLI long long int
#define fst first
#define snd second

#ifdef DEBUG
#include <misc/C++/Debug.cpp>
#else
#define dump(x)
#endif

#define gcd __gcd

using namespace std;
template <class T> constexpr T lcm(T m, T n){return m/gcd(m,n)*n;}

template <typename I> void join(ostream &ost, I s, I t, string d=" "){for(auto i=s; i!=t; ++i){if(i!=s)ost<<d; ost<<*i;}ost<<endl;}
template <typename T> istream& operator>>(istream &is, vector<T> &v){for(auto &a : v) is >> a; return is;}
template <typename T, typename U> istream& operator>>(istream &is, pair<T,U> &p){is >> p.first >> p.second; return is;}

template <typename T, typename U> T& chmin(T &a, const U &b){return a = (a<=b?a:b);}
template <typename T, typename U> T& chmax(T &a, const U &b){return a = (a>=b?a:b);}
template <typename T, size_t N, typename U> void fill_array(T (&a)[N], const U &v){fill((U*)a, (U*)(a+N), v);}

int main(){
  cin.tie(0);
  ios::sync_with_stdio(false);

  LLI n,d;
  while(cin >> n >> d){
    LLI ans = 0;

    vector<int> count(1000000 + 2000*2000*2+100);
    FORE(x,1,n){
      FORE(y,1,n){
	++count[x*x+y*y];
      }
    }

    FORE(z,1,n){
      FORE(w,1,n){
	if(w*w-z*z+d > 0){
	  ans += count[w*w-z*z+d];
	}
      }
    }

    cout << ans << endl;
  }
  
  return 0;
}
0