【BZOJ3629】【JLOI2014】聪明的燕姿 dfs 素数筛

该博客探讨了如何使用DFS和素数筛方法解决BZOJ3629与JLOI2014比赛中的问题。博主提出将数分解质因数,根据质因数的个数计算约数和,以根号时间复杂度进行分解。通过暴力拆解每个数并递归搜索,实现问题的解决方案。
摘要由CSDN通过智能技术生成

链接:

#include <stdio.h>
int main()
{
    puts("转载请注明出处[vmurder]谢谢");
    puts("网址:blog.csdn.net/vmurder/article/details/44698555");
}

题解:

我们发现把一个数分解质因数以后然后可以根据每种质因数的个数算出这个数的约数和。所以我们可以暴力拆解每个数,根号时间复杂度分解。
就是枚举每种质数它用了多少,然后这个数除一下再往下一层深搜。

代码:

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 50100
using namespace std;
int prime[N],cnt;
bool vis[N];
void shake(int n)
{
    int i,j,k;
    for(i=2;i<=n;i++)
    {
        if(!vis[i])prime[++cnt]=i;
        for(j=1;j<=cnt&&i*prime[j]<=n;j++)
        {
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)break;
        }
    }
}
bool check(int x)
{
    bool flag=0;
    for(int i=1;(long long)prime[i]*prime[i]<=x;i++)
        if(x%prime[i]==0)return 0;
    return 1;
}
int ans[N],num;
void dfs(long long now,int remain,int t)
{
    if(remain==1){ans[++num]=now;return ;}
    if(remain==2)return ;
    if(remain-1>=prime[t]&&check(remain-1))ans[++num]=now*(remain-1);
    for(int i=t;(long long)prime[i]*prime[i]<=remain;i++)
    {
        long long sum=prime[i]+1,x=prime[i];
        while(sum<=remain)
        {
            if(remain%sum==0)
                dfs(now*x,remain/sum,i+1);
            x*=prime[i],sum+=x;
        }
    }
}
int n;
int main()
{
//  freopen("test.in","r",stdin);

    shake(50000);
    while(scanf("%d",&n)!=EOF)
    {
        num=0;
        dfs(1,n,1);
        sort(ans+1,ans+num+1),n=0;
        for(int i=1;i<=num;i++)if(ans[i]!=ans[i-1])ans[++n]=ans[i];
        printf("%d\n",n);
        if(n>=1)printf("%d",ans[1]);
        for(int i=2;i<=n;i++)printf(" %d",ans[i]);
        if(num)puts("");
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值