bzoj 1225: [HNOI2001] 求正整数

题意：

求出具有n个不同因子的最小正整数m。

题解:

显然可以枚举因数爆搜，然而不会比大小。
于是学了对数，很好用。
这篇不错
code：

#include<cmath>#include<cfloat>#include<cstdio>#include<cstdlib>#include<cstring>#include<iostream>using namespace std;int n,ans[100005],res[21],tmp[21],pri[]={0,2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53};double mn=DBL_MAX,lg[21];void print(){    ans[0]=ans[1]=1;    for(int i=1;i<=16;i++)        for(;res[i]>0;res[i]--)        {            for(int j=1;j<=ans[0];j++)                ans[j]*=pri[i];            for(int j=1;j<=ans[0];j++)                ans[j+1]+=ans[j]/10,ans[j]%=10;            if(ans[ans[0]+1]!=0)                ans[0]++;            while(ans[ans[0]]/10!=0)                ans[ans[0]+1]+=ans[ans[0]]/10,ans[ans[0]]%=10,++ans[0];        }    for(int i=ans[0];i>=1;i--) printf("%d",ans[i]);    printf("\n");}void dfs(double x,int y,int z){    if(x>mn) return;    if(y==1)    {        mn=x;        memset(res,0,sizeof(res));        for(int i=1;i<z;i++) res[i]=tmp[i];        return;    }    if(z>16) return;    for(int i=0;(i+1)*(i+1)<=y;i++)        if(y%(i+1)==0)        {            if(i!=0)            {                tmp[z]=i;                dfs(x+lg[z]*i,y/(i+1),z+1);            }            if((i+1)*(i+1)!=y)            {                tmp[z]=y/(i+1)-1;                dfs(x+lg[z]*(y/(i+1)-1),i+1,z+1);            }        }}int main(){    scanf("%d",&n);    for(int i=1;i<=16;i++)        lg[i]=log(pri[i]);    dfs(0,n,1);    print();}