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

　　#define N 100

　　#define HALF N/2

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

　　{

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

　　printf("c method:

");

　　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

　　 if(f[k]!=1)

　　 {

　　 i=k;

　　 break;

　　 }

　　}

　　 for(i=0;i

　　 if(f[i]==0)

　　 printf("%d,",i);

　　 printf("

");

　　}

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

　　{

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

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

　　 printf("asm method:

");

　　 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

　　 if(f[i]==0)

　　 printf("%d,",i);

　　 return;

　　}

void main()

　　{

　　 clrscr();

　　 c_method();

　　 asm_method();

　　 getch();

　　}