首页 > 其他 > 详细

hdu 2138 How many prime numbers_Miller-Rabin算法模板

时间:2014-04-24 13:05:45      阅读:524      评论:0      收藏:0      [点我收藏+]

题目链接

题意:给你n个数,判读有几个素数

思路:求一个很大的数是不是素数,要用Miller-Rabin算法

#include <iostream>
#include<cstdio>
#include<cmath>
using namespace std;
#define int64 long long
int64 qpow(int64 a,int64 b,int64 M){
	int64 ans=1;
	while(b){
		if(b&1)ans*=a,ans%=M;
		a*=a;
		a%=M;
		b>>=1;
	}
	return ans;
}
bool MillerRabinTest(int64 x,int64 n){
	int64 y=n-1;
	while(!(y&1))y>>=1;
	x=qpow(x,y,n);
	while(y<n-1&&x!=1&&x!=n-1)
		x=(x*x)%n,y<<=(int64)1;
	return x==n-1||y&1==1;
}
bool isprime32(int64 n){
	if(n==2||n==7||n==61)return 1;
	if(n==1||(n&1)==0)return 0;
	return MillerRabinTest(2,n)&&MillerRabinTest(7,n)&&MillerRabinTest(61,n);
}
int main(int argc, char** argv) {
	int n,sum,i;
	int64 a;
	while(scanf("%d",&n)!=EOF){
		sum=0;
		for(i=0;i<n;i++){
			scanf("%lld",&a);
			if(isprime32(a))
				sum++;
		}
		printf("%d\n",sum);
	}
	return 0;
}


hdu 2138 How many prime numbers_Miller-Rabin算法模板,布布扣,bubuko.com

hdu 2138 How many prime numbers_Miller-Rabin算法模板

原文:http://blog.csdn.net/neng18/article/details/24365579

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!