login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A056581
Nearest integer to 1/(A056580(n) - exp(sqrt(n)*Pi)).
7
-7, -51, 4, -2, -5, 110, 15, -3, 3, 5, -7, -3, 19, 4, 5, -3, 430, 141, 4, 4, -2, 574, 3, 7, 1518, -3, 62, 84, -2, -10, 11, -7, -13, -4, 4, -3, 45551, -5, 3, 3, 2, -33, 4494, -8, -5, -6, 3, -2, 7, 2, 9, -3, -4, -4, 3, -17, -2, 5624716, 147, -5, 4, 3, 3, 2, 6, -2, 747638
OFFSET
1,1
COMMENTS
A measure of how close e^(Pi*sqrt(n)) is to an integer (higher absolute value of a(n) means closer, negative value means the closest integer is smaller than it).
The sign convention is chosen so that most terms and in particular record values such as those occurring for the Heegner numbers A003173, are positive, so that A069014 lists record indices of this sequence (except for A069014(2)=2 instead of 3 for signed values). The sequence is not defined for n=0,-1 where e^(sqrt(n)*Pi) is an integer. - M. F. Hasler, Apr 15 2008
Negative resp. positive values of a(n) correspond to 2nd resp. 3rd term of the continued fraction expansion of exp(sqrt(n)*Pi), up to a difference of -1 or -2 depending on the direction of rounding. - M. F. Hasler, Apr 15 2008
REFERENCES
For links, references and more information see A019296 and other cross-referenced sequences.
FORMULA
a(n) = 1/(A056580(n) - e^(sqrt(n)*Pi)).
A019296 ={-1, 0} U { n | abs(A056581(n)) > 100} U { some n for which abs(A056581(n)) = 100 }. - M. F. Hasler, Apr 15 2008
EXAMPLE
a(6)=110, since e^(Pi*sqrt(6)) = 2197.9908695... and 1/(2198-2197.9908695...) = 109.52... which rounds to 110.
e^(Pi*sqrt(163)) = 262537412640768743.9999999999992500725971981... (the Ramanujan number) and so a(163)=1333462407513.
PROG
(PARI) default(realprecision, 100); dZ(x)=round(x)-x
A056581(n)=round(1/dZ(exp(sqrt(n)*Pi))
KEYWORD
sign
AUTHOR
Henry Bottomley, Jun 30 2000
EXTENSIONS
Definition, formulas and values corrected and extended by M. F. Hasler, Apr 15 2008
STATUS
approved