1. 程式人生 > >HDU 4548 美素數 素數題解

HDU 4548 美素數 素數題解

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 美素數 素數題解