OFFSET
0,3
COMMENTS
The smallest nonnegative reduced residue system modulo N is the ordered set RRS(N) (written as a list) with integers k from {0, 1, ..., N-1} satisfying gcd(k, N) = 1, for N >= 1. See A038566 (with A038566(1) = 0).
If only odd members of RRS(N) are considered, name this list RRSodd(N), e.g., RRSodd(1) = [], the empty list, RRSodd(2) = [1], etc. See A216319 (but there A216319(1) = 1). The number of elements of RRSodd(N) is delta(N) = A055034(N), for N >= 2, and 0 for N = 1.
Here only numbers N = 2*n + 1 >= 1 are considered, and for the empty list RRSodd(1) a(0) is set to 0.
a(n) gives for n >= 1 also the sum of the numbers of the primitive period of the unsigned Schick sequences SBB(2*n+1, q0 = 1) (BB for Brändli and Beyne), for which 2*n + 1 satisfies A135303(n) = 1 (in Schick's notation B(2*n+1) = 1, implying initial value q0 = 1). The numbers n satisfying A135303(n) = 1 are given in A333854.
REFERENCES
Carl Schick, Trigonometrie und unterhaltsame Zahlentheorie, Bokos Druck, Zürich, 2003 (ISBN 3-9522917-0-6). Tables 3.1 to 3.10, for odd p = 3..113 (with gaps), pp. 158-166.
LINKS
Michael De Vlieger, Table of n, a(n) for n = 0..10000
Gerold Brändli and Tim Beyne, Modified Congruence Modulo n with Half the Amount of Residues, arXiv:1504.02757 [math.NT], 2016.
Wolfdieter Lang, On the Equivalence of Three Complete Cyclic Systems of Integers, arXiv:2008.04300 [math.NT], 2020.
FORMULA
a(n) = Sum_{j=1..delta(2*n+1)} RRSodd(2*n+1)_j, for n >= 1, with delta(k) = A055034(k). a(0) = 0 (undefined case).
EXAMPLE
n = 4: RRSodd(9) = {1, 5, 7} with sum a(4) = 13. Schick's unsigned cycle is SBB(9, 1) = (1, 7, 5). Because A135303(4) = B(9) = 1 there is only this cycle for n = 9.
MATHEMATICA
{0}~Join~Table[Total@ Select[Range[1, m, 2], GCD[#, m] == 1 &], {m, Array[2 # + 1 &, 50]}] (* Michael De Vlieger, Oct 15 2020 *)
PROG
(PARI) a(n) = if (n==0, 0, my(m=2*n+1); vecsum(select(x->((gcd(m, x)==1) && (x%2)), [1..m]))); \\ Michel Marcus, May 05 2020
(PARI) apply( {A333848(n)=vecsum([2*m-1|m<-[1..n], gcd(m*2-1, n*2+1)==1])}, [0..50]) \\ M. F. Hasler, Jun 04 2020
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Wolfdieter Lang, May 01 2020
STATUS
approved