login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A085918
Primes p such that for some k the number of terms > 0 and < 1 in the Farey sequence of order k is p.
1
3, 5, 11, 17, 31, 41, 71, 79, 101, 127, 139, 149, 179, 199, 211, 229, 241, 269, 277, 307, 359, 383, 431, 449, 541, 773, 829, 881, 1259, 1307, 1327, 1493, 1831, 1933, 2141, 2551, 3373, 3947, 4127, 4831, 4957, 5021, 5153, 5323, 5431, 5569, 5813, 6091, 6329
OFFSET
1,1
COMMENTS
Or, prime numbers of the form Sum_{j=2..k} phi(j). - Jorge Coveiro, Dec 22 2004. Examples: phi(2)+phi(3) = 3; phi(2)+phi(3)+phi(4) = 5; phi(2)+phi(3)+phi(4)+phi(5)+phi(6) = 11; phi(2)+phi(3)+phi(4)+phi(5)+phi(6)+phi(7) = 17.
Does this sequence have an infinite number of terms?
LINKS
EXAMPLE
The Farey sequence of order 4 is {0, 1/4, 1/3, 1/2, 2/3, 3/4, 1}. The number of terms > 0 and < 1 is 5, which is prime, so 5 is a term.
MATHEMATICA
Select[Accumulate[Table[EulerPhi[k], {k, 2, 150}]], PrimeQ] (* Amiram Eldar, Jul 06 2024 *)
PROG
(PARI) /* Farey sequence of order n */ fareycountp(n) = { for(x=2, n, y = farey(x); if(isprime(y), print1(y", ")); ) }
farey(n) = { c=1; m=n*(n-2)+2; a=vector(m); for(x=1, n, for(y=x, n, v = x/y; if(v<1, c++; a[c]=v; ) ) ); a = vecsort(a); c=0; for(x=2, m, if(a[x]<>a[x-1] & a[x]<>0, \ print1(a[x]", "); c++; ) ); return(c) }
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Cino Hilliard, Aug 16 2003
EXTENSIONS
Definition corrected by Jonathan Sondow, Apr 21 2005
STATUS
approved