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.

A104901 Numbers n such that sigma(n) = 8*phi(n).

This page as a plain text file.
%I A104901 #32 Dec 04 2019 18:18:51
%S A104901 42,594,744,1254,7668,8680,10788,11868,12192,14630,15642,16188,25908,
%T A104901 28458,49842,60078,70122,77142,105222,124968,125860,138460,142240,
%U A104901 165462,168402,169608,188860,201924,242316,259160,302260,553000,561906,700910,726440
%N A104901 Numbers n such that sigma(n) = 8*phi(n).
%C A104901 If p>3 and 2^p-1 is prime (a Mersenne prime) then 35*2^(p-2)*(2^p-1) is in the sequence. So 35*2^(A000043-2)*(2^A000043-1) is a subsequence of this sequence.
%C A104901 If p>2 and 2^p-1 is prime (a Mersenne prime) then 3*2^(p-2)*(2^p-1) is in the sequence (the proof is easy). - _Farideh Firoozbakht_, Dec 23 2007
%H A104901 Amiram Eldar, <a href="/A104901/b104901.txt">Table of n, a(n) for n = 1..10000</a> (calculated using data from Jud McCranie, terms 1..1000 from Donovan Johnson)
%H A104901 Kevin A. Broughan and Daniel Delbourgo, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL16/Broughan/broughan26.html">On the Ratio of the Sum of Divisors and Euler’s Totient Function I</a>, Journal of Integer Sequences, Vol. 16 (2013), Article 13.8.8.
%H A104901 Kevin A. Broughan and Qizhi Zhou, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL17/Broughan/bro32.html">On the Ratio of the Sum of Divisors and Euler's Totient Function II</a>, Journal of Integer Sequences, Vol. 17 (2014), Article 14.9.2.
%e A104901 p>3, q=2^p-1(q is prime); m=35*2^(p-2)*q so sigma(m)=48*(2^(p-1)-1)*2^p=8*(24*2^(p-3)*(2^p-2))=8*phi(m) hence m is in the sequence.
%e A104901 sigma(553000) = 1497600 = 8*187200 = 8*phi(553000) so 553000 is in the sequence but 553000 is not of the form 35*2^(p-2)*(2^p-1).
%t A104901 Do[If[DivisorSigma[1, m] == 8*EulerPhi[m], Print[m]], {m, 1000000}]
%t A104901 Select[Range[800000],DivisorSigma[1,#]==8*EulerPhi[#]&] (* _Harvey P. Dale_, Sep 12 2018 *)
%o A104901 (PARI) is(n)=sigma(n)==8*eulerphi(n) \\ _Charles R Greathouse IV_, May 09 2013
%Y A104901 Cf. A000043, A062699, A068390, A104900, A104902.
%K A104901 easy,nonn
%O A104901 1,1
%A A104901 _Farideh Firoozbakht_, Apr 01 2005