| FazBrowse GitHub Viewer | Trending | | Home |
| Tools: [Download Repo ZIP] [Original HTTPS Page] |
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,79 @@ | ||
| /** | ||
| * Repunit theorem helpers. | ||
| * | ||
| * A repunit of length n is: | ||
| * R_n = (10^n - 1) / 9 | ||
| * | ||
| * For a prime p (p != 2, 3, 5), p divides R_n iff ord_p(10) divides n. | ||
| * Reference: https://en.wikipedia.org/wiki/Repunit | ||
| */ | ||
|
|
||
| const gcd = (a, b) => { | ||
| let x = BigInt(a) | ||
| let y = BigInt(b) | ||
| while (y !== 0n) { | ||
| ;[x, y] = [y, x % y] | ||
| } | ||
| return x < 0n ? -x : x | ||
| } | ||
|
|
||
| const modPow = (base, exp, mod) => { | ||
| let result = 1n | ||
| let b = BigInt(base) % BigInt(mod) | ||
| let e = BigInt(exp) | ||
| const m = BigInt(mod) | ||
|
|
||
| while (e > 0n) { | ||
| if (e & 1n) result = (result * b) % m | ||
| b = (b * b) % m | ||
| e >>= 1n | ||
| } | ||
|
|
||
| return result | ||
| } | ||
|
|
||
| const multiplicativeOrder10 = (prime) => { | ||
| const p = BigInt(prime) | ||
| if (p <= 1n) throw new RangeError('prime must be > 1') | ||
| if (gcd(10n, p) !== 1n) throw new RangeError('10 and prime must be coprime') | ||
|
|
||
| // For prime p, ord_p(10) divides p-1. | ||
| const upper = p - 1n | ||
| for (let k = 1n; k <= upper; k++) { | ||
| if (upper % k === 0n && modPow(10n, k, p) === 1n) { | ||
| return k | ||
| } | ||
| } | ||
|
Comment thread
nickzerjeski marked this conversation as resolved.
|
||
|
|
||
| throw new Error('multiplicative order not found') | ||
| } | ||
|
|
||
| const repunitMod = (length, mod) => { | ||
| if (!Number.isInteger(length) || length < 1) { | ||
| throw new RangeError('length must be a positive integer') | ||
| } | ||
| const m = BigInt(mod) | ||
| if (m <= 0n) throw new RangeError('mod must be > 0') | ||
|
|
||
| let remainder = 0n | ||
| for (let i = 0; i < length; i++) { | ||
| remainder = (remainder * 10n + 1n) % m | ||
| } | ||
| return remainder | ||
| } | ||
|
|
||
| const isRepunitDivisibleByPrime = (length, prime) => { | ||
| if (!Number.isInteger(length) || length < 1) { | ||
| throw new RangeError('length must be a positive integer') | ||
| } | ||
|
|
||
| const p = BigInt(prime) | ||
| if (p === 2n || p === 5n) return false | ||
| if (p === 3n) return BigInt(length) % 3n === 0n | ||
| if (gcd(10n, p) !== 1n) return false | ||
|
|
||
| const order = multiplicativeOrder10(p) | ||
| return BigInt(length) % order === 0n | ||
|
Comment thread
nickzerjeski marked this conversation as resolved.
|
||
| } | ||
|
|
||
| export { multiplicativeOrder10, repunitMod, isRepunitDivisibleByPrime } | ||
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,50 @@ | ||
| import { | ||
| isRepunitDivisibleByPrime, | ||
| multiplicativeOrder10, | ||
| repunitMod | ||
| } from '../RepunitTheorem' | ||
|
|
||
| describe('RepunitTheorem', () => { | ||
| it('computes multiplicative order examples', () => { | ||
| expect(multiplicativeOrder10(11n)).toBe(2n) | ||
| expect(multiplicativeOrder10(37n)).toBe(3n) | ||
| expect(multiplicativeOrder10(7n)).toBe(6n) | ||
| }) | ||
|
|
||
| it('checks repunit divisibility using the theorem', () => { | ||
| // 111111 is divisible by 3, 7, 11, 13, 37 | ||
| expect(isRepunitDivisibleByPrime(6, 3n)).toBe(true) | ||
| expect(isRepunitDivisibleByPrime(6, 7n)).toBe(true) | ||
| expect(isRepunitDivisibleByPrime(6, 11n)).toBe(true) | ||
| expect(isRepunitDivisibleByPrime(6, 13n)).toBe(true) | ||
| expect(isRepunitDivisibleByPrime(6, 37n)).toBe(true) | ||
| }) | ||
|
|
||
| it('returns false when divisibility condition does not hold', () => { | ||
| expect(isRepunitDivisibleByPrime(1, 3n)).toBe(false) | ||
| expect(isRepunitDivisibleByPrime(3, 3n)).toBe(true) | ||
| expect(isRepunitDivisibleByPrime(6, 19n)).toBe(false) | ||
| expect(isRepunitDivisibleByPrime(9, 2n)).toBe(false) | ||
| expect(isRepunitDivisibleByPrime(9, 5n)).toBe(false) | ||
| }) | ||
|
Comment thread
nickzerjeski marked this conversation as resolved.
|
||
|
|
||
| it('computes repunit modulo without building huge integers', () => { | ||
| expect(repunitMod(6, 37n)).toBe(0n) | ||
| expect(repunitMod(6, 11n)).toBe(0n) | ||
| expect(repunitMod(7, 13n)).toBe(1n) | ||
| }) | ||
|
|
||
| it('validates multiplicative order input', () => { | ||
| expect(() => multiplicativeOrder10(1n)).toThrow(RangeError) | ||
| expect(() => multiplicativeOrder10(10n)).toThrow(RangeError) | ||
| }) | ||
|
|
||
| it('validates repunitMod input', () => { | ||
| expect(() => repunitMod(0, 7n)).toThrow(RangeError) | ||
| expect(() => repunitMod(5, 0n)).toThrow(RangeError) | ||
| }) | ||
|
|
||
| it('validates repunit divisibility input length', () => { | ||
| expect(() => isRepunitDivisibleByPrime(0, 7n)).toThrow(RangeError) | ||
| }) | ||
| }) | ||
| Back | FazBrowse Home | New Git URL |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Choose a reason Spam Abuse Off Topic Outdated Duplicate Resolved Low QualityThe file header comment says the theorem holds for primes p != 2,5, but it also needs p != 3 (equivalently gcd(p,9)=1). As written, it documents an incorrect condition and matches the buggy behavior for p=3.
Sorry, something went wrong.
Uh oh!
There was an error while loading. Please reload this page.