OFFSET
1,1
LINKS
Robert G. Wilson v, Table of n, a(n) for n = 1..100.
EXAMPLE
9 is a term since the divisors of 9 are [1, 3, 9] and 1^1+3^2+9^3 = 739, a prime.
MAPLE
isA180852 := proc(n) d := sort(convert(numtheory[divisors](n), list)) ; add( op(i, d)^i, i=1..nops(d)) ; isprime(%) ; end proc:
for n from 1 do if isA180852(n) then printf("%d, \n", n) ; end if; end do: # R. J. Mathar, Sep 23 2010
MATHEMATICA
fQ[n_] := Block[{d = Divisors@ n}, PrimeQ[ Plus @@ (d^Range@ Length@ d)]]; Select[ Range@ 4000000, fQ] (* Robert G. Wilson v, Sep 25 2010 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Sep 21 2010
EXTENSIONS
Slightly more precise definition by R. J. Mathar, Sep 23 2010
a(27) onwards from Robert G. Wilson v, Sep 25 2010
STATUS
approved