HDU 4548 美素數 素數題解
阿新 • • 發佈:2018-04-24
int 思想 csdn prime ali menu clu true max
本題就是能夠直接打表的,推斷能否夠打表也須要技巧的:
1 推斷最大的數值為1000000。百萬下面的數打表都是能夠的
2 能夠線性預處理好。使用素數篩子法是能夠接近線性預處理的。
故此能夠打表了。
須要熟悉的基本知識點:
1 素數篩子法 - 一兩分鐘之內寫出代碼
2 一般素數推斷法,由於位數相加之後的數值很小,故此一般素數推斷就能夠了,假設寫個primality test 算法會大材小用了。
3 然後是帶點動態規劃法的思想把前面的美素數疊加起來,方便查找。
算是基礎題目了,也是有人說的水題,我還是喜歡叫基礎題吧。
難在於推斷好,並運用好這些基礎知識。簡單優雅地解決,寫出代碼。
#include <stdio.h> #include <string.h> const int MAX_N = 1000001; bool isPrime(int n) { if (n == 2) return true; for (int r = 2; r * r <= n; r++) { if (n % r == 0) return false; } return true; } bool isMeiPrime(int n) { int d = 0; while (n) { d += n % 10; n /= 10; } return isPrime(d); } int primeNums[MAX_N]; bool primes[MAX_N]; void seive() { memset(primes, 0, MAX_N * sizeof(bool)); for (int i = 2; i < MAX_N; i++) { if (!primes[i]) { for (int j = i << 1; j < MAX_N; j += i) { primes[j] = true; } } } primeNums[0] = 0, primeNums[1] = 0; for (int i = 2; i < MAX_N; i++) { if (!primes[i] && isMeiPrime(i)) primeNums[i] = primeNums[i-1] + 1; else primeNums[i] = primeNums[i-1]; } } int main() { seive(); int T, a, b; scanf("%d", &T); for (int t = 1; t <= T; t++) { scanf("%d %d", &a, &b); printf("Case #%d: %d\n", t, primeNums[b] - primeNums[a-1]); } return 0; }
HDU 4548 美素數 素數題解