A026741 -id:A026741 - cp's OEIS frontend

A106614 a(n) = numerator of n/(n+13).

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 2, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 3, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 4, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 5, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75

Offset: 0

N. J. A. Sloane, May 15 2005

In general, the numerators of n/(n+p) for prime p and n >= 0, form a sequence with the g.f.: x/(1-x)^2 - (p-1)*x^p/(1-x^p)^2. - Paul D. Hanna, Jul 27 2005

a(n) <> n iff n = 13 * k, in this case, a(n) = k. - Bernard Schott, Feb 19 2019

G. C. Greubel, Table of n, a(n) for n = 0..10000
Index entries for linear recurrences with constant coefficients, signature (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,-1).

Cf. Sequences given by the formula numerator(n/(n + k)): A026741 (k = 2), A051176 (k = 3), A060819 (k = 4), A060791 (k = 5), A060789 (k = 6), A106608 thru A106612 (k = 7 thru 11), A051724 (k = 12), A106615 thru A106621 (k = 14 thru 20).

List([0..80],n->NumeratorRat(n/(n+13))); # Muniru A Asiru, Feb 19 2019

[Numerator(n/(n+13)): n in [0..100]]; // Vincenzo Librandi, Apr 18 2011

seq(numer(n/(n+13)),n=0..80); # Muniru A Asiru, Feb 19 2019

f[n_]:=Numerator[n/(n+13)];Array[f,100,0] (* Vladimir Joseph Stephan Orlovsky, Feb 17 2011 *)

vector(100, n, n--; numerator(n/(n+13))) \\ G. C. Greubel, Feb 19 2019

a(n)=if(n%13,n,n/13) \\ Charles R Greathouse IV, Jan 24 2022

[lcm(n,13)/13for n in range(0, 100)] # Zerinvary Lajos, Jun 09 2009

G.f.: x/(1-x)^2 - 12*x^13/(1-x^13)^2. - Paul D. Hanna, Jul 27 2005
Dirichlet g.f.: zeta(s-1)*(1-12/13^s). - R. J. Mathar, Apr 18 2011
a(n) = 2*a(n-13) - a(n-26). - G. C. Greubel, Feb 19 2019
From Amiram Eldar, Nov 25 2022: (Start)
Multiplicative with a(13^e) = 13^(e-1), and a(p^e) = p^e if p != 13.
Sum_{k=1..n} a(k) ~ (157/338) * n^2. (End)
Sum_{n>=1} (-1)^(n+1)/a(n) = 25*log(2)/13. - Amiram Eldar, Sep 08 2023

A162610 Triangle read by rows in which row n lists n terms, starting with 2n-1, with gaps = n-1 between successive terms.

Original entry on oeis.org

1, 3, 4, 5, 7, 9, 7, 10, 13, 16, 9, 13, 17, 21, 25, 11, 16, 21, 26, 31, 36, 13, 19, 25, 31, 37, 43, 49, 15, 22, 29, 36, 43, 50, 57, 64, 17, 25, 33, 41, 49, 57, 65, 73, 81, 19, 28, 37, 46, 55, 64, 73, 82, 91, 100, 21, 31, 41, 51, 61, 71, 81, 91, 101, 111, 121

Offset: 1

Omar E. Pol, Jul 09 2009

Note that the last term of the n-th row is the n-th square A000290(n).

Row sums are n*(n^2+2*n-1)/2, apparently in A127736. - R. J. Mathar, Jul 20 2009

			Triangle begins:
1
3, 4
5, 7, 9
7, 10, 13, 16
9, 13, 17, 21, 25
11, 16, 21, 26, 31, 36

Harvey P. Dale, Table of n, a(n) for n = 1..10000

Cf. A000027, A000290, A159797, A159798.

Cf. A209297; A005408 (left edge), A000290 (right edge), A127736 (row sums), A056220 (central terms), A026741 (number of odd terms per row), A142150 (number of even terms per row), A221491 (number of primes per row).

a162610 n k = k * n - k + n
a162610_row n = map (a162610 n) [1..n]
a162610_tabl = map a162610_row [1..]
-- Reinhard Zumkeller, Jan 19 2013

