codechef (大数据素数筛法）

题目：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;
}