A356518 Maximal numerators in approximations to the Aurifeuillian factors of p^p +- 1.
2, 28, 1706, 25082, 816634, 157704814
Offset: 1
Examples
The r(n) are 2/3, 28/45, 1706/2835, 25082/42525, 816634/1403325, ...
This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
Patrick A. Thomas's wiki page.
Patrick A. Thomas has authored 43 sequences. Here are the ten most recent ones:
The r(n) are 2/3, 28/45, 1706/2835, 25082/42525, 816634/1403325, ...
870542161121 is the larger Aurifeuillian factor of 19^19 + 1, and 19 is the 4th term of A002145, so a(4) = 870542161121.
105318477338066161993242388018074119617 is the smaller Aurifeuillian factor of 47^47 + 1, and 47 is the 8th term of A002145, so it is a(8).
813955076015309926319 is the larger Aurifeuillian factor of 29^29-1, and 29 is the 4th term of A002144, so a(4) = 813955076015309926319.
112663560435723374699 is the smaller Aurifeuillian factor of 29^29-1, and 29 is the 4th term of A002144, so a(4) = 112663560435723374699.
8^25 in base 20 ends in the digits 13, 18, 8 (or ...DI8 in extended hexadecimal notation).
x = ...65762C0520697E8CA1A31469. x^2 = ...8978C2E9CE8570624D4BDA86 = A331549. x^3 = ...65762C0520697E8CA1A31469 = x.
\\ See A331548 with initial b=9 instead of b=3.
Vecrev(digits(lift((9+O(15^99))^5^99),15)) \\ M. F. Hasler, Jan 26 2020
x^2 = ...8978C2E9CE8570624D4BDA86 = x.
\\ See A331548, with initial b=6 instead of b=3.
A331549_vec(n)=Vecrev(digits(lift(chinese(Mod(0,3^n),Mod(1,5^n))),15)) \\ or simpler but slower: Vecrev(digits(lift(Mod(3^4,15^n)^5^(n-1)),15)) \\ M. F. Hasler, Jan 26 2020
x equals the limit of the (n+1) trailing digits of 3^(5^n): 3^(5^0) = (3), 3^(5^1) = 1(13), 3^(5^2) = 1708EB01(D13), ... x = ...2AA66B44A40E43797853AD13. x^2 = ...65762C0520697E8CA1A31469 = A331550. x^3 = ...C44883AA4AE0AB75769B41DC = -x. x^4 = ...8978C2E9CE8570624D4BDA86 = A331549. x^5 = ...2AA66B44A40E43797853AD13 = x.
\\ after Paul D. Hanna's program in A120817 {a(n)=local(b=3, v=[]); for(k=1, n+1, b=b^5%15^k; v=concat(v, (15*b\15^k))); v[n+1]} for(k=0,80,print1(a(k),", ")) \\ Hugo Pfoertner, Jan 26 2020
(A331548_vec(n)=Vecrev(digits(lift(Mod(3,15^n)^5^(n-1)),15)))(99) \\ M. F. Hasler, Jan 26 2020
Comments