hdu2138 How many prime numbers(Miller_Rabbin随机素数测试)

给定几个数,问其中素数的个数。

由于数量很大,打表肯定过不了,所以必须用Miller_Rabbin随机素数测试的方法,下面是我的C++ Miller_Rabbin模板


#include <iostream>
#include <cstdio>
#include <algorithm>

using namespace std;

typedef long long ll;

ll exp_mod(ll a,ll b,ll mod)
{
	ll res = 1;
	while(b){
		if(b & 1) res = (res * a) % mod;
		a = (a * a) % mod;
		b >>= 1;
	}
	return res;
}

bool Miller_Rabbin(ll n)
{
	ll a,s = (n-1) >> 2,i = 1, t = n-1;
	s = min(s,(ll)50);
	do{
		a = rand() % (t) + 1;
		if(exp_mod(a,t,n) != 1) return false;
	}while(i++ < s);
	return true;
}

int main()
{
	int n,ans;
	ll num;
	while(scanf("%d",&n)!=EOF)
	{
		ans = 0;
		for(int i = 0;i < n;i++)
		{
			scanf("%I64d",&num);
			if(Miller_Rabbin(num))
				ans++;
		}
		printf("%d\n",ans);
	}
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值