[go: up one dir, main page]

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!)
A292303 a(1) = 1; a(n+1) = Sum_{k=1..n} lcm(a(k),n)/n. 0
1, 1, 2, 4, 4, 12, 9, 33, 50, 78, 99, 173, 264, 658, 570, 1056, 1099, 4113, 2443, 10129, 18866, 23226, 39775, 102665, 171529, 256039, 610467, 815809, 1795028, 3854202, 3044396, 10752800, 5509162, 22665306, 25847226, 66558954, 25219183, 167266731, 264535960, 163511658, 346473322, 1109093102 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,3
LINKS
FORMULA
a(1) = 1; a(n+1) = Sum_{k=1..n} a(k)/gcd(a(k),n).
EXAMPLE
a(1) = 1;
a(2) = lcm(a(1),1)/1 = lcm(1,1)/1 = 1;
a(3) = lcm(a(1),2)/2 + lcm(a(2),2)/2 = lcm(1,2)/2 + lcm(1,2)/2 = 2;
a(4) = lcm(a(1),3)/3 + lcm(a(2),3)/3 + lcm(a(3),3)/3 = lcm(1,3)/3 + lcm(1,3)/3 + lcm(2,3)/3 = 4, etc.
MATHEMATICA
a[1] = 1; a[n_] := a[n] = Sum[LCM[a[k - 1], n - 1]/(n - 1), {k, 2, n}]; Table[a[n], {n, 42}]
a[1] = 1; a[n_] := a[n] = Sum[a[k - 1]/GCD[a[k - 1], n - 1], {k, 2, n}]; Table[a[n], {n, 42}]
CROSSREFS
Sequence in context: A225232 A326824 A319210 * A000936 A319594 A065449
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Sep 14 2017
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 August 29 21:34 EDT 2024. Contains 375518 sequences. (Running on oeis4.)