首页 / 算法 / HDU 2138 How many prime numbers （判素数，米勒拉宾算法）

HDU 2138 How many prime numbers （判素数，米勒拉宾算法）

内容导读

互联网集市收集整理的这篇技术教程文章主要介绍了HDU 2138 How many prime numbers （判素数，米勒拉宾算法），小编现在分享给大家，供广大互联网技能从业者学习和参考。文章包含1192字，纯文字阅读大概需要2分钟。

内容图文

题意：给定一个数，判断是不是素数。

析：由于数太多，并且太大了，所以以前的方法都不适合，要用米勒拉宾算法。

代码如下：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <vector>
#include <cstring>
#include <map>
#include <cctype>

using namespace std;
typedef long long LL;
const int maxn = 1000 + 5;

LL qpow(int a, int b, int r){
    LL ans = 1;
    LL k = a % r;
    while(b){
        if(b & 1)  ans = (ans * k) % r;
        k = (k * k) % r;
        b >>= 1;
    }
    return ans;
}

bool miller_rabbin(int n, int a){
    int r = 0,s = n-1;
    if(!(n % a))  return false;
    while(!(s & 1)){  s >>= 1;  ++r; }

    LL k = qpow(a, s, n);
    if(1 == k)  return true;
    for(int j = 0; j < r; ++j, k = k * k % n)
        if(k == n-1)  return true;
    return false;
}

bool is_prime(int n){
    int tab[] = {2, 3, 5, 7};
    for(int i = 0; i < 4; ++i){
        if(n == tab[i])  return true;
        if(!miller_rabbin(n, tab[i])) return false;
    }
    return true;
}

int main(){
//    freopen("in.txt", "r", stdin);
    int n, x;
    while(~scanf("%d", &n)){
        int cnt = 0;
        for(int i = 0; i < n; ++i){
            scanf("%d", &x);
            if(is_prime(x))  ++cnt;
        }
        printf("%d\n", cnt);
    }
}

 

原文：http://www.cnblogs.com/dwtfukgv/p/5557082.html

内容总结

以上是互联网集市为您收集整理的HDU 2138 How many prime numbers （判素数，米勒拉宾算法）全部内容，希望文章能够帮你解决HDU 2138 How many prime numbers （判素数，米勒拉宾算法）所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错，欢迎将互联网集市网站推荐给程序员好友。

内容备注

版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 gblab@vip.qq.com 举报，一经查实，本站将立刻删除。

内容手机端

扫描二维码推送至手机访问。

本文链接：https://qyyshop.com/info/1288578.html

来源：【匿名】