OFFSET
1,1
COMMENTS
Conjectures: sigma(6n) > sigma(6n - 1) and sigma(6n) > sigma(6n + 1).
Conjectures are false. Try prime 73961483429 for n. One finds sigma(6*73961483429) < sigma(6*73961483429+1). The number n = 105851369791 provides a counterexample for the other case. - T. D. Noe, Apr 22 2013
Sum of the divisors of the numbers k which have the property that the width associated to the vertex of the first (also the last) valley of the smallest Dyck path of the symmetric representation of sigma(k) is equal to 2 (see example). Other positive integers have width 0 or 1 associated to the mentioned valley. - Omar E. Pol, Aug 11 2021
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = A000203(6n).
Sum_{k=1..n} a(k) = (55*Pi^2/72) * n^2 + O(n*log(n)). - Amiram Eldar, Dec 16 2022
EXAMPLE
From Omar E. Pol, Aug 11 2021: (Start)
Illustration of initial terms:
----------------------------------------------------------------------
n 6*n a(n) Diagram: 1 2 3 4
----------------------------------------------------------------------
_ _ _ _
| | | | | | | |
| | | | | | | |
* _ _| | | | | | | |
| _ _| | | | | | |
_ _ _| |_| | | | | | |
1 6 12 |_ _ _ _| * _ _ _| | | | | |
_| _ _ _| | | | |
* _| _| | | | | |
| _| _| * _ _ _ _| | | |
| |_ _| | _ _ _ _| | |
_ _ _ _ _ _| | _| | | | |
2 12 28 |_ _ _ _ _ _ _| _| _|_| * _ _ _ _ _| |
* _ _| _| | _ _ _ _ _|
| _ _| _ _| | |
| |_ _| _| _ _| |
| | _| _| _ _|
_ _ _ _ _ _ _ _ _| | | _| _|
3 18 39 |_ _ _ _ _ _ _ _ _ _| * _ _| | _|
| _ _| |
| |_ _ _|
| |
| |
_ _ _ _ _ _ _ _ _ _ _ _| |
4 24 60 |_ _ _ _ _ _ _ _ _ _ _ _ _|
.
Note that the mentioned vertices are aligned on two straight lines that meet at point (3,3).
a(n) equals the area (also the number of cells) in the n-th diagram. (End)
MATHEMATICA
DivisorSigma[1, 6*Range[60]] (* Harvey P. Dale, Apr 16 2016 *)
PROG
(PARI) a(n)=sigma(6*n) \\ Charles R Greathouse IV, Apr 22 2013
(Python)
from sympy import divisor_sigma
def a(n): return divisor_sigma(6*n)
print([a(n) for n in range(1, 54)]) # Michael S. Branicky, Dec 28 2021
(Python)
from math import prod
from collections import Counter
from sympy import factorint
def A224613(n): return prod((p**(e+1)-1)//(p-1) for p, e in (Counter(factorint(n))+Counter([2, 3])).items()) # Chai Wah Wu, Sep 07 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Apr 22 2013
EXTENSIONS
Corrected by Harvey P. Dale, Apr 16 2016
STATUS
approved