【代码超详解】LOJ #143 质数判定（Miller-Rabin算法 · 模板题）

最新推荐文章于 2024-02-01 15:09:49 发布
原创 最新推荐文章于 2024-02-01 15:09:49 发布 · 775 阅读 · 0 ·
CC 4.0 BY-SA版权
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

博客详细解析了LOJ #143质数判定问题，介绍了Miller-Rabin算法，并指出许多AC代码无法通过所有数据，强调了使用__int128类型的重要性，以及在不同平台的兼容性问题。

一、题目描述

#143. 质数判定

内存限制：256 MiB
时间限制：5000 ms
标准输入输出
题目类型：传统
评测方式：文本比较
上传者： fjzzq2002

题目描述

判定输入的数是不是质数。

输入格式

若干行，一行一个数x 。

行数不超过1e5 。

输出格式

对于输入的每一行，如果 x是质数输出一行

Y

否则输出一行

N

样例

样例输入

1
2
6
9
666623333

样例输出

N
Y
N
N
Y

数据范围与提示

1 ≤ x ≤1e18
欢迎hack（如果你不是管理员，可以在题目讨论区发帖）。

二、算法分析说明与代码编写指导


三、AC 代码

尼玛，这题巨几把毒瘤，搞了我整整3天。

快速幂要用__int128作为中间变量类型的版本，不然各种WA。不过VS是没有支持这个类型的，就没法本地调试了，除非临时改掉这个类型，提交的时候改回去。LOJ和洛谷都可以支持__int128，HDU不行。

这题好多网上AC代码根本就不能全部数据都过，估计数据是反复被加强过了。

老子就想找个说得清楚点的Miller-Rabin算法的解释和对应代码，结果这么多代码全都过不了。至于上面的算法概述，这里抄一点那里借来改一下再补充一些注意事项，再缩减一下，就成了现在这样子。整理过程真是吐血啊。

