pat 素数

最新推荐文章于 2025-05-16 14:00:49 发布

honortech

最新推荐文章于 2025-05-16 14:00:49 发布

阅读量383

点赞数 10

文章标签：算法 c++ 数据结构

本文链接：https://blog.youkuaiyun.com/honortech/article/details/135299997

版权

两段代码分别展示了如何用C++编程语言计算给定整数的素数因子以及生成一定范围内的素数。A1059PrimeFactors函数用于找出所有小于100010的素数并找到输入数N的素数因子；B1013数素则用于生成M到N范围内的素数序列。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

A1059 Prime Factors

#include <iostream>
#include <vector>
#include <cmath>

using namespace std;
vector<int> primes;
int index = 0;

bool is_prime(int num)
{
    for(int i = 2; i <= sqrt(num * 1.0); ++i)
    {
        if(num % i == 0)
            return false;
    }
    return true;
}

void find_prime()
{
    for(int i = 2; i < 100010; ++i)
    {
        if(is_prime(i))
        {
            primes.push_back(i);
            index++;
        }
    }
}

struct factor
{
    int fac; // 素数因子及其个数
    int num;
};
vector<factor> factors;
int f_index = 0;
int main()
{
    find_prime();
    int N;
    cin >> N;
    if(N == 1) // 1的特殊情况
    {
        cout << "1=1";
        return 0;
    }
    cout << N << "=";
    
    for(int i = 0; i < index; ++i)
    {
        if(N % primes[i] == 0)
        {
            factor f;
            f.fac = primes[i];
            f.num = 0;
            while(N % primes[i] == 0) // 可能有多个相同素数因子
            {
                f.num++;
                N /= primes[i];
            }
            factors.push_back(f);
            f_index++;
        }
        if(N == 1) break; // 已找出全部因子，break
    }
    
    for(int i = 0; i < f_index; ++i)
    {
        cout << factors[i].fac;
        if(factors[i].num > 1)
            cout << "^" << factors[i].num;
        if(i < f_index - 1)
            cout << "*";
    }

    return 0;
}

B1013 数素数

#include <iostream>
#include <vector>
#include <cmath>

using namespace std;
const int maxn = 1e8 + 5;
int primes[maxn];
int index = 0;
int M, N;
bool is_prime(int n)
{
    for(int i = 2; i <= sqrt(n * 1.0); ++i)
    {
        if(n % i == 0)
            return false;
    }
    return true;
}

void find_prime()
{
    for(int i = 2; i < maxn; ++i)
    {
        if(index >= N)
            break;
        if(is_prime(i))
            primes[index++] = i;
    }
}

int main()
{
    cin >> M >> N;
    find_prime();
    int count = 0;
    for(int i = M - 1; i < N; ++i)
    {
        cout << primes[i];
        count++;
        if(count % 10 != 0 && i != N - 1)
            cout << " ";
        else 
            cout << endl;
    }

    return 0;
}