login
This site is supported by donations to The OEIS Foundation.

 

Logo

Invitation: celebrating 50 years of OEIS, 250000 sequences, and Sloane's 75th, there will be a conference at DIMACS, Rutgers, Oct 9-10 2014.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A193383 G.f. is the imaginary part of the function C(x) that satisfies: C(x) = 1 + x/C(I*x). 2
-1, -1, 0, 1, 2, 2, 0, -3, -6, -6, 0, 10, 20, 20, 0, -35, -70, -70, 0, 126, 252, 252, 0, -462, -924, -924, 0, 1716, 3432, 3432, 0, -6435, -12870, -12870, 0, 24310, 48620, 48620, 0, -92378, -184756, -184756, 0, 352716, 705432, 705432, 0, -1352078, -2704156, -2704156, 0 (list; graph; refs; listen; history; text; internal format)
OFFSET

2,5

LINKS

Table of n, a(n) for n=2..52.

FORMULA

G.f.: imag(C(x)) = (x/2)*(1 - sqrt(1 + 4*x^4)/(1 - 2*x + 2*x^2)),

where C(x) = (1 + (1+I)*x)/2 + sqrt(1 + 4*x^4)/(2*(1 - (1-I)*x)),

and real(C(x)) = (1+x)/2 + (1-x)/2 * sqrt(1 + 4*x^4)/(1 - 2*x + 2*x^2).

EXAMPLE

C(x) = 1 + x - I*x^2 + (-1 - I)*x^3 - x^4 + (-1 + I)*x^5 + 2*I*x^6 + (2 + 2*I)*x^7 + 3*x^8 + (3 - 3*I)*x^9 - 6*I*x^10 + (-6 - 6*I)*x^11 - 10*x^12 +...

where

C(I*x)^-1 = 1 - I*x + (-1 - I)*x^2 - x^3 + (-1 + I)*x^4 + 2*I*x^5 +...

The real part of C(x) begins:

real(C(x)) = 1 + x - x^3 - x^4 - x^5 + 2*x^7 + 3*x^8 + 3*x^9 - 6*x^11 - 10*x^12 - 10*x^13 + 20*x^15 + 35*x^16 + 35*x^17 - 70*x^19 - 126*x^20 - 126*x^21 +...

The imaginary part of C(x) begins:

imag(C(x)) = -x^2 - x^3 + x^5 + 2*x^6 + 2*x^7 - 3*x^9 - 6*x^10 - 6*x^11 + 10*x^13 + 20*x^14 + 20*x^15 - 35*x^17 - 70*x^18 - 70*x^19 + 126*x^21 +...

PROG

(PARI) {a(n)=local(A=1+x); for(i=1, n, A=1+x/subst(A, x, I*x +x*O(x^n))); imag(polcoeff(A, n))}

(PARI) {a(n)=polcoeff((x/2)*(1-sqrt(1+4*x^4+x*O(x^n))/(1-2*x+2*x^2+x*O(x^n))), n)}

CROSSREFS

Cf. A193382 (imag), A193384 (norm).

Sequence in context: A214639 A216973 A061314 * A218033 A118262 A065484

Adjacent sequences:  A193380 A193381 A193382 * A193384 A193385 A193386

KEYWORD

sign

AUTHOR

Paul D. Hanna, Jul 24 2011

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Superseeker | Recent | More pages
The OEIS Community | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified September 21 10:11 EDT 2014. Contains 247025 sequences.