BZOJ 1041 [HAOI2008] 圆上的整点题解与分析

最新推荐文章于 2020-02-07 15:40:17 发布

csyzcyj

最新推荐文章于 2020-02-07 15:40:17 发布

阅读量6k

点赞数 4

分类专栏：题库 BZOJ(大视野测评)

本文链接：https://blog.csdn.net/csyzcyj/article/details/10044629

版权

题库同时被 2 个专栏收录

54 篇文章

订阅专栏

BZOJ(大视野测评)

19 篇文章

订阅专栏

该博客详细介绍了如何解决BZOJ 1041题目，即求解半径为n的圆周上整点的数量。通过数学分析，博主提出了一种优化算法，避免了暴力枚举，降低了时间复杂度。算法分为两种情况，分别针对d等于2R的约数的不同情况，并考虑了圆的对称性。最终的时间复杂度分析包括了枚举d和a的过程，以及求最大公约数的操作。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

[HAOI2008]圆上的整点

Time Limit: 10 Sec Memory Limit: 162 MB

Description

平面上有一个圆, 圆心坐标为(0,0),半径为n. 问圆周上有多少个整点. 整点的定义即x,y坐标均为整数的点.

input

输入文件只有一个正整数n.

Output

输出文件为一个正整数, 即圆周上有多少个整点.

Sample Input

Sample Output

HINT

20%的测试数据, n<=10000.
100%的测试数据, n<=2000 000 000, 保证答案在long long/int64 范围内.

样例说明：

4个整点分别为：(-4,0),(4,0),(0,4),(0,-4)

【分析】：

样例图示：

首先,最暴力的算法显而易见：枚举x轴上的每个点，带入圆的方程，检查是否算出的值是否为整点，这样的枚举量为2*N，显然过不了全点。

然后想数学方法。

有了上面的推理，那么实现的方法为：

枚举d∈[1,sqrt(2R)]，然后根据上述推理可知：必先判d是否为2R的一约数。

此时d为2R的约数有两种情况：d=d或d=2R/d。

第一种情况：d=2R/d。枚举a∈[1,sqrt(2R/2d)] <由2*a*a < 2*R/d转变来>，算出对应的b=sqrt(2R/d-a^2)，检查是否此时的A,B满足：A≠B且A,B互质 <根据上面的推理可知必需满足此条件>，若是就将答案加1

第二种情况：d=d。枚举a∈[1,sqrt(d/2)] <由2*a*a < d转变来>，算出对应的b=sqrt(d-a^2)，检查是否此时的A,B满足：A≠B且A,B互质 <根据上面的推理可知必需满足此条件>，若是就将答案加1

因为这样只算出了第一象限的情况<上面枚举时均是从1开始枚举>，根据圆的对称性，其他象限的整点数与第一象限中的整点数相同，最后，在象限轴上的4个整点未算，加上即可，那么最后答案为ans=4*第一象限整点数+4

【时间复杂度分析】：

枚举d：O(sqrt(2R)),然后两次枚举a：O(sqrt(d/2))+O(sqrt(R/d))，求最大公约数：O(logN)

【代码】：

/************************************************************** 
    Problem: 1041 
    Language: C++ 
    Result: Accepted 
    Time:188 ms 
    Memory:1284 kb 
****************************************************************/
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
#include<iostream>
#include<vector>
#include<stack>
#include<queue>
#include<cmath>
using namespace std;
long long R,ans=0;
long long gcd(long long x,long long y){return x%y==0 ? y : gcd(y,x%y);}
bool check(long long y,double x)
{
      if(x==floor(x))//判断整点 
      {
            long long x1=(long long)floor(x);
            if(gcd(x1*x1,y*y)==1 && x1*x1!=y*y)//gcd(A,B)=1并且A!=B 
                  return true;
      }
      return false;
}
int main()
{
      //freopen("input.in","r",stdin);
	  //freopen("output.out","w",stdout);
      scanf("%lld",&R);
      for(long long d=1;d<=(long long)sqrt(2*R);d++)
      {
            if((2*R)%d==0)
            {
                  for(long long a=1;a<=(long long)sqrt(2*R/(2*d));a++)//2*a^2<2*r/d
                  {
                        double b=sqrt(((2*R)/d)-a*a);
                        if(check(a,b))
                              ans++;
                  }
                  if(d!=(2*R)/d)
                  {
                        for(long long a=1;a<=(long long)sqrt(d/2);a++)//2*a^2<d
                        {
                              double b=sqrt(d-a*a);
                              if(check(a,b))
                                    ans++;
                        }
                  }
            }
      }    
      printf("%lld\n",ans*4+4); 
	  //system("pause");
      return 0;
}

转载注明出处：http://blog.csdn.net/u011400953