标签:style blog http color 2014 os
题意:给定a + bi,判断是否是高斯素数,i = sqrt(-2)。
思路:普通的高斯素数i = sqrt(-1),判断方法为:
1、如果a或b为0,判断另一个数为4 * n + 3形式的素数(用到费马平方和定理)
2、如果a、b都不为0,判断a ^ 2 + b ^ 2 是否为素数
那么这题,提取出sqrt(2)来,就和基本情况一样了。
对于2,变成: 如果a、b都不为0,判断a ^ 2 + 2 b ^ 2是否为素数
对于1有点不同,变成x ^ 2 + 2 y ^ 2形式,不在适用费马平方和定理,由于a,b最大才1W,完全可以去枚举,判断能不能拆成x ^ 2 + 2 * y ^ 2形式。这样一来问题就解决了
代码:
#include <stdio.h> #include <string.h> #include <math.h> const int N = 20005; int t, a, b, prime[N], pn = 0, vis[N]; bool judge() { if (a == 0) { int m = sqrt(m); for (int i = 0; i <= m; i++) { double t = sqrt((b - i * i) * 1.0 / 2); if (ceil(t) == floor(t)) return true; } return false; } int tmp = a * a + 2 * b * b; for (int i = 0; i < pn && prime[i] < tmp; i++) if (tmp % prime[i] == 0) return false; return true; } int main() { vis[1] = 1; for (int i = 2; i < N; i++) { if (vis[i]) continue; prime[pn++] = i; for (int j = i * i; j < N; j += i) { vis[j] = 1; } } scanf("%d", &t); while (t--) { scanf("%d%d", &a, &b); printf("%s\n", judge()?"Yes":"No"); } return 0; }
UVA 1415 - Gauss Prime(数论,高斯素数拓展),布布扣,bubuko.com
UVA 1415 - Gauss Prime(数论,高斯素数拓展)
标签:style blog http color 2014 os
原文地址:http://blog.csdn.net/accelerator_/article/details/36015937