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.

A318652 The 2-adic valuation of A046644(n)/A318512(n) (A318651).

Original entry on oeis.org

0, 1, 0, 3, 0, 2, 0, 4, 0, 2, 0, 4, 0, 2, 0, 7, 0, 2, 0, 4, 0, 2, 0, 5, 0, 2, 0, 4, 0, 2, 0, 8, 0, 2, 0, 4, 0, 2, 0, 5, 0, 2, 0, 4, 0, 2, 0, 8, 0, 2, 0, 4, 0, 2, 0, 5, 0, 2, 0, 4, 0, 2, 0, 10, 0, 2, 0, 4, 0, 2, 0, 6, 0, 2, 0, 4, 0, 2, 0, 8, 0, 2, 0, 4, 0, 2, 0, 5, 0, 2, 0, 4, 0, 2, 0, 9, 0, 2, 0, 4, 0, 2, 0, 5, 0
Offset: 1

Views

Author

Antti Karttunen, Aug 31 2018

Keywords

Crossrefs

Programs

  • PARI
    A318652(n) = valuation(A046644(n)/A318512(n),2); \\ Needs also code from those two respective entries.

Formula

a(n) = A007814(A318651(n)).
a(n) = A046645(n) - A318513(n).

A318512 Denominators (in their lowest terms) of the sequence whose Dirichlet convolution with itself yields squares (A000290), or equally A064549.

Original entry on oeis.org

1, 1, 2, 1, 2, 1, 2, 1, 8, 1, 2, 1, 2, 1, 4, 1, 2, 4, 2, 1, 4, 1, 2, 1, 8, 1, 16, 1, 2, 2, 2, 1, 4, 1, 4, 4, 2, 1, 4, 1, 2, 2, 2, 1, 16, 1, 2, 1, 8, 4, 4, 1, 2, 8, 4, 1, 4, 1, 2, 2, 2, 1, 16, 1, 4, 2, 2, 1, 4, 2, 2, 2, 2, 1, 16, 1, 4, 2, 2, 1, 128, 1, 2, 2, 4, 1, 4, 1, 2, 8, 4, 1, 4, 1, 4, 1, 2, 4, 16, 4, 2, 2, 2, 1, 8
Offset: 1

Views

Author

Antti Karttunen, Aug 30 2018

Keywords

Comments

