BZOJ 3629: [JLOI2014]聪明的燕姿

#include<cstdio>
#include<algorithm>
#define N 50005
using namespace std;

int t,S,cnt;
int mark
,prime
,ans
;

void init(int n)
{
for(int i=2;i<=n;++i)
{
if(!mark[i]) prime[++cnt]=i,mark[i]=1;
for(int j=1;j<=cnt&&prime[j]*i<=n;++j)
{
mark[prime[j]*i]=1;
if(i%prime[j]==0) break;
}
}
}

int check(int x)
{
if(x==1) return 0;
for(int i=1;prime[i]*prime[i]<=x;++i)
if(x%prime[i]==0) return 0;
return 1;
}

void dfs(int num,int pre,int now)
{
if(num==1) {ans[++ans[0]]=now;return;}
if(num-1>prime[pre]&&check(num-1)) ans[++ans[0]]=now*(num-1);
for(int i=pre+1;prime[i]*prime[i]<=num;++i)
{
int p=prime[i],sum=p+1;
while(sum<=num)
{
if(num%sum==0)
dfs(num/sum,i,now*p);
p*=prime[i];
sum+=p;
}
}
}

int main()
{
init(50000);
//  scanf("%d",&t);
//  while(t--)
while(scanf("%d",&S)!=EOF)
{
ans[0]=0;
dfs(S,0,1);
printf("%d\n",ans[0]);
sort(ans+1,ans+ans[0]+1);
for(int i=1;i<=ans[0];++i)
{
printf("%d",ans[i]);
if(i==ans[0]) printf("\n");
else printf(" ");
}
}
}