arrow-left

All pages
gitbookPowered by GitBook
1 of 1

Loading...

Discrete Log Problem

hashtag
Discrete log problem

Given any groupGGGand elementsa,ba,ba,bsuch that an=ba^n=ban=b, the problem of solving fornnnis known as the disctete log problem (DLP). In sage, this can be done for general groups by calling discrete_log

hashtag
Discrete log over

Typically, one considers the discrete log problem in , i.e. the multiplicative group of integers. Explicitly, the problem asks forgiven . This can be done by calling b.log(a) in sage:

This section is devoted to helping the reader understand which functions are called when for this specific instance of DLP.

Whenis composite and not a prime power, discrete_log() will be used, which uses generic algorithms to solve DLP (e.g. Pohlig-Hellman and baby-step giant-step).

When is a prime, Pari znlog will be used, which uses a linear sieve index calculus method, suitable for .

When , SageMath will fall back on the generic implementation discrete_log()which can be slow. However, Pari znlog can handle this as well, again using the linear sieve index calculus method. To call this within SageMath we can use either of the following (the first option being a tiny bit faster than the second)

hashtag
Example

Given a small prime, we can compare the Pari method with the Sage defaults

We can also solve this problem with even larger primes in a very short time

hashtag
Discrete log over

// elliptic curve discrete log functions

sage: G = DihedralGroup(99)
sage: g = G.random_element()
sage: discrete_log(g^9,g) # note that if the order of g is less than 9 we would get 9 mod g.order()
9
(Z/nZ)∗\left(\mathbb Z/n\mathbb Z\right)^*(Z/nZ)∗
(Z/nZ)∗\left(\mathbb Z/n\mathbb Z\right)^*(Z/nZ)∗
mod n\text{mod }nmod n
xxx
ax=b(modn)a^x=b\pmod nax=b(modn)
nnn
n=pn=pn=p
p<1050∼2166p < 10^{50} \sim 2 ^{166}p<1050∼2166
n=pkn = p^kn=pk
E(k)E(k)E(k)
sage: R = Integers(99)
sage: a = R(4)
sage: b = a^9
sage: b.log(a)
9
x = int(pari(f"znlog({int(b)},Mod({int(a)},{int(n)}))"))
x = gp.znlog(b, gp.Mod(a, n))
p = getPrime(36)
n = p^2
K = Zmod(n)
a = K.multiplicative_generator()
b = a^123456789

time int(pari(f"znlog({int(b)},Mod({int(a)},{int(n)}))")) 
# CPU times: user 879 µs, sys: 22 µs, total: 901 µs
# Wall time: 904 µs
# 123456789

time b.log(a)
# CPU times: user 458 ms, sys: 17 ms, total: 475 ms
# Wall time: 478 ms
# 123456789

time discrete_log(b,a)
# CPU times: user 512 ms, sys: 24.5 ms, total: 537 ms
# Wall time: 541 ms
# 123456789
p = getPrime(100)
n = p^2
K = Zmod(n)
a = K.multiplicative_generator()
b = a^123456789

time int(pari(f"znlog({int(b)},Mod({int(a)},{int(n)}))")) 
# CPU times: user 8.08 s, sys: 82.2 ms, total: 8.16 s
# Wall time: 8.22 s
# 123456789