|
|
A070554
|
|
Number of positive integers, k, where k <= 2n+1 and gcd(k, 2n+1) = gcd(k+1, 2n+1) = 1.
|
|
1
|
|
|
1, 1, 3, 5, 3, 9, 11, 3, 15, 17, 5, 21, 15, 9, 27, 29, 9, 15, 35, 11, 39, 41, 9, 45, 35, 15, 51, 27, 17, 57, 59, 15, 33, 65, 21, 69, 71, 15, 45, 77, 27, 81, 45, 27, 87, 55, 29, 51, 95, 27, 99, 101, 15, 105, 107, 35, 111, 63, 33, 75, 99, 39, 75, 125, 41, 129, 85, 27, 135, 137
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
LINKS
|
|
|
PROG
|
(PARI) a(n) = my(n = 2*n+1); n*prod(p=1, n, if (isprime(p) && !(n % p), (1-2/p), 1)); \\ Michel Marcus, Feb 02 2016
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|