zl_5.c

Eratosthenes筛法求素数

 
 

 



/*

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