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.
%I A078316 #14 Sep 08 2024 22:33:43 %S A078316 1,1,1,2,1,1,2,1,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,1,1,1,1, %T A078316 1,1,1,2,1,1,2,1,2,1,2,1,1,2,3,1,1,1,1,2,1,1,3,1,1,1,1,1,2,1,1,1,1,2, %U A078316 1,2,1,1,1,1,1,3,1,1,1,2,2,2,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,2,1,1,1,1,1,1 %N A078316 Maximum exponent in the prime factorization of n*rad(n)+1, where rad = A007947 (squarefree kernel). %H A078316 Reinhard Zumkeller, <a href="/A078316/b078316.txt">Table of n, a(n) for n = 1..10000</a> %F A078316 a(n) = A051903(A078310(n)). %t A078316 a[n_] := Max[FactorInteger[1 + n * Times @@ FactorInteger[n][[;;, 1]]][[;;, 2]]]; Array[a, 100] (* _Amiram Eldar_, Sep 07 2024 *) %o A078316 (Haskell) %o A078316 a078316 = a051903 . a078310 -- _Reinhard Zumkeller_, Jul 23 2013 %o A078316 (PARI) a(n) = vecmax(factor(1 + n * vecprod(factorint(n)[, 1]))[, 2]); \\ _Amiram Eldar_, Sep 07 2024 %Y A078316 Cf. A007947, A051903, A078310, A078313, A078315. %K A078316 nonn %O A078316 1,4 %A A078316 _Reinhard Zumkeller_, Nov 23 2002