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.

A368925 The minimal exponent in the prime factorization of the powerful numbers.

This page as a plain text file.
%I A368925 #11 Feb 16 2024 02:12:20
%S A368925 0,2,3,2,4,2,3,5,2,2,6,2,4,2,2,2,3,7,2,2,2,2,3,2,5,8,2,2,2,3,2,2,2,3,
%T A368925 2,2,2,9,2,2,4,3,2,2,6,2,2,2,3,2,2,2,2,3,10,2,2,2,2,2,4,2,3,2,2,2,2,2,
%U A368925 2,2,2,3,2,2,2,2,3,3,2,11,2,7,3,2,2,2,4
%N A368925 The minimal exponent in the prime factorization of the powerful numbers.
%H A368925 Amiram Eldar, <a href="/A368925/b368925.txt">Table of n, a(n) for n = 1..10000</a>
%F A368925 a(n) = A051904(A001694(n)).
%F A368925 a(n) >= 2 for n >= 2.
%F A368925 Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 2.
%t A368925 s[n_] := If[n == 1, 0, Min @@ Last /@ FactorInteger[n]]; s /@ Select[Range[3000], # == 1 || Min[FactorInteger[#][[;;, 2]]] > 1 &]
%t A368925 (* or *)
%t A368925 f[n_] := Module[{e = Min[FactorInteger[n][[;; , 2]]]}, If[n == 1, 0, If[e > 1, e, Nothing]]]; Array[f, 3000]
%o A368925 (PARI) lista(kmax) = {my(e); for(k = 1, kmax, if(k == 1, print1(0, ", "), e = vecmin(factor(k)[,2]); if(e > 1, print1(e, ", "))));}
%Y A368925 Cf. A001694, A051904, A368710.
%K A368925 nonn,easy
%O A368925 1,2
%A A368925 _Amiram Eldar_, Feb 15 2024