A303661 Powers of squarefree semiprimes that are not squarefree.

36, 100, 196, 216, 225, 441, 484, 676, 1000, 1089, 1156, 1225, 1296, 1444, 1521, 2116, 2601, 2744, 3025, 3249, 3364, 3375, 3844, 4225, 4761, 5476, 5929, 6724, 7225, 7396, 7569, 7776, 8281, 8649, 8836, 9025, 9261, 10000, 10648, 11236, 12321, 13225, 13924, 14161, 14884

Offset

1,1

Links

Amiram Eldar, Table of n, a(n) for n = 1..10000

Eric Weisstein's World of Mathematics, Squarefree.

Eric Weisstein's World of Mathematics, Semiprime.

Formula

Sum_{n>=1} 1/a(n) = Sum_{n>=1} 1/((A006881(n)-1)*A006881(n)) = Sum_{k>=2} (P(k)^2 - P(2*k))/2 = 0.07160601536406295068..., where P(k) is the prime zeta function. - Amiram Eldar, Feb 12 2021

Example

1089 is in the sequence because 1089 = 3^2*11^2.
1296 is in the sequence because 1296 = 2^4*3^4.

Mathematica

Select[Range[15000], Length[Union[FactorInteger[#][[All, 2]]]] == 1 && PrimeNu[#] == 2 && ! SquareFreeQ[#] &]
seq[max_] := Module[{sp = Select[Range[Floor@Sqrt[max]], SquareFreeQ[#] && PrimeNu[#] == 2 &], s = {}}, Do[s = Join[s, sp[[k]]^Range[2, Floor@Log[sp[[k]], max]]], {k, 1, Length[sp]}]; Union@s]; seq[10000] (* Amiram Eldar, Feb 12 2021 *)

Crossrefs

Cf. A006881, A013929, A072774, A072777, A085155, A123711, A200511, A246547, A303606.

Sequence in context: A044604 A340017 A303606 * A175391 A030627 A177492

Adjacent sequences: A303658 A303659 A303660 * A303662 A303663 A303664

Keyword

nonn

Author

Ilya Gutkovskiy, Apr 28 2018

Status

approved