Flatten[Table[NestList[#+n-1&,2n-1,n-1], {n,15}]] (* Harvey P. Dale, Oct 20 2011 *)

# From R. J. Mathar, Oct 20 2009
def A162610(n, k):
    return 2*n-1+(k-1)*(n-1)
print([A162610(n,k) for n in range(1,20) for k in range(1,n+1)])

T(n,k) = n+k*n-k, 1<=k<=n. - R. J. Mathar, Oct 20 2009

T(n,k) = (k+1)*(n-1)+1. - Reinhard Zumkeller, Jan 19 2013

More terms from R. J. Mathar, Oct 20 2009

A051715 Denominators of table a(n,k) read by antidiagonals: a(0,k) = 1/(k+1), a(n+1,k) = (k+1)(a(n,k)-a(n,k+1)), n >= 0, k >= 0.

Original entry on oeis.org

1, 2, 2, 3, 3, 6, 4, 4, 6, 1, 5, 5, 20, 30, 30, 6, 6, 15, 20, 30, 1, 7, 7, 42, 35, 140, 42, 42, 8, 8, 28, 84, 105, 28, 42, 1, 9, 9, 72, 84, 1, 105, 140, 30, 30, 10, 10, 45, 120, 140, 28, 105, 20, 30, 1, 11, 11, 110, 495, 3960, 924, 231, 165, 220, 66, 66, 12, 12, 66, 55, 495, 264, 308, 132, 165, 44, 66, 1

Offset: 0

N. J. A. Sloane

Leading column gives the Bernoulli numbers A027641/A027642.

			Table begins:
    1    1/2   1/3    1/4   1/5  1/6  1/7 ...
   1/2   1/3   1/4    1/5   1/6  1/7 ...
   1/6   1/6   3/20   2/15  5/42 ...
    0    1/30  1/20   2/35  5/84 ...
  -1/30 -1/30 -3/140 -1/105 ...

Alois P. Heinz, Antidiagonals n = 0..140, flattened
M. Kaneko, The Akiyama-Tanigawa algorithm for Bernoulli numbers, J. Integer Sequences, 3 (2000), #00.2.9.
Index entries for sequences related to Bernoulli numbers.

Rows 2, 3, 4 give A026741/A045896, A051712/A051713, A051722/A051723.
Columns 0, 1, 2, 3 give A000367/A002445, A051716/A051717, A051718/A051719, A051720/A051721.
Numerators are in A051714.

a:= proc(n,k) option remember;
      `if`(n=0, 1/(k+1), (k+1)*(a(n-1,k)-a(n-1,k+1)))
    end:
seq(seq(denom(a(n, d-n)), n=0..d), d=0..12); # Alois P. Heinz, Apr 17 2013

nmax = 12; a[0, k_] := 1/(k+1); a[n_, k_] := a[n, k] = (k+1)(a[n-1, k]-a[n-1, k+1]); Denominator[ Flatten[ Table[ a[n-k, k], {n, 0, nmax}, {k, n, 0, -1}]]](* Jean-François Alcover, Nov 28 2011 *)

a(n,k) = denominator(Sum_{j=0..n} (-1)^(n-j)*j!*Stirling2(n,j)/(j+k+1)). - Fabián Pereyra, Jan 14 2023

More terms from James Sellers, Dec 08 1999

A106621 a(n) = numerator of n/(n+20).

Original entry on oeis.org

0, 1, 1, 3, 1, 1, 3, 7, 2, 9, 1, 11, 3, 13, 7, 3, 4, 17, 9, 19, 1, 21, 11, 23, 6, 5, 13, 27, 7, 29, 3, 31, 8, 33, 17, 7, 9, 37, 19, 39, 2, 41, 21, 43, 11, 9, 23, 47, 12, 49, 5, 51, 13, 53, 27, 11, 14, 57, 29, 59, 3, 61, 31, 63, 16, 13, 33, 67, 17, 69, 7, 71, 18, 73, 37, 15, 19, 77, 39, 79

Offset: 0

N. J. A. Sloane, May 15 2005

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

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).

Cf. A000010, A010879.

List([0..80],n->NumeratorRat(n/(n+20))); # Muniru A Asiru, Feb 19 2019

[Numerator(n/(n+20)): n in [0..100]]; // Vincenzo Librandi, Mar 06 2018

seq(numer(n/(n+20)),n=0..80); # Muniru A Asiru, Feb 19 2019

f[n_]:=Numerator[n/(n+20)];Array[f,100,0] (* Vladimir Joseph Stephan Orlovsky, Feb 17 2011 *)

a(n) = numerator(n/(n+20)); \\ Michel Marcus, Mar 07 2018

[lcm(20,n)/20 for n in range(0, 80)] # Zerinvary Lajos, Jun 12 2009

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)

Keyword:mult added by Andrew Howroyd, Jul 25 2018

A195817 Multiples of 10 and odd numbers interleaved.

Original entry on oeis.org

0, 1, 10, 3, 20, 5, 30, 7, 40, 9, 50, 11, 60, 13, 70, 15, 80, 17, 90, 19, 100, 21, 110, 23, 120, 25, 130, 27, 140, 29, 150, 31, 160, 33, 170, 35, 180, 37, 190, 39, 200, 41, 210, 43, 220, 45, 230, 47, 240, 49, 250, 51, 260, 53, 270, 55, 280, 57, 290, 59, 300

Offset: 0

Omar E. Pol, Sep 29 2011

A008592 and A005408 interleaved.
Partial sums give the generalized 14-gonal (or tetradecagonal) numbers A195818.
a(n) is also the length of the n-th line segment of a rectangular spiral on the infinite square grid. The vertices of the spiral are the generalized 14-gonal numbers. - Omar E. Pol, Jul 27 2018

Vincenzo Librandi, Table of n, a(n) for n = 0..10000
Index entries for linear recurrences with constant coefficients, signature (0,2,0,-1).

Column 10 of A195151.
Sequences whose partial sums give the generalized n-gonal numbers, if n>=5: A026741, A001477, zero together with A080512, A022998, A195140, zero together with A165998, A195159, A195161, A195312, this sequence.
Cf. A091998, A195152.

[(2*(-1)^n+3)*n: n in [0..60]]; // Vincenzo Librandi, Sep 30 2011

With[{nn=30},Riffle[10*Range[0,nn],Range[1,2*nn+1,2]]] (* or *) LinearRecurrence[{0,2,0,-1},{0,1,10,3},70] (* Harvey P. Dale, Nov 24 2013 *)

a(n) = (2*(-1)^n+3)*n; \\ Andrew Howroyd, Jul 23 2018

a(n) = (2*(-1)^n+3)*n. - Vincenzo Librandi, Sep 30 2011
From Bruno Berselli, Sep 30 2011: (Start)
G.f.: x*(1+10*x+x^2)/((1-x)^2*(1+x)^2).
a(n) = -a(-n) = a(n-2)*n/(n-2) = 2*a(n-2)-a(n-4).
a(n) * a(n+1) = a(n(n+1)).
a(n) + a(n+1) = A091998(n+1). (End)
a(0)=0, a(1)=1, a(2)=10, a(3)=3, a(n)=2*a(n-2)-a(n-4). - Harvey P. Dale, Nov 24 2013
Multiplicative with a(2^e) = 5*2^e, a(p^e) = p^e for odd prime p. - Andrew Howroyd, Jul 23 2018
Dirichlet g.f.: zeta(s-1) * (1 + 2^(3-s)). - Amiram Eldar, Oct 25 2023

A005013 a(n) = 3*a(n-2) - a(n-4), a(0)=0, a(1)=1, a(2)=1, a(3)=4. Alternates Fibonacci (A000045) and Lucas (A000032) sequences for even and odd n.

Original entry on oeis.org

0, 1, 1, 4, 3, 11, 8, 29, 21, 76, 55, 199, 144, 521, 377, 1364, 987, 3571, 2584, 9349, 6765, 24476, 17711, 64079, 46368, 167761, 121393, 439204, 317811, 1149851, 832040, 3010349, 2178309, 7881196, 5702887, 20633239, 14930352, 54018521, 39088169, 141422324

Offset: 0

N. J. A. Sloane

S(n,sqrt(5)), with the Chebyshev polynomials A049310, is an integer sequence in the real quadratic number field Q(sqrt(5)) with basis numbers <1,phi>, phi:=(1+sqrt(5))/2. S(n,sqrt(5)) = A(n) + 2*B(n)*phi, with A(n)= a(n+1)*(-1)^n and B(n)= A147600(n-1), n>=0, with A147600(-1):=0.
a(n) = p(n+1) where p(x) is the unique degree-(n-1) polynomial such that p(k) = Fibonacci(k) for k = 1, ..., n. - Michael Somos, Jan 08 2012
Row sums of A227431. - Richard R. Forberg, Jul 29 2013
This is the sequence of Lehmer numbers u_n(sqrt(R),Q) with the parameters R = 5 and Q = 1. It is a strong divisibility sequence, that is, gcd(a(n), a(m)) = a(gcd(n,m)) for all natural numbers n and m. The sequence satisfies a linear recurrence of order four. - Peter Bala, Apr 18 2014
The sequence of convergents of the 2-periodic continued fraction [0; 1, -5, 1, -5, ...] = 1/(1 - 1/(5 - 1/(1 - 1/(5 - ...)))) = (1/2)*(5 - sqrt(5)) begins [0/1, 1/1, 5/4, 4/3, 15/11, 11/8, 40/29, ...]; the denominators give the present sequence. The sequence of numerators [0, 1, 5, 4, 15, 11, 40, ...] is A203976. Cf. A108412 and A026741. - Peter Bala, May 19 2014
Define a binary operation o on the real numbers by x o y = x*sqrt(1 + y^2) + y*sqrt(1 + x^2). The operation o is commutative and associative with identity 0. We have (1/2)*a(2*n + 1) = 1/2 o 1/2 o ... o 1/2 (2*n + 1 terms) and (1/2)*sqrt(5)* a(2*n) = 1/2 o 1/2 o ... o 1/2 (2*n terms). Cf. A084068 and A049629. - Peter Bala, Mar 23 2018

			G.f. = x + x^2 + 4*x^3 + 3*x^4 + 11*x^5 + 8*x^6 + 29*x^7 + 21*x^8 + 76*x^9 + ...
a(3) = 4 since p(x) = (x^2 - 3*x + 4) / 2 interpolates p(1) = 1, p(2) = 1, p(3) = 2, and p(4) = 4. - _Michael Somos_, Jan 08 2012

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

G. C. Greubel, Table of n, a(n) for n = 0..2000 (terms 0..500 from T. D. Noe)
A. F. Horadam, R. P. Loh and A. G. Shannon, Divisibility properties of some Fibonacci-type sequences, pp. 55-64 of Combinatorial Mathematics VI (Armidale 1978), Lect. Notes Math. 748, 1979. [Annotated scanned copy]
Seong Ju Kim, R. Stees and L. Taalman, Sequences of Spiral Knot Determinants, Journal of Integer Sequences, Vol. 19 (2016), #16.1.4.
D. H. Lehmer, An extended theory of Lucas' functions, Annals of Mathematics, Second Series, Vol. 31, No. 3 (Jul., 1930), pp. 419-448.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992.
E. L. Roettger and H. C. Williams, Appearance of Primes in Fourth-Order Odd Divisibility Sequences, J. Int. Seq., Vol. 24 (2021), Article 21.7.5.
E. W. Weisstein, MathWorld: Lehmer Number
H. C. Williams and R. K. Guy, Odd and even linear divisibility sequences of order 4, INTEGERS, 2015, #A33.
Index to divisibility sequences
Index entries for linear recurrences with constant coefficients, signature (0,3,0,-1).

Cf. A000032, A000045, A001906, A002878, A005247, A096140, A026741, A108412, A084068, A049629.

a:=[0,1,1,4];; for n in [5..40] do a[n]:=3*a[n-2]-a[n-4]; od; a; # Muniru A Asiru, Oct 21 2018

a005013 n = a005013_list !! n
a005013_list = alt a000045_list a000032_list where
   alt (f::fs) (:l:ls) = f : l : alt fs ls
-- Reinhard Zumkeller, Jan 10 2012

I:=[0,1,1,4]; [n le 4 select I[n]  else 3*Self(n-2) - Self(n-4): n in [1..40]]; // Vincenzo Librandi, Feb 09 2016

with(combinat): A005013 := n-> if n mod 2 = 0 then fibonacci(n) else fibonacci(n+1)+fibonacci(n-1); fi;
A005013:=z*(z**2+z+1)/((z**2+z-1)*(z**2-z-1)); # Simon Plouffe in his 1992 dissertation

CoefficientList[Series[(x + x^2 + x^3)/(1 - 3x^2 + x^4), {x, 0, 40}], x]
f[n_] = Product[(1 + 4*Sin[k*Pi/n]^2), {k, 1, Floor[(n - 1)/2]}]; a = Table[f[n], {n, 0, 30}]; Round[a]; FullSimplify[ExpandAll[a]] (* Roger L. Bagula and Gary W. Adamson, Nov 26 2008 *)
LinearRecurrence[{0, 3, 0, -1}, {0, 1, 1, 4}, 100] (* G. C. Greubel, Feb 08 2016 *)

{a(n) = if( n%2, fibonacci(n+1) + fibonacci(n-1), fibonacci(n))}; /* Michael Somos, Jan 08 2012 */

{a(n) = if( n<0, -a(-n), subst( polinterpolate( vector( n, k, fibonacci(k))), x, n+1))}; /* Michael Somos, Jan 08 2012 */

a(1) = a(2) = 1, a(3) = 4, a(n) = (a(n-1) * a(n-2) - 1) / a(n-3), unless n=3. a(-n) = -a(n).
a(2n) = A001906(n), a(2n+1) = A002878(n). a(n)=F(n+1)+(-1)^(n+1)F(n-1). - Mario Catalani (mario.catalani(AT)unito.it), Sep 20 2002
G.f.: x*(1+x+x^2)/((1-x-x^2)*(1+x-x^2)).
a(n) = Product_{k=1..floor((n-1)/2)} (1 + 4*sin(k*Pi/n)^2). - Roger L. Bagula and Gary W. Adamson, Nov 26 2008
Binomial transform is A096140. - Michael Somos, Apr 13 2012
From Peter Bala, Apr 18 2014: (Start)
a(n) = (alpha^n - beta^n)/(alpha - beta) for n odd, and a(n) = (alpha^n - beta^n)/(alpha^2 - beta^2) for n even, where alpha = (1/2)*(sqrt(5) + 1) and beta = (1/2)*(sqrt(5) - 1). Equivalently, a(n) = U(n-1, sqrt(5)/2) for n odd and a(n) = (1/sqrt(5))*U(n-1, sqrt(5)/2) for n even, where U(n,x) is the Chebyshev polynomial of the second kind. (End)
E.g.f.: (Phi/sqrt(5))*exp(-Phi*x)*(exp(x)-1)*(exp(sqrt(5)*x) - 1/(Phi)^2), where Phi = (1+sqrt(5))/2. - G. C. Greubel, Feb 08 2016
a(n) = (5^floor((n-1)/2)/2^(n-1))*Sum_{k=0..n-1} binomial(n-1,k)/5^floor(k/2). - Tony Foster III, Oct 21 2018
a(n) = hypergeom([(1 - n)/2, (n + 1) mod 2 - n/2], [1 - n], -4) for n >= 2. - Peter Luschny, Sep 03 2019

Additional comments from Michael Somos, Jun 01 2000

A064680 Halve every even number, double every odd number.

Original entry on oeis.org

0, 2, 1, 6, 2, 10, 3, 14, 4, 18, 5, 22, 6, 26, 7, 30, 8, 34, 9, 38, 10, 42, 11, 46, 12, 50, 13, 54, 14, 58, 15, 62, 16, 66, 17, 70, 18, 74, 19, 78, 20, 82, 21, 86, 22, 90, 23, 94, 24, 98, 25, 102, 26, 106, 27, 110, 28, 114, 29, 118, 30, 122, 31, 126, 32, 130, 33, 134, 34, 138, 35

Offset: 0

Reinhard Zumkeller, Oct 16 2001

a(a(m)) = m for all m with gcd(m, 4) <= 2.

Harry J. Smith, Table of n, a(n) for n = 0..1000, (corrected by _Peter Luschny_, Jan 19 2019)
Index entries for linear recurrences with constant coefficients, signature (0,2,0,-1).

Cf. A001477, A016825, A026741, A060819, A109043, A176447.

a064680 n = a064680_list !! n
a064680_list = zipWith ($) (cycle [(`div` 2), (* 2)]) [0..]
-- Reinhard Zumkeller, Jul 25 2012

[IsEven(n) select n/2 else 2*n: n in [0..70]]; // Bruno Berselli, Mar 09 2011

A064680:=n->n*(5-3*(-1)^n)/4: seq(A064680(n), n=0..100); # Wesley Ivan Hurt, Jul 24 2016

f[n_] := 2^(2 Mod[n, 2] - 1) n; Array[f, 70, 0] (* Or *)
f[n_] := If[ OddQ@ n, 2 n, n/2]; Array[f, 71, 0] (* Or *)
CoefficientList[ Series[x (2 + x + 2 x^2)/(1 - x^2)^2, {x, 0, 70}], x] (* Robert G. Wilson v *)
LinearRecurrence[{0,2,0,-1},{0,2,1,6},80] (* Harvey P. Dale, Mar 31 2025 *)

a(n) = if (n%2, 2*n, n/2); \\ Harry J. Smith, Sep 22 2009

a(n) = n * 2^(2 * (n mod 2) - 1).
G.f.: x*(2+x+2*x^2)/(1-x^2)^2.
a(n) = 2*a(n-2) - a(n-4) for n>3.
a(n)*a(n+3) = -2 + a(n+1)*a(n+2).
a(2*n) = A001477(n), a(1+2*n) = A016825(n). - Paul Curtz, Mar 09 2011
a(n) = n*(5-3*(-1)^n)/4. - Bruno Berselli, Mar 09 2011
a(n)= (period 4 sequence: repeat 2, 2, 1, 2) * (A060819(n)=0,1,1,3,1,5,...). - Paul Curtz, Mar 10 2011
E.g.f.: x*(sinh(x) + 4*cosh(x))/2. - Ilya Gutkovskiy, Jul 24 2016
a(n) = lcm(numerator(n/2), denominator(n/2)). - Wesley Ivan Hurt, Jul 24 2016
a(n) = A176447(n) + n. - Filip Zaludek, Dec 10 2016
From Amiram Eldar, Oct 07 2023: (Start)
a(n) = lcm(n,2) / gcd(n,2).
Sum_{k=1..n} a(k) ~ (5/8)*n^2. (End)

A384055 The number of integers k from 1 to n such that the greatest divisor of k that is a unitary divisor of n is odd.

Original entry on oeis.org

1, 1, 3, 3, 5, 3, 7, 7, 9, 5, 11, 9, 13, 7, 15, 15, 17, 9, 19, 15, 21, 11, 23, 21, 25, 13, 27, 21, 29, 15, 31, 31, 33, 17, 35, 27, 37, 19, 39, 35, 41, 21, 43, 33, 45, 23, 47, 45, 49, 25, 51, 39, 53, 27, 55, 49, 57, 29, 59, 45, 61, 31, 63, 63, 65, 33, 67, 51, 69

Offset: 1

Amiram Eldar, May 18 2025

Amiram Eldar, Table of n, a(n) for n = 1..10000

Unitary analog of A026741.
Cf. A000265, A006519, A384046, A384047.
The number of integers k from 1 to n such that the greatest divisor of k that is a unitary divisor of n is: A047994 (1), A384048 (squarefree), A384049 (cubefree), A384050 (powerful), A384051 (cubefull), A384052 (square), A384053 (cube), A384054 (exponentially odd), this sequence (odd), A384056 (power of 2), A384057 (3-smooth), A384058 (5-rough).

f[p_, e_] := p^e - If[p == 2, 1, 0]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a,100]

