cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A238229 Numbers n such that if x = sigma(n)-phi(n)+tau(n)-n then n = sigma(x)-phi(x)+tau(x)-x.

Original entry on oeis.org

2, 4, 20, 66, 342, 34092, 40842, 41922, 46242, 55422, 207624, 259448, 533172, 547300, 571992, 667408, 1531032, 1786288, 10983114, 114013752, 133506680, 323277822, 347360860, 386144360, 387415458, 459603716, 476991704, 1443279992, 1539484232, 15537978552
Offset: 1

Views

Author

Paolo P. Lava, Feb 20 2014

Keywords

Comments

The fixed points (terms with x = n) are 2, 4, 20, 66, 342, 41922, 10983114, ... - Amiram Eldar, Mar 31 2019

Examples

			Fixed points: 2, 4, 20, 66, 342, 41922, ...
Amicable pairs: (34092, 40842), (46242, 55422), (207624, 259448), ...
sigma(34092) = 86268, phi(34092) = 11352, tau(34092) = 18 and 86268 - 11352 + 18 - 34092 = 40842.
sigma(40842) = 88530, phi(40842) = 13608, tau(40842) = 12 and 88530 - 13608 + 12 - 40842 = 34092.
		

Crossrefs

Programs

  • Maple
    with(numtheory); P:=proc(q)local a,n;
    for n from 1 to q do a:=sigma(n)-phi(n)+tau(n)-n;
    if a>0 and sigma(a)-phi(a)+tau(a)-a=n then print(n);
    fi; od; end: P(10^6);
  • Mathematica
    f[n_] := If[n > 0, DivisorSigma[1, n] - EulerPhi[n] + DivisorSigma[0, n] - n, 0]; s={}; Do[ If[f[f[n]] == n, AppendTo[s, n]], {n, 1, 60000}]; s (* Amiram Eldar, Mar 31 2019 *)

Extensions

a(11)-a(29) from Amiram Eldar, Mar 31 2019
a(30) from Giovanni Resta, Apr 04 2019