离散对数&&大步小步算法及扩展
bsgs algorithm
ax≡b(mod n)
大步小步算法,这个算法有一定的局限性,只有当gcd(a,m)=1时才可以用
原理
此处讨论n为素数的时候。
ax≡b(mod n)(n为素数)
由费马小定理可知,只需要验证0,1,2...n-1是不是解即可,因为an-1 = 1mod(n)
算法过程
1、首先求出a0,a1,a2,...,am-1 模上n的值是否为b,存储在e[i]中,求出am的逆a-m
2、下面考虑am,am+1,...,a2m-1 模上n的值是否为b
此时不用一一检查,如果当中有解,相当于存在e[i],使得e[i] * am = b mod(n)
两边乘上a-m,e[i] = b * a-m mod(n),只需要检查存不存在这样的e[i]即可
3、同理,可以递推检查出a2m - a3m-1中解的情况
为了方便,把e[i]存储在map<int, int>x中,x[j]表示满足ei =j 的最小下标(因为可能有多个值相同)
map<ll, int>x;
ll log_mod(ll a, ll b, ll n)//n为素数
{
//if(b >= n)return -1;
a %= n;b %= n;//注意题目,如果b >= n是不存在解的
if(a == )
{
if(b == )return ;if(b == 1)return 0;
else return -;
}
ll m = ceil(sqrt(n + 0.5));
ll v = pow(a, n - - m, n);
x.clear();
x[] = m;
ll e = ;
for(int i = ; i < m; i++)//计算a的i次方mod n,并存下来
{
e = e * a % n;
if(!x[e])x[e] = i;
}
for(int i = ; i < m; i++)//计算a^(im), a^(im+1),...,a^(im+m-1)
{
int num = x[b];
if(num)return i * m + (num == m ? : num);
b = b * v % n;
}
return -;
}
扩展:
当n不为素数的时候,如何求解ax≡b(mod n) ?
转化成gcd(a, n) = 1即可
如何转化呢?
利用公式:ax≡b(mod n)⇔ax/d ≡ b/d (mod n/d),d=gcd(a,n)(这里是a乘上x,上述求解方程为ax)
每次用一个a和m和b消去gcd(a, n),消δ次,每次n /= g, b /= g,迭代更新下一个g,直到g=1
最后ax变成了ax - δ *a', 方程变成:ax - δ *a' = b' (mod n') 这里a' b' n'都是消因子之后的值
满足:a' * g = aδ b' * g = b n' * g = n g = gcd(aδ, n)
如果某一步g不整除b',直接返回-1
如果某一步b' = a',假设此时消因子次数达到cnt次,那就返回cnt
因为消因子次数cnt次等价于
原方程:ax-cntacnt = b (mod n)
消因子后:ax-cnta' = b' (mod n')
若此时b' = a' 那么等式就可以直接消去a' b',得到:ax-cnt = 1 (mod n),解就是cnt
下面求解:ax - δ *a' = b' (mod n')
可以求出a'的逆元,然后乘过去,用上述的大步小步算法求解
下面介绍一个技巧不用求逆元。
先把上述式子写成ax * tmp = b (mod n) 求出x之后加上δ就是解。
设m = sqrt(n +0.5), x = k * m - q (1 <= k <= m, q <= m)
上式可写成tmp * a k*m-q = b (mod n)
等价于 tmp * a k*m = b * aq(mod n)
可以从1-m枚举k,每次求出tmp * a k*m 判断是不是存在b*aq
最开始的大步小步算法也可以这样写
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll pow(ll a, ll b, ll m)
{
ll ans = ;
a %= m;
while(b)
{
if(b & )ans = (ans % m) * (a % m) % m;
b /= ;
a = (a % m) * (a % m) % m;
}
ans %= m;
return ans;
}
ll ext_log_mod(ll a, ll b, ll n)
{
if(b >= n)return -;//一些特殊情况的判断
a %= n;
if(b == )return ;
//if(n == 1)return -1;
ll cnt = ;//记录消因子次数
ll tmp = ;//存当前a'的值
for(ll g = __gcd(a, n); g != ; g = __gcd(a, n))
{
if(b % g)return -;//不能整除
b /= g; n /= g; tmp = tmp * a / g % n;
cnt++;
if(b == tmp)return cnt;
} ll m = sqrt(n + 0.5);
ll t = b;
map<ll, ll>Map;//记录b * a ^ i, i
Map[b] = ;
for(int i = ; i <= m; i++)
{
b = b * a % n;
Map[b] = i;
}
a = pow(a, m, n);
for(int k = ; k <= m; k++)//枚举k
{
tmp = tmp * a % n;//求出tmp*a^(k*m)
if(Map.count(tmp))return k * m - Map[tmp] + cnt;
}
return -;
}
int main()
{
ll a, b, p;
while(scanf("%lld%lld%lld", &a, &p, &b) != EOF)
{
ll ans = ext_log_mod(a, b, p);
if(ans == -)printf("Orz,I can’t find D!\n");
else printf("%lld\n", ans);
}
return ;
}
离散对数&&大步小步算法及扩展的更多相关文章
- 【题解】Matrix BZOJ 4128 矩阵求逆 离散对数 大步小步算法
传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=4128 大水题一道 使用大步小步算法,把数字的运算换成矩阵的运算就好了 矩阵求逆?这么基础的线 ...
- UVA 11916 Emoogle Grid 离散对数 大步小步算法
LRJ白书上的题 #include <stdio.h> #include <iostream> #include <vector> #include <mat ...
- 离散对数及其拓展 大步小步算法 BSGS
离散对数及其拓展 离散对数是在群Zp∗Z_{p}^{*}Zp∗而言的,其中ppp是素数.即在在群Zp∗Z_{p}^{*}Zp∗内,aaa是生成元,求关于xxx的方程ax=ba^x=bax=b的解, ...
- [模板]大步小步算法——BSGS算法
大步小步算法用于解决:已知A, B, C,求X使得 A^x = B (mod C) 成立. 我们令x = im - j | m = ceil(sqrt(C)), i = [1, m], j = [0, ...
- 大步小步算法模板题, poj2417
大步小步模板 (hash稍微有一点麻烦, poj不支持C++11略坑) #include <iostream> #include <vector> #include <c ...
- 洛谷 - P4861 - 按钮 - 扩展大步小步算法
https://www.luogu.org/problemnew/show/P4861 把好像把一开始b==1的特判去掉就可以AC了. #include<bits/stdc++.h> us ...
- BSGS-Junior·大步小步算法
本文原载于:http://www.orchidany.cf/2019/02/06/BSGS-junior/#more \(\rm{0x01}\) \(\mathcal{Preface}\) \(\rm ...
- [BSGS]大步小步算法
问题 BSGS被用于求解离散对数,即同余方程: \[ A^x\equiv B\pmod{P} \] 求\(x\)的最小非负整数解. 保证\(A\perp P\)(互质). 分析 首先,我们根据费马小定 ...
- BSGS算法(大步小步算法)
计算\(y^x ≡ z \ mod\ p\) 中 \(x\) 的解. 这个模板是最小化了\(x\) , 无解输出\(No \ Solution!\) map<ll,ll>data; ll ...
随机推荐
- HandlerInterceptor与MethodInterceptor
HandlerInterceptor是springMVC项目中的拦截器,它拦截的目标是请求的地址,比MethodInterceptor先执行.实现一个HandlerInterceptor拦截器可以直接 ...
- SVN使用指南
一:SVN服务器搭建和使用. 1. 首先来下载和搭建SVN服务器,下载地址如下: http://subversion.apache.org/packages.html,进入网址后,滚动到浏览器 ...
- Cheatsheet: 2018 04.01 ~ 04.30
Web Writing VS Code Extensions with TypeScript Docker How to write excellent Dockerfiles Raspberry P ...
- jsp技术知识点
1.jsp被Tomcat翻译成.java文件后,会被放在Tomcat安装目录下的\work\Catalina\localhost\station\org\apache\jsp文件夹下 2.El表达式表 ...
- win下gosublime配置ctag
ctags(Generate tag files for source code)是vim下方便代码阅读的工具.尽管ctags也可以支持其它编辑器,但是它正式支持的只有VIM.并且VIM中已经默认安装 ...
- Q:接口与抽象类
博文回答一下两个问题: 接口和抽象类的区别 选用接口和抽象类的依据 对于问题1: 从java语法的角度上来说,接口的所有成员和方法都是public的,且其方法均为abstract的.直到jdk1.8之 ...
- 5.Resource注解解析
Resource有两种使用场景 1.Resource 当Resource后面没带参数的时候是根据它所注释的属性名称到applicationContext.xml文件中查找是否有bean的id与之匹配, ...
- Oracle数据库采用数据泵方式导入导出数据
特别说明:Oralce的数据泵导入导出技术只能用在数据库服务器上,在只有客户端的机器上是无法使用数据泵技术的. 1.创建备份文件目录 mkdir d:\dmp 2.在Oralce中注册该目录,将目录 ...
- Code Signal_练习题_alternatingSums
Several people are standing in a row and need to be divided into two teams. The first person goes in ...
- 01:谁考了第k名 个人博客:doubleq.win
个人博客:doubleq.win 01:谁考了第k名 查看 提交 统计 提问 总时间限制: 1000ms 内存限制: 65536kB 描述 在一次考试中,每个学生的成绩都不相同,现知道了每个学生的 ...