|
MATHEMATICA
|
SemiprimeQ[n_Integer] := Plus @@ (Last /@ FactorInteger[n]) == 2 a = Table[ SeriesCoefficient[ Series[x/(x^3 - 2*x^2 - x + 1), {x, 0, 50}], n], {n, 0, 50}] f[n_] = If[SemiprimeQ[a[[n]]] == True, a[[n]], {}] Flatten[Table[f[n], {n, 1, Length[a]}]]
|