RosettaCodeData/Task/Attractive-numbers/Maxima/attractive-numbers-2.maxima

attractivep(n):=block(ifactors(n),apply("+",map(second,%%)),if primep(%%) then true)$
sublist(makelist(i,i,120),attractivep);