login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A331691 Resultant of the Shapiro polynomials P_n(x) and Q_n(x). 2
1, 2, -16, 2048, -67108864, 144115188075855872, -1329227995784915872903807060280344576, 226156424291633194186662080095093570025917938800079226639565593765455331328 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

COMMENTS

The Shapiro polynomials P_n(x) and Q_n(x) are defined by P_0(x) = Q_0(x) = 1 and then mutual recurrences P_{n+1}(x) = P_n(x) + x^(2^n)*Q_n(x) and Q_{n+1}(x) = P_n(x) - x^(2^n)*Q_n(x).  The coefficients of P are the Golay-Rudin-Shapiro sequence A020985.  a(n) is the polynomial resultant R(P_n(x),Q_n(x)) as considered by Brillhart and Carlitz.

LINKS

Table of n, a(n) for n=0..7.

John Brillhart and L. Carlitz, Note on the Shapiro Polynomials, Proceedings of the American Mathematical Society, volume 25, number 1, May 1970, pages 114-118.  Also at JSTOR.  See A001782 for a scanned copy.

Harold S. Shapiro, Extremal Problems for Polynomials and Power Series, Masters Thesis, Massachusetts Institute of Technology, 1951.  See pages 40-41.

FORMULA

a(n) = (-1)^(n-1) * 2^(2^(n+1) - n - 2) for n >= 1 [Brillhart and Carlitz theorem 2].

a(n) = (-1)^(n-1) * A016031(n+2) for n >= 1.

a(n) = - 2^(2^n-1) * a(n-1) for n >= 2 [Brillhart and Carlitz in proof of theorem 2].

PROG

(PARI) a(n) = if(n==0, 1, -(-2)^(2^(n+1) - n - 2));

(PARI) a(n) = my(P=1, Q=1); for(i=0, n-1, [P, Q]=[P+x^(2^i)*Q, P-x^(2^i)*Q]); polresultant(P, Q);

CROSSREFS

Cf. A016031 (absolute values), A001782 (discriminant).

Sequence in context: A060597 A091479 A016031 * A001309 A132569 A165644

Adjacent sequences:  A331688 A331689 A331690 * A331692 A331693 A331694

KEYWORD

sign,easy

AUTHOR

Kevin Ryde, Jan 24 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 | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 4 14:58 EST 2021. Contains 341796 sequences. (Running on oeis4.)