Project_Euler-07 题解

Project_Euler-07 题解

题目

思路

一个线性筛解决问题,当然也可以用埃式筛或者标准的暴力破解,这里选用最优秀的方式,顺便复习一下线性筛的内容:

c 复制代码
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<string.h>
#include<time.h>

// 返回一个指向存放n个数之前所有素数的数组
int *__PrimeGet(int n) {
	// arr以下标为数字,作为标记数组,0表示质数,1表示合数
    int *arr = (int *)malloc(sizeof(int) * n);
    // prime存放已经被筛出来的素数
    int *prime = (int *)malloc(sizeof(int) * n);
    // 全部初始化为0
    memset(arr, 0, sizeof(int) * n);
    memset(prime, 0, sizeof(int) * n);
    
    // pos表示prime中最后一个素数后一个地址,从1开始存放
    int pos = 1;
    //遍历arr中所有元素,都判断一遍,i表示正在判断的数
    for (int i = 2; i <= n; i++) {
    	// 如果为素数,存放到prime数组中。
        if (!arr[i]) prime[pos++] = i;
        // 遍历已经筛出的素数
        for (int j = 1; j < pos; j++) {
        	// 如果当前判断值i与当前遍历到的素数乘积大于要求范围n,停止判断
            if (i * prime[j] > n) break;
            // 将当前遍历到的素数的i倍的数标记为1
            arr[i * prime[j]] = 1;
            // 当前判断值是第一个可以除掉j值对应素数的,说明prime[j]是该数的最小质因数
            if (i % prime[j] == 0) break;
        }
    }
    free(arr);
    return prime;
}


// 这个程序是用来输出n之前所有素数的程序,如果你想要知道第10001个质数是多少,
// 你可以将结果重定向到一个带行号的文本文件中,查看第10001行的值,就是答案。
int main(){
    int n;
    scanf("%d", &n);
    int *prime = __PrimeGet(n);
    for (int i = 1; prime[i]; i++) printf("%d\n", prime[i]);
    free(prime);
    return 0;
}
```
相关推荐
ajassi20005 小时前
开源 C++ QT Widget 开发(七)线程--多线程及通讯
linux·c++·qt·开源
bdgtd881786 小时前
动态修补C扩展模块的函数指针有哪些风险?安全的修补方案是什么?
c语言·开发语言·安全
mit6.8246 小时前
8.27 网格memo
c++·算法
jeffery8926 小时前
4056:【GESP2403八级】接竹竿
数据结构·c++·算法
Forward♞6 小时前
Qt——界面美化 QSS
开发语言·c++·qt
Ghost-Face7 小时前
图论基础
算法
默归7 小时前
分治法——二分答案
python·算法
快乐的划水a8 小时前
解释器模式及优化
c++·设计模式·解释器模式
一枝小雨8 小时前
【数据结构】排序算法全解析
数据结构·算法·排序算法
略知java的景初8 小时前
深入解析十大经典排序算法原理与实现
数据结构·算法·排序算法