OFFSET
1,1
EXAMPLE
a(18) = 100 because the divisors of 100 are {1, 2, 4, 5, 10, 20, 25, 50, 100} with lambda(2+4+10+20+50+100) = lambda(186) = 30 and phi(1+5+25) = phi(31) = 30.
MAPLE
with(numtheory):nn:=800:
for n from 2 by 2 to nn do:
d:=divisors(n):n0:=nops(d):s0:=0:s1:=0:
for i from 1 to n0 do:
if irem(d[i], 2)=0
then
s0:=s0+d[i] else
s1:=s1+d[i]:
fi:
od:
if lambda(s0)=phi(s1) then
printf(`%d, `, n):else fi:
od:
MATHEMATICA
Select[Range[2, 1000, 2], EulerPhi[DivisorSigma[1, #/2^IntegerExponent[#, 2]]] == CarmichaelLambda[2*DivisorSigma[1, #/2]] &] (* Amiram Eldar, Jul 23 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Jul 23 2024
STATUS
approved