OFFSET
1,1
COMMENTS
First differs from A004083 at n = 19564195 where a(n) = 122925461 and A004083(n) = 122925462. - Martin Fuller, Apr 10 2025
LINKS
Martin Fuller, Proof for formula
Anthony G. Shannon, Bahar Kuloğlu, and Engin Özkan, Rhaly terraced sequences their generalizations, properties and applications, Comp. Appl. Math., Vol. 44 (2025), Article No. 226. See p. 4.
FORMULA
a(n) = ceiling(sqrt(x^2+1)) where x is the root of 1+sqrt(x^2+1)*cos(x) in the range [(2n-1)Pi, (2n-1/2)Pi]. - Martin Fuller, Apr 10 2025
EXAMPLE
For k = 0..4, there is 1 root; for k = 5..10, there are 5 roots; for k = 11..17, there are 9 roots; etc... .
PROG
(PARI) a(n)=my(x=solve(x=(2*n-1)*Pi, (2*n-1/2)*Pi, 1+sqrt(x^2+1)*cos(x))); ceil(sqrt(x^2+1)); \\ Martin Fuller, Apr 10 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Sébastien Dumortier, Sep 28 2008
EXTENSIONS
a(0)=0 removed, and a(16) onwards from Martin Fuller, Apr 10 2025
STATUS
approved
