結果
問題 | No.800 四平方定理 |
ユーザー |
|
提出日時 | 2019-03-17 22:21:02 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 45 ms / 2,000 ms |
コード長 | 2,427 bytes |
コンパイル時間 | 1,453 ms |
コンパイル使用メモリ | 168,256 KB |
実行使用メモリ | 42,352 KB |
最終ジャッジ日時 | 2024-07-08 00:04:28 |
合計ジャッジ時間 | 2,557 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
// template version 1.6// {{{ define basic macrousing 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#endiftypedef 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_FILEstd::ifstream in(infile);std::cin.rdbuf(in.rdbuf());#endifsolve();return 0;} //}}}