|
|
| Line 238: |
Line 238: |
| end | | end |
| return d | | return d |
| end
| |
|
| |
| -- stylua: ignore
| |
| p.primes = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97,
| |
| 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199,
| |
| 211, 223, 227, 229, 233, 239, 241, 251, 257, 263, 269, 271}
| |
|
| |
| -- get monzo of n/d
| |
| -- e.g. for 3/2: {[2] = -1, [3] = 1}
| |
| function p.get_monzo(n, d)
| |
| local n_pf = p.prime_factorization_raw(n)
| |
| local d_pf = p.prime_factorization_raw(d)
| |
| local result = {}
| |
| for i = 1, #p.primes do
| |
| local t = (n_pf[p.primes[i]] or 0) - (d_pf[p.primes[i]] or 0)
| |
| if t ~= 0 then
| |
| result[p.primes[i]] = t
| |
| end
| |
| end
| |
| return result
| |
| end | | end |
|
| |
|
| return p | | return p |