結果
| 問題 | 
                            No.800 四平方定理
                             | 
                    
| コンテスト | |
| ユーザー | 
                             Shibuyap
                         | 
                    
| 提出日時 | 2019-10-14 15:03:30 | 
| 言語 | C++14  (gcc 13.3.0 + boost 1.87.0)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 99 ms / 2,000 ms | 
| コード長 | 738 bytes | 
| コンパイル時間 | 1,617 ms | 
| コンパイル使用メモリ | 167,488 KB | 
| 実行使用メモリ | 53,504 KB | 
| 最終ジャッジ日時 | 2024-12-24 00:07:06 | 
| 合計ジャッジ時間 | 3,843 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge1 / judge5 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 3 | 
| other | AC * 30 | 
ソースコード
#include <bits/stdc++.h>
#define rep(i,n) for(int i = 0; i < (n); ++i)
#define srep(i,s,t) for (int i = s; i < t; ++i)
#define drep(i,n) for(int i = (n)-1; i >= 0; --i)
using namespace std;
typedef long long int ll;
typedef pair<int,int> P;
#define yn {puts("Yes");}else{puts("No");}
#define MAX_XY 8000005
int xy[MAX_XY];
int zw[MAX_XY];
int main() {
    int n, d;
    cin >> n >> d;
    
    srep(i,1,n+1){
        srep(j,1,n+1){
            xy[i*i + j*j]++;
        }
    }
    srep(i,1,n+1){
        srep(j,1,n+1){
            int tmp = d - (i + j) * (i - j);
            if(0 <= tmp)zw[tmp]++;
        }
    }
    int ans = 0;
    rep(i,MAX_XY){
        ans += xy[i] * zw[i];
    }
    cout << ans << endl;
    return 0;
}
 
 
 
            
            
            
        
            
Shibuyap