一个简单的素数产生函数

使用最常用的方法,判断从1到sqrt(n)的数,看是否i%这个数为0,为0,则不是素数...

 

ContractedBlock.gif ExpandedBlockStart.gif Code
#include <stdlib.h>
#include 
<iostream.h>
#include 
<math.h>

inline 
long sqrt(long i)
{
    
long r,rnew=1,rold=r;

    
do
    {
        rold
=r;
        r
=rnew;
        rnew 
= (r+(i/r));
        rnew 
>>= 1;
    }
    
while(rold != rnew);
    
return rnew;
}

inline 
int isprime(long i)
{
    
long si,j;

    si 
= sqrt(i);

    
for (j=2; (j<=si); j++)
    {
        
if (i%== 0)
            
return 0;
    }

    
return 1;

}

void main(int argc, char *argv[])
{
    
long i,i1,i2,f=0;

    
if (argc < 2)
    {
        cout 
<< endl << endl;
        cout 
<< "Prime number finder: Hugo Elias 1998" << endl << "http://freespace.virgin.net/hugo.elias" << endl << endl;
        cout 
<< "Usage:" << endl << "  primes a b > primes.txt" << endl << endl;
        cout 
<< "will find all primes between a and b" << endl << "and will write the results to the file PRIMES.TXT" << endl;
        
return;
    }

    i1 
= atol(argv[1]);
    i2 
= atol(argv[2]);

    
for (i=i1; i<i2; i++)
        
if (isprime(i))
        {
            f
++;
            
if (f==16)
            {
                cout 
<< endl;
                f
=0;
            }
            cout 
<< i << " ";
        }

}

 

posted on 2009-11-22 15:00 迈克老狼 阅读( ...) 评论( ...) 编辑 收藏

转载于:https://www.cnblogs.com/mikewolf2009/archive/2009/11/22/1608068.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值