0、参考链接:质数 | xiaoxiaoxy
注:一些我参考过的感觉好的文章链接都放在以上链接中了。
1、素数集:
素数集(A):
B = [(2*3)*n±1] n∈(1,2,3...)
C = [2,3]
D = A - C
E(2) = [ D(0)*D(0),D(0)*D(1),...,D(1)*D(1),D(1)*D(2),...,D(m)*D(n) ] (m,n∈N*,m≤n)
E(3) = [ D(0)*D(0)*D(0),...,D(x)*D(y)*(z) ] (x,y,z∈N*,x≤y≤z)
....
E(n) = [ ... ]
⇒ A = B∪C - E(2)∪E(3)∪...∪E(n)
稍微解释一下:在十进制下,所有的素数都包含于 [2,3,6n-1,6n+1] 这个集合之内,而这个集合之内的非素数又是由多个 素数 与 自身 或者 其他素数 多次相乘得来的。(个人总结的规律而已,也不会证明,所以这句话不一定对)
其他规律(未证明的规律,不一定对):
①:十进制下,除2、3、5外,素数以 1,3,7,9 结尾
②:素数模 6 的倍数得到的一定是素数 (1除外)
③:素数有无穷多个
2、素数的生成工具:
2.1 C代码
#include<stdio.h>
int main(){
long long int n,i,x,y,a[100000]; //arr[] 用于存储素数 每个元素64bit ,10万个大概占用 6M 的内存空间
FILE *fp=NULL;
int num=0,j,k=0;
printf("请输入一个大于 5 小于 2^63 的数:");
scanf("%ld",&n);
fp = fopen("./1.txt","w"); // 覆写当前目录下的 1.txt 文件,没有则先自动生成 1.txt
fprintf(fp,"2,3,"); // 先把不等于 6n-1,6n+1 的 2,3 写入
a[0] = 2;
a[1] = 3;
for(i = 1 ; i<=(n/6+1) ; i++){
x = 6*i-1;
y = 6*i+1;
if(num == 0)a[num] = x;
//---------------------- ↓x ↓----------------------------------------------------
for(j=0;a[j]*a[j]<=x;j++){ // 遍历素数集合就行,不需要从1开始
if( x % a[j] == 0 )break; //x能由素数相互乘出,说明 x 是素数
}
if(a[j]*a[j]>x){ //x不能由素数相互乘出,说明 x 是素数
num++;
a[num] = x; //素数存储至数组
fprintf(fp,"%ld,",x); //素数输出至文件
k++; //17个素数换行一次
if(k%17 == 0)fprintf(fp,"\n");
}
//---------------------- ↑x ↑----------------------------------------------------
//---------------------- ↓y ↓----------------------------------------------------
for(j=0;a[j]*a[j]<=y;j++){ // 遍历素数集合就行,不需要从1开始
if( y % a[j] == 0 )break; //y能由素数相互乘出,说明 y 是素数
}
if(a[j]*a[j]>y){ //y不能由素数相互乘出,说明 y 是素数
num++;
a[num] = y; //素数存储至数组
fprintf(fp,"%ld,",y); //素数输出至文件
k++; //17个素数换行一次
if(k%17 == 0)fprintf(fp,"\n");
}
//---------------------- ↑y ↑----------------------------------------------------
}
//for(i = 0;i<100;i++)printf("%d ",a[i]);
printf("输出成功,请在当前代码所在目录下查看1.txt文件\n");
printf("可能会多打印1-2项,例如输入范围上限:100,结果101和103也被打印出来了");
}
//方法解析:输入一个阈值,找出从 0 到阈值之内的所有素数
2.2 python代码
//方法解析:
// 1、输入一个阈值
// 2、生成 0-阈值 的 X = [ 2,3,6*n-1,6*n+1 ] 的列表
// 3、从列表中逐一剔除非素数(利用 X 中的非素数是由多个素数与自身或其他素数相乘得来这一性质)
//分析:
// 1、不好实现,即使实现了,似乎也不如上边那个C代码的思路
// 2、所以这里就不写了
2.3 bug&缺陷
本来想生成 0-2^30 范围内的所有素数的,结果因为数组满了导致输出失败,停在了第200000个素数左右,大概是 2^21 左右,输出不了 2^30 之内的素数,所以想要在不改变电脑配置的情况下生成更多的素数,就不得不拿时间来换空间了。这里的话就不能开数组记录出现过的素数了,而是要对每个新的6n-1与6n+1,用6m-1与6m+1来判断其是否能被整除。代码如下 ↓
2.4 优化1
#include<stdio.h>
int main(){
long long int n,i,x,y;
FILE *fp=NULL;
int j,k=0;
printf("请输入一个大于 5 小于 2^63 的数:");
scanf("%ld",&n);
fp = fopen("./1.txt","w"); // 覆写当前目录下的 1.txt 文件,没有则先自动生成 1.txt
fprintf(fp,"2,3,"); // 先把不等于 6n-1,6n+1 的 2,3 写入
for(i = 1 ; i<=(n/6+1) ; i++){
x = 6*i-1;
y = 6*i+1;
//---------------------- ↓x ↓----------------------------------------------------
for(j=6; (j-1)<=x ;j+=6){
if( (x%(j-1)==0) || (x%(j+1)==0) )break;
}
if(j>x){
fprintf(fp,"%ld,",x); //素数输出至文件
k++; //17个素数换行一次
if(k%17 == 0)fprintf(fp,"\n");
}
//---------------------- ↑x ↑----------------------------------------------------
//---------------------- ↓y ↓----------------------------------------------------
for(j=6; (j+1)<=y ;j+=6){
if( (y%(j-1)==0) || (y%(j+1)==0) )break;
}
if(j+2>=y){
fprintf(fp,"%ld,",y); //素数输出至文件
k++; //17个素数换行一次
if(k%17 == 0)fprintf(fp,"\n");
}
//---------------------- ↑y ↑----------------------------------------------------
}
printf("输出成功,请在当前代码所在目录下查看1.txt文件\n");
printf("可能会多打印1-2项,例如输入范围上限:100,结果101和103也被打印出来了");
}
这优化1确实很慢,2.1的代码都是秒算200w,但是空间不允许卡在了200w。而2 .4 优化1 的代码挂了将近一个小时才跑到600w( 2^22 - 2^23 之间) 。。。
2.5 优化2
学了一招,10分钟内可以打印出 0-2^30 以内的素数,总共5000w+个
#include<stdio.h>
int main(){
int n,i,x,y;
FILE *fp=NULL;
int *a = new int[1000000000]; //之前栈不够用,现在申请堆空间
int num=0,j,k=0;
printf("请输入一个大于 5 小于 2^31 的数:");
scanf("%d",&n);
fp = fopen("./1.txt","w"); // 覆写当前目录下的 1.txt 文件,没有则先自动生成 1.txt
fprintf(fp,"2,3,"); // 先把不等于 6n-1,6n+1 的 2,3 写入
a[0] = 2;
a[1] = 3;
for(i = 1 ; i<=(n/6+1) ; i++){
x = 6*i-1;
y = 6*i+1;
if(num == 0)a[num] = x;
//---------------------- ↓x ↓----------------------------------------------------
for(j=0;a[j]*a[j]<=x;j++){ // 遍历素数集合就行,不需要从1开始
if( x % a[j] == 0 )break; //x能由素数相互乘出,说明 x 是素数
}
if(a[j]*a[j]>x){ //x不能由素数相互乘出,说明 x 是素数
num++;
a[num] = x; //素数存储至数组
fprintf(fp,"%d,",x); //素数输出至文件
k++; //17个素数换行一次
if(k%17 == 0)fprintf(fp,"\n");
}
//---------------------- ↑x ↑----------------------------------------------------
//---------------------- ↓y ↓----------------------------------------------------
for(j=0;a[j]*a[j]<=y;j++){ // 遍历素数集合就行,不需要从1开始
if( y % a[j] == 0 )break; //x能由素数相互乘出,说明 x 是素数
}
if(a[j]*a[j]>y){ //x不能由素数相互乘出,说明 x 是素数
num++;
a[num] = y; //素数存储至数组
fprintf(fp,"%d,",y); //素数输出至文件
k++; //17个素数换行一次
if(k%17 == 0)fprintf(fp,"\n");
}
//---------------------- ↑y ↑----------------------------------------------------
}
//for(i = 0;i<100;i++)printf("%d ",a[i]);
printf("\nnum=%d\n",num);
printf("输出成功,请在当前代码所在目录下查看1.txt文件\n");
printf("可能会多打印1-2项,例如输入范围上限:100,结果101和103也被打印出来了");
}
3、10进制转6进制的实现
#include<stdio.h>
long long int six(long long int ten);
//--之后的代码只给出main函数,把后面给出的main函数替换进来就行--//
int main(){
printf("\n\t %ld \n",six(79));
}
//-- ↑ --//
long long int six(long long int ten){
long long int m,m1,n,x,num;
x = 6;
m = m1 = num = 0;
// 这里用 ten = 79 举例
while(ten!=0){ // ①ten = 79 ②ten = 78 ③ten = 72
n = ten%x; // ①n = 1 ②n = 78%36 = 6 ③n = 72%216 = 72
ten -= n; // ①ten = 78 ②ten = 78-6 = 72 ③ten = 72-72 = 0
n /= (x/6); // ①n = 1 ②n /= (36/6) = 1 ③n /= (216/6) = 2
while(m1--){ // ①未进入循环
n *= 10; // ②n = 10 ③n = 200
}
num += n; // ①num = 1 ②num = 11 ③num = 211
x *= 6; // ①x = 36 ②x = 216 ......
m++; // ①m = 1 ......
m1 = m; // ①m1 = 1 ......
}
return num; // 返回 211
}
注:从w中我们可以看出:6 似乎和素数有着某种联系,所以我们在 6 进制下对素数进行研究
3.1 简单查看素数对应的6进制
# define M 100000
int main(){
long long int a[M];
int i=0,j=0;
FILE *fp,*fp1;
fp = fopen("./1.txt","r"); //读取当前目录下的 1.txt 文件,存储着素数,由之前的代码生成
fp1 = fopen("./2.txt","w");
while(fscanf(fp,"%ld,",&a[i])!=EOF){ //读入1.txt中的素数集
fprintf(fp1,"%ld,",six(a[i]));
i++;
j++;
if(j%17==0)fprintf(fp1,"\n");
}
printf("输出成功,请查看当前代码目录下的 2.txt 文件");
}
//结果:200以内,除了末尾全是1、5外,似乎没什么规律可循
// 2,3,5,11,15,21,25,31,35,45,51,101,105,111,115,125,135,
// 141,151,155,201,211,215,225,241,245,251,255,301,305,331,335,345,351,
// 405,411,421,431,435,445,455,501,515,521,525,531,
// 5*5乘法表
// 1*1=1
// 1*2=2 2*2=4
// 1*3=3 2*3=10 3*3=13
// 1*4=4 2*4=12 3*4=20 4*4=24
// 1*5=5 2*5=14 3*5=23 4*5=32 5*5=41
4、素数的数量
代码不贴了,这是从 0、参考链接里搬过来的图,这里只是想表达一下随着整数集的指数倍增加,素数的数量也在以指数倍增长,只不过整数集的增长指数为2,素数的数量增长的指数逐渐趋近于1 (个人根据2^20以内的数字进行的猜测而已),总之就是素数有无穷多个。
5、尝试从不同角度看素数
5.1 素数包含在 [2,3,6n±1] 中,所以看看 6n±1 中的 n 与素数的关系。
尝试寻找1000以内的素数的规律,列表a可由 2.5优化2 中的代码生成后复制进去
##!/usr/bin/env python
# -*- coding: utf-8 -*-
import math
import string
a = [2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,
73,79,83,89,97,101,103,107,109,113,127,131,137,139,149,151,157,163,167,173,
179,181,191,193,197,199,211,223,227,229,233,239,241,251,257,263,269,271,277,281,
283,293,307,311,313,317,331,337,347,349,353,359,367,373,379,383,389,397,401,409,
419,421,431,433,439,443,449,457,461,463,467,479,487,491,499,503,509,521,523,541,
547,557,563,569,571,577,587,593,599,601,607,613,617,619,631,641,643,647,653,659,
661,673,677,683,691,701,709,719,727,733,739,743,751,757,761,769,773,787,797,809,
811,821,823,827,829,839,853,857,859,863,877,881,883,887,907,911,919,929,937,941,
947,953,967,971,977,983,991,997,1009,1013,1019,1021,1031,1033,1039,1049,1051,1061]
def s_m1(x): #分离出6n-1, m : minus
y = []
y.append(2)
y.append(3)
for i in x:
if ( (i+1)%6 == 0 ):
y.append(i)
return y
def s_a1(x): #分离出6n+1, a : add
y = []
y.append(2)
y.append(3)
for i in x:
if ( (i-1)%6 == 0 ):
y.append(i)
return y
def s_d6(x): #输出6n±1中的n, d : div
y = []
for i in x:
y.append((i+1)//6)
return y
def s_2(x): #输出平方
y = []
y.append(4)
y.append(9)
for i in x:
y.append(i*i)
return y
def s_3(x): #输出立方
y = []
y.append(8)
y.append(27)
for i in x:
y.append(i*i*i)
return y
def s_l(x): #输出缺少的整数
y = []
for i in range(0,x[-1]):
if i not in x:
y.append(i)
return y
b = s_m1(a) #输出[2,3,6n-1]
c = s_a1(a) #输出[2,3,6n+1]
x1 = s_d6(b)
x2 = s_d6(c)
x3 = s_d6(a)
'''
y1 = s_2(x1)
y2 = s_2(x2)
y3 = s_2(x3)
z1 = s_3(x1)
z2 = s_3(x2)
z3 = s_3(x3)
'''
a1 = s_l(x1) #输出 6n-1 中缺少的n
a2 = s_l(x2) #输出 6n+1 中缺少的n
a3 = s_l(x3) #输出 6n±1 中缺少的n
'''
y1 = s_2(a1)
y2 = s_2(a2)
y3 = s_2(a3)
z1 = s_3(a1)
z2 = s_3(a2)
z3 = s_3(a3)
'''
print ("6n-1:",b,'\n',"6n+1:",c,'\n',"6n±1:",a,'\n\n')
print ("x1:",x1,'\n',"x2:",x2,'\n',"x3:",x3,'\n\n')
#print ("y1:",y1,'\n',"y2:",y2,'\n',"y3:",y3,'\n\n')
#print ("z1:",z1,'\n',"z2:",z2,'\n',"z3:",z3,'\n\n')
print ("a1:",a1,'\n',"a2:",a2,'\n',"a3:",a3,'\n\n')
'''
x1: 稍微解释一下, x1 是由 [2,3,6n-1] 中的每个元素都 加 1 后 除以 6 得来的
[ 0, 0, 1, 2, 3, 4, 5, 7, 8, 9, 10, 12, 14, 15, 17, 18, 19, 22, 23, 25, 28, 29,
30, 32, 33, 38, 39, 40, 42, 43, 44, 45, 47, 49, 52, 53, 58, 59, 60, 64, 65,
67, 70, 72, 74, 75, 77, 78, 80, 82, 84, 85, 87, 93, 94, 95, 98, 99, 100, 103,
107, 108, 109, 110, 113, 114, 117, 120, 124, 127, 129, 133, 135, 137, 138,
140, 143, 144, 147, 148, 152, 155, 157, 158, 159, 162, 163, 164, 169, 170,
172, 175, 177]
x2:
[ 0, 0, 1, 2, 3, 5, 6, 7, 10, 11, 12, 13, 16, 17, 18, 21, 23, 25, 26, 27, 30,
32, 33, 35, 37, 38, 40, 45, 46, 47, 51, 52, 55, 56, 58, 61, 62, 63, 66, 68,
70, 72, 73, 76, 77, 81, 83, 87, 90, 91, 95, 96, 100, 101, 102, 103, 105, 107,
110, 112, 115, 118, 121, 122, 123, 125, 126, 128, 131, 135, 137, 138, 142,
143, 146, 147, 151, 153, 156, 161, 165, 166, 168, 170, 172, 173, 175 ]
a1: 6=2*3 11=... , a1 是记录 0 - x[-1] (也就是列表末尾元素) 中未出现的整数
[ 6, 11, 13, 16, 20, 21, 24, 26, 27, 31, 34, 35, 36, 37, 41, 46, 48, 50, 51,
54, 55, 56, 57, 61, 62, 63, 66, 68, 69, 71, 73, 76, 79, 81, 83, 86, 88, 89,
90, 91, 92, 96, 97, 101, 102, 104, 105, 106, 111, 112, 115, 116, 118, 119,
121, 122, 123, 125, 126, 128, 130, 131, 132, 134, 136, 139, 141, 142, 145,
146, 149, 150, 151, 153, 154, 156, 160, 161, 165, 166, 167, 168, 171, 173, 174, 176 ]
a2: 4=2*2 8=2*2*2=2*4 9=3*3 14=2*7 15=3*5 19=...
[ 4, 8, 9, 14, 15, 19, 20, 22, 24, 28, 29, 31, 34, 36, 39, 41, 42, 43, 44,
48, 49, 50, 53, 54, 57, 59, 60, 64, 65, 67, 69, 71, 74, 75, 78, 79, 80, 82,
84, 85, 86, 88, 89, 92, 93, 94, 97, 98, 99, 104, 106, 108, 109, 111, 113,
114, 116, 117, 119, 120, 124, 127, 129, 130, 132, 133, 134, 136, 139, 140,
141, 144, 145, 148, 149, 150, 152, 154, 155, 157, 158, 159, 160, 162, 163,
164, 167, 169, 171, 174 ]
a3: 20=4*5 24=4*6 31=... 似乎也没啥规律
[ 20, 24, 31, 34, 36, 41, 48, 50, 54, 57, 69, 71, 79, 86, 88, 89, 92, 97,
104, 106, 111, 116, 119, 130, 132, 134, 136, 139, 141, 145, 149, 150, 154,
160, 167, 171, 174, 176 ]
'''