\(\color{#0066ff}{ 題目描述 }\)

判定輸入的數是不是質數。

\(\color{#0066ff}{輸入格式}\)

若干行，一行一個數 \(x\)。

行數不超過 \(10^5\)。

\(\color{#0066ff}{輸出格式}\)

對于輸入的每一行，如果 \(x\) 是質數輸出一行 \(Y\)，否則輸出一行 \(N\)。

\(\color{#0066ff}{輸入樣例}\)

1 2 6 9 666623333

\(\color{#0066ff}{輸出樣例}\)

N Y N N Y

\(\color{#0066ff}{數據范圍與提示}\)

\(1≤x≤10^{18}\)。

\(\color{#0066ff}{ 題解 }\)

就是個Miller Rabbin的板子題

詳解

#include<bits/stdc++.h> #define LL long long LL in() {char ch; LL x = 0, f = 1;while(!isdigit(ch = getchar()))(ch == '-') && (f = -f);for(x = ch ^ 48; isdigit(ch = getchar()); x = (x << 1) + (x << 3) + (ch ^ 48));return x * f; } int prime[] = {2, 3, 5, 7, 11, 61, 24251}; LL msc(LL x, LL y, LL mod) {x %= mod;y %= mod;LL c = (long double)x / mod * y;LL d = x * y - c * mod;return ((d % mod) + mod) % mod; } LL ksm(LL x, LL y, LL mod) {LL re = 1LL;while(y) {if(y & 1) re = msc(re, x, mod);x = msc(x, x, mod); y >>= 1;}return (re + mod) % mod; } bool judge(LL a, LL p) {LL s = p - 1;while(!(s & 1)) s >>= 1;LL k = ksm(a, s, p);while (s != p - 1 && k != 1 && k != p - 1) k = msc(k, k, p), s <<= 1;return (k == p - 1) || ((s & 1)); } bool judge(LL n) {if(n == 1) return false;for(int i = 0; i < 7; i++) {if(n == prime[i]) return true;if(n % prime[i] == 0) return false;if(!judge(prime[i], n)) return false;}for(int i = 1; i <= 10; i++) if(!judge(2 + rand() % (n - 2), n)) return false;return true; }int main() {LL n;while(~scanf("%lld", &n)) printf(judge(n)? "Y\n" : "N\n");return 0; }

轉載于:https://www.cnblogs.com/olinr/p/10305930.html

總結

以上是生活随笔為你收集整理的loj #143. 质数判定的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。