码迷,mamicode.com
首页 > 其他好文 > 详细

HDU 2138

时间:2014-09-06 10:53:13      阅读:165      评论:0      收藏:0      [点我收藏+]

标签:blog   os   io   for   div   sp   log   amp   c   

这题用MILLER测试应该是不可避免的。

#include <iostream>
#include <cstdio>
#include <stdlib.h>
#include <time.h>
#define LL __int64
using namespace std;

LL random(LL n){
	return (LL)((double)rand()/RAND_MAX*n+0.5);
}

LL quick(LL a,LL k,LL m){
	LL ans=1;
	while(k){
		if(k&1){
			ans=ans*a%m;
		}
		k>>=1;
		a=a*a%m;
	}
	return ans;
}

bool judgep(LL p){
	for(int i=1;i<=15;i++){
		LL a=random(p-2)+1;
		if(quick(a,p,p)!=a)
		return false;
	}
	return true;
}

int main(){
	int n;
	LL p;
	srand(time(0));
	while(scanf("%d",&n)!=EOF){
		int flag=0;
		for(int i=1;i<=n;i++){
			scanf("%I64d",&p);
			if(judgep(p)){
				flag++;
			}
		}
		printf("%d\n",flag);
	}
	return 0;
}

  

HDU 2138

标签:blog   os   io   for   div   sp   log   amp   c   

原文地址:http://www.cnblogs.com/jie-dcai/p/3959160.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!