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 A085156 #10 Apr 19 2025 03:27:55 %S A085156 1,2,3,4,5,6,7,8,9,10,11,13,14,15,16,17,19,21,22,23,25,26,27,29,31,32, %T A085156 33,34,35,36,37,38,39,41,43,46,47,49,51,53,55,57,58,59,61,62,64,65,67, %U A085156 69,71,73,74,77,79,81,82,83,85,86,87,89,91,93,94,95,97,100,101 %N A085156 Powers of primes or powers of semiprimes. %C A085156 m is a term iff A067029(m) = A071178(m) and A001221(m) <= 2. %H A085156 Charles R Greathouse IV, <a href="/A085156/b085156.txt">Table of n, a(n) for n = 1..10000</a> %t A085156 q[n_] := Module[{e = FactorInteger[n][[;;, 2]]}, Length[e] == 1 || (Length[e] == 2 && Equal @@ e)]; Select[Range[100], q] (* _Amiram Eldar_, Apr 19 2025 *) %o A085156 (PARI) is(n)=my(f=factor(n)[, 2]); #f<2 || (#f==2 && f[1]==f[2]) \\ _Charles R Greathouse IV_, Oct 19 2015 %Y A085156 Union of A000961 and A085155. %Y A085156 Cf. A001221, A067029, A071178. %K A085156 nonn %O A085156 1,2 %A A085156 _Reinhard Zumkeller_, Jun 21 2003