BZOJ(2) 1041: [HAOI2008]圆上的整点

1041: [HAOI2008]圆上的整点

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 4966  Solved: 2258
[Submit][Status][Discuss]

Description

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

Input

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

Output

整点个数

Sample Input

4

Sample Output

4

HINT

科普视频

Source

思路：看一眼题目，感觉可做。然后开始想思路，想着想着发现不可做，所以就写了个暴力。

这是第一个：

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
long long n;
long long ans;
int main(){
    scanf("%lld",&n);
    int now=sqrt(n);
    long long maxn=n*n;
    for(int i=;i<=n;i++){
        long long k=i*i;
        long long z=maxn-k,s=sqrt(z);
        if(s*s!=z)    continue;
        if(k!=&&s!=)    ans+=;
        else ans+=;
    }
    cout<<ans;
}

后来看了一下题解里的视频，想了想，感觉很有道理，就开始了第二次奋斗历程。

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
long long n;
long long ans=,sum;
int p[],c[];
int main(){
    scanf("%lld",&n);
    int m=n*n;
    for(int i=;i<=n;i++)
        if(m%i==){
            p[++sum]=i;
            while(m%i==) m/=i,c[sum]++;
        }
    if(m>){ p[++sum]=m;c[sum]++; }
    for(int i=;i<=sum;i++)    cout<<p[i]<<" "<<c[i]<<endl;
    for(int i=;i<=sum;i++){
        if(p[i]%==)    continue;
        else if(p[i]%==)    ans*=c[i];
        else if(p[i]%==&&c[i]%!=)    ans=;
    }
    if(ans==){ printf("0\n");return ; }
    else if(ans==){ printf("4\n");return ; }
    else cout<<ans*+;
}

结果发现，调试忘记注释了。。。

只有某个奇怪的点得分了。。。

然后我就以为我把调试注释以后，就能AC4个点。然并卵，还是20分。

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
long long n;
long long ans=,sum;
int p[],c[];
int main(){
    scanf("%lld",&n);
    int m=n*n;
    for(long long i=;i<=n;i++)
        if(m%i==){
            p[++sum]=i;
            while(m%i==) m/=i,c[sum]++;
        }
    if(m>){ p[++sum]=m;c[sum]++; }
//    for(int i=1;i<=sum;i++)    cout<<p[i]<<" "<<c[i]<<endl;
    for(int i=;i<=sum;i++){
        if(p[i]%==)    continue;
        else if(p[i]%==)    ans*=c[i];
        else if(p[i]%==&&c[i]%!=)    ans=;
    }
    if(ans==){ printf("0\n");return ; }
    else if(ans==){ printf("4\n");return ; }
    else cout<<ans*+;
}

最后，终于胜利了！！

这是对上一个代码的优化。

首先 先把n的因数中的2都消去，因为2对答案并没有贡献。

其次 一个数 n*n=n^2 可以因式分解n^2=p1^a1*2*p2^a1*2*p3^a3*2*......

所以，因式分解后的质数，如果模4后余1那ans=ans*（ai*2+1）;最后ans*4。

还有自身是%4=1的质数的情况特判一下。ans=ans*（1*2+1）=ans*3;

然后输出就好了。

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
int ans=;
int n,sum=;
int main(){
    scanf("%d",&n);
    while((n&)^)    n>>=;
    while(sum*sum<=n)    sum++;
    for(int i=;i<=sum;i++)
        if(n%i==){
            int tmp=;
            while(n%i==){ n/=i;tmp+=; }
            if(i%==)    ans=ans*(tmp+);
        }
    if(n>&&n%==)    ans*=;
    printf("%d",ans);
}