BSGS_Baby steps giant steps算法

BSGS这个主要是用来解决这个题:

A^x=B(mod C)(C是质数),都是整数,已知A、B、C求x。

在具体的题目中,C一般是所有可能事件的总数。

解:

m = ceil(sqrt(C))(ceil为上取整), x = i * m + j

那么A^x = (A^m)^i * A^j, (0 <= i < m, 0 <= j < m)

然后可以枚举iO(sqrt(C))级别的枚举。

对于一个枚举出来的i,令D = (A^m)^i现在问题转化为求D * A^j ≡ B (mod C)

如果把A^j当作一个整体,那么套上exgcd就可以解出来了(而且因为C是质数,A是C的倍数的情况容易特判,除此之外必有

(D, C) = 1,所以一定有解):exgcd求逆元,可以参考:http://www.cnblogs.com/PJQOOO/p/3873654.html

求出了A^j,现在的问题就是我怎么知道j是多少?

先用O(sqrt(C))的时间,将A^j全部存进hash表里面。然后只要查表就在O(1)的时间内知道j是多少了。

时间: 2024-10-27 01:19:45

BSGS_Baby steps giant steps算法的相关文章

bsgs(Baby Steps Giant Steps)算法

BSGS算法(Baby Steps Giant Steps算法,大步小步算法,北上广深算法,拔山盖世算法) 适用问题 对于式子: $$x^y=z(mod_p)$$ 已知x,z,p,p为质数: 求解一个最小非负整数y: 存在一个y,属于[0,p-2](费马小定理) 于是有了一个笨拙的方法,枚举y 枚举y,期望效率:O(P) 寻求一种优化: 对式子变型: 设:$$y=i\sqrt{p}-j$$ 则$$x^{i\sqrt{p}-j}=z(mod_p)$$ ——这个变型的用意在于把y拆开 枚举y,变成枚

BSGS算法_Baby steps giant steps算法(无扩展)最强详解,你从未见过的详细

Baby Steps-Varsity Giant Step-Astronauts(May'n?椎名慶治) 阅读时可以听听这两首歌,加深对这个算法的理解.(Baby steps少女时代翻唱过,这个原唱反而不是很有名……Giant Step就比较碉,是一个假面骑士片的插曲,由超碉的May'n和一个人建立的临时组合唱的,怕不怕) 这个主要是用来解决这个题: A^x=B(mod C)(C是质数),都是整数,已知A.B.C求x. 我在网上看了好多介绍,觉得他们写得都不够碉,我看不懂…于是我也来写一发. 先

BSGS(Baby Steps,Giant Steps)

#include<bits/stdc++.h> using namespace std; typedef long long ll; ll n,m,y,z,p,x,ans,block; map<ll,int>mp; ll quick(ll a,ll b,ll p) { ll res=1%p; while (b) { if (b&1) { res=res*a%p; } a=a*a%p; b=b>>1; } return res; } void solve1() {

HDU 2815 扩展baby step giant step 算法

题目大意就是求 a^x = b(mod c) 中的x 用一般的baby step giant step 算法会超时 这里参考的是http://hi.baidu.com/aekdycoin/item/236937318413c680c2cf29d4 map平衡树查找值 1 #include <cstdio> 2 #include <cstring> 3 #include <iostream> 4 #include <cmath> 5 #include <

HDU 2815 Mod Tree 离散对数 扩展Baby Step Giant Step算法

链接:http://acm.hdu.edu.cn/showproblem.php?pid=2815 题意: 思路:与上题不同,这道题不要求m是素数,是利用扩展Baby Step Giant Step算法求离散对数. 以下转载自:AekdyCoin [扩展Baby Step Giant Step] [问题模型] 求解 A^x = B (mod C) 中 0 <= x < C 的解,C 无限制(当然大小有限制--) [写在前面] 这个问题比较麻烦,目前网络上流传许多版本的做法,不过大部分已近被证明

BSGS算法学习小记(大步小步算法)

简介 先看一个式子xy≡z(modp),z是质数 现在只知道x和z,要求y. 大步小步算法(BSGS,Baby Steps Giant Steps)就是解决这个问题. 算法流程 暴搜的枚举范围 根据费马小定理:xz?1≡1. 如果y已经枚举到了z-1了,继续枚举的话就会产生循环. 所以,在暴搜中y的枚举范围就是0--z-1. 如何优化暴搜 我们想一想可不可以用分块来解决枚举的y. 把y分成p?1????√分别枚举行不行? 设m=p?1????√,y=a?m+b,这样枚举a和b就相当于分块枚举了.

BSGS算法(大小步算法)

$BSGS$ 算法 $Baby\ Steps\ Giant\ Steps$. 致力于解决给定两个互质的数 $a,\ p$ 求一个最小的非负整数 $x$ 使得 $a^x\equiv b(mod\ p)$ 其中 $b$ 为任意正整数,$2≤a<p$,$2≤b<p$ 该算法使用的原理与欧拉定理有关,其中$a,\ p$互质 $a^{\phi (p)}\equiv 1(mod\ p)$ 又因为 $a^0\equiv 1(mod\ p)$ 所以$0到\phi p$是一个循环节,也就是说该算法最多查找$\p

【学术篇】SDOI2011 计算器

好一道三合一...(然而被我做成了四合一) 其实1 2 3是独立的OvO 然后就可以逐个分析了... 1 快速幂..就不说了..(我省选的时候有这么水的20pts部分分么←_← 2 两种做法(写在标题里面了).. 2.1(扩展欧几里得) \(xy\equiv z(mod\ p)\) 很显然可以写成\(xy=np+z\), 移项得\(xy-np=z\) 为保证\(y,p\)互质, 我们让两边同除\(gcd(y,p)\), 也就是让\(y,p,z\)分别除\(gcd(y,p)\) 由扩展欧几里得,

Steps UVA 846

说说:此题要求求出从整数x到达整数y所要经过的最短步数,且第一步和最后一步必须为一,同一时候每一步都比前一步多一步,少一步或一样.如果想搞清楚每一步详细是如何走的,那么这道题是相当麻烦的.考虑到前后两步之间最多差一,那么一到最大数之间的每个数从左到右以及从右到左的时候必然都出现.那么我们能够预想如果这样一个数列. 1,2,3,....n-2,n-1,n,n-1,n-2....3,2,1,若数列的和小于y和x的差肯定是不行的.所以要调整n至刚好大于或等于y和x的差.显然当y与x的差恰好为数列的和那