[bzoj1041][HAOI2008]圆上的整点

我能想得出怎么做才奇怪好吗

题解:http://blog.csdn.net/csyzcyj/article/details/10044629

#include<iostream>
#include<cstdio>
#include<cmath>
#define ll long long
using namespace std;
ll read()
{
    ll x=,f=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
    while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
    return x*f;
}

ll r;
ll ans=;

ll gcd(ll a,ll b) {return b==?a:gcd(b,a%b);}
inline ll sqr(ll x){return x*x;}

int main()
{
    r=read();ll to=sqrt(*r);
    for(ll i=;i<=to;i++)
        if(*r%i==)
        {
            for(ll j=;j<=sqrt((ll)*r/(*i));j++)
            {
                double b=sqrt((*r/i)-sqr(j));
                if(floor(b)==b&&sqr(j)!=sqr(floor(b))&&gcd(sqr(floor(b)),sqr(j))==)
                    ans++;
            }
            for(ll j=;j<=sqrt(i/);j++)
            {
                double b=sqrt(i-sqr(j));
                if(floor(b)==b&&sqr(j)!=sqr(floor(b))&&gcd(sqr(floor(b)),sqr(j))==)
                    ans++;
            }
        }
    printf("%lld",ans*+);
    return ;
}