OFFSET
1,2
COMMENTS
Definition of 'being coprime' and special-case conventions are as in Wikipedia. In particular, when m < 16 then floor(m/16) = 0, and zero is coprime only to 1. The complementary sequence is A248502.
The asymptotic density of this sequence is A250031(16)/A250033(16) = 280817/480480 = 0.58445... . - Amiram Eldar, Nov 30 2024
LINKS
Stanislav Sykora, Table of n, a(n) for n = 1..20000
Wikipedia, Coprime integers.
FORMULA
gcd(a(n),floor(a(n)/16)) = 1.
EXAMPLE
1 is a term because gcd(1,0) = 1.
2 is not a term because gcd(2,0) = 2.
129 is a term because 129 is coprime to floor(129/16) = 8.
MATHEMATICA
Select[Range[120], CoprimeQ[#, Floor[#/16]]&] (* Harvey P. Dale, Mar 12 2023 *)
PROG
(PARI) a=vector(20000);
i=n=0; while(i++, if(gcd(i, i\16)==1, a[n++]=i; if(n==#a, break))); a
CROSSREFS
KEYWORD
nonn,base,easy
AUTHOR
Stanislav Sykora, Oct 07 2014
STATUS
approved