2018.09.11 bzoj3629: [JLOI2014]聪明的燕姿（搜索）

传送门

一道神奇的搜索。

直接枚举每个质因数的次数，然后搜索就行了。

显然质因数k次数不超过logkn" role="presentation" style="position: relative;">logknlogkn，因此搜索很快。

注意，如果最后剩下的乘积-1是一个质数，那么这是一个可行解。

另外对于这道题，我们只需要筛出1e5的素数就够了，太大的直接枚举自己打的素数表判就行了。

代码：

#include<bits/stdc++.h>
#define N 100005
using namespace std;
int pri[N],tot=0,s,cnt=0,ans[N];
bool vis[N];
inline void init(int len){
    for(int i=2;i<=len;++i){
        if(!vis[i])pri[++cnt]=i;
        for(int j=1;j<=cnt;++j){
            if(i*pri[j]>len)break;
            vis[i*pri[j]]=true;
            if(i%pri[j]==0)break;
        }
    }
}
inline bool check(int x){
    if(x==1)return false;
    if(x<=100000)return !vis[x];
    for(int i=1;pri[i]*pri[i]<=x;++i)if(x%pri[i]==0)return false;
    return true;
}
inline void dfs(int lastpri,int mul,int lim){
    if(lim==1){ans[++tot]=mul;return;}
    if(lim-1>pri[lastpri]&&check(lim-1))ans[++tot]=mul*(lim-1);
    for(int i=lastpri+1;pri[i]*pri[i]<=lim;++i){
        for(int mult=pri[i]+1,tmp=pri[i];mult<=lim;mult+=(tmp*=pri[i]))
            if(lim%mult==0)dfs(i,mul*tmp,lim/mult);
    }
}
int main(){
    init(100000);
    while(~scanf("%d",&s)){
        tot=0,dfs(0,1,s),printf("%d\n",tot),sort(ans+1,ans+tot+1);
        for(int i=1;i<=tot;++i)printf("%d%c",ans[i],i==tot?'\n':' ');
    }
    return 0;
}