OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
PROG
(PARI) isok(n) = (prime(n) - sigma(n) - eulerphi(n) == prime(n+1) - sigma(n+1) - eulerphi(n+1)) \\ Michel Marcus, Jul 14 2013
(PARI) list(lim) = {my(v1 = -1, k = 0, f); forprime(p = 1, lim, k++; f = factor(k); v2 = p - sigma(f) - eulerphi(f); if(v2 == v1, print1(k-1, ", ")); v1 = v2); } \\ Amiram Eldar, Mar 06 2026
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
STATUS
approved
