Eratosthenes筛法求素数的方法技巧
|
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(); } |
本文地址：http://www.45fan.com/dnjc/70270.html