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.

A095952 Initial values for f(x)=phi(sigma(x)) such that iteration of f ends in cycle of length=1.

This page as a plain text file.
%I A095952 #7 Mar 28 2015 22:44:27
%S A095952 1,2,3,5,8,9,12,14,15,19,20,22,23,26,29,41,43,128,156,168,186,189,200,
%T A095952 201,217,231,237,240,248,254,260,266,271,285,291,297,303,304,313,314,
%U A095952 329,332,335,337,341,346,350,366,368,383,387,395,413,427,430,436,437
%N A095952 Initial values for f(x)=phi(sigma(x)) such that iteration of f ends in cycle of length=1.
%H A095952 Charles R Greathouse IV, <a href="/A095952/b095952.txt">Table of n, a(n) for n = 1..10000</a>
%e A095952 n=100: trajectory={980, 648, 880, 720, 720}
%t A095952 g[n_] := EulerPhi[ DivisorSigma[1, n]]; f[n_] := f[n] = Block[{lst = NestWhileList[g, n, UnsameQ, All ]}, -Subtract @@ Flatten[ Position[lst, lst[[ -1]]]]]; Select[ Range[1540], f[ # ] == 1 &] (* _Robert G. Wilson v_, Jul 14 2004 *)
%o A095952 (PARI) is(n)=my(t=f(n),h=f(t));while(t!=h,t=f(t);h=f(f(h)));t==f(t) \\ _Charles R Greathouse IV_, Nov 22 2013
%Y A095952 Cf. A000010, A000203, A096887, A095953, A096526, A095954, A096888, A096889, A096890, A095955, A095956.
%K A095952 nonn
%O A095952 1,2
%A A095952 _Labos Elemer_, Jul 13 2004