A168512 -id:A168512 - cp's OEIS frontend

A326712 Numbers with a record sum of divisors, weighted by divisor multiplicity (A168512).

1, 2, 3, 4, 6, 8, 12, 16, 18, 20, 24, 30, 32, 36, 48, 60, 72, 84, 90, 96, 108, 120, 144, 168, 180, 210, 216, 240, 288, 300, 324, 336, 360, 420, 480, 504, 540, 576, 600, 660, 672, 720, 840, 900, 960, 1008, 1080, 1200, 1260, 1440, 1560, 1620, 1680, 1800, 1920, 1980, 2016

Offset: 1

Amiram Eldar, Oct 18 2019

nonn

The corresponding record values are 1, 3, 4, 9, 12, 19, 30, 41, 42, 44, 64, 72, 75, 102, 134, 170, 208, 226, 237, 264, 294, ...

			The first values of A168512(n) for n=1..8 are {1, 3, 4, 9, 6, 12, 8, 19}. The record values are 1, 3, 4, 9, 12, 19 for 1, 2, 3, 4, 6, 8. Therefore this sequence begins with 1, 2, 3, 4, 6, 8.

Amiram Eldar, Table of n, a(n) for n = 1..600

Cf. A002093, A168512.

s[n_] := 1 + DivisorSum[n, #*IntegerExponent[n, #] &, # > 1 &]; seq = {}; sm = 0; Do[s1 = s[n]; If[s1 > sm, sm = s1; AppendTo[seq, n]], {n, 1, 100000}]; seq (* after Michael De Vlieger at A168512 *)

A326807 Numbers m such that s(m)/m > s(k)/k for all k < m, where s(m) = A168512(m) is the sum of divisors of m, weighted by divisor multiplicity.

Original entry on oeis.org

1, 2, 4, 8, 12, 16, 24, 36, 72, 120, 180, 240, 360, 720, 840, 1260, 1680, 2520, 5040, 10080, 15120, 25200, 27720, 55440, 110880, 166320, 277200, 332640, 554400, 665280, 720720, 1441440, 2162160, 3603600, 4324320, 7207200, 8648640, 10810800, 21621600, 36756720, 61261200, 73513440

Offset: 1

Amiram Eldar, Oct 19 2019

nonn

The least number m such that A168512(m)/m > k, for k = 2, 3, ... is 4, 120, 27720, 122522400, ...

Cf. A004394, A168512, A326677, A326712.

s[n_] := 1 + DivisorSum[n, #*IntegerExponent[n, #] &, # > 1 &]; seq = {}; sm = 0; Do[s1 = s[n]/n; If[s1 > sm, sm = s1; AppendTo[seq, n]], {n, 1, 100000000}]; seq (* after Michael De Vlieger at A168512 *)

A286561 Square array A(n,k): A(n,1) = 1, and for k > 1, A(n,k) = the highest exponent e such that k^e divides n, read by descending antidiagonals as A(1,1), A(1,2), A(2,1), etc.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 2, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 3, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 2, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 2, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 1

Offset: 1

Antti Karttunen, May 20 2017

			The top left 18 X 18 corner of the array:
  n \k 1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18
     .-----------------------------------------------------
   1 | 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
   2 | 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
   3 | 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
   4 | 1, 2, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
   5 | 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
   6 | 1, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
   7 | 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
   8 | 1, 3, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
   9 | 1, 0, 2, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0
  10 | 1, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0
  11 | 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0
  12 | 1, 2, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0
  13 | 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0
  14 | 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0
  15 | 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0
  16 | 1, 4, 0, 2, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0
  17 | 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0
  18 | 1, 1, 2, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1
  ---------------------------------------------------------
A(18,2) = 1, because 2^1 divides 18, but 2^2 does not. A(18,3) = 2, because 3^2 divides 18 (but 3^3 does not). A(18,4) = 0, because 4^0 (= 1) divides 18, but 4^1 does not. A(18,18) = 1, because 18^1 divides 18, but 18^2 does not.
A(2,18) = 0, because 18^0 divides 2, but 18^1 does not.

Antti Karttunen, Table of n, a(n) for n = 1..10585; the first 145 antidiagonals of array

Cf. A286562 (transpose), A286563 (lower triangular region), A286564 (lower triangular region reversed).
Cf. A169594 (row sums), also A168512, A178638, A186643.
Cf. also array A286156.

Table[Function[m, If[k == 1, 1, IntegerExponent[m, k]]][n - k + 1], {n, 15}, {k, n}] // TableForm (* Michael De Vlieger, May 20 2017 *)

A286561(n,k) = if(1==k, 1, valuation(n, k)); \\ Antti Karttunen, May 27 2017

def a(n, k):
    i=1
    if k==1: return 1
    while n%(k**i)==0:
        i+=1
    return i-1
for n in range(1, 21): print([a(k, n - k + 1) for k in range(1, n + 1)]) # Indranil Ghosh, May 20 2017

(define (A286561 n) (A286561bi (A002260 n) (A004736 n)))
(define (A286561bi row col) (if (= 1 col) 1 (let loop ((i 1)) (if (not (zero? (modulo row (expt col i)))) (- i 1) (loop (+ 1 i))))))

A169594 Number of divisors of n, counting divisor multiplicity in n.

Original entry on oeis.org

1, 2, 2, 4, 2, 4, 2, 6, 4, 4, 2, 7, 2, 4, 4, 9, 2, 7, 2, 7, 4, 4, 2, 10, 4, 4, 6, 7, 2, 8, 2, 11, 4, 4, 4, 12, 2, 4, 4, 10, 2, 8, 2, 7, 7, 4, 2, 14, 4, 7, 4, 7, 2, 10, 4, 10, 4, 4, 2, 13, 2, 4, 7, 15, 4, 8, 2, 7, 4, 8, 2, 16, 2, 4, 7, 7, 4, 8, 2, 14, 9, 4, 2, 13, 4, 4, 4, 10, 2, 13, 4, 7, 4, 4, 4, 17, 2, 7

Offset: 1

Joseph L. Pe, Dec 02 2009

The multiplicity of a divisor d > 1 in n is defined as the largest power i for which d^i divides n; and for d = 1 it is defined as 1.
a(n) is also the sum of the multiplicities of the divisors of n.
In other words, a(n) = 1 + sum of the highest exponents e_i for which each number k_i in range 2 .. n divide n, as {k_i}^{e_i} | n. For nondivisors of n this exponent e_i is 0, for n itself it is 1. - Antti Karttunen, May 20 2017
From Gus Wiseman, Mar 25 2021: (Start)
Also the number of strict chains of divisors ending with n and having constant (equal) first quotients. The case starting with 1 is A089723. For example, the a(1) = 1 through a(12) = 7 chains are:
  1  2    3    4      5    6    7    8        9      10    11    12
     1|2  1|3  1|4    1|5  1|6  1|7  1|8      1|9    1|10  1|11  1|12
               2|4         2|6       2|8      3|9    2|10        2|12
               1|2|4       3|6       4|8      1|3|9  5|10        3|12
                                     2|4|8                       4|12
                                     1|2|4|8                     6|12
                                                                 3|6|12
(End)
a(n) depends only on the prime signature of n. - David A. Corneth, Mar 28 2021

			The divisors of 8 are 1, 2, 4, 8 of multiplicity 1, 3, 1, 1, respectively. So a(8) = 1 + 3 + 1 + 1 = 6.

Antti Karttunen, Table of n, a(n) for n = 1..10000

Cf. A168512.
Row sums of A286561, A286563 and A286564.
A001055 counts factorizations (strict: A045778, ordered: A074206).
A057567 counts chains of divisors with weakly increasing first quotients.
A067824 counts strict chains of divisors ending with n.
A253249 counts strict chains of divisors.
A334997 counts chains of divisors of n by length.
A342086 counts chains of divisors with strictly increasing first quotients.
A342496 counts partitions with equal first quotients (strict: A342515, ranking: A342522, ordered: A342495).
A342530 counts chains of divisors with distinct first quotients.
Cf. A003238, A003242, A069916, A122651, A309891, A318991, A318992, A325545.
First differences of A078651.

a := n -> ifelse(n < 2, 1, 1 + add(padic:-ordp(n, k), k = 2..n)):
seq(a(n), n = 1..98);  # Peter Luschny, Apr 10 2025

divmult[d_, n_] := Module[{output, i}, If[d == 1, output = 1, If[d == n, output = 1, i = 0; While[Mod[n, d^(i + 1)] == 0, i = i + 1]; output = i]]; output]; dmt0[n_] := Module[{divs, l}, divs = Divisors[n]; l = Length[divs]; Sum[divmult[divs[[i]], n], {i, 1, l}]]; Table[dmt0[i], {i, 1, 40}]
Table[1 + DivisorSum[n, IntegerExponent[n, #] &, # > 1 &], {n, 98}] (* Michael De Vlieger, May 20 2017 *)

A286561(n,k) = { my(i=1); if(1==k, 1, while(!(n%(k^i)), i = i+1); (i-1)); };
A169594(n) = sumdiv(n,d,A286561(n,d)); \\ Antti Karttunen, May 20 2017

a(n) = { if(n == 1, return(1)); my(f = factor(n), u = vecmax(f[, 2]), cf = f, res = numdiv(f) - u + 1); for(i = 2, u, cf[, 2] = f[, 2]\i; res+=numdiv(factorback(cf)) ); res } \\ David A. Corneth, Mar 29 2021

A169594(n) = {my(s=0, k=2); while(k<=n, s+=valuation(n, k); k=k+1); s + 1} \\ Zhuorui He, Aug 28 2025

def a286561(n, k):
    i=1
    if k==1: return 1
    while n%(k**i)==0:
        i+=1
    return i-1
def a(n): return sum([a286561(n, d) for d in divisors(n)]) # Indranil Ghosh, May 20 2017

(define (A169594 n) (add (lambda (k) (A286561bi n k)) 1 n))
;; Implements sum_{i=lowlim..uplim} intfun(i)
(define (add intfun lowlim uplim) (let sumloop ((i lowlim) (res 0)) (cond ((> i uplim) res) (else (sumloop (1+ i) (+ res (intfun i)))))))
;; For A286561bi see A286561. - Antti Karttunen, May 20 2017

From Friedjof Tellkamp, Feb 29 2024: (Start)
a(n) = A309891(n) + 1.
G.f.: x/(1-x) + Sum_{k>=2, j>=1} x^(k^j)/(1-x^(k^j)).
Dirichlet g.f.: zeta(s) * (1 + Sum_{k>=1} (zeta(k*s) - 1)).
Sum_{n>=1} a(n)/n^2 = (7/24) * Pi^2. (End)

Extended by Ray Chandler, Dec 08 2009

A168654 The sum of the proper divisors of n, weighted by divisor multiplicity, equals n.

Original entry on oeis.org

6, 152, 656, 2888, 18632, 36224, 55328384, 1082574464

Offset: 1

Joseph L. Pe, Dec 01 2009

The multiplicity of a divisor d > 1 in n is defined as the largest power i for which d^i divides n; otherwise is defined as 1 if d = 1.
From Ray Chandler, Dec 08 2009: (Start)
Also in the sequence, but not necessarily the next terms,
2^k * p where p = A168512(2^k) is prime:
2^18 * 525529 = 137764274176,
2^25 * 67117859 = 2252101635801088,
2^26 * 134234921 = 9008353057439744,
2^30 * 2147551801 = 2305916187940225024,
2^40 * 2199025372073 = 2417853966368708281499648,
2^50 * 2251799880936649 = 2535301276174804923929356926976,
as well as k = 150, 348, 694, ... (End)

			The proper divisors of 152 are 1, 2, 4, 8, 19, 38, 76 of multiplicity 1, 3, 1, 1, 1, 1, 1 respectively. Since 1*1 + 3*2 + 1*4 + 1* 8 + 1*19 + 1*38 + 1*76 = 152, then 152 belongs to the sequence.

Cf. A168512

(*multiplicity of d in n*) divmult[d_, n_] := Module[{output, i}, If[d == 1, output = 1, If[d == n, output = 1, i = 0; While[Mod[n, d^(i + 1)] == 0, i = i + 1]; output = i]]; output]; (*sum of divisors weighted by divisor multiplicity*) dmt[n_] := Module[{divs, l}, divs = Divisors[n]; l = Length[divs]; Sum[divmult[divs[[i]], n]*divs[[i]], {i, 1, l}]]; (*search for sequence terms*) ls = {}; Do[If[dmt[n] == 2 n, ls = Append[ls, n]], {n, 2, 10^7}]; ls

a(7)-a(8) from Ray Chandler, Dec 08 2009

A171183 Numbers n such that sigmawt(n) = sigmawt(n+1), where sigmawt(n) is the sum of the divisors of n weighted by divisor multiplicity in n.

Original entry on oeis.org

14, 957, 1334, 1634, 2402, 2685, 20145, 33998, 42818, 74918, 79826, 79833, 84134, 111506, 122073, 138237, 147454, 166934, 201597, 274533, 289454, 347738, 383594, 416577, 440013, 544334, 605985, 649154, 655005, 1060802, 1642154, 1674513

Offset: 1

Joseph L. Pe, Dec 05 2009

nonn

Ray Chandler, Table of n, a(n) for n=1..200

See A168512 for definition of divisor multiplicity.

divmult[d_, n_] := Module[{output, i}, If[d == 1, output = 1, If[d == n, output = 1, i = 0; While[Mod[n, d^(i + 1)] == 0, i = i + 1]; output = i]]; output]; dmt[n_] := Module[{divs, l}, divs = Divisors[n]; l = Length[divs]; Sum[divmult[divs[[i]], n]*divs[[i]], {i, 1, l}]]; l = {}; Do[If[dmt[i] == dmt[i + 1], l = Append[l, i]], {i, 1, 10^6}]; l

Extended by Ray Chandler, Dec 08 2009

A344783 Numbers k such that 1 + Sum_{i=1..k} floor(k/i)*(2^i) is a prime number.

Original entry on oeis.org

1, 3, 4, 7, 18, 25, 26, 30, 40, 50, 95, 150, 348, 694, 1052, 1222, 1808, 2567, 4917, 5399, 7438, 10720, 12152, 30412, 38313, 53620, 121419, 123523

Offset: 1

Amiram Eldar, May 28 2021

Equivalently, numbers k such that A168512(2^k) is a prime number.
The corresponding primes are 3, 19, 41, 283, 525529, 67117859, 134234921, 2147551801, ...
If k is a term of this sequence then 2^k * A168512(2^k) is a term of A168654 (see Ray Chandler's comment in A168654).

			1 is a term since 1 + Sum_{i=1..1} floor(k/i)*(2^i) = 1 + 2 = 3 is a prime.
3 is a term since 1 + Sum_{i=1..3} floor(k/i)*(2^i) = 1 + 6 + 4 + 8 = 19 is a prime.

Cf. A168512, A168654.

Select[Range[100], PrimeQ[1 + Sum[Floor[#/i]*2^i,{i, 1, #}]] &]

a(27)-a(28) from Michael S. Branicky, Sep 23 2024

cp's OEIS Frontend

A326712 Numbers with a record sum of divisors, weighted by divisor multiplicity (A168512).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A326807 Numbers m such that s(m)/m > s(k)/k for all k < m, where s(m) = A168512(m) is the sum of divisors of m, weighted by divisor multiplicity.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

A286561 Square array A(n,k): A(n,1) = 1, and for k > 1, A(n,k) = the highest exponent e such that k^e divides n, read by descending antidiagonals as A(1,1), A(1,2), A(2,1), etc.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Scheme

A169594 Number of divisors of n, counting divisor multiplicity in n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

PARI

Python

Scheme

Formula

Extensions

A168654 The sum of the proper divisors of n, weighted by divisor multiplicity, equals n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Extensions

A171183 Numbers n such that sigmawt(n) = sigmawt(n+1), where sigmawt(n) is the sum of the divisors of n weighted by divisor multiplicity in n.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Extensions

A344783 Numbers k such that 1 + Sum_{i=1..k} floor(k/i)*(2^i) is a prime number.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Extensions