OFFSET
0,4
COMMENTS
Contains as subsequences A026741, A017281, A017305, A005408, A017353, and A017377. - Luce ETIENNE, Nov 04 2018
Multiplicative and also a strong divisibility sequence: gcd(a(n),a(m)) = a(gcd(n,m)) for n, m >= 1. - Peter Bala, Feb 24 2019
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..2000
Peter Bala, A note on the sequence of numerators of a rational function, 2019.
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-1).
FORMULA
a(n) = lcm(20, n)/20. - Zerinvary Lajos, Jun 12 2009
a(n) = n/gcd(n, 20). - Andrew Howroyd, Jul 25 2018
From Luce ETIENNE, Nov 04 2018: (Start)
a(n) = 9*a(n-20) - 36*a(n-40) + 84*a(n-60) - 126*a(n-80) + 126*a(n-100) - 84*a(n-120) + 36*a(n-140) - 9*a(n-160) + a(n-180).
a(n) = (5*(119*m^9 - 4923*m^8 + 86250*m^7 - 832230*m^6 + 4807887*m^5 - 16882299*m^4 + 34770400*m^3 - 37855620m^2 + 16581744*m + 54432)*floor(n/10) + 72*m*(3*m^8 - 120*m^7 + 2030*m^6 - 18900*m^5 + 105329*m^4 - 356580*m^3 + 706220*m^2 - 733200*m + 300258) + ((19*m^9 - 855*m^8 + 15810*m^7 - 154350*m^6 + 849387*m^5 - 2597175*m^4 + 4037840*m^3 - 2600100*m^2 + 540144*m - 90720)*floor(n/10) - 72*m*(m^7 - 35*m^6 + 490*m^5 - 3500*m^4 + 13489*m^3 - 27335*m^2 + 26340*m - 9450))*(-1)^floor(n/10))/362880 where m = (n mod 10). (End)
From Peter Bala, Feb 24 2019: (Start)
a(n) = n/gcd(n,20) is a quasi-polynomial in n since gcd(n,20) is a purely periodic sequence of period 20.
O.g.f.: F(x) - F(x^2) - F(x^4) - 4*F(x^5) + 4*F(x^10) + 4*F(x^20), where F(x) = x/(1 - x)^2.
O.g.f. for reciprocals: Sum_{n >= 1} x^n/a(n) = Sum_{d divides 20} (phi(d)/d) * log(1/(1 - x^d)) = log(1/(1 - x)) + (1/2)*log(1/(1 - x^2)) + (2/4)*log(1/(1 - x^4)) + (4/5)*log(1/(1 - x^5)) + (4/10)*log(1/(1 - x^10)) + (8/20)*log(1/(1 - x^20)), where phi(n) denotes the Euler totient function A000010. (End)
From Amiram Eldar, Nov 25 2022: (Start)
Multiplicative with a(2^e) = 2^max(0, e-2), a(5^e) = 5^max(0,e-1), and a(p^e) = p^e otherwise.
Dirichlet g.f.: zeta(s-1)*(1 - 1/2^s - 1/4^s - 4/5^s + 4/10^s + 4/20^s).
Sum_{k=1..n} a(k) ~ (231/800) * n^2. (End)
MAPLE
seq(numer(n/(n+20)), n=0..80); # Muniru A Asiru, Feb 19 2019
MATHEMATICA
f[n_]:=Numerator[n/(n+20)]; Array[f, 100, 0] (* Vladimir Joseph Stephan Orlovsky, Feb 17 2011 *)
PROG
(Sage) [lcm(20, n)/20 for n in range(0, 80)] # Zerinvary Lajos, Jun 12 2009
(Magma) [Numerator(n/(n+20)): n in [0..100]]; // Vincenzo Librandi, Mar 06 2018
(PARI) a(n) = numerator(n/(n+20)); \\ Michel Marcus, Mar 07 2018
(GAP) List([0..80], n->NumeratorRat(n/(n+20))); # Muniru A Asiru, Feb 19 2019
CROSSREFS
KEYWORD
nonn,easy,frac,mult
AUTHOR
N. J. A. Sloane, May 15 2005
EXTENSIONS
Keyword:mult added by Andrew Howroyd, Jul 25 2018
STATUS
approved