I can't help thinking that the answer ought to come out of a combination of standard techniques. First, indeed, reduce to fractional parts: but with a caveat, that it would be sensible to do some manipulation first. I'm thinking about "summation by parts". Then we know the fractional parts are uniformly distributed mod 1
Edit: There is a classic reference: G. On the faceH. Hardy & J. E. Littlewood, "Some problems of itDiophantine approximation: The lattice points of a right-angled triangle, though" (1st memoir), this won't be enoughProc. London Math. Soc. (2), iv.e 20, 1922, pp. 15-36. They consider the averagemodified fractional part will therefore be 0.5. But the problem issum, with {x} set up so that this "main term" contribution drops out? Otherwise why should the error term be small?
So the next recourse should be to exploit the Fourier seriesas x - [x] - 1/2, of the sawtooth function (not just its constant term). The Fourier series should be truncated at around${k\theta}$ up to n terms, and the error estimated away. Then there are geometric progressionswhere to sum. This is what brings inbe compatible with their notation $\theta$ would be the diophantine approximation propertiesreciprocal of $\tau$, in the denominators, which basically can be smallnot that this matters at Fibonacci numbers. We do know how small. We do know how many Fibonacci numbers are involved for the terms under considerationall.
Forgive me if I'm not writing down the many summations and estimates here The bound they get is O(log n) (Hardy's Works vol. I don't find MO well adapted to talking about analytic number theoryp. We're basically talking here about145), which depends only on the technology of Weyl's equidistribution theorem being applied to acontinued fraction having partial quotients. The particular case, relevant to say more than simply that$\tau$ is worked out in detail over the fractional parts are spread uniformly on [0, 1], but to dig down a bit further using exponential sumsnext few pages.
Edit: So I want to rely on The result is sharp. Off the abilitytop of my head this looks enough to estimate sums like
$\sum_{j=1}^{M} exp(2\pi i\tau rj)/j$
as bounded by
$|cosec(\pi \tau r)|$
andget the corresponding
$\sum_{j=M+1}^{infinity}$
error term O(nlog n) for the sum as posted, by log Mbreaking into at most n sums of this type.