結果

問題 No.800 四平方定理
ユーザー kjnh10kjnh10
提出日時 2019-03-17 22:21:02
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 61 ms / 2,000 ms
コード長 2,427 bytes
コンパイル時間 1,753 ms
コンパイル使用メモリ 166,060 KB
実行使用メモリ 42,168 KB
最終ジャッジ日時 2023-09-22 07:32:51
合計ジャッジ時間 3,873 ms
ジャッジサーバーID
(参考情報)
judge13 / judge15
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
4,380 KB
testcase_01 AC 2 ms
4,380 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 2 ms
4,376 KB
testcase_04 AC 2 ms
4,376 KB
testcase_05 AC 2 ms
4,376 KB
testcase_06 AC 2 ms
4,376 KB
testcase_07 AC 3 ms
4,376 KB
testcase_08 AC 4 ms
6,004 KB
testcase_09 AC 3 ms
4,380 KB
testcase_10 AC 27 ms
22,344 KB
testcase_11 AC 27 ms
22,364 KB
testcase_12 AC 28 ms
24,468 KB
testcase_13 AC 19 ms
19,248 KB
testcase_14 AC 22 ms
20,252 KB
testcase_15 AC 28 ms
24,364 KB
testcase_16 AC 22 ms
20,528 KB
testcase_17 AC 22 ms
20,512 KB
testcase_18 AC 31 ms
28,168 KB
testcase_19 AC 30 ms
28,484 KB
testcase_20 AC 1 ms
4,380 KB
testcase_21 AC 2 ms
4,376 KB
testcase_22 AC 31 ms
28,408 KB
testcase_23 AC 58 ms
42,168 KB
testcase_24 AC 50 ms
34,304 KB
testcase_25 AC 61 ms
42,156 KB
testcase_26 AC 1 ms
4,376 KB
testcase_27 AC 7 ms
10,828 KB
testcase_28 AC 41 ms
34,032 KB
testcase_29 AC 52 ms
38,176 KB
testcase_30 AC 48 ms
34,100 KB
testcase_31 AC 39 ms
32,120 KB
testcase_32 AC 40 ms
35,092 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

// template version 1.6
// {{{ define basic macro
using namespace std;
#include <iostream>
#include <bits/stdc++.h>
// varibable settings
// #define int long long
#define INF 1000000000000000000LL
// #define INF 2147483647
#define MOD 1000000007LL
// #define MOD 998244353LL
#define infile "../test/sample-1.in"

#define _overload3(_1,_2,_3,name,...) name
#define _rep(i,n) repi(i,0,n)
#define repi(i,a,b) for(int i=(int)(a);i<(int)(b);++i)
#define rep(...) _overload3(__VA_ARGS__,repi,_rep,)(__VA_ARGS__)
#define _rrep(i,n) rrepi(i,0,n)
#define rrepi(i,a,b) for(int i=(int)(b-1);i>=(int)(a);--i)
#define rrep(...) _overload3(__VA_ARGS__,rrepi,_rrep,)(__VA_ARGS__)
#define all(x) (x).begin(),(x).end()
#define sz(x) ((int)(x).size())
#define pb(a) push_back(a)
#define mp(a, b) make_pair(a, b)
#define uni(x) sort(all(x));x.erase(unique(all(x)),x.end())
#define ten(n) ((int)1e##n)
template<class T, size_t N> size_t sza(const T (&array)[N]) { return N; }
template<class T> inline void chmax(T &a, const T &b) { if(a < b) a = b; }
template<class T> inline void chmin(T &a, const T &b) { if(a > b) a = b; }
template<class T=int> T in() {T x; cin>>x; return (x);}
struct Fast { Fast(){ std::cin.tie(0); ios::sync_with_stdio(false); } } fast;

// dump macro
#ifdef PCM
  #include "dump.hpp"
#else
  #define dump(...) 42
  #define dump_1d(...) 42
  #define dump_2d(...) 42
#endif

typedef long long ll;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef long double ld;
typedef pair<int,int> pii;
typedef tuple<int,int,int> iii;
template<typename T> using PQ = priority_queue<T, vector<T>, greater<T>>;
int dx[]={1, -1, 0, 0};
int dy[]={0, 0, 1, -1};
#define fi first
#define se second

// }}}

int solve(){
  int n,d;cin>>n>>d;

  // map<int, int> d1;
  // map<int, int> d2;
  vi d1(n*n+d+1, 0);
  vi d2(n*n+d+1, 0);
  // dump("hello");
  int m = n*n - 1 + d;
  rep(x, 1, n+1){
    rep(y, 1, n+1){
      int t = x*x + y*y;
      if (t<=m) d1[t] += 1;
      if (t>=m) break;
    }
  }
  // dump("hello1");
  rep(w, 1, n+1){
    rep(z, 1, n+1){
      int t = w*w - z*z + d;
      if (t>=0) d2[t] += 1;
      else break;
    }
  }
  // dump("hello2");
  // dump(d1);
  // dump(d2);

  int res = 0;
  rep(i, 2, n*n+d+1){
    res += d1[i]*d2[i];
  }
  cout << res << endl;

  return 0;
}

signed main() { //{{{
#ifdef INPUT_FROM_FILE
  std::ifstream in(infile);
  std::cin.rdbuf(in.rdbuf());
#endif
  solve();
  return 0;
} //}}}
0