login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A107816
a(n) = sigma(A067651(n)).
2
1, 8, 20, 48, 144, 240, 720, 480, 756, 720, 864, 2304, 2880, 4080, 5040, 6240, 6048, 6720, 9984, 21600, 90720, 506880, 806400, 748800, 973440, 1153152
OFFSET
1,2
COMMENTS
Numbers n such that n = sigma(k) = phi(prime(k)-1) for some k.
Values of k are in A067651, values of prime(k) are in A107815.
EXAMPLE
a(3) = sigma(A067651(3)) = sigma(19) = 20.
PROG
(PARI) m=400000; for(n=1, m, s=sigma(n); if(s==eulerphi(prime(n)-1), print1(s, ", ")))
CROSSREFS
Sequence in context: A273266 A273304 A169878 * A361073 A361215 A205219
KEYWORD
nonn,more
AUTHOR
Klaus Brockhaus, May 24 2005
STATUS
approved