Anonymous user
P-Adic square roots: Difference between revisions
m
no edit summary
mNo edit summary |
|||
Line 96:
function padic.sqrt (byref g as ratio, byval sw as integer) as integer
dim as longint a = g.a, b = g.b
dim as longint q, x, pk, pm
dim as long f1, r, s, t
dim i as integer, f as double
Line 217:
pk = 1
for i = v + 1 to k - 1 + v
pm = pk: pk *= p
if pk \ pm - p then exit for
evalf(x)
d(i) = q * f1 mod p
| |||