BZOJ 1420: Discrete Root (原根+BSGS)
题意
已知kkk, aaa, ppp. 求 xk≡a (mod p)x^k\equiv a\ (mod\ p)xk≡a (mod p) 的所有根. 根的范围[0,p−1][0,p-1][0,p−1]. ppp为质数
分析
因为ppp是质数,那么一定有原根.设为ggg.
原根的性质如下:
- 对于[1,p−1][1,p-1][1,p−1]的所有iii,一定存在x∈[1,p−1]x\in[1,p-1]x∈[1,p−1]使得gx≡i (mod p)g^x\equiv i\ (mod\ p)gx≡i (mod p). 此时设xxx为I(i)I(i)I(i).
1.1.1.那么当aaa等于000
- 只有一个根就是000
2.2.2.当a∈[1,p−1]a\in[1,p-1]a∈[1,p−1]
- 画画柿子xk≡a (mod p)(gI(x))k≡a (mod p)(gk)I(x)≡a (mod p)\begin{aligned}x^k&\equiv a\ (mod\ p)\\(g^{I(x)})^k&\equiv a\ (mod\ p)\\(g^{k})^{I(x)}&\equiv a\ (mod\ p)\end{aligned}xk(gI(x))k(gk)I(x)≡a (mod p)≡a (mod p)≡a (mod p)
我们知道gkg^kgk,知道aaa,知道ppp.直接BSGSBSGSBSGS就行了.
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
int g, prime[100], cnt;
inline int qpow(int a, int b, int c) {
int re = 1;
while(b) {
if(b&1) re = 1ll * re * a % c;
a = 1ll * a * a % c; b >>= 1;
}
return re;
}
inline void Factor(int N) {
int tmp = N;
for(int i = 2; i*i <= N; ++i)
if(tmp % i == 0) {
prime[++cnt] = i;
while(tmp % i == 0) tmp /= i;
}
if(tmp > 1) prime[++cnt] = tmp;
}
inline int Get_g(int p) { //找原根
Factor(p-1);
for(int g = 2; ; ++g) {
bool flg = true;
for(int i = 1; i <= cnt; ++i)
if(qpow(g, (p-1)/prime[i], p) == 1)
{ flg = 0; break; }
if(flg) return g;
}
}
map<int, int>myhash;
vector<int>ans;
inline void Baby_Step_Giant_Step(int a, int b, int p) {
if(b == 1) ans.push_back(0);
myhash.clear();
int m = int(sqrt(p) + 1);
LL base = b;
for(int i = 0; i < m; ++i) {
myhash[base] = i;
base = 1ll * base * a % p;
}
base = qpow(a, m, p);
LL tmp = 1;
for(int i = 1; i <= m+1; ++i) {
tmp = 1ll * tmp * base % p;
if(myhash.count(tmp))
ans.push_back(i*m - myhash[tmp]);
}
}
int main() {
int p, k, a, g;
scanf("%d%d%d", &p, &k, &a);
if(!a) return puts("1"), puts("0"), 0;
//(g^I(x))^k = a (mod p)
//(g^k)^I(x) = a
Baby_Step_Giant_Step(qpow(g = Get_g(p), k, p), a, p);
for(int i = 0, siz = ans.size(); i < siz; ++i)
ans[i] = qpow(g, ans[i], p);
sort(ans.begin(), ans.end());
int siz = ans.size();
siz = unique(ans.begin(), ans.end()) - ans.begin();
printf("%d\n", siz);
for(int i = 0; i < siz; ++i)
printf("%d\n", ans[i]);
}
BZOJ 1420: Discrete Root (原根+BSGS)的更多相关文章
- bzoj 1420 Discrete Root - 原根 - exgcd - BSGS
题目传送门 戳我来传送 题目大意 给定$k, p, a$,求$x^{k}\equiv a \pmod{p}$在模$p$意义下的所有根. 考虑模$p$下的某个原根$g$. 那么$x = g^{ind_ ...
- BZOJ 1420 Discrete Root
思路:数学大汇总 提交:\(3\)次 错因:有一个\(j\)写成\(i\) 题解: 求:\(x^k \equiv a \mod p\) 我们先转化一下:求出\(p\)的原根\(g\) 然后我们用\(B ...
- BZOJ 3239 Discrete Logging(BSGS)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=3239 [题目大意] 计算满足 Y^x ≡ Z ( mod P) 的最小非负整数 [题解 ...
- BSGS 扩展大步小步法解决离散对数问题 (BZOJ 3239: Discrete Logging// 2480: Spoj3105 Mod)
我先转为敬? orz% miskcoo 贴板子 BZOJ 3239: Discrete Logging//2480: Spoj3105 Mod(两道题输入不同,我这里只贴了3239的代码) CODE ...
- 【BZOJ 1319】 Sgu261Discrete Rootsv (原根+BSGS+EXGCD)
1319: Sgu261Discrete Roots Time Limit: 1 Sec Memory Limit: 64 MBSubmit: 389 Solved: 172 Descriptio ...
- bzoj 3239: Discrete Logging && 2480: Spoj3105 Mod【BSGS】
都是BSGS的板子题 此时 \( 0 \leq x \leq p-1 \) 设 \( m=\left \lceil \sqrt{p} \right \rceil ,x=i*m-j \)这里-的作用是避 ...
- bzoj1420/1319 Discrete Root
传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=1420 http://www.lydsy.com/JudgeOnline/problem.ph ...
- POJ2417 Discrete Logging【BSGS】
Discrete Logging Time Limit: 5000MS Memory Limit: 65536K Total Submissions: 5577 Accepted: 2494 ...
- 51Nod1123 X^A Mod B 数论 中国剩余定理 原根 BSGS
原文链接https://www.cnblogs.com/zhouzhendong/p/51Nod1123.html 题目传送门 - 51Nod1123 题意 $T$ 组数据. 给定 $A,B,C$,求 ...
随机推荐
- vsftpd下载文件时内容乱码
windows客户端访问Linux服务端的ftp并下载文档时,内容会出现乱码,这是由于vsftpd文件服务器不支持转码功能 通过java FTPClient下载 方法为 OutputStream is ...
- 精通Java中的volatile关键字
在一些开源的框架的源码当中时不时都可以看到volatile这个关键字,最近特意学习一下volatile关键字的使用方法. 很多资料中是这样介绍volatile关键字的: volatile是轻量级的sy ...
- [转帖]关于Java中SQL语句的拼接规则
关于Java中SQL语句的拼接规则 自学demo 的时候遇到的问题 结果应该是 '"+e.getName()+"' 注意 一共有三组标点符号 (除去 方法函数后面的括号) 实现目标 ...
- 进行hcmcloud 数据库备份以及设置的处理过程.
导入数据库以及简单设置. 最近进行了一个数据库备份的简单工作: create database hcmcloud create user hcm with password 'Test6530' 执行 ...
- 获取web项目的绝对路径的方法总结
一.用Jsp获取 1.获取文件的绝对路径 String file="文件";(例如:data.mdb) String path=application.getRealPath(fi ...
- CrawlerRunner没有Log输出
官网log说明:https://docs.scrapy.org/en/latest/topics/logging.html#scrapy.utils.log.configure_logging 这里记 ...
- 201803-3 URL映射
问题描述 URL 映射是诸如 Django.Ruby on Rails 等网页框架 (web frameworks) 的一个重要组件.对于从浏览器发来的 HTTP 请求,URL 映射模块会解析请求中的 ...
- Hystrix服务容错保护
一.什么是灾难性雪崩效应? 造成灾难性雪崩效应的原因,可以简单归结为下述三种: 服务提供者不可用.如:硬件故障.程序BUG.缓存击穿.并发请求量过大等. 重试加大流量.如:用户重试.代码重试逻辑等. ...
- windows环境下安装selenium+python
selenium 是一个web的自动化测试工具,不少学习功能自动化的同学开始首选selenium ,相因为它相比QTP有诸多有点: * 免费,也不用再为破解QTP而大伤脑筋 * 小巧,对于不同的语 ...
- python基础知识0-2
# !/usr/bin/env python # 提示输入用户名和密码 # 验证用户名和密码# 如果错误,则输出用户名或密码错误# 如果成功,则输出 欢迎,XXX! import ge ...