BZOJ 1041 [HAOI2008]圆上的整点:数学
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1041
题意:
给定n(n <= 2*10^9),问你在圆x^2 + y^2 = n^2的圆周上,有多少个坐标为整数的点。
题解:
科普视频:http://www.bilibili.com/video/av12131743/
推导的大致思路:
推导:
一、17 = 4^2 + 1^2
求圆周上有多少个点,就是求有多少个整数对(a,b)满足a^2 + b^2 = R^2。
二、17 = (4+i)*(4-i)
变形:a^2 + b^2 = (a + b*i) * (a - b*i) = R^2。
其中,a + b*i 与 a - b*i 复共轭。
也就是将R^2分解成(a + b*i) * (a - b*i)。
有一个结论,对于整数a来说:
(1)如果a为4n + 1型的素数,则a可以被分解为两个不同的高斯素数。
(2)如果a为4n + 3型的素数,则不能被分解。因为它们不仅是普通素数,还是高斯素数。
(即费马平方和定理:只有4n+1型的素数,才能表示成两个数的平方和)
分解方法:
(1)首先将R^2分解质因数,R^2 = a1^p1 + a2^p2 +...
(2)然后将R^2继续分解成若干高斯素数之积。
(3)将这些高斯素数分成两组,如果这两组各自之积复共轭,则为一对合法的(a,b)。
其中,将高斯素数分组时,对于一个素因子ai,有pi+1中分组方法。
特别地,2^k对于最终答案没有影响。
根据乘法原理,在能够分组(分成复共轭数)的前提下,最终的分组方法数 = 4*∏(pi+1)。
(这就是本题的做法。分解质因数,复杂度O(sqrt(N)))
三、积性函数χ(n),求π的表达式(这部分跟此题无关)
对于函数χ(n),定义为:
(1)n = 4k + 1时,χ(n) = 1
(2)n = 4k + 3时,χ(n) = -1
(3)n为偶数时,χ(n) = 0
函数χ(n)对于任意整数满足性质:χ(ab) = χ(a)*χ(b),所以χ(n)为积性函数。
将圆上点的数量写成如下形式:
即:N = 4*∏(∑ χ(ki)),ki为R^2的因子。
将上式拆开,每一项χ(n)的n为R的因子:
圆内所有点的个数:
移动之后:
所以得到了圆内点的个数,也就是圆面积的另一种表达形式。
最终得到了一个π的表达式。
AC Code:
#include <iostream>
#include <stdio.h>
#include <string.h> using namespace std; long long n;
long long ans=; int main()
{
cin>>n;
n=n*n;
long long t=n;
while(!(t&)) t>>=;
for(int i=;i*i<=n && t>;i++)
{
int p=;
while(t%i==)
{
p++;
t/=i;
}
if(i%==) ans*=(p+);
else if(i%== && (p&))
{
ans=;
break;
}
}
if(t%==) ans=;
cout<<ans*<<endl;
}
BZOJ 1041 [HAOI2008]圆上的整点:数学的更多相关文章
- bzoj 1041: [HAOI2008]圆上的整点 数学
1041: [HAOI2008]圆上的整点 Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://www.lydsy.com/JudgeOnline/ ...
- BZOJ 1041: [HAOI2008]圆上的整点【数论,解方程】
1041: [HAOI2008]圆上的整点 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 4210 Solved: 1908[Submit][Sta ...
- BZOJ 1041: [HAOI2008]圆上的整点
1041: [HAOI2008]圆上的整点 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 3621 Solved: 1605[Submit][Sta ...
- bzoj 1041: [HAOI2008]圆上的整点 本原勾股數組
1041: [HAOI2008]圆上的整点 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 2027 Solved: 853[Submit][Stat ...
- BZOJ 1041 [HAOI2008]圆上的整点:数学【费马平方和定理】
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1041 题意: 给定n(n <= 2*10^9),问你在圆x^2 + y^2 = n^ ...
- BZOJ(2) 1041: [HAOI2008]圆上的整点
1041: [HAOI2008]圆上的整点 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 4966 Solved: 2258[Submit][Sta ...
- 1041: [HAOI2008]圆上的整点
1041: [HAOI2008]圆上的整点 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 4298 Solved: 1944[Submit][Sta ...
- 【BZOJ】1041: [HAOI2008]圆上的整点(几何)
http://www.lydsy.com:808/JudgeOnline/problem.php?id=1041 所谓的神题,我不会,直接题解..看了半天看懂题解了.详见hzwer博客 这题呢,我只能 ...
- 1041: [HAOI2008]圆上的整点 - BZOJ
Description 求一个给定的圆(x^2+y^2=r^2),在圆周上有多少个点的坐标是整数.Input rOutput 整点个数Sample Input4Sample Output4HINT n ...
随机推荐
- 图像处理:卷积模块FPGA 硬件加速
本文记录了利用FPGA加速图像处理中的卷积计算的设计与实现.实现环境为Altera公司的Cyclone IV型芯片,NIOS II软核+FPGA架构. 由于这是第一次设计硬件加速模块,设计中的瑕疵以及 ...
- C# yield return用法
本文实例讲述了C#中yield return用法,并且对比了使用yield return与不使用yield return的情况,以便读者更好的进行理解.具体如下: yield关键字用于遍历循环中,yi ...
- web项目的发布
1.我们将web项目布署到IIS上之前需要将我们的源码发布,然后再将发布后的代码布署到II上 2.首先在VS上打开我们的应用程序,在WEB层上右键->发布,会弹出一个"发布web&qu ...
- MySQL索引实战经验总结
MySQL索引对数据检索的性能至关重要,盲目的增加索引不仅不能带来性能的提升,反而会消耗更多的额外资源,本篇总结了一些MySQL索引实战经验. 索引是用于快速查找记录的一种数据结构.索引就像是数据库中 ...
- html标签缺省(自带)样式大全
html标签默认样式整理 作者:佚名 来源:互联网 时间:07-30 16:54:48 文为大家整理了html标签默认样式属性及浏览器默认样式等等,喜欢css布局的朋友们可以学下,希望对大家有所帮助 ...
- 【2017集美大学1412软工实践_助教博客】团队作业4——第一次项目冲刺(Alpha版本)小组 成绩
第四次团队作业成绩公布 题目 团队作业4: http://www.cnblogs.com/happyzm/p/6722264.html 团队成绩 成绩公示如下: 检查项 会议内容 代码签入 心得体会或 ...
- 结对编程1---四则运算---GUI
201421123003 黄建英 201421123004 黄美海 题目描述: 我们在个人作业1中,用各种语言实现了一个命令行的四则运算小程序.进一步,本次要求把这个程序做成GUI(可以是Window ...
- 团队作业8----第二次项目冲刺(beta阶段)5.24
Day6-05.24 1.每日会议 会议内容: 1.组长林乔桦对昨日的工作进行了总结并且安排今日的任务. 2.阶段进入尾声,大家再一次集中对软件进行了优化讨论. 3.今天主要大家的工作重心放在异常的测 ...
- 201521123056 《Java程序设计》第1周学习总结
1. 本周学习总结 java语言结构简单,但功能齐全,同时能够在不同系统平台上编译,但编译的前提是系统要有安装JVM(即java虚拟机),JVM是java实现跨平台的最核心部分.本周还学习了JDK.J ...
- JavaSE集合(十)之Map
前面给大家介绍了集合家族中的Collection家族,这一篇给大家分享的是集合中的另一个家族就是Map家族.以前的时候学习Map的时候没有很认真的去学习,我觉得很多东西还是不是很清楚. 这次我将总结的 ...