|
| |
|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
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.
|
|
|
LINKS
|
|
|
|
FORMULA
|
a(n) = 1/(A056580(n) - e^(sqrt(n)*Pi)).
|
|
|
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))
|
|
|
CROSSREFS
|
|
|
|
KEYWORD
|
sign
|
|
|
AUTHOR
|
|
|
|
EXTENSIONS
|
Definition, formulas and values corrected and extended by M. F. Hasler, Apr 15 2008
|
|
|
STATUS
|
approved
|
| |
|
|