#include<cstdio>
#include<algorithm>
#pragma warning(disable:4996)
using namespace std;
unsigned long long x, prime[] = {
   
    2,3,</
最低0.47元/天 解锁文章
【LOJ143】质数判定
cz_xuyixuan的博客
09-02 689
【题目链接】 点击打开链接 【思路要点】 MillerRabinMillerRabinMillerRabin 素数测试模板题。 时间复杂度 O(TLog2N)O(TLog2N) O(TLog^2N) 。 【代码】 #include&amp;lt;bits/stdc++.h&amp;gt; using namespace std; const int...
LOJ #143. 质数判定 题解
bifanwen的博客
07-19 414
博客园同步 原题链接 简要题意： 给定 TTT 个数 nnn，判素数。 T≤105,n≤1018T \leq 10^5 , n \leq 10^{18}T≤105,n≤1018. 可能你判一个都有困难是不是 ⋯⋯\cdots \cdots⋯⋯. 二次探测定理 若 x2≡1(modp),x<px^2 \equiv 1 \pmod p , x < px2≡1(modp),x<p，则 x=1x= 1x=1 或 x=p−1x = p-1x=p−1. 简单证明： x2≡1(modp)x^2 \eq
miller-rabin素性检测算法源码
06-03
miller-rabin素性检测算法的源代码 能够运行，很好的资源哦
Loj#143-[模板]质数判定【Miller-Rabin】
QuantAsk
01-08 444
正题 题目链接:https://loj.ac/p/143 题目大意 给出一个数ppp，让你判定是否为质数。 解题思路 Miller−RabinMiller-RabinMiller−Rabin是一种基于费马小定理和二次探测定理的具有较高正确性的高效质数判定算法。 首先讲一下两个定理 费马小定理：gcd(a,p)=1   ⇒   ap−1=1(mod p)gcd(a,p)=1\ \ \ \Rightarrow\ \ \ a^{p-1}=
LOJ143（Miller Rabin素数测试）
IQ^QI的博客
09-07 478
LOJ143 题意：1e5个(1e18)的数，分别判断它们是不是质素。 题解：Miller Rabin素数测试(费马小定理，二次探测） 当被测试数在某一个范围内时，我们可以通过选取适当的测试底数让 MillerRabin 素数测试对于这个范围内的每一个数都能够正确地得出结果。下面是节选自 维基百科 的选取 a的方式。 当 N<4,759,123,141选取 a=2,7,61即可确保算...
Miller–Rabin 素性测试
jxy250的博客
02-01 1616
Miller–Rabin 素性测试
【代码超详解】LOJ #132. 树状数组 3 ：区间修改，区间查询
COFACTOR
02-07 723
一、题目描述 样例 样例输入 5 10 2 6 6 1 1 2 1 4 1 2 5 10 2 1 3 2 2 3 1 2 2 8 1 2 3 7 1 4 4 10 2 1 2 1 4 5 6 2 3 4 样例输出 15 34 32 33 50 二、算法分析说明与代码编写指导 如果不了解树状数组，建议先阅读之前的题解： https://blog.youkuaiyun.com/COFACTOR/arti...
【学习笔记】Miller–Rabin素数测试
cz_xuyixuan的博客
09-02 4574
【算法简介】 MillerRabinMillerRabinMillerRabin 素数测试是一种判断一个数是否是质数的方式。 其单次测试的时间复杂度不会超过 O(Log2N)O(Log2N) O(Log^2N) ，期望为 O(LogN)O(LogN) O(LogN) ，几乎不需要额外的空间。 MillerRabinMillerRabinMillerRabin 素数测试不是...
LOJ #143. 质数判定
weixin_30508241的博客
09-11 479
题目描述 判定输入的数是不是质数。 输入格式 若干行，一行一个数x。 行数不超过1.5×104​​。 输出格式 对于输入的每一行，如果x是质数输出一行Y，否则输出一行N。 样例 样例输入 1 2 6 9 666623333 样例输出 ...
loj #143. 质数判定
weixin_30416871的博客
01-22 321
\(\color{#0066ff}{ 题目描述 }\) 判定输入的数是不是质数。 \(\color{#0066ff}{输入格式}\) 若干行，一行一个数 \(x\)。 行数不超过 \(10^5\)。 \(\color{#0066ff}{输出格式}\) 对于输入的每一行，如果 \(x\) 是质数输出一行 \(Y\)，否则输出一行 \(N\)。 \(\color{#0066ff}{输入样例}\) 1 ...
Miller_Rabin 米勒拉宾概率筛【模板】
jisuanji2606414的博客
07-30 363
Miller_Rabin筛法是一种正确率极高，复杂度优秀的概率筛法，将质数分为2和奇质数，对于奇质数n，将n-1其分解为2^p+m，然后，取若干[1,n-1]的随机数a,本模板取了20次，每次判断是否有。a^m三1（modn)或者是否有[0,p-1]的j,满足a^(2jm)三n-1(modn)验证正确性如下，取2-1e6全部数字，进行正确性判断，正确个数1e6-1，全部正确。甚至可以将判断次数缩小至5，也能保证完全正确。......
Miller-Rabin测试素数模板
bahuan1974的博客
07-14 171
ll qpow(ll a,ll b,ll M) //快速幂 { ll ans =1; while(b) { if(b&1) ans*=a;ans%=M; a*=a;a%=M;b>>=1; } return ans; } bool Mille...
HDU 2138 Miller-Rabin 模板题
weixin_30405421的博客
09-19 174
求素数个数。 /** @Date : 2017-09-18 23:05:15 * @FileName: HDU 2138 miller-rabin 模板.cpp * @Platform: Windows * @Author : Lweleth (SoungEarlf@gmail.com) * @Link : https://github.com/ ...
miller_rabin 模板
lifehappy
07-22 400
miller_rabin 模板 ll quick_mult(ll a, ll b, ll mod) { ll ans = 0; while(b) { if(b & 1) ans = (ans + a) % mod; a = (a + a) % mod; b >>= 1; } return ans; } ll quick_pow(ll a, ll n, ll mod) { ll ans = 1;
ComSec作业一：Miller-Rabin算法---编程题
HHHHUIIII的博客
10-02 225
Miller-Rabin算法 算法的理论基础： Fermat定理：若n是奇素数，a是任意正整数(1≤ a≤ n−1)，则 a^(n-1) ≡ 1 mod n。 推演自Fermat定理, 如果n是一个奇素数，将n−1表示成2^s*r 的形式，r是奇数，a与n是互素的任何随机整数，那么a^r ≡ 1 mod n或者对某个j (0 ≤ j≤ s−1, j∈Z) 等式a^(2jr) ≡ −1 mod n 成立。 #include <iostream> #include <cstdio>
质数判定
Dh_virgil的博客
10-13 130
这道题呢数据量显然是非常大的，所以容易超时，普通的优化做不到AC，所以要优化算法。对于输入的每一行，如果 x是质数输出一行Y ，否则输出一行 N。Miller-Rabin素性测试算法。判定输入的数是不是质数。若干行，一行一个数 x。行数不超过 10^5。
Miller_Rabin和Pollard_rho素数模板——POJ-2429的题解
AlphaWA的博客
03-14 347
网络上有很多讲解，不再赘述。此处记载此题题解以供复习之用。#include&lt;cstdio&gt; #include&lt;iostream&gt; #include&lt;cstdlib&gt; #include&lt;cmath&gt; #include&lt;set&gt; #include&lt;algorithm&gt; using namespace std; typedef l...
miller_rabin 证明与实现
lifehappy
07-22 526
143. 质数判定__模板题链接 前置知识 费马小定理：ppp是质数，则对于任意的aaa，aaa与ppp互质，则有ap−1≡1(modp)a ^ {p - 1} \equiv 1 \pmod {p}ap−1≡1(modp)。 二次探测定理：如果ppp是一个质数，x2≡1(modp)x ^ 2 \equiv 1 \pmod px2≡1(modp)，则有解为x1=1,x2=p−1x_1 = 1, x_2 = p - 1x1​=1,x2​=p−1，推理如下 ⇒x2−1≡0(modp)\Rightarrow x
【LOJ143】素数判定
weixin_30362233的博客
08-20 141
题面 https://loj.ac/problem/143 题解 直接$Miller$就好。 #pragma GCC optimize(2) #include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #define LLL __i...
loj#P188 可并堆
最新发布
11-15
可并堆是一种支持合并操作的堆数据结构，常见的可并堆有左偏树、斜堆、二项堆等。对于 LOJ#P188 可并堆的问题，下面以左偏树为例给出解题思路和代码实现。 ### 解题思路 1. **左偏树的性质**： - 左偏树是一种可并堆，它满足堆性质（小根堆或大根堆），即每个节点的值小于（或大于）其子节点的值。 - 左偏树还满足左偏性质，即每个节点的左子树的距离（到最近的叶子节点的距离）不小于右子树的距离。 2. **合并操作**： - 合并两个左偏树时，比较两个根节点的值，将值较大的根节点的树合并到值较小的根节点的右子树中。 - 合并后，检查右子树的距离是否大于左子树的距离，如果是，则交换左右子树，以维护左偏性质。 3. **插入操作**： - 插入一个新节点可以看作是合并一个只有一个节点的左偏树和原左偏树。 4. **删除操作**： - 删除根节点后，将其左右子树合并成一个新的左偏树。 ### 代码实现 ```python class Node: def __init__(self, val): self.val = val self.left = None self.right = None self.dist = 0 def merge(x, y): if not x: return y if not y: return x if x.val > y.val: x, y = y, x x.right = merge(x.right, y) if not x.left or (x.right and x.left.dist < x.right.dist): x.left, x.right = x.right, x.left x.dist = (x.right.dist + 1) if x.right else 0 return x def insert(root, val): new_node = Node(val) return merge(root, new_node) def delete(root): return merge(root.left, root.right) # 示例使用 root = None root = insert(root, 3) root = insert(root, 1) root = insert(root, 5) print(root.val) # 输出堆顶元素 root = delete(root) print(root.val) # 输出删除堆顶元素后的堆顶元素 ```
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
表情包 代码片
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值