OFFSET
0,3
COMMENTS
a(n) = smallest m such that A007953(n^m) > n, or 0 if no such m exists.
If n is a term of A178501, then a(n) = 0. - Felix Fröhlich, Jun 10 2018
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..10000
MATHEMATICA
myLen::usage = "myLen[n] gives the smallest m > 1 such that Total[IntegerDigits[n^m]] > n.";
myLen[n_Integer] := Module[{m = 2}, If[Total[IntegerDigits[n]] < 2, 0, While[Total[IntegerDigits[n^m]] <= n, m++]; m]];
Table[myLen[n], {n, 0, 100}] (* Vincenzo Librandi, Nov 19 2025 *)
PROG
(PARI) a(n) = if(sumdigits(n) < 2, return(0), my(m=2); while(1, if(sumdigits(n^m) > n, return(m)); m++)) \\ Felix Fröhlich, Jun 10 2018
(Magma) a := function(n) if &+Intseq(n) lt 2 then return 0; end if; m := 2; while &+Intseq(n^m) le n do m +:= 1; end while; return m; end function;
[ a(n) : n in [0..100] ]; // Vincenzo Librandi, Nov 18 2025
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Max Alekseyev, Jun 08 2018
STATUS
approved
