51nod网站链接:士兵的数字游戏
看懂题目后就知道,其实这道题就是求区间[1,a]内的所有数的质因子个数之和。
方法利用线性筛素数的方法,可以以O(n)的复杂度求出所有素数。
同时对程序稍作变动,就可以在求素数的同时求质因子个数,并累加
#include<stdio.h>
//using namespace std;
const int MaxPrimNum = 183072, MaxSize = 5000001; //, MaxSize1 = 31251
int MinPrime[MaxSize]; //分别记录可以整除MaxSize的最小与最大prime的index //MaxPrime可以改成bit类型数组来记录是否是prime或者是否是Prime幂次方
int prime[MaxPrimNum],result[MaxSize];
int j, len, num_prime = 1, T, tmpresult;
__int64 i, tmp, tmp1;
int tmpsumresult = 0;
int main()
{
// result[1] = 0;
// tmp1 = (MaxSize >> 1) + 1;
for (i = 2; i < MaxSize/*tmp1*/; ++i)
{
if (MinPrime[i])
{
len = MinPrime[i];
tmpresult = result[i] + 1;
tmpsumresult = result[i] += tmpsumresult;
}
else
{
prime[len = MinPrime[i] = num_prime++] = i;
tmpresult = 2;
result[i] = ++tmpsumresult;
}
for (j = 1; j <= len && (tmp = i * prime[j]) < MaxSize; ++j)
{
MinPrime[tmp] = j;
result[tmp] = tmpresult;
}
}
scanf("%d", &T);
for (i = 0; i < T; ++i)
{
scanf("%d%d", &tmp, &tmp1);
printf("%d\n", result[tmp]-result[tmp1]);
}
return 0;
}