結果
問題 | No.800 四平方定理 |
ユーザー |
![]() |
提出日時 | 2020-05-07 08:42:41 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 305 ms / 2,000 ms |
コード長 | 912 bytes |
コンパイル時間 | 2,739 ms |
コンパイル使用メモリ | 193,824 KB |
最終ジャッジ日時 | 2025-01-10 07:27:07 |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
#include<bits/stdc++.h> using namespace std; #define ALL(x) x.begin(),x.end() #define rep(i,n) for(int i=0;i<(n);i++) #define debug(v) cout<<#v<<":";for(auto x:v){cout<<x<<' ';}cout<<endl; #define INF 1000000000 #define mod 1000000007 using ll=long long; const ll LINF=1001002003004005006ll; int dx[]={1,0,-1,0}; int dy[]={0,1,0,-1}; // ll gcd(ll a,ll b){return b?gcd(b,a%b):a;} template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return true;}return false;} template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return true;}return false;} signed main(){ cin.tie(0); ios::sync_with_stdio(0); ll n,d;cin>>n>>d; vector<ll> x(10000000,0),y(10000000,0); for(ll i=1;i<=n;i++){ for(ll j=1;j<=n;j++){ x[i*i+j*j]++; if(d+i*i-j*j>=0) y[d+i*i-j*j]++; } } ll ans=0; rep(i,10000000) ans+=x[i]*y[i]; cout<<ans<<endl; return 0; }