OFFSET
1,1
COMMENTS
Subsequence of A002808. - Altug Alkan, Oct 07 2015
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..108 (terms < 10^13)
EXAMPLE
For n = 1372 we have phi(n) = 588, n' = 1960 and 1960 = 1372 + 588.
MAPLE
with(numtheory); P:= proc(q) local a1, a2, n, p;
for n from 1 to q do a1:=n*add(op(2, p)/op(1, p), p=ifactors(n)[2]);
if a1=n+phi(n) then print(n); fi; od; end: P(10^6);
PROG
(PARI) for(n=2, 10^10, if((k = n + eulerphi(n)) && (d(n) = local(fac); if(n<1, 0, fac=factor(n); sum(i=1, matsize(fac)[1], n*fac[i, 2]/fac[i, 1]))) && k==d(n), print1(n", "))) \\ Altug Alkan, Oct 06 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Oct 25 2013
EXTENSIONS
a(21)-a(33) from Giovanni Resta, Oct 25 2013
STATUS
approved