A196439 a(n) = the sum of numbers k <= n such that GCQ_A(n, k) = LCQ_A(n, k) = 0 (see definition in comments).
1, 3, 3, 6, 3, 10, 3, 6, 7, 6, 3, 15, 3, 6, 7, 6, 3, 10, 3, 12, 7, 6, 3, 15, 3, 6, 7, 6, 3, 10, 3, 6, 7, 6, 3, 15, 3, 6, 7, 12, 3, 10, 3, 6, 7, 6, 3, 15, 3, 6, 7, 6, 3, 10, 3, 6, 7, 6, 3, 28, 3, 6, 7, 6, 3, 10, 3, 6, 7, 6, 3, 15, 3, 6, 7, 6, 3, 10, 3, 12, 7, 6, 3, 15, 3, 6, 7, 6, 3, 10, 3, 6, 7, 6, 3, 15, 3, 6, 7, 12, 3, 10, 3, 6, 7
Offset: 1
Keywords
Examples
For n = 6, a(6) = 10 because there are 4 cases k (k = 1, 2, 3, 4) with GCQ_A(6, k) = 0: GCQ_A(6, 1) = 0, GCQ_A(6, 2) = 0, GCQ_A(6, 3) = 0, GCQ_A(6, 4) = 0, GCQ_A(6, 5) = 4, GCQ_A(6, 6) = 5. Sum of such numbers k is 10. Also there are 4 same cases k with LCQ_A(6, k) = 0: LCQ_A(6, 1) = 0, LCQ_A(6, 2) = 0, LCQ_A(6, 3) = 0, LCQ_A(6, 4) = 0, LCQ_A(6, 5) = 4, LCQ_A(6, 6) = 4.
Links
- Antti Karttunen, Table of n, a(n) for n = 1..65537
Programs
-
PARI
GCQ_A(a, b) = { forstep(m=min(a, b)-1, 2, -1, if(a%m && b%m, return(m))); 0; }; \\ From PARI-program in A196438. A196440(n) = sum(k=1,n,(2<=GCQ_A(n,k))*k); A196439(n) = (((n*(n+1))/2) - A196440(n)); \\ Antti Karttunen, Jun 12 2018
Extensions
More terms from Antti Karttunen, Jun 12 2018
Comments