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”).

A236866
Positions of primes in A007775 (numbers not divisible by 2, 3 or 5).
1
2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 15, 16, 17, 18, 19, 20, 22, 23, 24, 26, 27, 28, 29, 30, 31, 34, 35, 37, 38, 40, 41, 42, 44, 45, 47, 48, 49, 51, 52, 53, 54, 57, 60, 61, 62, 63, 64, 65, 67, 69, 71, 72, 73, 74, 75, 76, 79, 82, 83, 84, 85, 89, 90, 93, 94, 95
OFFSET
1,1
COMMENTS
From Antti Karttunen, Feb 01 2014: (Start)
Positions of primes among natural numbers coprime to 30.
Term 1 is missing from the sequence, because A007775(1)=1 is not considered a prime, terms 2 - 13 are all present, and 14 is the next term missing from here, as A007775(14)=49 is the first composite in that sequence.
(End)
LINKS
PROG
(Python)
import sympy
from sympy import isprime
i=0
for n in range(1000):
if n%2 and n%3 and n%5:
i+=1 # A007775(i)=n
if isprime(n): print str(i)+', ',
(Scheme, with Antti Karttunen's IntSeq-library):
(define A236866 (MATCHING-POS 1 1 (lambda (n) (prime? (A007775 n)))))
;; Where a slow version of A007775 can be defined for example like this:
(define A007775 (MATCHING-POS 1 1 (lambda (n) (= 1 (gcd n 30)))))
;; from Antti Karttunen, Feb 01 2014
CROSSREFS
Sequence in context: A004437 A298534 A291442 * A122526 A120401 A285356
KEYWORD
nonn,easy
AUTHOR
Alex Ratushnyak, Jan 31 2014
STATUS
approved