bzoj 3629 [JLOI2014]聪明的燕姿（约数和，搜索）

【题目链接】

　　http://www.lydsy.com/JudgeOnline/problem.php?id=3629

【题意】

　　给定S，找出所有约数和为S的数。

【思路】

　　若n=p1^a1*p2^a2*...*pk^ak

　　则约数和f(n)为(p1^0+p1+p1^2+...+p1^a1)*(p2^0+p2+p2^2+...+p2^a2)*...*(pk^0+pk+pk^2+...+pk^ak)

　　考虑搜索，使得和为S。至于这个搜索怎么写的，我能说我看不懂吗=_=

【代码】

 #include<cmath>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 using namespace std;

 typedef long long ll;
 const int N = 1e5+;

 ll su[N],vis[N]; ll ans[N];

 void get_prime()
 {
     vis[]=;
     for(int i=;i<N;i++) {
         if(!vis[i]) su[++su[]]=i;
         for(int j=;j<=su[]&&su[j]*i<N;j++) {
             vis[i*su[j]]=;
             if(i%su[j]==) break;
         }
     }
 }

 bool is_prime(ll x)
 {
     if(x<N) return !vis[x];
     for(ll i=;su[i]*su[i]<=x;i++)
         if(x%su[i]==) return ;
     return ;
 }

 void dfs(ll now,int pos,ll left)
 {
     if(left==) {
         ans[++ans[]]=now;
         return ;
     }
     if( left->=su[pos] && is_prime(left-) )
         ans[++ans[]]=(left-)*now;
     for(int i=pos; su[i]*su[i]<=left ;i++) {
         ll sum=su[i]+,pow=su[i];
         for(;sum<=left;pow*=su[i],sum+=pow)
             if(left%sum==)
                 dfs(now*pow,i+,left/sum);
     }
 }

 ll S;

 int main()
 {
     get_prime();
     while(scanf("%lld",&S)!=EOF) {
         ans[]=;
         dfs(,,S);
         sort(ans+,ans+ans[]+);
         printf("%d\n",ans[]);
         for(int i=;i<=ans[];i++)
             printf("%lld%c",ans[i],i==ans[]?'\n':' ');
     }
     return ;
 }