login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A333849 a(n) = gcd(A333848(n), 2*(2*n+1)), for n >= 0. 3
2, 1, 2, 1, 1, 1, 2, 2, 2, 1, 6, 1, 2, 1, 2, 1, 2, 2, 2, 6, 2, 1, 2, 1, 1, 2, 2, 10, 6, 1, 2, 6, 2, 1, 2, 1, 2, 2, 2, 1, 1, 1, 2, 2, 2, 2, 6, 2, 2, 2, 2, 1, 6, 1, 2, 6, 2, 2, 6, 2, 1, 2, 2, 1, 6, 1, 2, 2, 2, 1, 2, 2, 2, 6, 2, 1, 2, 10, 2, 2, 2, 1, 10, 1, 2, 18, 2, 2, 2, 1, 2 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,1
COMMENTS
For n >= 1, a(n) enters the formula for the length L(2*n+1) = A332441(n) of the directed Euler tour ET(2*n+1, q0 = 1) based on the unsigned Schick sequence for 2*n+1, namely L(2*n+1) = A003558(n)*2*(2*n+1)/a(n). For Schick sequences and references see A332439.
LINKS
Wolfdieter Lang, On the Equivalence of Three Complete Cyclic Systems of Integers, arXiv:2008.04300 [math.NT], 2020.
FORMULA
a(n) = gcd(A333848(n), 2*(2*n+1)), for n >= 0.
MATHEMATICA
{2}~Join~Table[GCD[Total@ Select[Range[1, m, 2], GCD[#, m] == 1 &], 2 m], {m, Array[2 # + 1 &, 85]}] (* Michael De Vlieger, Oct 15 2020 *)
PROG
(PARI) f(n) = if (n==0, 0, my(m=2*n+1); vecsum(select(x->((gcd(m, x)==1) && (x%2)), [1..m]))); \\ A333848
a(n) = gcd(f(n), 2*(2*n+1)); \\ Michel Marcus, May 05 2020
CROSSREFS
Sequence in context: A081652 A302790 A363455 * A140816 A029433 A337144
KEYWORD
nonn,easy
AUTHOR
Wolfdieter Lang, May 01 2020
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 20:15 EDT 2024. Contains 371281 sequences. (Running on oeis4.)