OFFSET
1,1
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 A248501.
The asymptotic density of this sequence is 1 - A250031(16)/A250033(16) = 199663/480480 = 0.415549... . - 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
2 is a term because gcd(2,0) = 2 > 1.
21 is not a term because floor(21/16) = 1 and 1 is coprime to any number.
200 is a term because floor(200/16) = 12 and gcd(200,12) = 4 > 1.
MATHEMATICA
Select[Range[150], !CoprimeQ[#, Floor[#/16]] &] (* Amiram Eldar, Nov 30 2024 *)
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