TC SRM 552 DIV1 100PT(数论)

转载请注明出处，谢谢 http://blog.csdn.net/ACM_cxlove?viewmode=contents           by---cxlove

题目：对于一个数，如果所有的素数要么不是他的因子，要么就得小于某个上界并且这个数拥有奇数个这个素因子。

比如16就不满足，2^4,4为偶数。问1---UPTO范围里有多少个这样的数。

首先把范围内的素数表打出。

然后开始搜索素数，对于某个素数，只能取奇数个，或者不取。

看代码注释

bool flag[1000005]={0};
vector<int>v;
class HolyNumbers{
    public:
    //预处理素数表
    void Init(int n){
        for(int i=2;i<=sqrt(n+1.0);i++){
            if(flag[i]) continue;
            for(int j=2;j*i<=n;j++)
                flag[i*j]=true;
        }
        for(int i=2;i<=n;i++)
            if(!flag[i])
               v.push_back(i);
    }
    //搜索到第idx个素数,这时候的范围是 1--up
    LL slove(int idx,LL up){
        //如果素数没了，而且当前及以后的素数都超过范围，则返回1
        //这个1表示之前搜的素数组成的数刚好就是要求的数
        //初始的话就是返回1就代表1是满足条件的
        if(idx>=v.size()||v[idx]>up)  return 1;
        if((LL)v[idx]*v[idx]>=up){
            //这里说明这个素因子以及后面的素因子都不可能取3个
            //那么求出有多少个素因子可以取，都是取1个
            int k=upper_bound(v.begin()+idx,v.end(),up)-v.begin();
            return k-idx+1;
        }
        LL p=v[idx];
        //不取这个素因子
        LL ans=slove(idx+1,up);
        //取奇数个
        for(LL num=p;num<=up&&num>0;num*=p*p)
            ans+=slove(idx+1,up/num);
        return ans;
    }
    LL count(LL upTo, int maximalPrime){
        v.clear();
        Init(maximalPrime);
        return slove(0,upTo);
    }
};