OFFSET
1,2
COMMENTS
Inspired by A134922.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
a[n_] =
If[EulerPhi[2*n - 1] - EulerPhi[2*n] == 0, n,
If[DivisorSigma[1, 2*n - 1] - DivisorSigma[1, 2*n] == 0, n,
If[Length[Divisors[2*n - 1]] - Length[Divisors[2*n]] == 0, n, {}]]];
Flatten[Table[a[n], {n, 1, 400}]]
PROG
(PARI) isok(k) = {my(f1 = factor(2*k-1), f2 = factor(2*k)); eulerphi(f1) == eulerphi(f2) || sigma(f1) == sigma(f2) || numdiv(f1) == numdiv(f2); } \\ Amiram Eldar, Jun 07 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Roger L. Bagula, Oct 17 2008
EXTENSIONS
Name edited by Amiram Eldar, Jun 07 2025
STATUS
approved
