cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A100047 A Chebyshev transform of the Fibonacci numbers.

Original entry on oeis.org

0, 1, 1, -1, -1, 0, -1, -1, 1, 1, 0, 1, 1, -1, -1, 0, -1, -1, 1, 1, 0, 1, 1, -1, -1, 0, -1, -1, 1, 1, 0, 1, 1, -1, -1, 0, -1, -1, 1, 1, 0, 1, 1, -1, -1, 0, -1, -1, 1, 1, 0, 1, 1, -1, -1, 0, -1, -1, 1, 1, 0, 1, 1, -1, -1, 0, -1, -1, 1, 1, 0, 1, 1, -1, -1, 0, -1, -1, 1, 1, 0, 1, 1, -1, -1, 0, -1, -1, 1, 1
Offset: 0

Views

Author

Paul Barry, Oct 31 2004

Keywords

Comments

Multiplicative with a(p^e) = (-1)^(e+1) if p = 2, 0 if p = 5, 1 if p == 1 or 9 (mod 10), (-1)^e if p == 3 or 7 (mod 10). - David W. Wilson, Jun 10 2005
This sequence is a divisibility sequence, i.e., a(n) divides a(m) whenever n divides m. Case P1 = 1, P2 = -1, Q = 1 of the 3 parameter family of 4th-order linear divisibility sequences found by Williams and Guy. - Peter Bala, Mar 24 2014
From Peter Bala, Mar 24 2014: (Start)
This is the particular case P1 = 1, P2 = -1, Q = 1 of the following results:
Let P1, P2 and Q be integers. Let alpha and beta denote the roots of the quadratic equation x^2 - 1/2*P1*x + 1/4*P2 = 0. Let T(n,x;Q) denote the bivariate Chebyshev polynomial of the first kind defined by T(n,x;Q) = 1/2*( (x + sqrt(x^2 - Q))^n + (x - sqrt(x^2 - Q))^n ) (when Q = 1, T(n,x;Q) reduces to the ordinary Chebyshev polynomial of the first kind T(n,x)). Then we have
1) The sequence A(n) := ( T(n,alpha;Q) - T(n,beta;Q) )/(alpha - beta) is a linear divisibility sequence of the fourth order.
2) A(n) belongs to the 3-parameter family of fourth-order divisibility sequences found by Williams and Guy.
3) The o.g.f. of the sequence A(n) is the rational function x*(1 - Q*x^2)/(1 - P1*x + (P2 + 2*Q)*x^2 - P1*Q*x^3 + Q^2*x^4).
4) The o.g.f. is the Chebyshev transform of the rational function x/(1 - P1*x + P2*x^2), where the Chebyshev transform takes the function A(x) to the function (1 - Q*x^2)/(1 + Q*x^2)*A(x/(1 + Q*x^2)).
5) Let q = sqrt(Q) and set a = sqrt( q + (P2)/(4*q) + (P1)/2 ) and b = sqrt( q + (P2)/(4*q) - (P1)/2 ). Then the o.g.f. of the sequence A(n) is the Hadamard product of the rational functions x/(1 - (a + b)*x + q*x^2) and x/(1 - (a - b)*x + q*x^2). Thus A(n) is the product of two (usually, non-integer) Lucas-type sequences.
6) A(n) = the bottom left entry of the 2 X 2 matrix 2*T(n,1/2*M;Q), where M is the 2 X 2 matrix [0, -P2; 1, P1].
For examples of the above see A006238, A054493, A078070, A092184, A098306, A100048, A108196, A138573, A152090 and A218134. (End)

Examples

			A Chebyshev transform of the Fibonacci numbers A000045: if A(x) is the g.f. of a sequence, map it to ((1-x^2)/(1+x^2))A(x/(1+x^2)).
The denominator is the 10th cyclotomic polynomial.
G.f. = x + x^2 - x^3 - x^4 - x^6 - x^7 + x^8 + x^9 + x^11 + x^12 - x^13 + ...
		

Crossrefs

Programs

  • Maple
    a := n -> (-1)^iquo(n, 5)*signum(mods(n, 5)):
    seq(a(n), n=0..89); # after Michael Somos, Peter Luschny, Dec 30 2018
  • Mathematica
    a[ n_] := {1, 1, -1, -1, 0, -1, -1, 1, 1, 0}[[Mod[ n, 10, 1]]]; (* Michael Somos, May 24 2015 *)
    a[ n_] := (-1)^Quotient[ n, 5] Sign[ Mod[ n, 5, -2]]; (* Michael Somos, May 24 2015 *)
    a[ n_] := (-1)^Quotient[n, 5] {1, 1, -1, -1, 0}[[Mod[ n, 5, 1]]]; (* Michael Somos, May 24 2015 *)
    LinearRecurrence[{1, -1, 1, -1}, {0, 1, 1, -1}, 90] (* Jean-François Alcover, Jun 11 2019 *)
  • PARI
    {a(n) = (-1)^(n\5) * [0, 1, 1, -1, -1][n%5+1]}; /* Michael Somos, May 24 2015 */
    
  • PARI
    {a(n) = (-1)^(n\5) * sign( centerlift( Mod(n, 5)))}; /* Michael Somos, May 24 2015 */

