题目链接: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]圆上的整点:数学的更多相关文章

  1. bzoj 1041: [HAOI2008]圆上的整点 数学

    1041: [HAOI2008]圆上的整点 Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://www.lydsy.com/JudgeOnline/ ...

  2. BZOJ 1041: [HAOI2008]圆上的整点【数论,解方程】

    1041: [HAOI2008]圆上的整点 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 4210  Solved: 1908[Submit][Sta ...

  3. BZOJ 1041: [HAOI2008]圆上的整点

    1041: [HAOI2008]圆上的整点 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 3621  Solved: 1605[Submit][Sta ...

  4. bzoj 1041: [HAOI2008]圆上的整点 本原勾股數組

    1041: [HAOI2008]圆上的整点 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 2027  Solved: 853[Submit][Stat ...

  5. BZOJ 1041 [HAOI2008]圆上的整点:数学【费马平方和定理】

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1041 题意: 给定n(n <= 2*10^9),问你在圆x^2 + y^2 = n^ ...

  6. BZOJ(2) 1041: [HAOI2008]圆上的整点

    1041: [HAOI2008]圆上的整点 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 4966  Solved: 2258[Submit][Sta ...

  7. 1041: [HAOI2008]圆上的整点

    1041: [HAOI2008]圆上的整点 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 4298  Solved: 1944[Submit][Sta ...

  8. 【BZOJ】1041: [HAOI2008]圆上的整点(几何)

    http://www.lydsy.com:808/JudgeOnline/problem.php?id=1041 所谓的神题,我不会,直接题解..看了半天看懂题解了.详见hzwer博客 这题呢,我只能 ...

  9. 1041: [HAOI2008]圆上的整点 - BZOJ

    Description 求一个给定的圆(x^2+y^2=r^2),在圆周上有多少个点的坐标是整数.Input rOutput 整点个数Sample Input4Sample Output4HINT n ...

随机推荐

  1. JavaScript学习日志(六):事件

    这篇随笔,深恶痛绝,敲到快结束的时候,凌晨00:19,突然闪退,也不知道是Mac的原因还是chrome的原因,重新打开的时候,以为自动保存有效果,心想没关系,结果他么的只保存了四分之一,WTF?!!! ...

  2. 【Java线程】Java内存模型总结

    学习资料:http://www.infoq.com/cn/articles/Java-memory-model-1 Java的并发采用的是共享内存模型(而非消息传递模型),线程之间共享程序的公共状态, ...

  3. js Web存储方式

    JSON是数据交互中最常用的一种数据格式. 由于各种语言的语法都不同,在传递数据时,可以将自己语言中的数组.对象等转换为JSON字符串> 传递之后,可以讲JSON字符串,在解析为JSON对象. ...

  4. 集美大学网络1413第十二次作业成绩(个人作业3) -- Alpha阶段个人总结

    题目 个人作业3--个人总结(Alpha阶段) 优秀作业链接:**068 未交:**087 个人作业3成绩 学号 姓名 总结(4) 5个问题(2.5) 自我评价(2.5) 博客互动 (1) 总分(10 ...

  5. 201521123040《Java程序设计》第8周学习总结

    1. 本周学习总结 1.1 以你喜欢的方式(思维导图或其他)归纳总结集合与泛型相关内容. 2. 书面作业 本次作业题集集合 1.List中指定元素的删除(题目4-1) 1.1 实验总结 我们每次删除元 ...

  6. 201521123059 《Java程序设计》第二周学习总结

    1.本周总结 本周老师讲了和自己掌握了以下内容: (1).三元条件运算符 表达式1?表达式2:表达式3: (2). 字符串String类 String的不可变优点:编译器可以让字符串共享,效率高.但是 ...

  7. 201521123055《Java程序设计》第1周学习总结

     1. 本章学习总结 (1)JAVA环境配置(JDK,JVM) (2)编写简易程序熟练代码结构  2. 书面作业 1.为什么java程序可以跨平台运行?执行java程序的步骤是什么? JAVA程序需要 ...

  8. POJ-2299 Ultra-QuickSort (树状数组,离散化,C++)

    Problem Description In this problem, you have to analyze a particular sorting algorithm. The algorit ...

  9. linux(4) vi编辑/删除、复制、粘贴 /bash shell 环境变量设置/数据流重定向 | 的用法

    一.vi文字处理器1.vi与vimvi:文字处理器vim:程序开发工具2.vi介绍三种模式:一般模式(vi刚进入的,不可编辑),编辑模式(按i后,左下方是insert)和命令行模式(按esc退出,:w ...

  10. 国际化:Java平台下的Locale类

    Locale类有如下几个构造器: Locale(String language) Locale(String language, String country) Locale(String langu ...