题目:http://www.codechef.com/NOV13/problems/PRETNUM/
Yet Another Cute Girl
思路:
一个数的因子的个数为素数,即这个数等于一个素数的n次方,并且n+1为素数。L和R的最大值为10^12,要判断大数为素数,也可以用筛法,先把10^6内的素数打表,然后分别把这些素数累乘,设10^6<L<=R ,则此时n>=2,找到L到R之间的符合条件的数,并用数组d[ ]标记,再用筛法把L到R之间的数作上标记,最后找到那些未被标记的即为大素数。
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<stdlib.h>
#include<stack>
#include<iostream>
using namespace std;
typedef long long LL;
LL prime[100000],vis[1000001],c;
void pri()
{
LL i,j;
prime[0]=2;
vis[0]=1;vis[1]=1;
for(j=4;j<=1000000;j+=2)
vis[j]=1;
c=1;
for(i=3;i<=1000000;i+=2)
{
if(vis[i]==0)
{
prime[c++]=i;
for(j=i*i;j<=1000000;j+=i)
vis[j]=1;
}
}
}
LL d[1000009];
int main()
{
LL t,l,r,i,j,s;
scanf("%lld",&t);
pri();
while(t--)
{
for(i=0;i<=r-l;i++)d[i]=0;
s=0;
scanf("%lld%lld",&l,&r);
for(i=0;i<c;i++)
{
int time=0;
j=1;
while(j<=r)
{
if(j>=l)
{
if(vis[time+1]==0)s++;
d[j-l]=1;
}
j*=prime[i];
time++;
}
}
i=1000001;
if(i<l)i=l;
LL k;
if(i<=r)
for(j=0;j<c;j++)
{
if(i%prime[j]==0)k=i;
else k=(i/prime[j]+1)*prime[j];
while(k<=r)
{
d[k-l]=1;
k+=prime[j];
}
}
for(j=i;j<=r;j++)if(d[j-l]==0)s++;
printf("%lld\n",s);
}
return 0;
}