A064368 Number of 2 X 2 symmetric singular matrices with entries from {0,...,n}.

1, 4, 7, 10, 15, 18, 21, 24, 29, 36, 39, 42, 47, 50, 53, 56, 65, 68, 75, 78, 83, 86, 89, 92, 97, 108, 111, 118, 123, 126, 129, 132, 141, 144, 147, 150, 163, 166, 169, 172, 177, 180, 183, 186, 191, 198, 201, 204, 213, 228, 239, 242, 247, 250, 257, 260, 265, 268

Offset

0,2

Links

Vaclav Kotesovec, Table of n, a(n) for n = 0..10000

Formula

a(n) = n + 1 + 2*Sum_{k=1..n} Sum_{d^2|k} phi(d), where phi = Euler totient function A000010.

a(n) ~ (n/zeta(2)) * (log(n) + 3*gamma - 1 + zeta(2) - 2*zeta'(2)/zeta(2)), where gamma is Euler's constant (A001620). - Amiram Eldar, Nov 07 2024

Mathematica

f[p_, e_] := p^Floor[e/2]; a[0] = 0; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; With[{max = 100}, 1 + Range[0, max] + 2 * Accumulate[Array[a, max + 1, 0]]] (* Amiram Eldar, Nov 07 2024 *)

Prog

(PARI) a(n) = n + 1 + 2*sum(k=1, n, sumdiv(k, d, issquare(d)*eulerphi(sqrtint(d)))) \\ Michel Marcus, Jun 17 2013

Crossrefs

Cf. A000010, A000188, A064276, A059306, A062801, A059976, A039623.

Cf. A001620, A306016.

Sequence in context: A310702 A310703 A298789 * A026372 A014690 A274383

Adjacent sequences: A064365 A064366 A064367 * A064369 A064370 A064371

Keyword

nonn

Author

Vladeta Jovovic, Sep 27 2001

Status

approved