Eratosthenes筛法求素数

王朝other·作者佚名  2006-01-09
窄屏简体版  字體: |||超大  

Eratosthenes筛法求素数

作者:未知 来源:未知

我用了c的方法和在tc下混合编程的方法,如果用win-tc就可以直接编译此程序,如果是用tc2.0就必须用tcc编辑器自己编译。或许很多人用tc嵌入汇编的时候会觉得编译不通过,编译的时候需要将tc所在目录里的Lib目录里的所有文件复制到tc当前目录下,然后用在当前目录用 tcc -Iinclude -Llib 文件名字既可。比如以下程序在文件su.c中,那么你只需要在tc目录里用tcc -Iinclude -Llib su的命令既可编译此程序。当然有win-tc更为方便,直接编译就可以了。

/*

Description:

Eratosthenes筛法

创建一个数组,以下面的方式将某些位置以1标记:从位置2开始,将所有2的倍数的位置

标记为1,然后对2之后的素数3进行同样的操作这样就可以找到3之后未被标记的5,

然后再将所有5的倍数记1,如此重复便可以找到数组中的所有素数。

Author:跳蚤侦探

Date: 22/08/2004

Copyright:www.8623.com

*/

#include <dos.h>

#define N 100

#define HALF N/2

void c_method()/*c语言的方法*/

{

int f[N+1]={0},i,j,k;

printf("c method:\n");

f[0]=f[1]=1;

i=2;

while(i<=N/2)

{

j=2;

while(1)

{

f[j*i]=1;

j++;

if((j*i)>=N)

break;

}

for(k=i+1;k<N;k++)

if(f[k]!=1)

{

i=k;

break;

}

}

for(i=0;i<N;i++)

if(f[i]==0)

printf("%d,",i);

printf("\n");

}

void asm_method()/*混合汇编的方法*/

{

unsigned f[N+1]={0},i;

f[0]=f[1]=1;

printf("asm method:\n");

asm mov di,2

asm clc

PROCESS:

asm mov cx,2

PROSON:

asm mov ax,di

asm mov bx,ax

asm mul cx

asm cmp ax,N

asm jae PROSONEND

f[_AX]=1;

asm mov ax,bx

asm inc cx

asm jmp PROSON

PROSONEND:

asm mov si,di

asm inc si

PROSON1:

if(f[_SI]==0)

{

asm je PROSON1END

}

asm inc si

asm jmp PROSON1

PROSON1END:

asm mov di,si

asm cmp di,HALF

asm jae PROEND

asm jmp PROCESS

PROEND:

for(i=0;i<N;i++)

if(f[i]==0)

printf("%d,",i);

return;

}

void main()

{

clrscr();

c_method();

asm_method();

getch();

}

 
 
 
免责声明:本文为网络用户发布,其观点仅代表作者个人观点,与本站无关,本站仅提供信息存储服务。文中陈述内容未经本站证实,其真实性、完整性、及时性本站不作任何保证或承诺,请读者仅作参考,并请自行核实相关内容。
 
 
© 2005- 王朝網路 版權所有 導航