MatsGranvik

Dirichlet divisor problem and sum of roots of unity

Jul 21st, 2022 (edited)
100
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.21 KB | None | 0 0
  1. (*start*)
  2. nn = 18;
  3. Accumulate[Table[DivisorSigma[0, n], {n, 1, nn - 1}]]
  4. Monitor[Table[-1/2*n + n*HarmonicNumber[n] +
  5. Sum[(I*Log[-Exp[(2 I n Pi )/k]])/(2*Pi), {k, 1, n}], {n, 2, nn}], n]
  6. FullSimplify[%]
  7. %%% - %
  8. (*end*)
  9.  
  10. {1, 3, 5, 8, 10, 14, 16, 20, 23, 27, 29, 35, 37, 41, 45, 50, 52}
  11.  
  12. (*expanded form start*)
  13. TableForm[
  14. T = Table[
  15. Table[Chop[N[-1/2 + n/k + I*Log[-E^(2*I*n*Pi/k)]/(2*Pi)]], {k, 1,
  16. n}], {n, 1, 12}]]
  17. Flatten[Round[T]]
  18. (*end*)
  19.  
  20. (*start*)
  21. T(n,k) = -1/2 + n/k + I*Log[-E^(2*I*n*Pi/k)]/(2*Pi)
  22. (*end*)
  23.  
  24.  
  25. (*more*)
  26. (*expanded form start*)
  27. nn = 12; TableForm[
  28. T = Table[
  29. Table[k*FullSimplify[
  30. 0*(-1/2 + n/k) + I*Log[-E^(2*I*n*Pi/k)]/(2*Pi)], {k, 1,
  31. nn}], {n, 1, nn}]]
  32. TableForm[2*(T[[2 ;; nn, 1 ;; nn]] - Accumulate[Differences[T]])]
  33. TableForm[Accumulate[Differences[T]]]
  34. (*end*)
  35.  
  36. (*expanded form start*)
  37. nn = 12; TableForm[
  38. T = Table[
  39. Table[FullSimplify[(-1/2 + n/k) +
  40. I*Log[-E^(2*I*n*Pi/k)]/(2*Pi)], {k, 1, nn}], {n, 1, nn}]]
  41. TableForm[
  42. Table[Table[(-1/2 + n/k) + (k - 2 - 2*Mod[n - 1, k])/(2*k), {k, 1,
  43. nn}], {n, 1, nn}]]
  44. Clear[n, k, x]
  45. FullSimplify[(-1/2 + n/k) + (k - 2 - 2*x)/(2*k)]
  46. FullSimplify[(-1/2 + (n + 1)/k) + (k - 2 - 2*x)/(2*k)]
Add Comment
Please, Sign In to add comment