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();
}