数论渣表示已经报警了。
对于任意一个数x,我们可以将其进行质因数分解。
x=p1^a1*p2^a2*…*pn^an.
那么显然x的约数个数为(a1+1)(a2+1)…*(an+1)。
那么x的约数个数和为(p1^0+p1^1+…p1^a1)(p2^0+p2^1+…p2^a2)…(pn^0+pn^1+…pn^an).
那么我们就可以进行搜索了。枚举每个质数pi和合法的ai。
left=n/当前数字约数和。
当left==1时表示当前数字约数和==n也就是一个合法的解。
当left-1是一个素数时,且这个素数没有被枚举过,now*(left-1)也是一个合法的解(将这个质因子直接加入答案)。
#include<iostream>
#include<cstdio>
#include<algorithm>
#define MAXN 100005
using namespace std;
int n,prime[MAXN],ans[MAXN],tot;
bool flag[MAXN];
inline void get_prime()
{
for (int i=2;i<=MAXN;i++)
{
if (!flag[i]) prime[++prime[0]]=i;
for (int j=1;i*prime[j]<=MAXN&&j<=prime[0];j++)
{
flag[i*prime[j]]=1;
if (i%prime[j]==0) break;
}
}
}
bool is_prime(int x)
{
if (x==1) return 0;
if (x<=MAXN) return !flag[x];
for (int i=1;prime[i]*prime[i]<=x;i++)
if (x%prime[i]==0) return 0;
return 1;
}
void dfs(int now,int pos,int left)
{
if (left==1) {ans[++ans[0]]=now; return;}
if (left-1>=prime[pos]&&is_prime(left-1))
ans[++ans[0]]=(left-1)*now;
for (int i=pos;prime[i]*prime[i]<=left;i++)
for (int power_sum=prime[i]+1,power=prime[i];power_sum<=left;power*=prime[i],power_sum+=power)
if (left%power_sum==0) dfs(now*power,i+1,left/power_sum);
}
int main()
{
get_prime();
while (scanf("%d",&n)!=EOF)
{
ans[0]=0;
dfs(1,1,n);
sort(ans+1,ans+ans[0]+1);
cout << ans[0] << endl;
for (int i=1;i<=ans[0];i++)
printf("%d%c",ans[i],i==ans[0]?'\n':' ');
}
return 0;
}