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.

Previous Showing 11-20 of 38 results. Next

A322327 a(n) = A005361(n) * A034444(n).

Original entry on oeis.org

1, 2, 2, 4, 2, 4, 2, 6, 4, 4, 2, 8, 2, 4, 4, 8, 2, 8, 2, 8, 4, 4, 2, 12, 4, 4, 6, 8, 2, 8, 2, 10, 4, 4, 4, 16, 2, 4, 4, 12, 2, 8, 2, 8, 8, 4, 2, 16, 4, 8, 4, 8, 2, 12, 4, 12, 4, 4, 2, 16, 2, 4, 8, 12, 4, 8, 2, 8, 4, 8, 2, 24, 2, 4, 8, 8, 4, 8, 2, 16, 8, 4, 2, 16, 4, 4, 4, 12, 2, 16, 4, 8, 4, 4, 4, 20, 2, 8, 8, 16
Offset: 1

Views

Author

Werner Schulte, Dec 03 2018

Keywords

Comments

Conjecture: Let k be some fixed integer and a_k(n) = A005361(n) * k^A001221(n) for n > 0 with 0^0 = 1. Then a_k(n) is multiplicative with a_k(p^e) = k*e for prime p and e > 0. For k = 0 see A000007 (offset 1), for k = 1 see A005361, for k = 2 see this sequence, for k = 3 see A226602 (offset 1), and for k = 4 see A322328.
Dirichlet inverse b(n) [= A355837(n)] is multiplicative with b(p^e) = 2 * (e mod 2) * (-1)^((e+1)/2) for prime p and e > 0.

Crossrefs

Programs

  • Mathematica
    a[n_] := If[n==1, 1, Module[{f = FactorInteger[n]}, 2^Length[f] * Times@@f[[;;,2]]]]; Array[a, 100] (* Amiram Eldar, Dec 03 2018 *)
  • PARI
    a(n) = my(f=factor(n)); vecprod(f[,2])*2^omega(n); \\ Michel Marcus, Dec 04 2018
    
  • PARI
    A322327(n) = factorback(apply(e -> e+e, factor(n)[, 2])); \\ Antti Karttunen, Jul 18 2022
    
  • Python
    from math import prod
    from sympy import factorint
    def A322327(n): return prod(e<<1 for e in factorint(n).values()) # Chai Wah Wu, Dec 26 2022

Formula

Multiplicative with a(p^e) = 2*e for prime p and e > 0.
Dirichlet g. f.: (zeta(s))^2 * zeta(2*s) / zeta(4*s).
Equals Dirichlet convolution of A000005 and A227291.
Sum_{k=1..n} a(k) ~ 15*(log(n) + 2*gamma - 1 + 12*zeta'(2)/Pi^2 - 360*zeta'(4)/Pi^4) * n / Pi^2 + 6*zeta(1/2)^2 * sqrt(n) / Pi^2, where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, Sep 20 2020
a(n) = A005361(n^2) = A286324(n^2). - Amiram Eldar, Dec 09 2023

Extensions

Data section extended up to a(100) by Antti Karttunen, Jul 18 2022

A306071 Decimal expansion of Sum_{n>=1} (-1)^omega(n) phi(n)^2/n^4, where omega(n) is the number of distinct prime factors of n (A001221) and phi is Euler's totient function (A000010).

Original entry on oeis.org

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

Views

Author

Amiram Eldar, Jun 19 2018

Keywords

Comments

The constant A that appears in the asymptotic formulae for the sums of the bi-unitary divisor function (A306069) and the bi-unitary totient function (A306070).
The product in Suryanarayana's 1972 paper has a error that was corrected in his 1975 paper.
The probability that 2 randomly selected numbers will be unitary coprime (i.e. their largest common unitary divisor is 1). - Amiram Eldar, Aug 27 2019

Examples

			0.80733082163620503914...
		

References

  • József Sándor, Dragoslav S. Mitrinovic, Borislav Crstici, Handbook of Number Theory I, Springer Science & Business Media, 2005, page 72.

Crossrefs

