bzoj千题计划127：bzoj1041: [HAOI2008]圆上的整点

标签：center   span   void   计划   data   continue   bsp   php   a*   

http://www.lydsy.com/JudgeOnline/problem.php?id=1041

设 X>0 ,Y>0

X^2 + Y^2 = R^2

X^2 = R^2-Y^2 = (R+Y)(R-Y)

令  d=gcd（R+Y,R-Y），A=(R+Y)/d，B=(R-Y)/d

则 gcd（A，B）=1，且A != B                                                                                                                                                                              

 X^2= d^2 *A * B

所以 A * B 为 完全平方数

又因为 gcd（A，B）=1 ，A!=B，所以 A，B 都是 完全平方数

令 a= 根号A，b=根号B

a^2 + b ^2 = 2*R / d

所以 d 必须是 2*R 的 约数   

根号(2*R) 枚举 约数 d

1、a^2 + b^2 = 2*R / d

2、a^2 + b^2 = d

对于 每一种 情况 分别 根号复杂度 枚举 a，计算b

判断相应的 A ，B 是否满足  gcd=1 且 A!=B

满足则 ans+1

这只算出了第一象限的情况

根据园的对称性，ans*4 可得 所有 象限内的点

最后在加上4个在 坐标轴上的点即可

#include<cmath>
#include<cstdio>

using namespace std;

typedef long long LL;

LL R;

int ans=0;

int gcd(int A,int B) { return !B ? A : gcd(B,A%B); }

void solve(int t,int d)
{
    int n=sqrt(t*1.0);
    int A,B,b;
    for(LL a=1;a<=n;++a)
    {
        B=t-a*a; b=sqrt(B);
        if(b*b!=B || !B) continue;
        A=a*a; 
        if(gcd(A,B)==1 && A!=B) ans++;
    }
}

int main()
{
    scanf("%lld",&R);
    int n=sqrt(R*2.0);
    for(int d=1;d<=n;++d)
    {
        if(R*2%d==0) 
        {
            solve(2*R/d,d);
            if(d*d!=n) solve(d,2*R/d);
        }
    }
    ans/=2;
    ans=ans*4+4;
    printf("%d",ans);
}

Description

Input

Output

Sample Input

Sample Output

bzoj千题计划127：bzoj1041: [HAOI2008]圆上的整点

标签：center   span   void   计划   data   continue   bsp   php   a*   

原文地址：http://www.cnblogs.com/TheRoadToTheGold/p/7931736.html