bzoj 3629

给出数 $n$
记 $f(x)$ 表示 $x$ 的因子和
求出所有 $x$ 使得 $f(x) = n$
考虑 $x = p_1 ^{a_1} * p_2 ^ {a_2} * \cdots * p_k ^ {a_k}$
那么 $f(x) = (1 + p_1 + p_1 ^ 2 + \cdots + p_1 ^ {a_1}) * (1 + p_2 + p_2 ^ 2 + \cdots + p_2 ^ {a_2}) * \cdots * (1 + p_k + p_k

^ 2 + \cdots + p_k ^ {a_k})$
因此可以爆搜枚举 $p, a$
Dfs(Now_result, prime_pos, x_remind) {}
分别表示
1.当前结果，即枚举到的素数的指数次幂的乘积，即 $p_i ^ {a_i}$ 的乘积　
2.当前枚举到的素数位置 首先要求出 $1e5, \sqrt(limit)$ 内的素数
3.给出的 $x$ 在枚举了之前的数后还剩多少
对于答案的录入
1.如果 x_remind = 1 ,相当于枚举到了这样一种形式 $() * () * () * p_i ^ 0$, 显然当前 Now_result 可以录入.
2.如果 x_remind - 1 是一个 大于等于 Prime[Prime_pos] 的素数，显然 (x_remind - 1) * Now_result 可以录入.
考虑这样的话我们已经枚举到了这样的一种形式 $() * () * () * (p_i + 1)$, 所以还原之前的数就是 p_i * Now_result.

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>

using namespace std;

#define gc getchar()

inline int read() {
    int x = ; char c = gc;
    while(c < '' || c > '') c = gc;
    while(c >= '' && c <= '') x = x *  + c - '', c = gc;
    return x;
}

const int N = 1e5 + ;

#define LL long long

LL Prime[N], Mark[N], Ans_js;
LL Answer[N], tot, n;

void Get_prime() {
    for(int i = ; i <= N - ; i ++) {
        if(!Mark[i]) Prime[++ tot] = i;
        for(int j = ; j <= tot && Prime[j] * i < N; j ++) {
            Mark[Prime[j] * i] = ;
            if(i % Prime[j] == ) break;
        }
    }
}

bool Is_prime(LL x) {
    if(x == ) return ;
    for(int i = ; Prime[i] * Prime[i] <= x; i ++) {
        if(x % Prime[i] == ) return ;
    }
    return ;
}

void Dfs(LL Now_result, int prime_pos, LL x_remind) {
    if(x_remind == ) {
        Answer[++ Ans_js] = Now_result;
        return ;
    }
    if(x_remind -  >= Prime[prime_pos] && Is_prime(x_remind - )) {
        Answer[++ Ans_js] = (x_remind - ) * Now_result;
    }
    for(int i = prime_pos; Prime[i] * Prime[i] <= x_remind; i ++) {
        for(LL tmp = Prime[i] + , imp = Prime[i]; tmp <= x_remind; imp *= Prime[i], tmp += imp) {
            if(x_remind % tmp == ) Dfs(Now_result * imp, i + , x_remind / tmp);
        }
    }
}

int main() {
    Get_prime();
    while(scanf("%lld", &n) == ) {
        Ans_js = ;
        Dfs(, , n);
        sort(Answer + , Answer + Ans_js + );
        cout << Ans_js << "\n";
        for(int i = ; i < Ans_js; i ++) cout << Answer[i] << " ";
        if(Ans_js) cout << Answer[Ans_js] << "\n";
    }
    return ;
}