[BZOJ1041] [HAOI2008]圆上的整点

Description

　　求一个给定的圆(x^2+y^2=r^2)，在圆周上有多少个点的坐标是整数。

Input

　　只有一个正整数r,r<=2000 000 000
　
Output

　　整点个数
　
Sample Input

4

Sample Output

4

题解：
很简单的数学题（看起来）。
1.我们可以发现当x或者y等于0的时候，就有4个整点了（easy）。
2.我们把圆方程变式一下： y^2 = r^2 - x^2 = (r+x)(r-x) 。 那么y = sqrt( (r+x)(r-x) )
3.如果y要是整数的话，那么（r+x）*（r-x）要为完全平方数。
4.我们设i=gcd(r+x,r-x) ,那么 gcd((r+x) / i,(r-x) / i)==1
5.设 整数 s =( r - x ) / i , x = ( r + x ) / i 则 s + x = ( 2 * r ) / i ， 那么因为 s + x 为整数 ， 那么 i 要是 2 * r 的因数
6.那么我们可知y * y = s * z * i * i ，由于 i * i 和 y * y是完全平方数，s 与 z 是互质的，又因为 s * z 是完全平方数，那么s和z都是完全平方数。
7.由于s+z=（2 * r）/ i , 我们知道s+z是整数，那么i是2 * r 的因数。那么我们只要枚举 i 和 sqrt(s)就可以了。记得判断x，y大小与判重。（对于每一个2*r 的因数要计算两次，因为可能是 i 或者 (2 * r )/ i）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<cstdlib>
#include<string>
#include<cmath>
#include<algorithm>
#define LiangJiaJun main
#define ll long long
using namespace std;
ll gcd(ll a,ll b){return (b==0)?a:gcd(b,a%b);}
ll R,ans = 0;

bool pd(ll grep,double lp){
     if(lp != floor(lp))return 0;
     ll qs = floor(lp);
     if(gcd( grep*grep , qs*qs ) == 1 && grep * grep != qs*qs)return 1;
     return 0;
}
int LiangJiaJun (){
    cin>>R;
    ans = 4 ;
    for(ll i = 1; i <= (ll)sqrt(2ll*R) ; i++){
        if(2ll * R % i != 0) continue;
        for(ll s = 1 ; s <= (ll)sqrt(R/i) ; s++ ){
            double z = sqrt(((2 * R )/ i ) - s * s) ;
            if(pd(s,z)) ans += 4ll;
        }
        if(i!=2ll*R/i){
            for(ll s = 1; s <= (ll)sqrt(i/2ll); s++){
                double z = sqrt(i - s*s);
                if(pd(s,z)) ans += 4ll;
            }
        }
    }
    cout<<ans<<endl;
    return 0;
}