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);    }};