a(n) = {my(f = factor(n)); prod(i = 1, #f~, f[i,1]^f[i,2] - if(f[i,1] == 2, 1, 0));}

Multiplicative with a(2^e) = 2^e-1, and p^e if p is an odd prime.
a(n) = n * A047994(n) / A384056(n).
a(n) = A047994(A006519(n)) * A000265(n).
Dirichlet g.f.: zeta(s-1) * (1 - 1/2^(s-1) + 1/2^(2*s-1))/(1 - 1/2^s).
Sum_{k=1..n} a(k) ~ (5/12) * n^2.

A154955 a(1) = 1, a(2) = -1, followed by 0, 0, 0, ... .

Original entry on oeis.org

1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0

Offset: 1

Mats Granvik, Jan 18 2009

Matrix inverse of A000012.
Moebius transform of the sequence A000035. Dirichlet inverse of A209229. Partial sums of a(n) is characteristic function of 1 (A063524). a(n)=(-1)^(n+1)*A019590(n). a(n) for n >= 1 is Dirichlet convolution of following functions b(n), c(n), a(n) = Sum_{d|n} b(d)*c(n/d): a(n) = A000012(n) * A092673(n). Examples of Dirichlet convolutions with function a(n), i.e. b(n) = Sum_{d|n} a(d)*c(n/d): a(n) * A000012(n) = A000035(n), a(n) * A000027(n) = A026741(n), a(n) * A008683(n) = A092673(n), a(n) * A036987(n-1) = A063524(n), a(n) * A000005(n) = A001227(n). - Jaroslav Krizek, Mar 21 2009
The Kn21 sums, see A180662, of triangle A108299 equal the terms of this sequence. - Johannes W. Meijer, Aug 14 2011
{a(n-1)}A132393.%20-%20_Wolfdieter%20Lang">{n>=1}, gives the alternating row sums of A132393. - _Wolfdieter Lang, May 09 2017
With offset 0 the alternating row sums of A097805. - Peter Luschny, Sep 07 2017

Antti Karttunen, Table of n, a(n) for n = 1..10000
Index entries for linear recurrences with constant coefficients, signature (1).

Cf. A019590 (absolute values).
Cf. A000005, A000012, A000027, A000035, A001227, A008683, A019590, A026741, A036987, A063524, A092673.
Cf. A097805, A108299, A132393, A180662.

Join[{1, -1}, 0 Range[3, 100]] (* Wesley Ivan Hurt, Jun 22 2024 *)
PadRight[{1,-1},120,0] (* Harvey P. Dale, Dec 10 2024 *)

A154955(n)=(n==1)-(n==2) \\ M. F. Hasler, Jan 13 2012

G.f.: A(x) = x - x^2 = x / (1 + x / (1 - x)). - Michael Somos, Jan 03 2013
a(n) = (2/sqrt(3))*sin((2*Pi/3)*n!). - Lorenzo Pinlac, Jan 16 2022
a(n) = [n = 1] - [n = 2], where [] is the Iverson bracket. - Wesley Ivan Hurt, Jun 22 2024
Multiplicative with a(2) = -1, a(2^e) = 0 if e > 1, a(p^e) = 0 if p > 2. - Antti Karttunen, Dec 17 2024

Keyword:mult added by Antti Karttunen, Dec 17 2024

A106609 Numerator of n/(n+8).

Original entry on oeis.org

0, 1, 1, 3, 1, 5, 3, 7, 1, 9, 5, 11, 3, 13, 7, 15, 2, 17, 9, 19, 5, 21, 11, 23, 3, 25, 13, 27, 7, 29, 15, 31, 4, 33, 17, 35, 9, 37, 19, 39, 5, 41, 21, 43, 11, 45, 23, 47, 6, 49, 25, 51, 13, 53, 27, 55, 7, 57, 29, 59, 15, 61, 31, 63, 8, 65, 33, 67, 17, 69, 35, 71, 9, 73, 37, 75, 19, 77, 39, 79

Offset: 0

N. J. A. Sloane, May 15 2005

The graph of this sequence is made up of four linear functions: a(n_odd)=n, a(n=2+4i)=n/2, a(4+8i)=n/4, a(8i)=n/8. - Zak Seidov, Oct 30 2006. [In general, f(n) = numerator of n/(n+m) consists of linear functions n/d_i, where d_i are divisors of m (including 1 and m).]
a(n+2), n>=0, is the denominator of the harmonic mean H(n,2) = 4*n/(n+2). a(n+2) = (n+2)/gcd(n+2,8). a(n+5) = A227042(n+2, 2), n >= 0. - Wolfdieter Lang, Jul 04 2013
The sequence p(n) = a(n-4), n>=1, with a(-3) = a(3) = 3, a(-2) = a(2) = 1 and a(-1) = a(1) = 1, appears in the problem of writing 2*sin(2*Pi/n) as an integer in the algebraic number field Q(rho(q(n))), where rho(k) = 2*cos(Pi/k) and q(n) = A225975(n). One has 2*sin(2*Pi/n) = R(p(n), x) modulo C(q(n), x), with x = rho(q(n)) and the integer polynomials R and C given in A127672 and A187360, respectively. See a comment on A225975. - Wolfdieter Lang, Dec 04 2013
A204455(n) divides a(n) for n>=1. - Alexander R. Povolotsky, Apr 06 2015
A multiplicative sequence. Also, a(n) is a strong divisibility sequence, that is, gcd(a(n),a(m)) = a(gcd(n,m)) for n >= 1, m >= 1. In particular, a(n) is a divisibility sequence: if n divides m then a(n) divides a(m). - Peter Bala, Feb 20 2019

N. J. A. Sloane, Table of n, a(n) for n = 0..10000
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,0,2,0,0,0,0,0,0,0,-1).

Cf. A109049, A204455, A225975, A227042 (second column, starting with a(5)).

List([0..80],n->NumeratorRat(n/(n+8))); # Muniru A Asiru, Feb 19 2019

[Numerator(n/(n+8)): n in [0..100]]; // Vincenzo Librandi, Apr 18 2011

a := n -> iquo(n, [8, 1, 2, 1, 4, 1, 2, 1][1 + modp(n, 8)]):
seq(a(n), n=0..79); # using Wolfdieter Lang's formula, Peter Luschny, Feb 22 2019

f[n_]:=Numerator[n/(n+8)];Array[f,100,0] (* Vladimir Joseph Stephan Orlovsky, Feb 16 2011 *)
LinearRecurrence[{0,0,0,0,0,0,0,2,0,0,0,0,0,0,0,-1},{0,1,1,3,1,5,3,7,1,9,5,11,3,13,7,15},100] (* Harvey P. Dale, Sep 27 2019 *)

vector(100, n, n--; numerator(n/(n+8))) \\ G. C. Greubel, Feb 19 2019

[lcm(n,8)/8 for n in range(0, 100)] # Zerinvary Lajos, Jun 09 2009

a(n) = 2*a(n-8) - a(n-16).
G.f.: x* (x^2-x+1) * (x^12 +2*x^11 +4*x^10 +3*x^9 +4*x^8 +4*x^7 +7*x^6 +4*x^5 +4*x^4 +3*x^3 +4*x^2 +2*x +1) / ( (x-1)^2 *(x+1)^2 *(x^2+1)^2 *(x^4+1)^2 ). - R. J. Mathar, Dec 02 2010
From R. J. Mathar, Apr 18 2011: (Start)
a(n) = A109049(n)/8.
Dirichlet g.f. zeta(s-1)*(1-1/2^s-1/2^(2s)-1/2^(3s)).
Multiplicative with a(2^e) = 2^max(0,e-3). a(p^e) = p^e if p>2. (End)
a(n) = n/gcd(n,8), n >= 0. See the harmonic mean comment above. - Wolfdieter Lang, Jul 04 2013
a(n) = n if n is odd; for n == 0 (mod 8) it is n/8, for n == 2 or 6 (mod 8) it is n/2 and for n == 4 (mod 8) it is n/4. - Wolfdieter Lang, Dec 04 2013
From Peter Bala, Feb 20 2019: (Start)
O.g.f.: Sum_{n >= 0} a(n)*x^n = F(x) - F(x^2) - F(x^4) - F(x^8), where F(x) = x/(1 - x)^2.
More generally, for m >= 1, Sum_{n >= 0} (a(n)^m)*x^n = F(m,x) + (1 - 2^m)*( F(m,x^2) + F(m,x^4) + F(m,x^8) ), where F(m,x) = A(m,x)/(1 - x)^(m+1) with A(m,x) the m-th Eulerian polynomial: A(1,x) = x, A(2,x) = x*(1 + x), A(3,x) = x*(1 + 4*x + x^2) - see A008292.
Sum_{n >= 1} (1/n)*a(n)*x^n = G(x) - (1/2)*G(x^2) - (1/4)*G(x^4) - (1/8)*G(x^8), where G(x) = x/(1 - x).
Sum_{n >= 1} (1/n^2)*a(n)*x^n = L(x) - (1/2^2)*L(x^2) - (1/4)^2*L(x^4) - (1/8)^2*L(x^8), where L(x) = Log(1/(1 - x)).
Sum_{n >= 1} (1/a(n))*x^n = L(x) + (1/2)*L(x^2) + (1/2)*L(x^4) + (1/2)*L(x^8). (End)
Sum_{k=1..n} a(k) ~ (43/128) * n^2. - Amiram Eldar, Nov 25 2022

cp's OEIS Frontend

A106614 a(n) = numerator of n/(n+13).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

PARI

Sage

Formula

A162610 Triangle read by rows in which row n lists n terms, starting with 2n-1, with gaps = n-1 between successive terms.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

Python

Formula

Extensions

A051715 Denominators of table a(n,k) read by antidiagonals: a(0,k) = 1/(k+1), a(n+1,k) = (k+1)(a(n,k)-a(n,k+1)), n >= 0, k >= 0.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A106621 a(n) = numerator of n/(n+20).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

Extensions

A195817 Multiples of 10 and odd numbers interleaved.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A005013 a(n) = 3*a(n-2) - a(n-4), a(0)=0, a(1)=1, a(2)=1, a(3)=4. Alternates Fibonacci (A000045) and Lucas (A000032) sequences for even and odd n.

Views

Author

Keywords

Comments

Examples

References

Links