首页 > 代码库 > UVa 10539 - Almost Prime Numbers

UVa 10539 - Almost Prime Numbers

题目:Almost Prime Numbers是只有一个素数因子的数,统计一个区间中的Almost Prime Numbers。

分析:数论,分治。根据定义可知Almost Prime Numbers就是素数的幂(大于1次)。

            首先,利用筛法计算出1000000内的素数(至少的素数的平方);

            然后,计算所有的素数的小于1000000000000的超过2次的幂,排序;

            最后利用二分查找,找到不超过当前数的下标,做差即可。

说明:打表计算输出,直接计算会TLE。

#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>

using namespace std;

long long ampr[100001];
long long prim[100001];
int 	  used[1000001];

int bs(int r, long long key)
{
	int l = 0,mid;
	while (l < r) {
		mid = (l+r+1)>>1;
		if (ampr[mid] > key)
			r = mid-1;
		else l = mid;
	}
	return l;
}

bool cmp(long long a, long long b)
{
	return a < b;
}

int main()
{
	memset(used, 0, sizeof(used));
	int count = 0;
	for (int i = 2 ; i < 1000000 ; ++ i)
		if (!used[i]) {
			prim[count ++] = i;
			for (int j = i<<1 ; j < 1000000 ; j += i)
				used[j] = 1;
		}
	
	int number = 0;
	for (int i = 0 ; i < count ; ++ i) {
		long long base = prim[i]*prim[i]+0LL;
		while (base <= 1000000000000LL) {
			ampr[number ++] = base;
			base = base*prim[i];
		}
	}
	
	sort(ampr, ampr+number, cmp);
	
	int n,x,y;
	while (cin >> n)
	while (n --) {
		long long a,b;
		cin >> a >> b;
		x = bs(number-1, a);
		y = bs(number-1, b);
		if (ampr[x] == a) x --;
		if (ampr[0] > a) x --;
		if (ampr[0] > b) y --;
		cout << y-x << endl;
	}
	return 0;
}

UVa 10539 - Almost Prime Numbers