HDU2138 & 米勒拉宾模板
Posted YCuangWhen
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了HDU2138 & 米勒拉宾模板相关的知识,希望对你有一定的参考价值。
题意:
给出n个数,判断它是不是素数.
SOL:
米勒拉宾裸题,思想方法略懂,并不能完全理解,所以实现只能靠背模板....
好在不是很长...
Code:
/*========================================================================== # Last modified: 2016-03-21 10:09 # Filename: miller-rabin.cpp # Description: ==========================================================================*/ #define me AcrossTheSky #include <cstdio> #include <cmath> #include <ctime> #include <string> #include <cstring> #include <cstdlib> #include <iostream> #include <algorithm> #include <set> #include <map> #include <stack> #include <queue> #include <vector> #define lowbit(x) (x)&(-x) #define FOR(i,a,b) for((i)=(a);(i)<=(b);(i)++) #define FORP(i,a,b) for(int i=(a);i<=(b);i++) #define FORM(i,a,b) for(int i=(a);i>=(b);i--) #define ls(a,b) (((a)+(b)) << 1) #define rs(a,b) (((a)+(b)) >> 1) #define getlc(a) ch[(a)][0] #define getrc(a) ch[(a)][1] #define maxn 100000 #define maxm 100000 #define pi 3.1415926535898 #define _e 2.718281828459 #define INF 1070000000 using namespace std; typedef long long ll; typedef unsigned long long ull; template<class T> inline void read(T& num) { bool start=false,neg=false; char c; num=0; while((c=getchar())!=EOF) { if(c==‘-‘) start=neg=true; else if(c>=‘0‘ && c<=‘9‘) { start=true; num=num*10+c-‘0‘; } else if(start) break; } if(neg) num=-num; } /*==================split line==================*/ int c[7]={2,3,5,7,11,13}; ll pow(ll a,ll k,ll mod){ ll s=1; for (;k;a=(a*a)%mod,k>>=1) if (k&1) s*=a,s%=mod; return s; } bool check(int x,int a){ if (!x%a) return false; int r=0,s=x-1; while (!(s&1)) s>>=1,r++; ll k=pow(a,s,x); if (k==1) return true; for (int j=0;j<r;j++,k=k*k%x) if (k==x-1) return true; return false; } bool miller_rabin(ll x){ FORP(i,0,5) if (x==c[i]) return true; FORP(i,0,5){ if (!check(x,c[i])) return false; } return true; } int main(){ int n; while (scanf("%d",&n)!=EOF){ int ans=0; FORP(i,1,n) { ll x; read(x); if (miller_rabin(x)) ans++; } printf("%d\n",ans); } }
以上是关于HDU2138 & 米勒拉宾模板的主要内容,如果未能解决你的问题,请参考以下文章