筛选法求素数

 
 C++ Code 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
 
筛法素数产生器

语法:result=prime( int  a[], int  n);

参数:

a[]:
用于返回素数的数组

n:
产生n以内的素数,按升序放入a[]中

返回值:
n以内素数的个数

注意:



其中W[],B[]已知,W[i]>0且W[i]与W[j]互质, 求a

源程序:



int  prime( int  a[], int  n)
{
    
int  i,j,k,x,num,*b;
    n++;
    n/=
2 ;
    b=(
int  *)malloc( sizeof ( int )*(n+ 1 )* 2 );
    a[
0 ]= 2 ;
    a[
1 ]= 3 ;
    num=
2 ;
    
for (i= 1 ; i<= 2 *n; i++)
        b[i]=
0 ;
    
for (i= 3 ; i<=n; i+= 3 )
        
for (j= 0 ; j< 2 ; j++)
        {
            x=
2 *(i+j)- 1 ;
            
while (b[x]== 0 )
            {
                a[num++]=x;
                
for (k=x; k<= 2 *n; k+=x)
                    b[k]=
1 ;
            }
        }
    
return  num;
}

转载于:https://www.cnblogs.com/Skyxj/p/3221881.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值