博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
【BZOJ】 1041: [HAOI2008]圆上的整点
阅读量:4548 次
发布时间:2019-06-08

本文共 1489 字,大约阅读时间需要 4 分钟。

题目链接:


${x^{2}+y^{2}=r^{2} }$

${\Rightarrow y^{2}=(r-x)(r+x)}$

令${d=gcd(r-x,r+x)}$

则${y^{2}=d^{2}*\frac{r+x}{d}*\frac{r-x}{d}}$

再令${A=\frac{r+x}{d}}$,${B=\frac{r-x}{d}}$

则${y^{2}=d^{2}*A*B}$

考虑${y^{2}}$是完全平方数,${d^{2}}$是完全平方数,又${gcd(A,B)=1}$那么${A,B}$都是完全平方数。

设${A=a^{2}}$,${B=b^{2}}$

${A+B=a^{2}+b^{2}}$

${\Rightarrow \frac{2*r}{d}=a^{2}+b^{2}}$

 

  考虑枚举${\frac{2*r}{d}}$,这一步的复杂度是${O(\sqrt{r})}$的,然后再在${\left [ 1,\sqrt{2*\frac{r}{d}}/2 \right ]}$的范围内枚举${a}$,进而可以算出${A,b,B}$,然后判断${A,B}$是否互质,$B$是否为完全平方数,这样子就算出了第一象限的答案,然后将$ans*4+4$,算是统计了每一个象限的并且加上了坐标轴上的四个点。

 

1 #include
2 #include
3 #include
4 #include
5 #include
6 #include
7 #include
8 #include
9 #include
10 using namespace std;11 #define llg long long12 #define maxn 10001013 #define yyj(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);14 llg ans,n;15 inline llg getint()16 {17 llg w=0,q=0; char c=getchar();18 while((c<'0' || c>'9') && c!='-') c=getchar();19 if (c=='-') q=1, c=getchar(); while (c>='0' && c<='9') w=w*10+c-'0', c=getchar();20 return q ? -w : w;21 }22 23 void calc(llg d)24 {25 for (llg a=1;a<=sqrt(d/2);a++)26 {27 llg A=a*a,B=d-A,b=sqrt(B);28 if (b*b==B && __gcd(A,B)==1 && A!=B) ans++;29 }30 }31 32 int main()33 {34 yyj("circle");35 cin>>n;36 for (llg i=1;i<=sqrt(n*2);i++)37 if ((2*n%i)==0)38 {39 calc(i);40 calc(2*n/i);41 }42 cout<

 

转载于:https://www.cnblogs.com/Dragon-Light/p/6744132.html

你可能感兴趣的文章
挂载光盘与rpm安装
查看>>
[Android学习系列18]线程,进程,异步的一些事
查看>>
腾讯 AI Lab 计算机视觉中心人脸 & OCR团队近期成果介绍(3)
查看>>
课堂练习-增加信息
查看>>
A little issue in Mathematical Thought from Ancient to Modern Times, Vol. 3
查看>>
Zabbix对接AD域
查看>>
django 将view视图中的对象传入forms表单验证模块中
查看>>
log4net配置
查看>>
中文词频统计及词云制作
查看>>
python面试题No5
查看>>
BeanShell PreProcessor数据base64加密
查看>>
10条建议帮助你创建更好的jQuery插件
查看>>
setPreferredSize和setSize的区别及用法
查看>>
Python简介及编码
查看>>
[转]Android:Layout_weight的深刻理解
查看>>
监听键盘弹出 隐藏
查看>>
iOS开发 - NSBundle, NSDevice, NSLocale
查看>>
innerHtml安全问题
查看>>
UVA 11992,。。。伪-二维线段树
查看>>
[原创]通过函数指针实现事件消息处理
查看>>