cp's OEIS Frontend

From Farideh Firoozbakht, Aug 28 2006: (Start)

For n < 6, the product of the first n Fermat primes is in the sequence because if m = 2^(2^n)-1 and n < 6 then d(m) = d(phi(m)) = 2^n.

(1). If p is a Sophie Germain prime greater than 3 then m = 69615*(2p+1) (A005385) is in the sequence because d(m) = d(phi(m)) = 96. 765765, 1601145, 3271905, 4107285, 5778045, 7448805, ... is the related subsequence.

(2). If p is a prime greater than 3 such that 4p+1 is prime then m = 700245*(4p+1) (A090866) is in the sequence because d(m) = d(phi(m)) = 160. 20307105, 37112985, 104336505, 121142385, ... is the related subsequence. (End)

It is an open question whether this sequence contains infinitely many terms; see Bellaouar et al., 2023. - Allen Stenger, Feb 16 2024

Least k such that A000005(A000010(k)) = n.

From Jon E. Schoenfield, Nov 13 2016: (Start)

For every n > 0, phi(2^n) = 2^(n-1) has exactly n divisors, so a(n) <= 2^n.

For every prime p, since phi(a(p)) has exactly p divisors, phi(a(p)) must be of the form q^(p-1), where q is a prime number. If q >= 3, we would have phi(a(p)) >= 3^(p-1), and since k > phi(k) for every k > 1, we would have a(p) >= 3^(p-1)+1, which would be contradicted by the upper bound a(p) <= 2^p (see above) unless 3^(p-1)+1 <= 2^p, which is true only for p = 2. Thus, for every prime p > 2, phi(a(p)) = 2^(p-1), so a(p) > 2^(p-1). In summary, we can state that, for every prime p > 2:

(1) a(p) is the least k such that phi(k) = 2^(p-1), and

(2) 2^(p-1) < a(p) <= 2^p.

After a(36)=1333, the next few known terms are a(38)=786433, a(39)=42653, a(40)=1763, and a(42)=2993; as shown above, known bounds on a(37) and a(41) are 2^36 < a(37) <= 2^37 and 2^40 < a(41) <= 2^41.

For prime p < 37, a(p) = A001317(p-1).

Observation: for prime p < 37, a(p) is the product of distinct Fermat primes 2^(2^j)+1 for j=0..4, i.e., 3, 5, 17, 257, and 65537 (see A019434), according to the locations of the 1-bits in p-1:

. p-1 in

p a(p) prime factorization of a(p) binary

== ========== =========================== ======

2 3 = 3 1

3 5 = 5 10

5 17 = 17 100

7 85 = 17 * 5 110

11 1285 = 257 * 5 1010

13 4369 = 257 * 17 1100

17 65537 = 65537 10000

19 327685 = 65537 * 5 10010

23 5570645 = 65537 * 17 * 5 10110

29 286331153 = 65537 * 257 * 17 11100

31 1431655765 = 65537 * 257 * 17 * 5 11110

.

This pattern does not continue to p=37, since 2^(2^5)+1 is not prime. (See also A038183 and the observation there from Arkadiusz Wesolowski.) (End)

As noted, for every prime p, phi(a(p))=2^(p-1), decompose a(p) = p_1^(e_1) *...* p_m^(e_m), then phi(a(p)) = p_1^(e_1-1)*(p_1 - 1) * ... * p_m^(e_m-1)*(p_m - 1). Thus a(p) is of the form 2^e * F_(a_1) *...* F_(a_l), where F_(a_i) = 2^(a_i) + 1 denote distinct Fermat primes. If e = 0, a_1 + ... + a_l = p - 1, while if e > 0, e + a_1 + ... + a_l = p. It can be deduced that a(p) = 2^p unless p-1 can be written as a_1 + ... a_l where 2^(a_i) + 1 are distinct Fermat primes. The only Fermat primes known have a_i in {1,2,4,8,16} and it is known that 2^a + 1 is composite for 16 < a < 2^33 (cf. A019434). It follows from the fact that 1 + 2 + 4 + 8 + 16 = 31 that a(p) = 2^p for primes p with 32 < p <= 2^33. - Pjotr Buys, Sep 18 2019

First negative value is a(120) = -2.

tau(k) is the number of divisors of k (A000005); phi(k) is the Euler totient function (A000010); and sopf(k) is the sum of the distinct primes dividing k without repetition (A008472).