Saturday, 14 September 2013

convergence divergence - Why does this sum converge: sumn1k=2frac1kfracnmodkn



I was playing around with numbers and wanted to create a function that somehow indicates if a number could be a prime. So I came up with this, with the intention that it should make small jumps if x is a prime number:




f(n)=n1k=21knmodkn



It inneed does small jumps at small primes, but in general the term seems to converge to something close around 0.422780.



Why is this? I would expect this sum to grow without bounds, as for a large number N, there should be a lot of numbers which do not divide N without remainder.



f(n) in the range 0 to 1000


Answer



We have n(modk)=k{n/k} where {x} is the fractional part of x.
Hence we are interested in the behaviour of the sort-of-average-value
f(n)=1nn1k=2{n/k}



which as n, by Riemann sums, converges to
I=10{1x}dxx1/x=+1{x}x2dx=m110x(x+m)2dx

which equals
m1(log(m+1)log(m)1m+1)=m1(1m1m+1)1m1(1mlog(1+1m))γ

i.e. 1γ0.422784335 as claimed by Daniel Fischer.


No comments:

Post a Comment

real analysis - How to find limhrightarrow0fracsin(ha)h

How to find limh0sin(ha)h without lhopital rule? I know when I use lhopital I easy get $$ \lim_{h\rightarrow 0}...