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.

A134618 Numbers such that the sum of cubes of their prime factors (taken with multiplicity) is a prime.

This page as a plain text file.
%I A134618 #19 Mar 25 2025 10:49:00
%S A134618 12,28,40,45,48,52,54,56,63,75,80,96,104,108,117,136,152,153,165,175,
%T A134618 210,224,232,245,250,261,268,300,320,325,333,344,350,363,384,387,390,
%U A134618 399,405,416,432,462,464,468,475,477,504,507,531,536,539,561,570,584
%N A134618 Numbers such that the sum of cubes of their prime factors (taken with multiplicity) is a prime.
%H A134618 Harvey P. Dale and Hieronymus Fischer, <a href="/A134618/b134618.txt">Table of n, a(n) for n = 1..10000</a> (first 1000 terms from _Harvey P. Dale_)
%F A134618 {k: A224787(k) in A000040}. - _R. J. Mathar_, Mar 25 2025
%e A134618 a(2) = 28, since 28 = 2*2*7 and 2^3 + 2^3 + 7^3 = 359 which is prime.
%t A134618 Select[Range[600],PrimeQ[Total[Flatten[Table[#[[1]],{#[[2]]}]&/@ FactorInteger[#]]^3]]&] (* _Harvey P. Dale_, Feb 01 2013 *)
%o A134618 (Python)
%o A134618 from sympy import factorint, isprime
%o A134618 def ok(n): return isprime(sum(p**3 for p in factorint(n, multiple=True)))
%o A134618 print([k for k in range(585) if ok(k)]) # _Michael S. Branicky_, Dec 28 2021
%Y A134618 Cf. A001597, A025475, A134333, A134344, A134376.
%Y A134618 Cf. A134600, A134602, A134605, A134608, A134612, A134616, A134620.
%K A134618 nonn
%O A134618 1,1
%A A134618 _Hieronymus Fischer_, Nov 11 2007
%E A134618 Example clarified by _Harvey P. Dale_, Feb 01 2013
%E A134618 Minor edits by _Hieronymus Fischer_, May 06 2013