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.

A107651 Numbers n such that phi(sigma(n)) + phi(phi(n)) = n.

This page as a plain text file.
%I A107651 #10 May 15 2013 01:31:05
%S A107651 3,28,108,2352,2544,7936,13632,26736,209904,256608,1394112,2052864,
%T A107651 2169456,2490864,11942400,18884416,258072480,415272960,2064579840,
%U A107651 3737456640,3963371520,4672512000
%N A107651 Numbers n such that phi(sigma(n)) + phi(phi(n)) = n.
%C A107651 If both (3^n-1)/2 and 2*3^n-1 are prime then 48*(2*3^n-1) is in the sequence (the proof is easy). So if n is in the intersection of A028491 and A003307 then 48*(2*3^n-1) is in this sequence. Conjecture: There exist only two such terms, namely 2544 and 209904.
%C A107651 If both (3^n*31-1)/2 and 2*3^n*31-1 are prime then 48*(2*3^n*31-1) is in the sequence (the proof is easy). Conjecture: There exist only three such terms, namely 26736, 2169456, and 26376103844085843261484656.
%e A107651 18884416 is in the sequence because phi(sigma(18884416)) + phi(phi(18884416)) = 18884416.
%t A107651 Do[If[n == EulerPhi[DivisorSigma[1, n]] + EulerPhi[EulerPhi[n]], Print[n] ], {n, 10000000}]
%o A107651 (PARI) is(n)=eulerphi(sigma(n))+eulerphi(eulerphi(n))==n \\ _Charles R Greathouse IV_, Mar 06 2013
%Y A107651 Cf. A003307, A028491, A107652, A107653.
%K A107651 nonn
%O A107651 1,1
%A A107651 _Farideh Firoozbakht_, May 26 2005
%E A107651 a(17)-a(22) from _Donovan Johnson_, Mar 06 2013