Flooring Tiles 反素数

http://acm.hdu.edu.cn/showproblem.php?pid=4228

#include
#include
#include
using namespace std;

typedef __int64 lld;
lld p[1000];
int cnt;
lld prime[25]={2,3,5,7,11,13,17,19,23,29,31,37,41,53,59,61};

void Gao(int fac,int limit,lld s,int k)
{
    if(fac>150 || s>((lld)1<<60)) return ;
    if(p[fac]!=-1 && p[fac]>s) p[fac]=s;
    if(p[fac]==-1)p[fac]=s;
    if(k>16)return ;
    lld temp=s;
    for(int i=1;i<=limit;i++)
    {
        temp*=prime[k];
        if(temp>((lld)1<<60))break;
        Gao(fac*(i+1),i,temp,k+1);
    }
}

int main()
{
    int i;
    int n;
    cnt=0;
    memset(p,-1,sizeof(p));
    p[1]=1;
    Gao(1,75,1,0);
    while(scanf("%d",&n)!=EOF && n)
    {
        lld ans=((lld)1<<60);
        for(i=1;i<=150;i++)
        {
        //    printf("%d  %I64d\n",i,p[i]);
            if((i+1)/2==n && p[i]!=-1)
            {
                
                if(ans>p[i])
                {
                    ans=p[i];
                }
            }
        }
        printf("%I64d\n",ans);
    }
    return 0;
}

（转）对于任何正整数x,起约数的个数记做g(x).例如g(1)=1,g(6)=4.
如果某个正整数x满足:对于任意i(0 现在给一个N,求出不超过N的最大的反素数.
比如:输入1000  输出 840
思维过程:
求[1..N]中约数在大的反素数-->求约数最多的数
如果求约数的个数 756=2^2*3^3*7^1
(2+1)*(3+1)*(1+1)=24
基于上述结论,给出算法:按照质因数大小递增顺序搜索每一个质因子,枚举每一个质因子
为了剪枝:
性质一:一个反素数的质因子必然是从2开始连续的质数.
因为最多只需要10个素数构造:2,3,5,7,11,13,17,19,23,29
性质二:p=2^t1*3^t2*5^t3*7^t4.....必然t1>=t2>=t3>=....
http://apps.hi.baidu.com/share/detail/16212584