码迷,mamicode.com
首页 > 其他好文 > 详细

bzoj千题计划297:bzoj3629: [JLOI2014]聪明的燕姿

时间:2018-03-20 20:42:21      阅读:169      评论:0      收藏:0      [点我收藏+]

标签:amp   prim   .com   body   分解   sig   +=   style   tar   

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

 

约数和定理:

若n的标准分解式为 p1^k1 * p2^k2 ……

那么n的约数和= π (Σ pi^xi ) xi∈[0,ki]

原本枚举小于S的质数,通过先判断S-1是不是质数 就可以 枚举根号S内的质数

 

#include<cstdio>
#include<algorithm>

using namespace std;

#define N 1000000

int prime[N+1],cnt;
bool vis[N+1];

int ans[N],tot;

void pre()
{
    vis[1]=true;
    for(int i=2;i<=N;++i)
    {
        if(!vis[i]) prime[++cnt]=i;
        for(int j=1;j<=cnt;++j)
        {
            if(i*prime[j]>N) break;
            vis[i*prime[j]]=true;
            if(!(i%prime[j])) break;
        }
    }
}

bool isprime(int x)
{
    if(x<=N) return !vis[x];
    for(int i=1;prime[i]*prime[i]<=x;++i)
        if(!(x%prime[i])) return false;
    return true;
}

void dfs(int last,int num,int rest)
{
    if(rest==1) { ans[++tot]=num; return; }
    if(rest-1>prime[last] && isprime(rest-1)) ans[++tot]=num*(rest-1);
    for(int i=last+1;prime[i]*prime[i]<=rest;++i)
        for(int sum=prime[i]+1,nnum=prime[i];sum<=rest;nnum*=prime[i],sum+=nnum)
            if(!(rest%sum)) dfs(i,num*nnum,rest/sum);
}

int main()
{
    pre();
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        tot=0;
        dfs(0,1,n);
        sort(ans+1,ans+tot+1);
        printf("%d\n",tot);
        for(int i=1;i<=tot;++i) printf("%d ",ans[i]);
        if(tot) printf("\n");
    }
}

 

bzoj千题计划297:bzoj3629: [JLOI2014]聪明的燕姿

标签:amp   prim   .com   body   分解   sig   +=   style   tar   

原文地址:https://www.cnblogs.com/TheRoadToTheGold/p/8611933.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!