Formula

G.f.: x*(1 - x^2)/(1 - x + x^2 - x^3 + x^4).
a(n) = a(n-1) - a(n-2) + a(n-3) - a(n-4).
a(n) = n*Sum_{k=0..floor(n/2)} (-1)^k *binomial(n-k, k)*A000045(n-2*k)/(n -k).
From Peter Bala, Mar 24 2014: (Start)
a(n) = (T(n,alpha) - T(n,beta))/(alpha - beta), where alpha = (1 + sqrt(5))/4 and beta = (1 - sqrt(5))/4 and T(n,x) denotes the Chebyshev polynomial of the first kind.
a(n) = bottom left entry of the matrix T(n, M), where M is the 2 X 2 matrix [0, 1/4; 1, 1/2].
The o.g.f. is the Hadamard product of the rational functions x/(1 - 1/2*(sqrt(5) + 1)*x + x^2) and x/(1 - 1/2*(sqrt(5) - 1)*x + x^2). (End)
Euler transform of length 10 sequence [ 1, -2, 0, 0, -1, 0, 0, 0, 0, 1]. - Michael Somos, May 24 2015
a(n) = a(-n) = -a(n + 5) for all n in Z. - Michael Somos, May 24 2015
|A011558(n)| = |A080891(n)| = |a(n)| = A244895(n). - Michael Somos, May 24 2015

A091394 a(n) = Product_{ p | n } (1 + Legendre(-5,p) ).

Original entry on oeis.org

1, 0, 2, 0, 1, 0, 2, 0, 2, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 4, 0, 2, 0, 1, 0, 2, 0, 2, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 2, 0, 2, 0, 2, 0, 2, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 4, 0, 0, 0, 2, 0, 4, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 2, 0, 2, 0, 0, 0, 4, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 2, 0, 4
Offset: 1

Views

Author

N. J. A. Sloane, Mar 02 2004

Keywords

Crossrefs

Programs

  • Maple
    with(numtheory); L := proc(n,N) local i,t1,t2; t1 := ifactors(n)[2]; t2 := mul((1+legendre(N,t1[i][1])),i=1..nops(t1)); end; [seq(L(n,-5),n=1..120)];
  • Mathematica
    a[n_] := Times@@ (1+KroneckerSymbol[-5, #]& /@ FactorInteger[n][[All, 1]]);
    Array[a, 105] (* Jean-François Alcover, Apr 08 2020 *)
  • PARI
    a(n)={my(f=factor(n)[, 1]); prod(i=1, #f, 1 + kronecker(-5, f[i]))} \\ Andrew Howroyd, Jul 25 2018

Formula

Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 5*sqrt(5)/(6*Pi) = 0.593135 . - Amiram Eldar, Oct 17 2022

A323378 Square array read by antidiagonals: T(n,k) = Kronecker symbol (-n/k), n >= 1, k >= 1.

Original entry on oeis.org

1, 1, 1, 1, 0, -1, 1, -1, 1, 1, 1, 0, 0, 0, 1, 1, -1, -1, 1, -1, -1, 1, 0, 1, 0, -1, 0, -1, 1, 1, 0, 1, 1, 0, -1, 1, 1, 0, -1, 0, 0, 0, 1, 0, 1, 1, 1, 1, 1, 1, -1, -1, -1, 1, 1, 1, 0, 0, 0, -1, 0, 1, 0, 0, 0, -1, 1, -1, -1, 1, -1, -1, 1, -1, 1, 1, 1, -1
Offset: 1

Views

Author

Jianing Song, Jan 12 2019

Keywords

Comments

If A215200 is arranged into a square array A215200(n,k) = kronecker symbol(n/k) with n >= 0, k >= 1, then this sequence gives the other half of the array.
Note that there is no such n such that the n-th row and the n-th column are the same.

Examples

			Table begins
  1,  1, -1,  1,  1, -1, -1,  1,  1,  1, ... ((-1/k) = A034947)
  1,  0,  1,  0, -1,  0, -1,  0,  1,  0, ... ((-2/k) = A188510)
  1, -1,  0,  1, -1,  0,  1, -1,  0,  1, ... ((-3/k) = A102283)
  1,  0, -1,  0,  1,  0, -1,  0,  1,  0, ... ((-4/k) = A101455)
  1, -1,  1,  1,  0, -1,  1, -1,  1,  0, ... ((-5/k) = A226162)
  1,  0,  0,  0,  1,  0,  1,  0,  0,  0, ... ((-6/k) = A109017)
  1,  1, -1,  1, -1, -1,  0,  1,  1, -1, ... ((-7/k) = A175629)
  1,  0,  1,  0, -1,  0, -1,  0,  1,  0, ... ((-8/k) = A188510)
  ...
		

Crossrefs

Cf. A215200.
The first rows are listed in A034947, A188510, A102283, A101455, A226162, A109017, A175629, A188510, ...

Programs

  • PARI
    T(n,k) = kronecker(-n, k)
Showing 1-3 of 3 results.