素数判定，miller rabin

最新推荐文章于 2022-12-17 17:06:44 发布

Thomas_ZQQ@Runespoor

最新推荐文章于 2022-12-17 17:06:44 发布

阅读量128

点赞数

CC 4.0 BY-SA版权

分类专栏：数论 ACM 文章标签：数学

本文链接：https://blog.youkuaiyun.com/weixin_42484877/article/details/80880889

ACM 同时被 2 个专栏收录

95 篇文章

订阅专栏

数论

12 篇文章

订阅专栏

本文提供了一种高效的素性检验算法实现，用于解决HDU2138竞赛题目。通过快速幂取模计算和Miller-Rabin素性测试，该算法能够快速判断给定数值是否为素数，并支持多组数据输入。

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

点击打开链接

裸题 hdu 2138 注意多组数据

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define maxn 10020

typedef long long ll;
int n,num,p;

ll power(ll x,int y){
	ll res = 1;
	while ( y ){
		if ( y & 1 ) res = res * x % p;
		x = x * x % p;
		y >>= 1;
	}
	return res;
}
bool cal(ll x,ll d,ll r){
	x = power(x,d);
	if ( x == 1 || x == p - 1 ) return 1;
	for (int i = 1 ; i <= r ; i++){
	   	x = x * x % p;
		if ( x == p - 1 ) return 1;
	}
	return 0;
}
bool judge(){
	int a = p - 1,cnt = 0;
	if ( p == 2 ) return 1;
	while ( (a & 1) == 0 ) a >>= 1 , cnt++;
	for (int i = 1 ; i <= 10 ; i++){
		int cur = rand() % (p - 1) + 1;
		if ( !cal(cur,a,cnt) ) return 0;	
	}
	return 1;
}
int main(){
	while (~scanf("%d",&n)){
		num = 0;
		while ( n-- ){
			scanf("%d",&p);
			if ( judge() ) num++;
		}
		printf("%d\n",num);
	}
	return 0;
}