[bzoj3629][JLOI2014]聪明的燕姿（数论+dfs）

题目：

我是超链接

题解：

有一个东西称为约数和定理
分解质因数：$n=\prod_{i=1} p_i^{k_i}$
约数和：

$$d(n)=\prod_{i=1} (\sum_{j=0}^{k_i}p_i^{k_j})$$
基本到这里就可以运用搜索了。。。
①若当前数可表示成一个并未搜索过的质数与1的和，则之前搜索过的数与这个质数的乘积符合题意。
②对于每一个未被搜索过且平方小于当前数的质数，则枚举所有可能符合题意的素数进行递归搜索。
Emmm…输出还是要乖乖的，不然容易PE（比如最后多个空格之类的

代码：

#include <cstdio>
#include <algorithm>
using namespace std;
const int N=1e5;
const int Q=1e6;
int p[N+5],ans[Q],cnt;
bool ss[N+5];
void pre()
{
    int ll=0;
    for (int i=2;i<=N;i++)
    {
        if (!ss[i]) p[++ll]=i;
        for (int j=1;j<=ll && p[j]*i<=N;j++)
        {
            ss[p[j]*i]=1;
            if(i%p[j]==0) break;  
        }  
    }
}
bool Isprime(int n)
{
    if (n==1) return 0;
    if (n<=N) return !ss[n];
    for (int i=1;p[i]*p[i]<=n;i++)
      if (n%p[i]==0) return 0;
    return 1; 
}
void dfs(int last,int now,int tot)
{
    if (tot==1) {ans[++cnt]=now; return;}
    if (tot-1>p[last] && Isprime(tot-1)) ans[++cnt]=now*(tot-1);
    for (int i=last+1;p[i]*p[i]<=tot;i++)
      for (int sum=p[i]+1,a=p[i];sum<=tot;a=a*p[i],sum+=a)
        if (tot%sum==0) dfs(i,now*a,tot/sum);
}
int main()
{
    int s;
    pre();
    while (scanf("%d",&s)!=EOF)
    {
        cnt=0;
        dfs(0,1,s);
        sort(ans+1,ans+cnt+1);
        printf("%d\n",cnt);
        for (int i=1;i<=cnt;i++) printf("%d%c",ans[i],i==cnt?'\n':' ');
    }
}