OFFSET
1,1
LINKS
Bhavik Mehta, Table of n, a(n) for n = 1..10000
Thomas Bloom, Problem 290
Wouter van Doorn, On the non-monotonicity of the denominator of generalized harmonic sums, arXiv:2411.03073 [math.NT], 2024.
FORMULA
a(n) < 4.374*n for all n > 1. - Wouter van Doorn, Nov 06 2024
EXAMPLE
1/3+1/4+1/5=47/60 and 1/3+1/4+1/5+1/6=19/20, and 60>20, so a(3)=5.
PROG
(PARI) a(n) = for(k=0, oo, my(s=sum(n=n, n+k, 1/n)); if(denominator(s)>denominator(s+1/(n+k+1)), return(n+k); break))
(Python)
from fractions import Fraction
from itertools import count
def A375081(n):
a = Fraction((n<<1)+1, n*(n+1))
for k in count(n+1):
if a.denominator > (a:=a+Fraction(1, k+1)).denominator:
return k # Chai Wah Wu, Jul 30 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Ralf Stephan, Jul 29 2024
EXTENSIONS
a(56) onwards from Bhavik Mehta, Jul 31 2024
STATUS
approved