OFFSET
1,1
LINKS
Scott R. Shannon, Table of n, a(n) for n = 1..10000 (1..500 from Harvey P. Dale).
EXAMPLE
250 is a term because 250 = 2*5^3 and 52 = 2^2*13 and 2+5+5+5 = 2+2+13 = 17.
MATHEMATICA
spf[n_]:=Total[Flatten[Table[#[[1]], #[[2]]]&/@FactorInteger[n]]];
spffQ[ n_]:=!PalindromeQ[n]&&spf[n]==spf[IntegerReverse[n]];
Select[Range[ 20000], spffQ] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, May 19 2017 *)
PROG
(Python)
from sympy import factorint
def sopfr(n): return sum(p*e for p, e in factorint(n).items())
def ok(n): return n != (r:=int(str(n)[::-1])) and sopfr(n) == sopfr(r)
print([k for k in range(20000) if ok(k)]) # Michael S. Branicky, Mar 05 2026
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Jason Earls, Jul 08 2003
EXTENSIONS
Corrected by T. D. Noe, Oct 25 2006
STATUS
approved
