求一个给定的圆(x^2+y^2=r^2),在圆周上有多少个点的坐标是整数。
只有一个正整数n,n<=2000 000 000
整点个数
#include<bits/stdc++.h>
using namespace std;
#define maxn 2000000000
typedef long long ll;
int ans;
int gcd(int a,int b){
if(b==0) return a;
else return gcd(b,a%b);
}
void solve(ll rr)
{
for(ll i=1;i*i<=rr;i++){
ll t=i*i;
ll j= sqrt(rr-t);
if(i>=j) break;
if((j*j==rr-t)&&(gcd(i,j)==1))
ans++;
}
}
int main()
{
ll r;
while(~scanf("%lld",&r)){
ans=0;
r<<=1;
for(ll i=1;i*i<=r;i++){
if(r%i!=0) continue;
ll d=i;
solve(r/d);
if(d*d==r) break;
solve(d);
}
cout<<ans*4+4<<endl;
}
}
4
4
思路:比赛时候没有想出来,感觉自己还是比较笨;看了题解才做出来;下面说一下详细思路;
r^2=x^2+y^2;所以y=sqrt(r^2-x^2)=sqrt((r+x)*(r-x));如果暴力枚举x,y,10秒的时间也会爆,所以我们要尽可能的缩小数据范围;首先求出d=gcd(r+x,r-x);我们可以写出y^2/d^2=((r-x)/d)*((r+x)/d);
又因为(r-x)/d和(r+x)/互质并且等式左边是平方的形式,所以(r-x)/d和(r+x)/d一定也是平方的形式;设(r-x)/d=u^2,(r+x)/d=v^2;容易得出u,v也是互质的。那么u^2+v^2=((r-x)/d)+((r+x)/d)=2r/d;我们现在可以枚举2r的约数,对于每个约数x,我们有u^2<=x并且v^2<=x;所以我们可以在x范围枚举u,再求v,然后判断是否互质,如果是就算找到了一个答案。