hdu 2138 How many prime numbers(数论:素数判定)

来源：互联网发布：mac粉尘壁纸编辑：程序博客网时间：2024/06/10 06:58

因为给出的数据是32位,所以可以直接对每个数暴力判定

当然也可以用大素数判定Miller Rabin算法

两份代码如下:

#include <cmath>#include <cstdio>#include <iostream>#include <algorithm>#define LL long longusing namespace std;bool judge(int n) {    int i;    int m = (int)sqrt(n+0.5);    for(i=2; i<=m; ++i)        if(n % i == 0)            return false;    return true;}int main(void) {    int n, i, cnt, x;    while(scanf("%d", &n) != EOF) {        cnt = 0;        for(i=0; i<n; ++i) {            scanf("%d", &x);            if(judge(x))                cnt++;        }        cout << cnt << endl;    }    return 0;}

#include<stdio.h>#include<stdlib.h>#include<cmath>bool witness(__int64 a,__int64 n) {    __int64 t, d, x;    d = 1;    int i=ceil(log(n-1.0)/log(2.0)) - 1;    for(; i>=0; i--) {        x = d;  d = (d*d)%n;        if(d==1 && x!=1 && x!=n-1) return true;        if( ((n-1) & (1<<i)) > 0)            d = (d*a)%n;    }    return d==1? false : true;}bool miller_rabin(__int64 n) {    if(n==2)    return true;    if(n==1 || ((n&1)==0))    return false;    for(int i=0;i<50;i++){        __int64 a=rand()*(n-2)/RAND_MAX +1;        if(witness(a, n))    return false;    }    return true;}int main() {    int n,cnt;    __int64 a;    while(scanf("%d",&n)!=EOF) {        cnt=0;        while(n--) {            scanf("%I64d",&a);            if(miller_rabin(a))                cnt++;        }        printf("%d\n",cnt);    }    return 0;}

0 0