bzoj1041题解
来源:互联网 发布:个人软件定制 编辑:程序博客网 时间:2024/05/22 01:46
求一个给定的圆(x^2+y^2=r^2),在圆周上有多少个点的坐标是整数。r<=2000 000 000
这道题刚看时,就明白暴力不能解决一切。要是r^2<=20亿,还可以sqrt循环,这样的话只能用数学方法了。
(嘘!本数学方法的关键一步是借鉴来的)x^2+y^2=r^2
x^2=r^2-y^2
x^2=(r+y)(r-y)
我们设gcd(r+y,r-y)=d //不知道什么是GCD?就是最大公因数
再设m=(r+y)/d,n=(r-y)/d,易得gcd(m,n)=1
因为m*n=x^2/d/d,即m*n是平方数。
又因为m和n互质,所以它们分别是平方数,我们还可以设m=u^2,n=v^2 。
带回去,r+y=d*u^2,r-y=d*v^2。
两式相加,2r=d(u^2+v^2)
加到这里,这些字母都表示整数。
现在我们的目标达成了:即把r^2降成了r。
我们可以枚举sqrt(2r),再依次判断是否有整数的u,v,且它们是否互质 。
最后还要注意以下:首先,所得结果要乘4 ,比如(3,4),还有(3,-4),(-3,4),(-3,-4)。
当然,如果x=0或y=0是一定可以的。最后还要加上4,即(0,r),(r,0),(0,-r),(-r,0)。(呵呵,r不等于0,O(∩_∩)O~~)
以下附代码:
#include<cstdio>#include<cmath>using namespace std;const double jing=0.00000001;long gcd(long a,long b){if (a%b==0) return b;return gcd(b,a%b);}long long ans=0,r,temp;long i,j;double t;int main(){ scanf("%lld",&r);r*=2; for (i=1;i<=trunc(sqrt(r));i++) if (r%i==0) { temp=r/i; for (j=1;j<=trunc(sqrt(temp));j++) { t=trunc(sqrt(temp-j*j));if (t<=j) break; if (t*t+j*j==temp){if (gcd(j,t)==1) ans++;} } if (r/i==i) continue;temp=i; for (j=1;j<=trunc(sqrt(temp));j++) { t=trunc(sqrt(temp-j*j));if (t<=j) break; if (t*t+j*j==temp){if (gcd(j,t)==1) ans++;} } } ans*=4;ans+=4;printf("%ld",ans); return 0;}
0 0
- bzoj1041题解
- bzoj1041
- bzoj1041
- 【数学分析】BZOJ1041(HAOI2008)[圆上的整点]题解
- bzoj1041未完成
- bzoj1041 圆上的整点
- BZOJ1041 圆上的整点 Solution
- 【BZOJ1041】[HAOI2008]圆上的整点 数学
- [BZOJ1041]HAOI2008圆上的整点
- 【HAOI2008】【BZOJ1041】圆上的整点
- 【BZOJ1041】[HAOI2008]圆上的整点【数学】
- [BZOJ1041][HAOI2008]圆上的整点
- bzoj1041: [HAOI2008]圆上的整点
- [枚举] BZOJ1041: [HAOI2008]圆上的整点
- bzoj1041: [HAOI2008]圆上的整点
- bzoj1041[HAOI2008]圆上的整点
- [BZOJ1041] [HAOI2008]圆上的整点
- 【bzoj1041】[HAOI2008]圆上的整点
- 《java编程思想系列之7》 内部类
- python for android : 显示图片
- Linux中ifcfg-eth0配置参数说明
- 射线与平面的相交检测(Ray-Plane intersection test)
- 第一周工作心情
- bzoj1041题解
- Android内核详解之Low memory killer
- Java关于String十个问题
- (转)在endnote中制作GB/T7714《文后参考文献著录规则》的输出格式
- *bzoj1083题解
- 国内也有不错的开源项目,只是感觉后劲不足,亟须自然科学基金赞助
- mysql错误代码
- Reorder List
- poj 2739 Sum of Consecutive Prime Numbers 小结