bzoj1041: [HAOI2008]圆上的整点

题面在这里
emmmm先来吐槽一下下~
%一发dsy上0ms的大佬们。。。
我就比较弱辣要200+ms，，QwQ
做法：
来做一些鬼畜的变形= =
$r^2=x^2+y^2$
$y^2=r^2-x^2$
$y^2=(r+x)(r-x)$
$y=\sqrt{(r+x)(r-x)}$
设$d=gcd(r+x,r-x)$
$a=\frac{r+x}{d}$，$b=\frac{r-x}{d}$
则$y=\sqrt{d^2\times{a}\times{b}}=d\sqrt{a\times{b}}$
因为$a$和$b$互质，所以$a$和$b$是完全平方数
设$a=A^2$,$b=B^2$
则$A^2=a=\frac{r+x}{d}$,$B^2=b=\frac{r-x}{d}$
两式相加得$\frac{2r}{d}=A^2+B^2$
所以我们可以考虑枚举$d$，再枚举$A$，计算出$B$，然后判断是否合法，即满足：1) $B$是整数，2) $gcd(A,B)=1$。
至此，我们求出某一象限内的个数，因为要四个象限，于是答案$\times{4}$。
再加上坐标轴的$4$个点，于是答案再$+4$。

/*************************************************************
    Problem: bzoj 1041 [HAOI2008]圆上的整点
    User: fengyuan
    Language: C++
    Result: Accepted
    Time: 276 ms
    Memory: 1300 kb
    Submit_Time: 2017-12-14 19:52:42
*************************************************************/

#include<bits/stdc++.h>
#define rep(i, x, y) for (int i = (x); i <= (y); i ++)
#define down(i, x, y) for (int i = (x); i >= (y); i --)
#define mid ((l+r)/2)
#define lc (o<<1)
#define rc (o<<1|1)
#define pb push_back
#define mp make_pair
#define PII pair<int, int>
#define F first
#define S second
#define B begin()
#define E end()
using namespace std;
typedef long long LL;
//head

LL r;

inline LL gcd(LL a, LL b){ return !b ? a : gcd(b, a%b); }

int main()
{
    scanf("%lld", &r); LL ans = 0;
    rep(i, 1, floor(sqrt(2*r))) if (2*r % i == 0){
        LL d = i;
        rep(j, 1, floor(sqrt(2*r/d))){
            LL a = j;
            LL b = sqrt(2*r/d-a*a);
            if (a*a + b*b != 2*r/d || gcd(a, b) != 1 || b <= a) continue;
            ans ++;
        }
        if (d*d == 2*r || d == 1) continue;
        d = 2*r/d;
        rep(j, 1, floor(sqrt(2*r/d))){
            LL a = j;
            LL b = sqrt(2*r/d-a*a);
            if (a*a + b*b != 2*r/d || gcd(a, b) != 1 || b <= a) continue;
            ans ++;
        }
    }
    printf("%lld\n", ans*4+4);
    return 0;
}