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
Antti Karttunen, Table of n, a(n) for n = 1..10000
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):
;; 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
KEYWORD
nonn,easy
AUTHOR
Alex Ratushnyak, Jan 31 2014
STATUS
approved