标签：limits 高斯 P2508 luogu 质数 整数 4k 圆上 prod
圆上的整点
题目链接：luogu P2508
题目大意
给你一个圆，问你圆周上有多少个点的坐标是整点。
思路
考虑一个东西叫做高斯整数。
其实它是复数，是 \(a+bi\) 中 \(a,b\) 都是整数的复数。
那它跟它共轭的乘积其实就是 \(a^2+b^2\)，所以我们可以把它转化成 \(a^2+b^2=N\) 这个东西，满足条件的高斯整数个数。
那既然是这样我们就考虑这个 \(N\) 要怎么在高斯整数上分解。
那在整数范围有唯一的分解，就是质因数分解嘛。
不过准确一点其实可以说是不唯一，因为你可以同时给两个质因数取反，得到的还是可以，毕竟负数也是整数。
那在高斯划分中，你也可以划分成一些高斯整数，它在高斯整数上不能再分，那这些就是高斯质数。
在弄 \(-1\) 自然可以，你会发现你甚至可以弄 \(i\)（准确来讲是一个 \(i\) 一个 \(-i\)，因为 \(i(-i)=-(-1)=1\)）。
那么就是一个问题了，怎么知道一个高斯整数是不是高斯质数。
我们可以用一个叫做费马平方和的定理。
费马平方和定理：奇素数 \(p\) 可以表示乘两个数的平方和当且仅当 \(p\) 是形如 \(4k+1\)，其中 \(k\) 为整数。
不考虑两个正整数顺序的时候方法唯一。
（有人不会证，不过想看证明的可以去看百度百科，看着好像不难懂）
那我们就分类讨论一下，对于质数分成 \(2,4k+1,4k+3\)。
- \(4k+3\) 是高斯质数，那就对应的圆上没有整点。
- \(4k+1\) 是恰好可以被分成一对共轭负数的乘积，那就对应的圆上有整点（至于有多少个我们后面再看）
- \(2\) 是可以分成 \((1+i)(1-i)\)，而且特殊的是这两个成 \(90°\)（为啥特殊后面会用到）
然后考虑一个数质因数分解之后要怎么弄，我们考虑把三种质数分开来：
\(N=2^p\prod\limits_{a_i=4k+3}a_i^{m_i}\prod\limits_{b_i=4k+1}b_i^{n_i}\)
其中 \(2,b_i\) 是可以分解的，那问题是我们要把 \(N\) 分成一对共轭质数。
我们思考有怎样的分配方式，一种是能分解成共轭的，一边一个，要么是两个一样的，一遍给一个（因为同乘也是可以的，相当于把圆放大罢了，对于的位置还是在整点）
对于每个分类讨论：
- \(4k+1\)
我们可以把它分解两种高斯质数，那如果是 \(n_i\) 个这样的，那我们可以选择给左边 \(0\sim n_i\) 个，所以方案是 \(n_i+1\)。 - \(4k+2\)
这种不能拆，就只能给两边，所以如果 \(m_i\) 是奇数就不行整个答案是 \(0\)，否则就不变。 - \(2\)
那看起来也是跟 \(4k+1\) 一样？
然后发现我们漏了一个问题，你上面这个拆只是在正整数的位置，你总的位置还有四个呢！
那答案就要乘上 \(4\)，如果在二维平面上形象地表示一下的话，就是每次转 \(90\) 度的感觉。
但是这个时候 \(2\) 就有问题了，你这样搞就会重复，具体来讲只有四个不一样的，所以你可以相当于 \(2\) 不用管，最后答案乘 \(4\) 即可。
然后注意到题目是 \(x^2+y^2=r^2\)。
所以 \(N=r^2\)，当然直接单纯的这样是不行的。
考虑看看式子：
\(r=2^p\prod\limits_{a_i=4k+3}a_i^{m_i}\prod\limits_{b_i=4k+1}b_i^{n_i}\)
\(N=2^r=2^{2p}\prod\limits_{a_i=4k+3}a_i^{2m_i}\prod\limits_{b_i=4k+1}b_i^{2n_i}\)
改一改就好了。
代码
#include<cstdio>
#define ll long long
using namespace std;
int n;
ll ans;
void slove(int x, int num) {
if (x == 2) {
return ;
}
if (x % 4 == 3) {
if ((2 * num) & 1) ans = 0; return ;
}
if (x % 4 == 1) {
ans *= (2 * num + 1); return ;
}
}
int main() {
scanf("%d", &n); ans = 1;
for (int i = 2; i * i <= n; i++) {
if (n % i != 0) continue;
int num = 0; while (n % i == 0) n /= i, num++;
slove(i, num);
}
if (n > 1) slove(n, 1);
printf("%lld", ans * 4);
return 0;
}
标签：limits,高斯,P2508,luogu,质数,整数,4k,圆上,prod 来源： https://www.cnblogs.com/Sakura-TJH/p/luogu_P2508.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。