A325989 -id:A325989 - cp's OEIS frontend

A359411 a(n) is the number of divisors of n that are both infinitary and exponential.

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 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, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1

Offset: 1

Amiram Eldar, Dec 30 2022

First differs from A318672 and A325989 at n = 32.
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 an infinitary and an exponential divisor, the exponent of the highest power of p dividing d is a number k such that k | e and the bitwise AND of e and k is equal to k.
The least term that is higher than 2 is a(216) = 4.
The position of the first appearance of a prime p in this sequence is 2^A359081(p), if A359081(p) > -1. E.g., 2^39 = 549755813888 for p = 3, 2^175 = 4.789...*10^52 for p = 5, and 2^1275 = 6.504...*10^383 for p = 7.
This sequence is unbounded since A246600 is unbounded (see A359082).

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

Cf. A037445, A049419, A077609, A080948, A138302, A246600, A322791, A359081, A359082, A359412.

Cf. A318672, A325989.

s[n_] := DivisorSum[n, 1 &, BitAnd[n, #] == # &]; f[p_, e_] := s[e]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]

s(n) = sumdiv(n, d, bitand(d, n)==d);
a(n) = {my(f = factor(n)); prod(i = 1, #f~, s(f[i,2]));}

from math import prod
from sympy import divisors, factorint
def A359411(n): return prod(sum(1 for d in divisors(e,generator=True) if e|d == e) for e in factorint(n).values()) # Chai Wah Wu, Sep 01 2023

Multiplicative with a(p^e) = A246600(e).
a(n) = 1 if and only if n is in A138302.
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p prime} (1 + Sum_{k>=1} A246600(k)/p^k) = 1.135514937... .

A325988 Number of covering (or complete) factorizations of n.

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, 4, 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, 5, 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

Gus Wiseman, May 30 2019

nonn

First differs from A072911 at a(64) = 5, A072911(64) = 4.

A covering factorization of n is an orderless factorization of n into factors > 1 such that every divisor of n is the product of some submultiset of the factors.

			The a(64) = 5 factorizations:
  (2*2*2*2*2*2)
  (2*2*2*2*4)
  (2*2*2*8)
  (2*2*4*4)
  (2*4*8)
The a(96) = 4 factorizations:
  (2*2*2*2*2*3)
  (2*2*2*3*4)
  (2*2*3*8)
  (2*3*4*4)

Cf. A002033, A103295, A126796, A188431.

Cf. A325684, A325781, A325786, A325788, A325986, A325989, A325990.

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Length[Select[facs[n],Union[Times@@@Subsets[#]]==Divisors[n]&]],{n,100}]

a(2^n) = A126796(n).

A325990 Numbers with more than one perfect factorization.

Original entry on oeis.org

8, 24, 27, 32, 40, 54, 56, 72, 88, 96, 104, 108, 120, 125, 128, 135, 136, 152, 160, 168, 184, 189, 200, 216, 224, 232, 243, 248, 250, 256, 264, 270, 280, 288, 296, 297, 312, 328, 343, 344, 351, 352, 360, 375, 376, 378, 384, 392, 408, 416, 424, 432, 440, 456

Offset: 1

Gus Wiseman, May 30 2019

nonn

First differs from A060476 in lacking 1 and having 432.

A perfect factorization of n is an orderless factorization of n into factors > 1 such that every divisor of n is the product of exactly one submultiset of the factors. This is the intersection of covering (or complete) factorizations (A325988) and knapsack factorizations (A292886).

Positions of terms > 1 in A325989.

Cf. A002033, A292886, A325780, A325781, A325782, A325787, A325789, A325988.

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Select[Range[100],Function[n,Length[Select[facs[n],Sort[Times@@@Union[Subsets[#]]]==Divisors[n]&]]>1]]

cp's OEIS Frontend

A359411 a(n) is the number of divisors of n that are both infinitary and exponential.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A325988 Number of covering (or complete) factorizations of n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A325990 Numbers with more than one perfect factorization.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica