%I #19 Sep 08 2022 08:45:56
%S 1,4,6,16,36,55,78,105,124,144,171,200,253,325,406,465,666,689,715,
%T 741,915,930,990,1027,1081,1136,1240,1421,1448,1610,1653,1711,1752,
%U 1764,1800,1827,2211,2352,2448,2667,2800,2835,3403,3600,3619,3620,3660,3900,4840,4970,5253,5264,5513,5671,5886,6100,6328,8001,8112
%N Numbers n such that lambda(n) = lambda(n + lambda(n)).
%C Lambda is the function (A002322). If there are infinitely many Sophie Germain primes (conjecture), then this sequence is infinite. Proof: The numbers of the form p(2p+1) are in a subsequence if p and 2p+1 are both prime with p > 3, because from the property that lambda(p(2p+1)) = p(p-1), if m = p(2p+1) then lambda(m+phi(m)) = lambda (p(2p+1) + p(p-1)) = lambda(3p^2) = p(p-1) = lambda(m).
%H Vincenzo Librandi, <a href="/A188466/b188466.txt">Table of n, a(n) for n = 1..1000</a>
%e 36 is in the sequence because lambda(36) = 6, and lambda(36 + 6) = lambda(42) = 6.
%t Select[Range[20000], CarmichaelLambda[ #] == CarmichaelLambda[ # + CarmichaelLambda[#] ] &]
%o (Magma) [1] cat [n: n in [2..8140] | CarmichaelLambda(n) eq CarmichaelLambda(n+CarmichaelLambda(n))]; // _Bruno Berselli_, Apr 10 2011
%o (PARI) lambda(n) = lcm(znstar(n)[2]);
%o isok(n) = lambda(n) == lambda(n+lambda(n)); \\ _Michel Marcus_, May 12 2018
%Y Cf. A002322, A005384.
%Y Cf. A185165: Numbers n such that lambda(n)= lambda(n - lambda(n)).
%Y Cf. A051487: Numbers n such that phi(n) = phi(n - phi(n)).
%Y Cf. A108569: Numbers n such that phi(n) = phi(n + phi(n)).
%K nonn
%O 1,2
%A _Michel Lagneau_, Apr 01 2011