7
$\begingroup$

I have to solve find the value of $$\sum_{k=1}^{n/2} k\log k$$ as a part of question.

How should I proceed on this ?

  • 0
    @user8250: It's submission or summation.2011-03-14
  • 0
    I may be asking a silly question, but what is a "submission"?2011-03-14
  • 0
    @DJC srry :) it was a typo .... i meant sigma and i am reading more on how to use tags to write equations here....2011-03-14
  • 0
    Just checking! :)2011-03-14
  • 3
    I doubt you will find a "closed form" formula. Maybe an asymptotic approximation?2011-03-14
  • 5
    $log(1^1) + log(2^2) + log(3^3) + ... + log((n/2)^{n/2}) = log(1*4*27*...)$?2011-03-14
  • 0
    Can I create an upper bound of it ?? I reached the same stage as @The Chaz2011-03-14
  • 0
    Just an idea if anybody could expand on it .... We can write k as n/2 in each of them .... then it would be that O(n/2 * log((n/2)factorial)) ..... Can anybody concise this expression more....2011-03-14

2 Answers 2

7

Got it. The constant in Moron's answer is $C = \log A$, where $A$ is the Glaisher-Kinkelin constant. Thus $C = \frac{1}{12} - \zeta'(-1)$.

The expression $H(n) = \prod_{k=1}^n k^k$ is called the hyperfactorial, and it has the known asymptotic expansion

$$H(n) = A e^{-n^2/4} n^{n(n+1)/2+1/12} \left(1 + \frac{1}{720n^2} - \frac{1433}{7257600n^4} + \cdots \right).$$ Taking logs and using the fact that $\log (1 + x) = O(x)$ yields an asymptotic expression for the OP's sum $$\sum_{k=1}^n k \log k = C - \frac{n^2}{4} + \frac{n(n+1)}{2} \log n + \frac{\log n}{12} + O \left(\frac{1}{n^2}\right),$$ the same as the one Aryabhata obtained with Euler-Maclaurin summation.


Added: Finding an asymptotic formula for the hyperfactorial is Problem 9.28 in Concrete Mathematics (2nd ed.). The answer they give uses Euler-Maclaurin, just as Aryabhata's answer does. They also mention that a derivation of the value of $C$ is in N. G. de Bruijn's Asymptotic Methods in Analysis, $\S$3.7.

  • 0
    +1: Very nice ! :-) I am surprised the Plouffe inverter does not have this, though.2011-03-14
  • 1
    @Moron: Thanks! It's satisfying to have finally figured out what that constant is. :)2011-03-14
  • 0
    @Moron: It looks like Plouffe does have it (or a scaled version of it) after all. The numerical estimate I was using before wasn't precise enough.2011-03-14
  • 0
    I see. I was about to suggest maybe you should add it there :-)2011-03-14
  • 1
    I suppose one proof of this would be to start from $\zeta'(s) = -\sum \log k/k^s$ for $Re(s) \gt 1$ and extend it for all $s$ (which can be done using Euler McLaurin I believe), similar to the Riemann Zeta function.2011-03-14
  • 0
    @Moron: That sounds like it would work.2011-03-14
6

Here is an asymptotic expression using EulerMcLaurin Summation.

$$ \sum _{k=1}^{n} k \log k = \int_{1}^{n} x \log x\ \text{d}x + (n\log n)/2 + C' + (\log n + 1)/12+ \mathcal{O}(1/n^2)$$

$$ = n^2(2 \log n - 1)/4 + (n\log n)/2 + (\log n)/12 + C + \mathcal{O}(1/n^2)$$

for some constant $C$.

  • 0
    +1, although I get $C = \frac{1}{4}$, $(\log n)/12$ rather than $(\log n)/18$, and $O(\frac{1}{n^2})$. (And I have verified this numerically.)2011-03-14
  • 0
    @MIke: You are correct. I have edited the answer. Thanks. That C=1/4 would need a proof, but that is a neat value for that constant.2011-03-14
  • 0
    @Moron: You're right to be suspicious of $C = 1/4$. It's not correct. The value of $C$ to six decimal places appears to be $0.248755$ - close to $1/4$ but not quite $1/4$. I'm not sure how to prove that or get an explicit expression, though.2011-03-14
  • 0
    @Mike: Perhaps the Plouffe inverter has something. We would probably need more than 6 digits though.2011-03-14
  • 0
    I tried Plouffe, with the more precise estimate $0.24875448644$, but I didn't get anything exact. Just so you know, this number comes from using Mathematica to find the difference between the OP's sum and your estimate for large values of $n$.2011-03-14
  • 0
    @Mike: There is an infinite series representation using the $B_{2k}$ which comes from the Euler McLaurin formula. It might be easier to use that to calculate the expression and might also lead us to a "closed" form expression. It is very similar (but not quite the same) to the series we get for $\sum \log k$ which gives us the Stirling's constant $\sqrt{2\pi}$, I believe.2011-03-14
  • 0
    I tried that, but I got $1/4 + \sum_{k=2}^{\infty} \frac{B_{2k}}{(2k)!} (2k-3)!$, which doesn't converge. I could have made a mistake, though.2011-03-14
  • 0
    @Mike: You are right that it does not converge! I didn't actually try to work it out. Sorry about that.2011-03-14
  • 1
    But you're right that the same procedure for finding Stirling's constant ought to work here. There's got to be a way to get this! :)2011-03-14