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 A342090 #21 Feb 17 2024 11:27:28
%S A342090 4,12,16,20,27,28,36,44,48,52,54,60,64,68,76,80,84,92,100,108,112,116,
%T A342090 124,132,135,140,144,148,156,164,172,176,180,188,189,192,196,204,208,
%U A342090 212,216,220,228,236,240,244,252,256,260,268,270,272,276,284,292,297,300
%N A342090 Numbers with at least one prime power p^e in their prime factorization such that p|e.
%C A342090 Numbers with a unitary divisor of the form p^(m*p) where p is a prime and m > 0.
%C A342090 The numbers of terms that do not exceed 10^k, for k = 1, 2, ..., are 1, 19, 188, 1883, 18825, 188244, 1882429, 18824297, 188242957, 1882429628, ...
%C A342090 The asymptotic density of this sequence is 1 - Product_{p prime} 1 - (p - 1)/(p*(p^p - 1)) = 0.18824296270011399086...
%H A342090 Amiram Eldar, <a href="/A342090/b342090.txt">Table of n, a(n) for n = 1..10000</a>
%e A342090 4 = 2^2 is a term since 2 divides 2.
%e A342090 8 = 2^3 is not a term since 2 does not divide 3.
%t A342090 q[n_] := AnyTrue[FactorInteger[n], Divisible[Last[#], First[#]] &]; Select[Range[2, 300], q]
%Y A342090 Subsequence of A013929.
%Y A342090 Subsequences: A000302, A000312, A009971, A017113, A051674.
%Y A342090 Cf. A072873, A369070 (characteristic function).
%K A342090 nonn
%O A342090 1,1
%A A342090 _Amiram Eldar_, Feb 27 2021
%E A342090 Wrong term 1 removed by _Amiram Eldar_, Jan 16 2024