cp's OEIS Frontend

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.

A109381 Maximum digit of n^2 written in factorial base.

Original entry on oeis.org

0, 1, 2, 1, 2, 1, 2, 2, 2, 3, 4, 1, 1, 2, 3, 4, 2, 2, 3, 3, 3, 3, 4, 4, 4, 5, 5, 1, 2, 1, 2, 2, 2, 3, 3, 4, 4, 5, 2, 3, 2, 2, 3, 3, 4, 4, 5, 3, 3, 3, 4, 3, 4, 5, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 1, 1, 2, 3, 4, 2, 2, 3, 4, 5, 3, 2, 3, 4, 3, 3, 3, 4, 4, 4, 4, 4, 5, 5, 5, 5, 6, 6, 6, 6, 2, 3, 4, 2
Offset: 0

Views

Author

Keywords

Comments

Conjecture: lim_{n->infinity} a(n) = infinity. If true, convergence is very slow, since a(1183893) = 1. Sequence is certainly unbounded, since for n >= 4, there is always a square between n*n! and (n+1)!.

Examples

			4^2 = 16 = 2*6+2*2 = 220(factorial base), so a(4) = max(2,2,0) = 2.
		

Crossrefs

Indices of 1's in this sequence are A014597.

Programs

  • Maple
    b:= proc(n, i) local r; `if`(n b(n^2, 2):
    seq(a(n), n=0..188);  # Alois P. Heinz, Mar 28 2018
  • Mathematica
    Block[{nn = 105^2, r}, r = Reverse@ Most@ NestWhileList[# + 1 &, 2, #! < nn &]; Array[Max@ IntegerDigits[#^2, MixedRadix@ r] &, Sqrt@ nn, 0]] (* Michael De Vlieger, Jan 01 2019 *)
  • PARI
    f(n, p=2) = if( nA007623
    a(n) = my(dfb=digits(f(n^2))); if (#dfb, vecmax(dfb), 0); \\ Michel Marcus, Mar 28 2018