聪明的燕姿

最新推荐文章于 2024-03-13 22:19:44 发布

胡耀文

最新推荐文章于 2024-03-13 22:19:44 发布

阅读量113

点赞数

文章标签：算法

本文链接：https://blog.csdn.net/supreme567/article/details/115710680

版权

该博客主要探讨了两个C++程序，用于寻找给定整数的质数因子和生成所有小于给定数的素数。第一个程序通过枚举和质数筛法找出质数因子，第二个程序使用深度优先搜索和前缀积方法生成所有可能的素数组合。内容涉及到质数检测、算法实现和代码优化。

摘要由CSDN通过智能技术生成

在这里插入图片描述

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
int s;
int tot;
int a[10000001];
int ysh(int x){
    int ans=1+x;
    
    for(int i=2;i<x;i++){
        if(x%i==0)ans+=i;
        
    }
    return ans;
}
int main(){
    while(cin>>s){
        tot=0;
        memset(a,0,sizeof a);
        for(int i=2;i<s;i++){
        if(s==ysh(i))a[tot++]=i;
        }
        
        cout<<tot<<endl;
        for(int i=0;i<tot;i++)printf("%d ",a[i]);
        cout<<endl;
    }
    
    
}

#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;

const int N = 50000;

int primes[N], cnt;
bool st[N];

int ans[N], len;

void get_primes(int n)
{
    for (int i = 2; i <= n; i ++ )
    {
        if (!st[i]) primes[cnt ++ ] = i;
        for (int j = 0; primes[j] * i <= n; j ++ )
        {
            st[primes[j] * i] = true;
            if (i % primes[j] == 0) break;
        }
    }
}

bool is_prime(int x)
{
    if (x < N) return !st[x];
    for (int i = 0; primes[i] <= x / primes[i]; i ++ )
        if (x % primes[i] == 0)
            return false;
    return true;
}

void dfs(int last, int prod, int s)
{
    if (s == 1)
    {
        ans[len ++ ] = prod;
        return;
    }

    if (s - 1 > (last < 0 ? 1 : primes[last]) && is_prime(s - 1))
        ans[len ++ ] = prod * (s - 1);

    for (int i = last + 1; primes[i] <= s / primes[i]; i ++ )
    {
        int p = primes[i];
        for (int j = 1 + p, t = p; j <= s; t *= p, j += t)
            if (s % j == 0)
                dfs(i, prod * t, s / j);
    }
}

int main()
{
    get_primes(N - 1);

    int s;
    while (cin >> s)
    {
        len = 0;
        dfs(-1, 1, s);

        cout << len << endl;
        if (len)
        {
            sort(ans, ans + len);
            for (int i = 0; i < len; i ++ ) cout << ans[i] << ' ';
            cout << endl;
        }
    }

    return 0;
}

作者：yxc
链接：https://www.acwing.com/activity/content/code/content/197156/
来源：AcWing
著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。