bzoj1041 [HAOI2008]圆上的整点 gcd
这个题非常恶心,因为条件很少,要求也很少,看起来没有任何特殊的性质
所以只能往约数、gcd上靠
然后就是
x^2=r^2-y^2=(r-y)(r+y)
设A=r-y, B=r+y
由于A*B是完全平方数, 设d=gcd(A,B)
则 gcd(A/d,B/d)=1
设A=u*d B=v*d
所以u=(r-y)/d , v=(r+y)/d
则v+u=2r/d
则u和v都是完全平方数
再设 u=a*a , v=b*b
枚举2r的因子d , 再枚举a, 统计答案
码:
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
long long n,i,j,a,b,ans;
long long gcd(long long a,long long b)
{if(!b)return a;return gcd(b,a%b);
}
int main()
{scanf("%lld",&n);for(i=1;i*i<=2*n;i++){if((2*n)%i==0){for(j=1;j*j<=2*n/i;j++)//d==i {a=j*j;b=2*n/i-a;if(b<=0)break;long long p=sqrt(b);if(gcd(a,b)==1&&(p*p)<=(b)&&(p*p)>=(b)){ans+=1; } }if(i*i!=2*n){for(j=1;j*j<=i;j++)//d==2*n/i {a=j*j;b=2*n/(2*n/i)-a;if(b<=0)break;long long p=sqrt(b);if(gcd(a,b)==1&&(p*p)<=(b)&&(p*p)>=(b)){ans++; } }} } }printf("%lld",(ans*4+4)/2);
}
bzoj1041 [HAOI2008]圆上的整点 gcd相关推荐
- bzoj千题计划127:bzoj1041: [HAOI2008]圆上的整点
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 ...
- BZOJ1041:[HAOI2008]圆上的整点(数论)
Description 求一个给定的圆(x^2+y^2=r^2),在圆周上有多少个点的坐标是整数. Input 只有一个正整数n,n<=2000 000 000 Output 整点个数 Samp ...
- bzoj 1041: [HAOI2008]圆上的整点
1041: [HAOI2008]圆上的整点 Time Limit: 10 Sec Memory Limit: 162 MB Submit: 2838 Solved: 1238 [Submit][S ...
- 【HAOI2008】【BZOJ1041】圆上的整点
Description 求一个给定的圆(x^2+y^2=r^2),在圆周上有多少个点的坐标是整数. Input r Output 整点个数 Sample Input 4 Sample Output 4 ...
- 【codevs1867】【Tyvj3508】【BZOJ1041】圆上的整点,数学乱搞
传送门1 传送门2 传送门3 写在前面:乱搞数学-- 思路: 原博文 有了上面的推理,那么实现的方法为: 枚举d∈[1,sqrt(2R)],然后根据上述推理可知:必先判d是否为2R的一约数. 此时d为 ...
- [BZOJ 1041] 圆上的整点
Description 求一个给定的圆( x2+y2=r2 x 2 + y 2 = r 2 x^2+y^2=r^2),在圆周上有多少个点的坐标是整数 Input 只有一个正整数 n n n,n≤200 ...
- BZOJ 1041 圆上的整点 数学
题目链接: https://www.lydsy.com/JudgeOnline/problem.php?id=1041 题目大意: 求一个给定的圆(x^2+y^2=r^2),在圆周上有多少个点的坐标是 ...
- 2019 ICPC 上海站网络赛 K.Peekaboo (圆上整点)
https://nanti.jisuanke.com/t/41421 题意:给定平面上两个点到原点之间的距离aaa和bbb以及这两点之间的距离ccc,且这两点的坐标均是整数,求该两点可能的坐标,按字典 ...
- BZOJ 4544: 椭圆上的整点
Sol 数学. 跟圆上的整点一样...TA写了个积性函数的算法...以后再说吧... \(x^2+3y^2=r^2\) \(3y^2=r^2-x^2\) \(3y^2=(r-x)(r+x)\) \(y ...
最新文章
- linux怎么装谷歌浏览器,如何在linux CentOS 上安装chrome 谷歌浏览器?
- UA MATH567 高维统计II 随机向量2 各向同性的随机向量
- python守护多线程_Python多线程Threading、子线程与守护线程实例详解
- 怎么用计算机知道别人手机密码,怎样才能知道别人家的wifi密码
- [活动 3.30]MAUI 跨平台应用开发实战
- C# 8 新特性 - 异步流 Asynchronous Streams
- 前端学习(702):while循环
- 如何动态创建asp.net 用户控件
- OpenInfra Days China 2020 官网正式上线,含免费注册通道!
- scala中的if esle使用
- 拷贝网页内容增加版权信息的 JavaScript 代码示例
- 使用STM8S003F3P6的硬件I2C读写AT24C16的EEPROM的经验心得
- html5%3chr%3e的样式,Vbs脚本编程简明教程
- 浅谈ARP攻击现象和处理办法
- Informatic学习总结_day02_增量抽取
- 如何打造一个能自动回复的钉钉机器人
- 医学图像中的窗宽(Window Width,WW)和窗位(Window Level,WL)
- 计算机时间小于会计期间错误,会计期间手工记账和电脑记账不一样问题
- 阿里云的ACA认证到底是个啥?有用吗?
- 人类已经无法阻止苹果了——吐槽PC厂商
热门文章
- 红橙Darren视频笔记 view的绘制流程(上) onMeasure测量代码分析 基于API27
- 【李宏毅机器学习】Unsupervised Learning - Word Embedding 无监督学习 - 词嵌入(p22) 学习笔记
- 【idea基础知识】常用快捷键整理
- navicate将远程数据库导入到本地数据库
- latex参考文献BibTeX的使用方法
- day16 java中的package和import
- centos7重启网卡命令_重启 CentOS 7 系统后的 IP 地址问题
- 上项线体表位置_体表定位
- ip动态分配痕迹会保留多久_段王爷《新国潮七剑》为你解剖新国潮还能潮多久?...
- pyqt tableview大数据自动退出_JAVAFX之tableview界面实时刷新导致的内存溢出(自己挖的坑,爬着也要出来啊0.0)...