cp's OEIS Frontend

All terms are primes. Note a connection to the Wieferich primes A001220: a(2) = a(3) = A001220(1), a(3) = a(4) = A001220(2).

From John Blythe Dobson, Jan 14 2017: (Start)

All Wieferich primes p will belong to this sequence twice, because if H([p/k]) denotes the harmonic number with index floor(p/k), then p divides all of H([p/4]), H([p/2]), and H(p-1). The first two of these elements gives one solution, and the second and third another. This property of the Wieferich primes predates their name, and was apparently first proved by Glaisher in "On the residues of r^(p-1) to modulus p^2, p^3, etc.," pp. 21-22, 23 (see References).

Note also a connection to the Mirimanoff primes A014127: a(1) = A014127(1), a(8) = A014127(2). All Mirimanoff primes p will belong to this sequence, because p divides both H([p/3]) and H([2p/3]). This property of the Mirimanoff primes likewise predates their name, and was apparently first proved by Glaisher in "A general congruence theorem relating to the Bernoullian function," p. 50 (see Links).

The Wieferich primes and Mirimanoff primes would seem to be the only cases for which the value of n in A126196(n) is predictable from knowledge of p. It is not obvious that all members of the present sequence are prime; however, by definition all their divisors must be non-harmonic primes A092102. Furthermore, it is clear from the cited literature under that entry that H([n/2]) == H(n) == 0 (mod p) is only possible when n < p. Thus, all divisors of the present sequence must belong to the harmonic irregular primes A092194.

One possible reason for interest in this sequence is a 1995 result of Dilcher and Skula (see Links) which among other things shows that if a prime p were an exception to the first case of Fermat's Last Theorem, then p would divide both H([p/k]) and H([2p/k]) for every value of k from 2 to 46. To date, the only values for which such coincidences have been found have k = 2, 3, or 4. For k = 6 to hold, p would have to be simultaneously a Wieferich prime and a Mirimanoff prime, while for k = 5 to hold, p would have to be simultaneously a Wall-Sun-Sun prime and a member of A123692. The sparse numerical results for the present sequence suggest that even the more relaxed condition H([n/2]) == H(n) == 0 (mod p) is rarely satisfied. (End)

For p = prime(n), Boyd defines J_p to be the set of numbers k such that p divides A001008(k), the numerator of the harmonic number H(k). For harmonic primes, J_p contains only the three numbers p-1, (p-1)p and (p-1)(p+1). It has been conjectured that there are an infinite number of these primes and that their density in the primes is 1/e.

Prime p=A000040(n) is in this sequence iff neither H(k) == 0 (mod p), nor H(k) == -A177783(n) (mod p) have solutions for 1 <= k <= p-2. - Max Alekseyev, May 13 2010

a(n)<=n for n =2,5,14,18,25,29,33,46,49,...

For p = prime(n), Boyd defines J_p to be the set of numbers k such that p divides A001008(k). This sequence gives the smallest elements of J_p. The largest elements of J_p are given by A177734. The sizes of J_p are given by A092103.

For p = prime(n), Boyd defines J_p to be the set of numbers k such that p divides A001008(k). This sequence gives the size of J_p. A072984 and A177734 give the smallest and largest elements of J_p, respectively.

A092101 gives primes prime(n) such that a(n) = 3 (i.e., a(A000720(A092101(m))) = 3 for all m). A092102 gives primes prime(n) such that a(n) > 3.

From Carlo Sanna, Apr 06 2016: (Start)

Eswarathasan and Levine conjectured that for any prime number p the set J_p is finite.

I proved that if J_p(x) is the number of integers in J_p that are less than x > 1, then J_p(x) < 129 p^(2/3) x^0.765 for any prime p. In particular, J_p has asymptotic density zero. (End)

Bing-Ling Wu and Yong-Gao Chen improved Sanna's (see previous comment) result showing that J_p(x) <= 3 x^(2/3 + 1/(25 log p)) for any prime p and any x > 1. - Carlo Sanna, Jan 12 2017

For any prime p, generation m consists of the numbers p^(m-1) <= k < p^m. The zeroth generation consists of just the number 0. When there is a k in generation m such that p divides A001008(k), then that k may generate solutions in generation m+1. It is conjectured that for all primes there are solutions for only a finite number of generations.

Boyd's table 3 states incorrectly that harmonic primes have 2 generations; harmonic primes have 3 generations.

These primes are a subset of the non-harmonic primes A092102. Because these primes are analogous to the irregular primes A000928 that divide the numerators of Bernoulli numbers, they might be called H-irregular primes. The density of these primes is about 0.4 -- very close to the density of irregular primes.

These primes are called Harmonic irregular primes in the Wikipedia entry for "Regular prime" (see links). It may be noted that if p is known to be of this type and H(k) is the smallest Harmonic number divisible by p, then not only does k < p-1 hold, but k <= (p-1)/2. This is because, by symmetry, H(p-1-n) == H(n) (mod p), so that any eligible k lying between (p+1)/2 and p-1 would have a counterpart in the range between 1 and (p-1)/2. Furthermore, the minimal k cannot be exactly equal to (p-1)/2, because then p would be a Wieferich prime (A001220) and would also divide H(Int(p/4)). Thus k <= (p-3)/2, and this inequality is sharp because exact equality holds for p = 29, 37, 3373 (see A072984). - John Blythe Dobson, Apr 09 2015