These are also denominators (in their lowest terms) for the sequence whose Dirichlet convolution with itself yields A064549, n * Product_{primes p|n} p.
From Antti Karttunen, Sep 02 2018: (Start)
Proof for the above claim:
This sequence is defined as the denominator (given in the lowest terms) of rational valued function r(1) = 1, r(n) = (1/2) * (A000290(n) - Sum_{d|n, d>1, d 1. Define sequence Ay(n) as the denominator of function s(n), with otherwise similar definition, but with A064549 in place of A000290. Let Ay(n) be the denominator of s(n), reduced also into the lowest terms. (Corresponding numerators are A318649 and A318511 respectively. Note that the denominators in both cases must always be of the form 2^k, with k >= 0).
By applying the distributive property of Dirichlet Convolution [which says that for any completely multiplicative function f, it doesn't matter whether one multiplies the result of convolution afterwards, or whether one multiplies the operands separately before convolution: f(g * g) = (fg) * (fg)], with A000027 in the role of f in both cases, one obtains a pair of equations:
A318649(n) A318681(n) n*A299149(n)
---------- = ---------- = ------------
and
A318511(n) A318680(n) n*A318653(n)
---------- = ---------- = ------------
Ay(n) A299150(n) A299150(n)
where the leftmost ratios are reduced into their lowest terms.
Sequence A318656 gives the 2-adic valuation of ratio A318649(n)/A318512(n), and because there are no even terms neither in A299149 nor in A318653, it also gives the 2-adic valuation of the latter ratio. As A318511/Ay is given in the lowest terms (not both of A318511(n) and Ay(n) can be even at same n), this implies that Ay must indeed be identical to A318512, and furthermore that A318655(n) = A007814(A318649(n)) = A007814(A318511(n)).
(End)

Crossrefs

Programs

  • Mathematica
    f[1] = 1; f[n_] := f[n] = 1/2 (n*Times @@ FactorInteger[n][[All, 1]] - Sum[f[d] f[n/d], {d, Divisors[n][[2 ;; -2]]}]); Table[Denominator[f[n]], {n, 1, 100}] (* Vaclav Kotesovec, May 10 2025 *)
  • PARI
    up_to = 65537;
    A064549(n) = { my(f=factor(n)); for (i=1, #f~, f[i, 2]++); factorback(f); };
    DirSqrt(v) = {my(n=#v, u=vector(n)); u[1]=1; for(n=2, n, u[n]=(v[n]/v[1] - sumdiv(n, d, if(d>1&&dA064549(n)));
    A318512(n) = denominator(v318511_12[n]);
    
  • PARI
    for(n=1, 100, print1(denominator(direuler(p=2, n, 1/(1-p^2*X)^(1/2))[n]), ", ")) \\ Vaclav Kotesovec, May 09 2025

Formula

a(n) = denominator of f(n), where f(1) = 1, f(n) = (1/2) * (A000290(n) - Sum_{d|n, d>1, d 1. [Equally, one could use A064549 in place of A000290.]
a(n) = 2^A318513(n).
a(n) = A046644(n)/A318651(n).
a(2n-1) = A046644(2n-1) = A318658(2n-1), for all n >= 1.

Extensions

The main definition changed, more formulas added by Antti Karttunen, Aug 31 2018

A318649 Numerators of the sequence whose Dirichlet convolution with itself yields squares, A000290.

Original entry on oeis.org

1, 2, 9, 6, 25, 9, 49, 20, 243, 25, 121, 27, 169, 49, 225, 70, 289, 243, 361, 75, 441, 121, 529, 90, 1875, 169, 3645, 147, 841, 225, 961, 252, 1089, 289, 1225, 729, 1369, 361, 1521, 250, 1681, 441, 1849, 363, 6075, 529, 2209, 315, 7203, 1875, 2601, 507, 2809, 3645, 3025, 490, 3249, 841, 3481, 675, 3721, 961, 11907, 924, 4225, 1089
Offset: 1

Views

Author

Antti Karttunen, Aug 31 2018

Keywords

Crossrefs

Cf. A000290, A318512 (denominators).
Cf. also A046643, A299149, A318511, A318651, A318654 (gives the positions of even terms), A318655 (the 2-adic valuation).

Programs

  • PARI
    up_to = 65537;
    DirSqrt(v) = {my(n=#v, u=vector(n)); u[1]=1; for(n=2, n, u[n]=(v[n]/v[1] - sumdiv(n, d, if(d>1&&dA318649(n) = numerator(v318649_aux[n]);
    
  • PARI
    for(n=1, 100, print1(numerator(direuler(p=2, n, 1/(1-p^2*X)^(1/2))[n]), ", ")) \\ Vaclav Kotesovec, May 09 2025

Formula

a(n) = numerator of f(n), where f(1) = 1, f(n) = (1/2) * ((n^2) - Sum_{d|n, d>1, d 1.
a(n) = n*A318512(n)*A299149(n)/A299150(n).
Sum_{k=1..n} A318649(k) / A318512(k) ~ n^3/(3*sqrt(Pi*log(n))) * (1 + (1 - 3*gamma/2) / (6*log(n))), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, May 09 2025

A318655 The 2-adic valuation of A318649, the numerators of "Dirichlet Square Root" of squares.

Original entry on oeis.org

0, 1, 0, 1, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0
Offset: 1

Views

Author

Antti Karttunen, Sep 01 2018

Keywords

Comments

Probably also the 2-adic valuation of A318511.

Crossrefs

Cf. A318511, A318649, A318651, A318652, A318654 (the positions of nonzero terms).

Programs

Formula

a(n) = A007814(A318649(n)).
It seems that for all n >= 1, a(n) <= A007814(A064549(n)) <= A007814(A000290(n)).
Showing 1-4 of 4 results.