Programs

  • Mathematica
    cc = CoefficientList[Series[Log[1 - (p - 1)/(p^2*(p + 1))] /. p -> 1/x, {x, 0, 36}], x]; f = FindSequenceFunction[cc]; digits = 20; A = Exp[NSum[f[n + 1 // Floor]*(PrimeZetaP[n]), {n, 2, Infinity}, NSumTerms -> 16 digits, WorkingPrecision -> 16 digits]]; RealDigits[A, 10, digits][[1]] (* Jean-François Alcover, Jun 19 2018 *)
    $MaxExtraPrecision = 1000; Do[Print[Zeta[2] * Exp[-N[Sum[q = Expand[(2*p^2 - 2*p^3 + p^4)^j]; Sum[PrimeZetaP[Exponent[q[[k]], p]] * Coefficient[q[[k]], p^Exponent[q[[k]], p]], {k, 1, Length[q]}]/j, {j, 1, t}], 120]]], {t, 300, 1000, 100}] (* Vaclav Kotesovec, May 29 2020 *)
  • PARI
    prodeulerrat(1 - (p-1)/(p^2 * (p+1))) \\ Amiram Eldar, Mar 18 2021

Formula

Equals Product_{p prime} (1 - (p-1)/(p^2 * (p+1))).
Equals zeta(2) * Product_{p prime} (1 - 2/p^2 + 2/p^3 - 1/p^4).

Extensions

a(1)-a(20) from Jean-François Alcover, Jun 19 2018
a(20)-a(24) from Jon E. Schoenfield, May 27 2019
More terms from Vaclav Kotesovec, May 29 2020

A348341 a(n) is the number of noninfinitary divisors of n.

Original entry on oeis.org

0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 2, 0, 0, 0, 3, 0, 2, 0, 2, 0, 0, 0, 0, 1, 0, 0, 2, 0, 0, 0, 2, 0, 0, 0, 5, 0, 0, 0, 0, 0, 0, 0, 2, 2, 0, 0, 6, 1, 2, 0, 2, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 2, 3, 0, 0, 0, 2, 0, 0, 0, 4, 0, 0, 2, 2, 0, 0, 0, 6, 3, 0, 0, 4, 0, 0, 0
Offset: 1

Views

Author

Amiram Eldar, Oct 13 2021

Keywords

Examples

			a(4) = 1 since 4 has one noninfinitary divisor, 2.
		

Crossrefs

Programs

  • Mathematica
    a[1] = 0; a[n_] := DivisorSigma[0, n] - Times @@ Flatten[2^DigitCount[#, 2, 1] & /@ FactorInteger[n][[;; , 2]]]; Array[a, 100]
  • PARI
    A348341(n) = (numdiv(n)-factorback(apply(a -> 2^hammingweight(a), factorint(n)[, 2]))); \\ Antti Karttunen, Oct 13 2021

Formula

a(n) = A000005(n) - A037445(n).
a(n) = 0 if and only if the number of divisors of n is a power of 2, (i.e., n is in A036537).
a(n) > 0 if and only if the number of divisors of n is not a power of 2, (i.e., n is in A162643).
Sum_{k=1..n} a(k) ~ c * n * log(n), where c = (1 - 2 * A327576) = 0.266749... . - Amiram Eldar, Dec 09 2022

A365498 Dirichlet g.f.: zeta(s) * Product_{p prime} (1 + 1/p^s - 1/p^(3*s)).

Original entry on oeis.org

1, 2, 2, 2, 2, 4, 2, 1, 2, 4, 2, 4, 2, 4, 4, 1, 2, 4, 2, 4, 4, 4, 2, 2, 2, 4, 1, 4, 2, 8, 2, 1, 4, 4, 4, 4, 2, 4, 4, 2, 2, 8, 2, 4, 4, 4, 2, 2, 2, 4, 4, 4, 2, 2, 4, 2, 4, 4, 2, 8, 2, 4, 4, 1, 4, 8, 2, 4, 4, 8, 2, 2, 2, 4, 4, 4, 4, 8, 2, 2, 1, 4, 2, 8, 4, 4, 4
Offset: 1

Views

Author

Vaclav Kotesovec, Sep 06 2023

Keywords

Comments

The number of unitary divisors of n that are cubefree numbers (A004709). - Amiram Eldar, Sep 06 2023

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := If[e <= 2, 2, 1]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 06 2023 *)
  • PARI
    for(n=1, 100, print1(direuler(p=2, n, 1/(1-X) * (1 + X - X^3))[n], ", "))

Formula

Dirichlet g.f.: zeta(s)^2 * Product_{p prime} (1 - 1/p^(2*s) - 1/p^(3*s) + 1/p^(4*s)).
Let f(s) = Product_{p prime} (1 - 1/p^(2*s) - 1/p^(3*s) + 1/p^(4*s)).
Sum_{k=1..n} a(k) ~ f(1) * n * (log(n) + 2*gamma - 1 + f'(1)/f(1)), where
f(1) = Product_{p prime} (1 - 1/p^2 - 1/p^3 + 1/p^4) = 0.5358961538283379998085026313185459506482223745141452711510108346133288...,
f'(1) = f(1) * Sum_{p prime} (-4 + 3*p + 2*p^2) * log(p) / (1 - p - p^2 + p^4) = f(1) * 1.4525924794451595590371439593828547341482465114411929136723476679...
and gamma is the Euler-Mascheroni constant A001620.
Multiplicative with a(p^e) = 2 if e <= 2, and 1 otherwise. - Amiram Eldar, Sep 06 2023
From Vaclav Kotesovec, Jan 27 2025: (Start)
Following formulas have been conjectured for this sequence by Sequence Machine, with each one giving the first 1000000 terms correctly:
a(n) = A056671(n) * A368885(n).
a(n) = A034444(n) / A368248(n).
a(n) = A158522(n) / A307428(n).
a(n) = A369310(n) / A190867(n).
a(n) = A286324(n) / A368172(n). (End)

A353898 a(n) is the number of divisors of n whose exponents in their prime factorizations are all powers of 2 (A138302).

Original entry on oeis.org

1, 2, 2, 3, 2, 4, 2, 3, 3, 4, 2, 6, 2, 4, 4, 4, 2, 6, 2, 6, 4, 4, 2, 6, 3, 4, 3, 6, 2, 8, 2, 4, 4, 4, 4, 9, 2, 4, 4, 6, 2, 8, 2, 6, 6, 4, 2, 8, 3, 6, 4, 6, 2, 6, 4, 6, 4, 4, 2, 12, 2, 4, 6, 4, 4, 8, 2, 6, 4, 8, 2, 9, 2, 4, 6, 6, 4, 8, 2, 8, 4, 4, 2, 12, 4, 4, 4
Offset: 1

Views

Author

Amiram Eldar, May 10 2022

Keywords

Comments

First differs from A049599 and A282446 at n=32.

Examples

			The divisors of 8 are 1, 2 = 2^1, 4 = 2^2 and 8 = 2^3. 3 of these divisors, 1, 2 and 4, are in A138302. Therefore, a(8) = 3.
		

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := Floor[Log2[e]] + 2; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]

Formula

Multiplicative with a(p^e) = floor(log_2(e)) + 2.
a(n) > 1 for n > 1 and a(n) = 2 if and only if n is a prime.
a(n) = A000005(n) if and only if n is cubefree (A004709).

A331109 The number of dual-Zeckendorf-infinitary divisors of n = Product_{i} p(i)^r(i): divisors d = Product_{i} p(i)^s(i), such that the dual Zeckendorf expansion (A104326) of each s(i) contains only terms that are in the dual Zeckendorf expansion of r(i).

Original entry on oeis.org

1, 2, 2, 2, 2, 4, 2, 4, 2, 4, 2, 4, 2, 4, 4, 4, 2, 4, 2, 4, 4, 4, 2, 8, 2, 4, 4, 4, 2, 8, 2, 4, 4, 4, 4, 4, 2, 4, 4, 8, 2, 8, 2, 4, 4, 4, 2, 8, 2, 4, 4, 4, 2, 8, 4, 8, 4, 4, 2, 8, 2, 4, 4, 8, 4, 8, 2, 4, 4, 8, 2, 8, 2, 4, 4, 4, 4, 8, 2, 8, 4, 4, 2, 8, 4, 4, 4
Offset: 1

Views

Author

Amiram Eldar, Jan 09 2020

Keywords

Comments

Dual-Zeckendorf-infinitary divisors are analogous to infinitary divisors (A077609) with dual Zeckendorf expansion instead of binary expansion.
First differs from A286324 at n = 32.

Examples

			a(32) = 4 since 32 = 2^5 and the dual Zeckendorf expansion of 5 is 110, i.e., its dual Zeckendorf representation is a set with 2 terms: {2, 3}. There are 4 possible exponents of 2: 0, 2, 3 and 5, corresponding to the subsets {}, {2}, {3} and {2, 3}. Thus 32 has 4 dual-Zeckendorf-infinitary divisors: 2^0 = 1, 2^2 = 4, 2^3 = 8, and 2^5 = 32.
		

Crossrefs

Programs

  • Mathematica
    fibTerms[n_] := Module[{k = Ceiling[Log[GoldenRatio, n*Sqrt[5]]], t = n, fr = {}}, While[k > 1, If[t >= Fibonacci[k], AppendTo[fr, 1]; t = t - Fibonacci[k], AppendTo[fr, 0]]; k--]; fr];
    dualZeck[n_] := Module[{v = fibTerms[n]}, nv = Length[v]; i = 1; While[i <= nv - 2, If[v[[i]] == 1 && v[[i + 1]] == 0 && v[[i + 2]] == 0, v[[i]] = 0; v[[i + 1]] = 1; v[[i + 2]] = 1; If[i > 2, i -= 3]]; i++]; i = Position[v, _?(# > 0 &)]; If[i == {}, 1, 2^Total[v[[i[[1, 1]] ;; -1]]]]];
    f[p_, e_] := dualZeck[e]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]

Formula

Multiplicative with a(p^e) = 2^A112310(e).

A362852 The number of divisors of n that are both bi-unitary and exponential.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Amiram Eldar, May 05 2023

Keywords

Comments

First differs from A061704 at n = 128, and from A304327 and abs(A307428) at n = 64.
If e > 0 is the exponent of the highest power of p dividing n (where p is a prime), then for each divisor d of n that is both a bi-unitary and an exponential divisor, the exponent of the highest power of p dividing d is a number k such that k | e but k != e/2.
The least term that is higher than 2 is a(64) = 3.
This sequence is unbounded. E.g., a(2^(2^prime(n))) = prime(n).

Examples

			a(8) = 2 since 8 has 2 divisors that are both bi-unitary and exponential: 2 and 8.
		

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := DivisorSigma[0, e] - If[OddQ[e], 0, 1]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    a(n) = {my(f = factor(n)); prod(i = 1, #f~, numdiv(f[i, 2]) - !(f[i, 2] % 2));}

Formula

Multiplicative with a(p^e) = d(e) if e is odd, and d(e)-1 if e is even, where d(k) is the number of divisors of k (A000005).
a(n) = 1 if and only if n is cubefree (A004709).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p prime} (1 + Sum_{k>=1} (d(k)+(k mod 2)-1)/p^k) = 1.1951330849... .

A340232 a(n) is the least number with exactly 2*n bi-unitary divisors.

Original entry on oeis.org

2, 6, 32, 24, 512, 96, 8192, 120, 131072, 1536, 2097152, 480, 33554432, 24576, 536870912, 840, 8589934592, 7776, 137438953472, 7680, 2199023255552, 6291456, 35184372088832, 3360, 562949953421312, 100663296, 9007199254740992, 122880, 144115188075855872, 124416
Offset: 1

Views

Author

Amiram Eldar, Jan 01 2021

Keywords

Comments

Every integer except 1 has an even number of bi-unitary divisors.

Examples

			a(1) = 2 since 2 is the least number with 2*1 = 2 bi-unitary divisors, 1 and 2.
a(2) = 6 since 6 is the least number with 2*2 = 4 bi-unitary divisors, 1, 2, 3 and 6.
		

Crossrefs

Subsequence of A025487.
Similar sequences: A005179 (all divisors), A038547 (odd divisors), A085629 (coreful divisors), A309181 (nonunitary), A340233 (exponential).

Programs

  • Mathematica
    f[p_, e_] := If[OddQ[e], e + 1, e]; d[1] = 1; d[n_] := Times @@ (f @@@ FactorInteger[n]);  max = 10; s = Table[0, {max}]; c = 0; n = 2;  While[c < max, i = d[n]/2; If[i <= max && s[[i]] == 0, c++; s[[i]] = n]; n++]; s

Formula

A286324(a(n)) = 2*n and A286324(k) != 2*n for all k < a(n).

A361782 Numerators of the harmonic means of the bi-unitary divisors of the positive integers.

Original entry on oeis.org

1, 4, 3, 8, 5, 2, 7, 32, 9, 20, 11, 12, 13, 7, 5, 64, 17, 12, 19, 8, 21, 22, 23, 16, 25, 52, 27, 14, 29, 10, 31, 64, 11, 68, 35, 72, 37, 38, 39, 32, 41, 7, 43, 44, 3, 23, 47, 32, 49, 100, 17, 104, 53, 18, 55, 56, 57, 116, 59, 4, 61, 31, 63, 384, 65, 11, 67, 136
Offset: 1

Views

Author

Amiram Eldar, Mar 24 2023

Keywords

Examples

			Fractions begin with 1, 4/3, 3/2, 8/5, 5/3, 2, 7/4, 32/15, 9/5, 20/9, 11/6, 12/5, ...
		

Crossrefs

Cf. A188999, A222266, A286324, A361783 (denominators).
Similar sequences: A099377, A103339, A361316.

Programs

  • Mathematica
    f[p_, e_] := p^e * If[OddQ[e], (e + 1)*(p - 1)/(p^(e + 1) - 1), e/((p^(e + 1) - 1)/(p - 1) - p^(e/2))]; a[1] = 1; a[n_] := Numerator[Times @@ f @@@ FactorInteger[n]]; Array[a, 100]
  • PARI
    a(n) = {my(f = factor(n), p, e); numerator(n * prod(i = 1, #f~, p = f[i, 1]; e = f[i, 2];  if(e%2, (e + 1)*(p - 1)/(p^(e + 1) - 1), e/((p^(e + 1) - 1)/(p - 1) - p^(e/2))))); }

Formula

a(n) = numerator(n*A286324(n)/A188999(n)).
f(n) = a(n)/A361783(n) is multiplicative with f(p^e) = (e+1)*(p-1)/(p^(e+1)-1) if e is odd, and e/((p^(e+1)-1)/(p-1) - p^(e/2)) if e is even.

A361783 Denominators of the harmonic means of the bi-unitary divisors of the positive integers.

Original entry on oeis.org

1, 3, 2, 5, 3, 1, 4, 15, 5, 9, 6, 5, 7, 3, 2, 27, 9, 5, 10, 3, 8, 9, 12, 5, 13, 21, 10, 5, 15, 3, 16, 21, 4, 27, 12, 25, 19, 15, 14, 9, 21, 2, 22, 15, 1, 9, 24, 9, 25, 39, 6, 35, 27, 5, 18, 15, 20, 45, 30, 1, 31, 12, 20, 119, 21, 3, 34, 45, 8, 9, 36, 25, 37, 57
Offset: 1

Views

Author

Amiram Eldar, Mar 24 2023

Keywords

Crossrefs

Cf. A188999, A222266, A286324, A286325 (positions of 1's), A361782 (numerators).
Similar sequences: A099378, A103340, A361317.

Programs

  • Mathematica
    f[p_, e_] := p^e * If[OddQ[e], (e + 1)*(p - 1)/(p^(e + 1) - 1), e/((p^(e + 1) - 1)/(p - 1) - p^(e/2))]; a[1] = 1; a[n_] := Denominator[Times @@ f @@@ FactorInteger[n]]; Array[a, 100]
  • PARI
    a(n) = {my(f = factor(n), p, e); denominator(n * prod(i = 1, #f~, p = f[i, 1]; e = f[i, 2];  if(e%2, (e + 1)*(p - 1)/(p^(e + 1) - 1), e/((p^(e + 1) - 1)/(p - 1) - p^(e/2))))); }

Formula

a(n) = denominator(n*A286324(n)/A188999(n)).
Previous Showing 11-20 of 38 results. Next