HDU2138 How many prime numbers

最新推荐文章于 2023-08-03 10:52:13 发布

原创最新推荐文章于 2023-08-03 10:52:13 发布 · 1.2k 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#算法

ACM_数学专项系列专栏收录该内容

23 篇文章

订阅专栏

本文介绍了一种基于Miller-Rabin素性测试的改进算法实现，通过特定的测试基数判断一个数是否为素数，并展示了完整的C++代码实现。该算法适用于需要高效素数判定的场景。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

PS: 素数判定问题。

#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;

long long pow_mod(long long a, long long i, long long n) {
    if(i==0) return 1%n;
    long long tmp = pow_mod(a, i>>1, n); // int tmp. WA.
    tmp = tmp*tmp%n;
    if(i&1) tmp = (long long)tmp*a%n;
    return tmp;
}

bool test(int n, int a, int d) {
    if(n==2) return true;
    if(n==a) return true;
    if((n&1)==0) return false;
    while(!(d&1)) d = d>>1;
    int t = pow_mod(a, d, n);
    while((d!=n-1) && (t!=1) && (t!=n-1)) {
        t = (long long)t*t%n;
        d = d << 1;
    }
    return (t==n-1 || (d&1)==1);
}

bool isPrime(int n) {
    if(n<2) return false;
   // int a[] = {2, 3 , 5, 7, 11, 13, 17};
    int a[] = {2, 3, 5};
    for(int i = 0; i <= 2; i++) {
        if(!test(n, a[i], n-1)) return false;
    }
    return true;
}
int main()
{
    int n, t, cnt;
    while(scanf("%d", &n)!=EOF) {
        int *p = new int[n+5];
        for(int i = 0; i < n; i++) {
            scanf("%d", &t);
            p[i] = t;
        }
        cnt = 0;
        for(int i = 0; i < n; i++) {
            if(isPrime(p[i])) cnt++;
        }
        printf("%d\n", cnt);
    }
    return 0;
}