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.

Showing 1-2 of 2 results.

A020491 Numbers k such that sigma_0(k) divides phi(k).

Original entry on oeis.org

1, 3, 5, 7, 8, 9, 10, 11, 13, 15, 17, 18, 19, 21, 23, 24, 26, 28, 29, 30, 31, 33, 34, 35, 37, 39, 40, 41, 43, 45, 47, 49, 51, 52, 53, 55, 56, 57, 58, 59, 61, 63, 65, 67, 69, 70, 71, 72, 73, 74, 76, 77, 78, 79, 82, 83, 84, 85, 87, 88, 89, 90, 91, 93, 95, 97, 98, 99, 101, 102, 103, 104
Offset: 1

Views

Author

Keywords

Comments

In other words, numbers k such that d(k) divides phi(k).
From Enrique Pérez Herrero, Aug 11 2010: (Start)
sigma_0(k) divides phi(k) when:
k is an odd prime: A065091;
k is an odd squarefree number: A056911;
k = 2^m, where m <> 1 is a Mersenne number (A000225).
If d divides (p-1), with p prime, then p^(d-1) is in this sequence, as are p^(p-1), p^(p-2) and p^(-1+p^n).
(End)
phi(n) and d(n) are multiplicative functions, so if m and n are coprime and both of them are in this sequence then m*n is also in this sequence. - Enrique Pérez Herrero, Sep 05 2010
From Bernard Schott, Aug 14 2020: (Start)
The corresponding quotients are in A289585.
About the 3rd case of Enrique Pérez Herrero's comment: if k = 2^M_m, where M_m = 2^m - 1 is a Mersenne number >= 3 (A000225), then the corresponding quotient phi(k)/d(k) is the integer 2^(2^m-m-2) = A076688(m); hence, these numbers k, A058891 \ {2}, form a subsequence. (End)

Crossrefs

Complement of A015733. [Enrique Pérez Herrero, Aug 11 2010]

Programs

  • Mathematica
    Select[ Range[ 105 ], IntegerQ[ EulerPhi[ # ]/DivisorSigma[ 0, # ] ]& ]
  • PARI
    isok(k) = !(eulerphi(k) % numdiv(k)); \\ Michel Marcus, Aug 10 2020

A289585 Quotients as they appear as k increases when tau(k) divides phi(k).

Original entry on oeis.org

1, 1, 2, 3, 1, 2, 1, 5, 6, 2, 8, 1, 9, 3, 11, 1, 3, 2, 14, 1, 15, 5, 4, 6, 18, 6, 2, 20, 21, 4, 23, 14, 8, 4, 26, 10, 3, 9, 7, 29, 30, 6, 12, 33, 11, 3, 35, 2, 36, 9, 6, 15, 3, 39, 10, 41, 2, 16, 14, 5, 44, 2, 18, 15, 18, 48, 7, 10, 50, 4, 51, 6, 6, 13, 53, 3, 54, 5, 18, 56, 22, 12, 24, 2
Offset: 1

Views

Author

Bernard Schott, Jul 08 2017

Keywords

Comments

Numbers k such that tau(k) divides phi(k) are in A020491.
Only for seven integers which are in A020488, we have a(n) = 1.
The integers such that a(n) = 2, 3, 4 are respectively in A062516, A063469, A063470.
When p is an odd prime then phi(p) = p-1, tau(p) = 2, so phi(p)/tau(p) = (p-1)/2 and A005097 is an infinite subsequence.
For k = A058891(m+1), that is 2^A000225(m), with m>=2, the corresponding quotient phi(k)/tau(k) is the integer A076688(m). - Bernard Schott, Aug 15 2020

Examples

			a(10) = 2 because A020491(10) = 15 and phi(15)/tau(15) = 8/4 = 2.
		

Crossrefs

Programs

  • Maple
    for n from 1 to 50 do q:=phi(n)/tau(n);
    if q=floor(q) then print(n,q,phi(n),tau(n)) else fi; od:
  • Mathematica
    f[n_] := Block[{d = EulerPhi[n]/DivisorSigma[0, n]}, If[ IntegerQ@d, d, Nothing]]; Array[f, 120] (* Robert G. Wilson v, Jul 09 2017 *)
  • PARI
    lista(nn) = {for (n=1, nn, q = eulerphi(n)/numdiv(n); if (denominator(q)==1, print1(q, ", ")););} \\ Michel Marcus, Jul 10 2017

Formula

a(n) = A000010(A020491(n)) / A000005(A020491(n)). - David A. Corneth, Jul 09 2017
Showing 1-2 of 2 results.