cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-4 of 4 results.

A347956 Dirichlet convolution of A003602 with A069359.

Original entry on oeis.org

0, 1, 1, 3, 1, 8, 1, 7, 5, 11, 1, 22, 1, 14, 13, 15, 1, 35, 1, 31, 16, 20, 1, 50, 8, 23, 20, 40, 1, 81, 1, 31, 22, 29, 19, 95, 1, 32, 25, 71, 1, 106, 1, 58, 62, 38, 1, 106, 11, 77, 31, 67, 1, 134, 25, 92, 34, 47, 1, 217, 1, 50, 78, 63, 28, 156, 1, 85, 40, 151, 1, 215, 1, 59, 95, 94, 28, 181, 1, 151, 74, 65, 1, 286, 34
Offset: 1

Views

Author

Antti Karttunen, Sep 20 2021

Keywords

Crossrefs

Programs

Formula

a(n) = Sum_{d|n} A003602(d) * A069359(n/d).

A328203 Expansion of Sum_{k>=1} k * x^k / (1 - x^(2*k))^2.

Original entry on oeis.org

1, 2, 5, 4, 8, 10, 11, 8, 20, 16, 17, 20, 20, 22, 42, 16, 26, 40, 29, 32, 58, 34, 35, 40, 53, 40, 74, 44, 44, 84, 47, 32, 90, 52, 94, 80, 56, 58, 106, 64, 62, 116, 65, 68, 174, 70, 71, 80, 102, 106, 138, 80, 80, 148, 146, 88, 154, 88, 89, 168, 92, 94, 241, 64, 172
Offset: 1

Views

Author

Ilya Gutkovskiy, Oct 07 2019

Keywords

Crossrefs

Programs

  • Magma
    a:=[]; for k in [1..65] do if IsOdd(k) then a[k]:=(k * #Divisors(k) + DivisorSigma(1,k)) / 2; else a[k]:=(k * (#Divisors(k) - #Divisors(k div 2)) + DivisorSigma(1,k) - DivisorSigma(1,k div 2)) / 2;  end if; end for; a; // Marius A. Burtea, Oct 07 2019
    
  • Mathematica
    nmax = 65; CoefficientList[Series[Sum[k x^k/(1 - x^(2 k))^2, {k, 1, nmax}], {x, 0, nmax}], x] // Rest
    a[n_] := DivisorSum[n, (n Mod[#, 2] + Boole[OddQ[n/#]] #)/2 &]; Table[a[n], {n, 1, 65}]
  • PARI
    A328203(n) = if(n%2,(1/2)*(sigma(n)+(n*numdiv(n))),2*A328203(n/2)); \\ Antti Karttunen, Nov 13 2021

Formula

a(n) = (n * d(n) + sigma(n)) / 2 if n odd, (n * (d(n) - d(n/2)) + sigma(n) - sigma(n/2)) / 2 if n even.
a(n) = (n * A001227(n) + A002131(n)) / 2.
a(2*n) = 2 * a(n).
From Antti Karttunen, Nov 13 2021: (Start)
The following two convolutions were found by Jon Maiga's Sequence Machine search algorithm. Both are easy to prove:
a(n) = Sum_{d|n} A003602(d) * A026741(n/d).
a(n) = Sum_{d|n} A109168(d) * A193356(n/d), where A109168(d) = A140472(d) = (d+A006519(d))/2.
(End)

A349391 Dirichlet convolution of A126760 with omega.

Original entry on oeis.org

0, 1, 1, 2, 1, 4, 1, 3, 2, 5, 1, 7, 1, 6, 5, 4, 1, 7, 1, 9, 6, 7, 1, 10, 3, 8, 3, 11, 1, 16, 1, 5, 7, 9, 7, 12, 1, 10, 8, 13, 1, 20, 1, 13, 9, 11, 1, 13, 4, 18, 9, 15, 1, 10, 8, 16, 10, 13, 1, 27, 1, 14, 11, 6, 9, 24, 1, 17, 11, 32, 1, 17, 1, 16, 18, 19, 9, 28, 1, 17, 4, 17, 1, 34, 10, 18, 13, 19, 1, 27, 10, 21
Offset: 1

Views

Author

Antti Karttunen, Nov 15 2021

Keywords

Crossrefs

Cf. A347233, A347234, A349390, A349392, A349393, A349395 for other Dirichlet convolutions of A126760. And also A347957.

Programs

  • Mathematica
    f[n_] := 2 * Floor[(m = n/2^IntegerExponent[n, 2]/3^IntegerExponent[n, 3])/6] + Mod[m, 3]; a[n_] := DivisorSum[n, f[#] * PrimeNu[n/#] &]; Array[a, 100] (* Amiram Eldar, Nov 16 2021 *)
  • PARI
    A126760(n) = {n&&n\=3^valuation(n, 3)<A126760
    A349391(n) = sumdiv(n,d,A126760(n/d)*omega(d));

Formula

a(n) = Sum_{d|n} A126760(n/d) * A001221(d).

A349346 Dirichlet inverse of A181988, where A181988(n) = A001511(n)*A003602(n).

Original entry on oeis.org

1, -2, -2, 1, -3, 4, -4, 0, -1, 6, -6, -2, -7, 8, 4, 0, -9, 2, -10, -3, 5, 12, -12, 0, -4, 14, -2, -4, -15, -8, -16, 0, 7, 18, 6, -1, -19, 20, 8, 0, -21, -10, -22, -6, 3, 24, -24, 0, -9, 8, 10, -7, -27, 4, 8, 0, 11, 30, -30, 4, -31, 32, 4, 0, 9, -14, -34, -9, 13, -12, -36, 0, -37, 38, 8, -10, 9, -16, -40, 0, -4, 42
Offset: 1

Views

Author

Antti Karttunen, Nov 15 2021

Keywords

Crossrefs

Programs

  • PARI
    up_to = 20000;
    DirInverseCorrect(v) = { my(u=vector(#v)); u[1] = (1/v[1]); for(n=2, #v, u[n] = (-u[1]*sumdiv(n, d, if(dA001511(n) = 1+valuation(n,2);
    A003602(n) = (1+(n>>valuation(n,2)))/2;
    A181988(n) = (A001511(n)*A003602(n));
    v349346 = DirInverseCorrect(vector(up_to,n,A181988(n)));
    A349346(n) = v349346[n];

Formula

a(1) = 1; a(n) = -Sum_{d|n, d < n} A181988(n/d) * a(d).
a(n) = A349347(n) - A181988(n).
Showing 1-4 of 4 results.