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.

A039766 Numbers k such that gcd(phi(k), k-1) = number of divisors of k.

This page as a plain text file.
%I A039766 #16 Mar 15 2019 21:59:38
%S A039766 1,3,21,33,57,69,77,85,93,105,125,129,141,161,175,177,201,205,209,213,
%T A039766 221,237,249,253,309,321,329,345,365,381,393,413,417,437,445,453,473,
%U A039766 475,485,489,493,497,501,517,533,537,565,573,581,597,629,633,649,665
%N A039766 Numbers k such that gcd(phi(k), k-1) = number of divisors of k.
%H A039766 Robert Israel, <a href="/A039766/b039766.txt">Table of n, a(n) for n = 1..10000</a>
%e A039766 phi(21)=12, gcd(12,20)=4, 21 is divisible by {1,3,7,21}.
%p A039766 filter:= t -> igcd(numtheory:-phi(t),t-1) = numtheory:-tau(t):
%p A039766 select(filter,[$1..1000]); # _Robert Israel_, Mar 15 2019
%o A039766 (PARI) isok(n) = gcd(eulerphi(n), n-1) == numdiv(n); \\ _Michel Marcus_, May 30 2014
%Y A039766 Cf. A000005, A000010, A049559.
%K A039766 nonn,easy
%O A039766 1,2
%A A039766 _Olivier Gérard_
%E A039766 Term 1 prepended by _Michel Marcus_, May 30 2014