結果

問題 No.800 四平方定理
ユーザー kjnh10kjnh10
提出日時 2019-03-17 22:15:31
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 2,415 bytes
コンパイル時間 1,722 ms
コンパイル使用メモリ 166,868 KB
実行使用メモリ 42,368 KB
最終ジャッジ日時 2024-07-07 23:48:21
合計ジャッジ時間 4,612 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 16 RE * 14
権限があれば一括ダウンロードができます

ソースコード

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;
      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;
      if (t<=0) 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