A027750 -id:A027750 - cp's OEIS frontend

A037278 Replace n with concatenation of its divisors.

1, 12, 13, 124, 15, 1236, 17, 1248, 139, 12510, 111, 1234612, 113, 12714, 13515, 124816, 117, 1236918, 119, 12451020, 13721, 121122, 123, 1234681224, 1525, 121326, 13927, 12471428, 129, 12356101530, 131, 12481632, 131133, 121734, 15735, 123469121836, 137

Offset: 1

N. J. A. Sloane

a(n) is the union of A176555(n) for n >= 1 and A176556(n) for n >= 2. See A176553 (numbers m such that concatenations of divisors of m are noncomposites) and A176554 (numbers m such that concatenations of divisors of m are nonprimes). [Jaroslav Krizek, Apr 21 2010]

a(n) is the concatenation of n-th row of the triangle in A027750.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A027750, A037283, A037277, A106708.

a037278 = read . concatMap show . a027750_row :: Integer -> Integer
-- Reinhard Zumkeller, Jul 13 2013, May 01 2012, Aug 07 2011

m=1;
for u=1:34 div=divisors(u); conc=str2num(strrep(num2str(div), ' ', ''));
   sol(m)=conc; m=m+1;
end
sol % Marius A. Burtea, Jun 01 2019

k:=1; sol:=[];
for u in [1..34] do D:=Divisors(u); conc:=D[1];
    for u1 in [2..#D] do a:=#Intseq(conc); a1:=#Intseq(D[u1]); conc:=10^a1*conc+D[u1];end for;
     sol[u]:=conc; k:=k+1;
end for;
sol; // Marius A. Burtea, Jun 01 2019

a[n_] := ToExpression[ StringJoin[ ToString /@ Divisors[n] ] ]; Table[ a[n], {n, 1, 34}] (* Jean-François Alcover, Dec 01 2011 *)
FromDigits[Flatten[IntegerDigits/@Divisors[#]]]&/@Range[40] (* Harvey P. Dale, Nov 09 2012 *)

a(n) = my(s=""); fordiv(n, d, s = concat(s, Str(d))); eval(s); \\ Michel Marcus, Jun 01 2019 and Sep 22 2022

from sympy import divisors
def a(n): return int("".join(str(d) for d in divisors(n)))
print([a(n) for n in range(1, 35)]) # Michael S. Branicky, Dec 31 2020

A134681(n) = A055642(a(n)). - Reinhard Zumkeller, Nov 06 2007

More terms from Erich Friedman

A054973 Number of numbers whose divisors sum to n.

Original entry on oeis.org

1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 2, 1, 1, 1, 0, 0, 2, 0, 1, 0, 0, 0, 3, 0, 0, 0, 1, 0, 1, 2, 2, 0, 0, 0, 1, 0, 1, 1, 1, 0, 3, 0, 1, 0, 0, 0, 3, 0, 0, 0, 0, 0, 2, 0, 2, 1, 0, 0, 3, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 5, 0, 1, 0, 0, 0, 1, 0, 2, 0, 0, 0, 3, 0, 0, 0, 0, 0, 3, 1, 0, 1, 0, 0, 4, 0

Offset: 1

Henry Bottomley, May 16 2000

nonn

a(n) = frequency of values n in A000203(m), where A000203(m) = sum of divisors of m. a(n) >= 1 for such n that A175192(n) = 1, a(n) >= 1 if A000203(m) = n for any m. a(n) = 0 for such n that A175192(n) = 0, a(n) = 0 if A000203(m) = n has no solution. - Jaroslav Krizek, Mar 01 2010
First occurrence of k: 2, 1, 12, 24, 96, 72, ..., = A007368. - Robert G. Wilson v, May 14 2014
a(n) is also the number of positive terms in the n-th row of triangle A299762. - Omar E. Pol, Mar 14 2018
Also the number of integer partitions of n whose parts form the set of divisors of some number (necessarily the greatest part). The Heinz numbers of these partitions are given by A371283. For example, the a(24) = 3 partitions are: (23,1), (15,5,3,1), (14,7,2,1). - Gus Wiseman, Mar 22 2024

			a(12) = 2 since 11 has factors 1 and 11 with 1 + 11 = 12 and 6 has factors 1, 2, 3 and 6 with 1 + 2 + 3 + 6 = 12.

T. D. Noe, Table of n, a(n) for n = 1..10000
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems (invphi.gp).

Cf. A000203 (sum-of-divisors function).
For partial sums see A074753.
Cf. A002191, A007609, A308039.
The non-strict version is A371284, ranks A371288.
These partitions have ranks A371283, unsorted version A275700.
A000005 counts divisors, row-lengths of A027750.
A000041 counts integer partitions, strict A000009.
Cf. A001221, A002865, A008289, A371286, A371285.

nn = 105; t = Table[0, {nn}]; k = 1; While[k < 6 nn^(3/2)/Pi^2, d = DivisorSigma[1, k]; If[d < nn + 1, t[[d]]++]; k++]; t (* Robert G. Wilson v, May 14 2014 *)
Table[Length[Select[IntegerPartitions[n],#==Reverse[Divisors[Max@@#]]&]],{n,30}] (* Gus Wiseman, Mar 22 2024 *)

a(n)=v = vector(0); for (i = 1, n, if (sigma(i) == n, v = concat(v, i));); #v; \\ Michel Marcus, Oct 22 2013

a(n)=sum(k=1,n,sigma(k)==n) \\ Charles R Greathouse IV, Nov 12 2013

first(n)=my(v=vector(n),t); for(k=1,n, t=sigma(n); if(t<=n, v[t]++)); v \\ Charles R Greathouse IV, Mar 08 2017

A054973(n)=#invsigma(n) \\ See Alekseyev link for invsigma(). - M. F. Hasler, Nov 21 2019

Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = A308039. - Amiram Eldar, Dec 23 2024

Incorrect comment deleted by M. F. Hasler, Nov 21 2019

A336811 Irregular triangle read by rows T(n,k) in which the length of row n equals the partition number A000041(n-1) and every column k gives the positive integers A000027, with n >= 1 and k >= 1.

Original entry on oeis.org

1, 2, 3, 1, 4, 2, 1, 5, 3, 2, 1, 1, 6, 4, 3, 2, 2, 1, 1, 7, 5, 4, 3, 3, 2, 2, 1, 1, 1, 1, 8, 6, 5, 4, 4, 3, 3, 2, 2, 2, 2, 1, 1, 1, 1, 9, 7, 6, 5, 5, 4, 4, 3, 3, 3, 3, 2, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1, 10, 8, 7, 6, 6, 5, 5, 4, 4, 4, 4, 3, 3, 3, 3, 2, 2, 2, 2, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1

Offset: 1

Omar E. Pol, Nov 20 2020

In other words: row n lists A028310(n-1) blocks where the m-th block consists of A187219(m) copies of n - m + [m=1], with n >= 1 and m >= 1, where [] is the Iverson bracket. [Corrected by Paolo Xausa, Feb 10 2023]
All divisors of all terms in row n are also all parts in the last section of the set of partitions of n.
Thus all divisors of all terms of the first n rows of triangle are also all parts of all partitions of n. In other words: all divisors of the first A000070(n-1) terms of the sequence are also all parts of all partitions of n. - Omar E. Pol, Jun 19 2021
From Omar E. Pol, Jul 31 2021: (Start)
The number of k's in row n is equal to A002865(n-k), 1 <= k <= n.
The number of terms >= k in row n is equal to A000041(n-k), 1 <= k <= n.
The number of k's in the first n rows (or in the first A000070(n-1) terms of the sequence) is equal to A000041(n-k), 1 <= k <= n.
The number of terms >= k in the first n rows (or in the first A000070(n-1) terms of the sequence) is equal to A000070(n-k), 1 <= k <= n.
First n rows of triangle (or first A000070(n-1) terms of the sequence) in nonincreasing order give the n-th row of A176206. (End)

			Triangle begins:
1;
2;
3, 1;
4, 2, 1;
5, 3, 2, 1, 1;
6, 4, 3, 2, 2, 1, 1;
7, 5, 4, 3, 3, 2, 2, 1, 1, 1, 1;
8, 6, 5, 4, 4, 3, 3, 2, 2, 2, 2, 1, 1, 1, 1;
9, 7, 6, 5, 5, 4, 4, 3, 3, 3, 3, 2, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1;
...
For n = 6, by definition the length of row 6 is A000041(6-1) = A000041(5) = 7, so the row 6 of triangle has seven terms. Since every column lists the positive integers A000027 so the row 6 is [6, 4, 3, 2, 2, 1, 1].
Then we have that the divisors of the numbers of the 6th row are:
.
6th row of the triangle ---------->   6 4 3 2 2 1 1
                                      3 2 1 1 1
                                      2 1
                                      1
.
There are seven 1's, four 2's, two 3's, one 4 and one 6.
In total there are 7 + 4 + 2 + 1 + 1 = 15 divisors.
On the other hand the last section of the set of the partitions of 6 can be represented in several ways, five of them as shown below:
._ _ _ _ _ _
|_ _ _      |       6    6                  6                       6
|_ _ _|_    |     3 3    3 3              3   3                     3   3
|_ _    |   |     4 2    4 2            4       2                     4     2
|_ _|_ _|_  |   2 2 2    2 2 2        2   2       2                 2 2   2
          | |       1      1                        1                           1
          | |       1        1                        1                       1
          | |       1        1                          1                   1
          | |       1          1                          1               1
          | |       1          1                            1           1
          | |       1            1                            1       1
          |_|       1              1                            1   1
.
   Figure 1.  Figure 2.  Figure 3.        Figure 4.                   Figure 5.
.
In every figure there are seven 1's, four 2's, two 3's, one 4 and one 6, as shown also the 6th row of A182703.
In total there are 7 + 4 + 2 + 1 + 1 = A138137(6) = 15 parts in every figure.
Figure 5 is an arrangement that shows the correspondence between divisors and parts since the columns give the divisors of the terms of 6th row of triangle.
Finally we can see that all divisors of all numbers in the 6th row of the triangle are the same positive integers as all parts in the last section of the set of the partitions of 6.
Example edited by _Omar E. Pol_, Aug 10 2021

Paolo Xausa, Table of n, a(n) for n = 1..11732 (rows 1..27 of triangle, flattened).

Row sums give A000070.
Row n has length A000041(n-1).
Every column k gives A000027.
Companion of A176206.
Cf. A000007, A000041, A027750, A028310, A002865, A133735, A135010, A138121, A138137, A182703, A187219, A207378, A221529, A336812, A339278, A340035, A340061, A346741.

A336811[row_]:=Flatten[Table[ConstantArray[row-m,PartitionsP[m]-PartitionsP[m-1]],{m,0,row-1}]];
Array[A336811,10] (* Generates 10 rows *) (* Paolo Xausa, Feb 10 2023 *)

f(n) = numbpart(n-1);
T(n, k) = {if (k > f(n), error("invalid k")); if (k==1, return (n)); my(s=0); while (k <= f(n-1), s++; n--;); 1+s;}
tabf(nn) = {for (n=1, nn, for (k=1, f(n), print1(T(n,k), ", ");); print;);} \\ Michel Marcus, Jan 13 2021

A000013 Definition (1): Number of n-bead binary necklaces with beads of 2 colors where the colors may be swapped but turning over is not allowed.

Original entry on oeis.org

1, 1, 2, 2, 4, 4, 8, 10, 20, 30, 56, 94, 180, 316, 596, 1096, 2068, 3856, 7316, 13798, 26272, 49940, 95420, 182362, 349716, 671092, 1290872, 2485534, 4794088, 9256396, 17896832, 34636834, 67110932, 130150588, 252648992, 490853416, 954444608, 1857283156, 3616828364

Offset: 0

N. J. A. Sloane

Definition (2): Equivalently, number of different output sequences from an n-stage pure cycling shift register when 2 sequences are considered the same if one is the complement of the other.
Definition (3): Also number of different output sequences from an n-stage pure cycling shift register constrained so contents have even weight.
Definition (4): Also number of output sequences from (n-1)-stage shift register which feeds back the mod 2 sum of the contents of the register.
The equivalence of definitions (1) and (2) follows at once from the definitions.
If u is an output sequence of type (2) then its derivative is of type (3) - so (2) and (3) count the same things.
If we have a shift register of type (4), append a new cell which contains the mod 2 sum of the contents to get a shift register of type (3). So (3) and (4) count the same things.
If n is even, a(n) = A000116(n/2). If 2^(n+1)-1 is prime, then a(n) = A128976(n+1), the number of cycles in the digraph of the Lucas-Lehmer operator LL(x) = x^2 - 2 acting on Z/(2^(n+1)-1). - M. F. Hasler, May 19 2007
Also number of 2n-bead balanced binary necklaces that are equivalent to their complements. - Andrew Howroyd, Sep 29 2017

			G.f. = 1 + x + 2*x^2 + 2*x^3 + 4*x^4 + 4*x^5 + 8*x^6 + 10*x^7 + 20*x^8 + ...

S. W. Golomb, Shift-Register Sequences, Holden-Day, San Francisco, 1967, p. 172.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Seiichi Manyama, Table of n, a(n) for n = 0..3334 (first 201 terms from T. D. Noe)
Nicolás Álvarez, Victória Becher, Martín Mereb, Ivo Pajor, and Carlos Miguel Soto, On extremal factors of de Bruijn-like graphs, Univ. Buenos Aires (Argentina 2023). See also arXiv:2308.16257 [math.CO], 2023. See references.
Joerg Arndt, Matters Computational (The Fxtbook), p. 151, p. 408.
Raghav Bhutani and Frederick Saia, Replacement dynamics of binary quadratic forms, arXiv:2508.05816 [math.NT], 2025. See p. 6.
Henry Bottomley, Initial terms of A000011 and A000013
Zachary E. Chin and Isaac L. Chuang, The quantum trajectory sensing problem and its solution, arXiv:2410.00893 [quant-ph], 2024. See p. 19.
N. J. Fine, Classes of periodic sequences, Illinois J. Math., 2 (1958), 285-302.
E. N. Gilbert and J. Riordan, Symmetry types of periodic sequences, Illinois J. Math., 5 (1961), 657-665.
Darij Grinberg and Peter Mao, Necklaces over a group with identity product, arXiv:2405.08937 [math.CO], 2024. See pp. 15, 22.
Karyn McLellan, Periodic coefficients and random Fibonacci sequences, Electronic Journal of Combinatorics, 20(4), 2013, #P32.
Frank Ruskey, Necklaces, Lyndon words, De Bruijn sequences, etc.
Frank Ruskey, Necklaces, Lyndon words, De Bruijn sequences, etc. [Cached copy, with permission, pdf format only]
N. J. A. Sloane, On single-deletion-correcting codes
N. J. A. Sloane, On single-deletion-correcting codes, arXiv:math/0207197 [math.CO], 2002; in Codes and Designs (Columbus, OH, 2000), 273-291, Ohio State Univ. Math. Res. Inst. Publ., 10, de Gruyter, Berlin, 2002.
N. J. A. Sloane, Maple code for this and related sequences
Index entries for sequences related to necklaces

Cf. A000031, A000016, A000116.
Cf. A128976.
Cf. A000010, A027750.
Cf. A385665, A000048.

a000013 0 = 1
a000013 n = sum (zipWith (*)
   (map (a000010 . (* 2)) ds) (map (2 ^) $ reverse ds)) `div` (2 * n)
   where ds = a027750_row n
-- Reinhard Zumkeller, Jul 08 2013

with(numtheory): A000013 := proc(n) local s,d; if n = 0 then RETURN(1) else s := 0; for d in divisors(n) do s := s+(phi(2*d)*2^(n/d))/(2*n); od; RETURN(s); fi; end;

a[n_] := Fold[ #1 + EulerPhi[2#2]2^(n/#2)/(2n) &, 0, Divisors[n]]
a[ n_] := If[ n < 1, Boole[n == 0], DivisorSum[ n, EulerPhi[2 #] 2^(n/#) &] / (2 n)]; (* Michael Somos, Dec 19 2014 *)
mx=40;CoefficientList[Series[1-Sum[EulerPhi[2i] Log[1-2*x^i]/(2i),{i,1,mx}],{x,0,mx}],x] (* Herbert Kociemba, Nov 01 2016 *)

{a(n) = if( n<1, n==0, sumdiv(n, k, eulerphi(2*k) * 2^(n/k)) / (2*n))}; /* Michael Somos, Oct 20 1999 */

from sympy import divisors, totient
def a(n): return 1 if n<1 else sum([totient(2*d)*2**(n//d) for d in divisors(n)])//(2*n) # Indranil Ghosh, Apr 28 2017

a(n) = Sum_{ d divides n } (phi(2*d)*2^(n/d))/(2*n) for n>0. - Michael Somos, Oct 20 1999
G.f.: 1 - Sum_{i>=1} phi(2*i)*log(1-2*x^i)/(2*i). - Herbert Kociemba, Nov 01 2016
From Richard L. Ollerton, May 11 2021: (Start)
For n >= 1:
a(n) = (1/(2*n))*Sum_{k=1..n} phi(2*gcd(n,k))*2^(n/gcd(n,k))/phi(n/gcd(n,k)), where phi = A000010.
a(n) = (1/(2*n))*Sum_{k=1..n} phi(2*n/gcd(n,k))*2^gcd(n,k)/phi(n/gcd(n,k)). (End)
a(n) ~ 2^(n-1)/n. - Cedric Lorand, Apr 24 2022
a(n) = Sum_{k=1..n} A385665(n,k) = Sum_{d|n} A000048(d). - Tilman Piesk, Jul 31 2025

A161906 Triangle read by rows in which row n lists the divisors of n that are <= sqrt(n).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 3, 1, 2, 1, 1, 2, 3, 1, 1, 2, 1, 3, 1, 2, 4, 1, 1, 2, 3, 1, 1, 2, 4, 1, 3, 1, 2, 1, 1, 2, 3, 4, 1, 5, 1, 2, 1, 3, 1, 2, 4, 1, 1, 2, 3, 5, 1, 1, 2, 4, 1, 3, 1, 2, 1, 5, 1, 2, 3, 4, 6, 1, 1, 2, 1, 3, 1, 2, 4, 5, 1, 1, 2, 3, 6, 1, 1, 2, 4, 1, 3, 5, 1, 2, 1, 1, 2, 3

Offset: 1

Omar E. Pol, Jun 27 2009

If we define a divisor d|n to be inferior if d <= n/d, then inferior divisors are counted by A038548 and listed by this sequence. - Gus Wiseman, Mar 08 2021

			Triangle begins:
   1....... 1;
   2....... 1;
   3....... 1;
   4..... 1,2;
   5....... 1;
   6..... 1,2;
   7....... 1;
   8..... 1,2;
   9..... 1,3;
  10..... 1,2;
  11....... 1;
  12... 1,2,3;
  13....... 1;
  14..... 1,2;
  15..... 1,3;
  16... 1,2,4;

Reinhard Zumkeller, Rows n = 1..1000 of triangle, flattened

Initial terms are A000012.
Final terms are A033676.
Row lengths are A038548 (number of inferior divisors).
Row sums are A066839 (sum of inferior divisors).
The prime terms are counted by A063962.
The odd terms are counted by A069288.
Row products are A072499.
Row LCMs are A072504.
The superior version is A161908.
The squarefree terms are counted by A333749.
The prime-power terms are counted by A333750.
The strictly superior version is A341673.
The strictly inferior version is A341674.
A001221 counts prime divisors, with sum A001414.
A000005 counts divisors, listed by A027750 with sum A000203.
A056924 count strictly superior (or strictly inferior divisors).
A207375 lists central divisors.
- Inferior: A217581.
- Superior: A033677, A051283, A059172, A063538, A063539, A070038, A116882, A116883, A341591, A341592, A341593, A341675, A341676.
- Strictly Inferior: A060775, A070039, A333805, A333806, A341596, A341677.
- Strictly Superior: A048098, A064052, A140271, A238535, A341594, A341595, A341642, A341643, A341644, A341645, A341646.
Cf. A000196, A001055, A001248, A006530, A020639, A050320, A068101, A161901.

a161906 n k = a161906_tabf !! (n-1) !! (k-1)
a161906_row n = a161906_tabf !! (n-1)
a161906_tabf = zipWith (\m ds -> takeWhile ((<= m) . (^ 2)) ds)
                       [1..] a027750_tabf'
-- Reinhard Zumkeller, Jun 24 2015, Mar 08 2013

div[n_] := Select[Divisors[n], # <= Sqrt[n] &]; div /@ Range[48] // Flatten (* Amiram Eldar, Nov 13 2020 *)

row(n) = select(x->(x<=sqrt(n)), divisors(n)); \\ Michel Marcus, Nov 13 2020

More terms from Sean A. Irvine, Nov 29 2010

A063655 Smallest semiperimeter of integral rectangle with area n.

Original entry on oeis.org

2, 3, 4, 4, 6, 5, 8, 6, 6, 7, 12, 7, 14, 9, 8, 8, 18, 9, 20, 9, 10, 13, 24, 10, 10, 15, 12, 11, 30, 11, 32, 12, 14, 19, 12, 12, 38, 21, 16, 13, 42, 13, 44, 15, 14, 25, 48, 14, 14, 15, 20, 17, 54, 15, 16, 15, 22, 31, 60, 16, 62, 33, 16, 16, 18, 17, 68, 21, 26

Offset: 1

Floor van Lamoen, Jul 24 2001

nonn

Similar to A027709, which is minimal perimeter of polyomino of n cells, or equivalently, minimal perimeter of rectangle of area at least n and with integer sides. Present sequence is minimal semiperimeter of rectangle with area exactly n and with integer sides. - Winston C. Yang (winston(AT)cs.wisc.edu), Feb 03 2002
Semiperimeter b+d, d >= b, of squarest (smallest d-b) integral rectangle with area bd = n. That is, b = largest divisor of n <= sqrt(n), d = smallest divisor of n >= sqrt(n). a(n) = n+1 iff n is noncomposite (1 or prime). - Daniel Forgues, Nov 22 2009
From Juhani Heino, Feb 05 2019: (Start)
Basis for any thickness "frames" around the minimal area. Perimeter can be thought as the 0-thick frame, it is obviously 2a(n). Thickness 1 is achieved by laying unit tiles around the area, there are 2(a(n)+2) of them. Thickness 2 comes from the second such layer, now there are 4(a(n)+4) and so on. They all depend only on a(n), so they share this structure:
Every n > 1 is included. (For different thicknesses, every integer that can be derived from these with the respective formula. So, the perimeter has every even n > 2.)
For each square n > 1, a(n) = a(n-1).
a(1), a(2) and a(6) are the only unique values - the others appear multiple times.
(End)
Gives a discrete Uncertainty Principle. A complex function on an abelian group of order n and its Discrete Fourier Transform must have at least a(n) nonzero entries between them. This bound is achieved by the indicator function on a subgroup of size closest to sqrt(n). - Oscar Cunningham, Oct 10 2021
Also two times the median divisor of n, where the median of a multiset is either the middle part (for odd length), or the average of the two middle parts (for even length). The version for mean instead of median is A057020/A057021. Other doubled medians of multisets are: A360005 (prime indices), A360457 (distinct prime indices), A360458 (distinct prime factors), A360459 (prime factors), A360460 (prime multiplicities), A360555 (0-prepended differences). - Gus Wiseman, Mar 18 2023

			Since 15 = 1*15 = 3*5 and the 3*5 rectangle gives smallest semiperimeter 8, we have a(15)=8.

T. D. Noe, Table of n, a(n) for n = 1..10000
Roy Meshulam, An uncertainty inequality for finite abelian groups, arXiv:math/0312407 [math.CO], 2003.
Roy Meshulam, An uncertainty inequality for finite abelian groups, European Journal of Combinatorics, 27 (2006) 63-67.

Cf. A027709, A033676, A033677, A092510, A162348.
Positions of odd terms are A139710.
Positions of even terms are A139711.
A000005 counts divisors, listed by A027750.
A000975 counts subsets with integer median.
Cf. A003601, A026424, A057020/A057021, A325347, A359893, A359901.

A063655 := proc(n)
    local i,j;
    for i from floor(sqrt(n)) to 1 by -1 do
        j := floor(n/i) ;
        if i*j = n then
            return i+j;
        end if;
    end do:
end proc:
seq(A063655(n), n=1..80); # Winston C. Yang, Feb 03 2002

Table[d = Divisors[n]; len = Length[d]; If[OddQ[len], 2*Sqrt[n], d[[len/2]] + d[[1 + len/2]]], {n, 100}] (* T. D. Noe, Mar 06 2012 *)
Table[2*Median[Divisors[n]],{n,100}] (* Gus Wiseman, Mar 18 2023 *)

A063655(n) = { my(c=1); fordiv(n,d,if((d*d)>=n,if((d*d)==n,return(2*d),return(c+d))); c=d); (0); }; \\ Antti Karttunen, Oct 20 2017

from sympy import divisors
def A063655(n):
    d = divisors(n)
    l = len(d)
    return d[(l-1)//2] + d[l//2] # Chai Wah Wu, Jun 14 2019

a(n) = A033676(n) + A033677(n).
a(n) = A162348(2n-1) + A162348(2n). - Daniel Forgues, Sep 29 2014
a(n) = Min_{d|n} (n/d + d). - Ridouane Oudra, Mar 17 2024

Corrected and extended by Larry Reeves (larryr(AT)acm.org) and Dean Hickerson, Jul 26 2001

A073093 Number of prime power divisors of n.

Original entry on oeis.org

1, 2, 2, 3, 2, 3, 2, 4, 3, 3, 2, 4, 2, 3, 3, 5, 2, 4, 2, 4, 3, 3, 2, 5, 3, 3, 4, 4, 2, 4, 2, 6, 3, 3, 3, 5, 2, 3, 3, 5, 2, 4, 2, 4, 4, 3, 2, 6, 3, 4, 3, 4, 2, 5, 3, 5, 3, 3, 2, 5, 2, 3, 4, 7, 3, 4, 2, 4, 3, 4, 2, 6, 2, 3, 4, 4, 3, 4, 2, 6, 5, 3, 2, 5, 3, 3, 3, 5, 2, 5, 3, 4, 3, 3, 3, 7, 2, 4, 4, 5, 2, 4, 2, 5, 4

Offset: 1

Reinhard Zumkeller, Aug 24 2002

Also, number of prime divisors of 2n (counted with multiplicity).
A001221(n) < a(n) <= A000005(n) for all n; a(n)=A001221(n)+1 iff n is squarefree (A005117); a(n)=A000005(n) iff n is a prime power (A000961).
a(n) is also the number of kBenoit Cloitre, Oct 13 2002
a(n) is also 1 + the number of divisors of n with omega(d)=1, where omega is A001221. - Enrique Pérez Herrero, Nov 05 2009
Length of n-th row of triangle A210208. - Reinhard Zumkeller, Mar 18 2012
a(n) depends only on the prime signature of n with a(A025487(n)) = 1, 2, 3, 3, 4, 4, 5, 5, 4, 6, 5, 6, 5, 7, 6, 7 ,.. = A036041(n)+1; (n>=1). - R. J. Mathar, May 28 2017

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
T. M. Apostol, Resultants of Cyclotomic Polynomials, Proc. Amer. Math. Soc. 24, 457-462, 1970.
T. M. Apostol, The Resultant of the Cyclotomic Polynomials Fm(ax) and Fn(bx), Math. Comput. 29, 1-6, 1975.
Eric Weisstein's World of Mathematics, Cyclotomic Polynomial

Cf. A000961, A023888, A054372. Bisection of A001222.

a073093 = length . a210208_row  -- Reinhard Zumkeller, Mar 18 2012

[n eq 1 select 1 else &+[p[2]: p in Factorization(n)]+1: n in [1..100]]; // Vincenzo Librandi, Jan 06 2017

seq(numtheory:-bigomega(n)+1, n=1..1000); # Robert Israel, Sep 06 2015

f[n_] := Plus @@ Flatten[ Table[1, {#[[2]]}] & /@ FactorInteger[n]]; Table[ f[2n], {n, 105}] (* Robert G. Wilson v, Dec 23 2004 *)
A001221[n_] := (Length[ FactorInteger[n]]); SetAttributes[A001221, Listable]; A073093[n_]:=Length[Select[A001221[Divisors[n]], # == 1 &]]; (* Enrique Pérez Herrero, Nov 05 2009 *)
PrimeOmega[Range[100]] + 1 (* Paolo Xausa, Nov 23 2024 *)

numlib::Omega (2*n)$ n=1..105 // Zerinvary Lajos, May 13 2008

a(n)=sum(k=1,n,if(1-polresultant(polcyclo(n),polcyclo(k)),1,0))

A073093(n)=bigomega(n)+1   \\ M. F. Hasler, Dec 08 2010

If n = Product (p_j^k_j), a(n) = 1 + Sum (k_j).
a(n) = bigomega(n)+1 = A001222(n)+1 = A001222(2*n).
a(n) = if n=1 then 1 else a(A032742(n)) + 1. - Reinhard Zumkeller, Sep 24 2009
a(n) = max { a(d) ; d 1. - David W. Wilson, Dec 08 2010
a(n) = Sum_{k = 1 .. A001221(n)} A010055(A027750(n,k)). - Reinhard Zumkeller, Mar 18 2012
G.f.: x/(1 - x) + Sum_{k>=2} floor(1/omega(k))*x^k/(1 - x^k), where omega(k) is the number of distinct prime factors (A001221). - Ilya Gutkovskiy, Jan 04 2017

A207375 Irregular array read by rows in which row n lists the (one or two) central divisors of n in increasing order.

Original entry on oeis.org

1, 1, 2, 1, 3, 2, 1, 5, 2, 3, 1, 7, 2, 4, 3, 2, 5, 1, 11, 3, 4, 1, 13, 2, 7, 3, 5, 4, 1, 17, 3, 6, 1, 19, 4, 5, 3, 7, 2, 11, 1, 23, 4, 6, 5, 2, 13, 3, 9, 4, 7, 1, 29, 5, 6, 1, 31, 4, 8, 3, 11, 2, 17, 5, 7, 6, 1, 37, 2, 19, 3, 13, 5, 8, 1, 41, 6, 7, 1, 43

Offset: 1

Omar E. Pol, Feb 23 2012

If n is a square then row n lists only the square root of n because the squares (A000290) have only one central divisor.
If n is not a square then row n lists the pair (j, k) of divisors of n, nearest to the square root of n, such that j*k = n.
Conjecture 1: It appears that the n-th record in this sequence is the last member of row A008578(n).
Column 1 gives A033676. Right border gives A033677. - Omar E. Pol, Feb 26 2019
The conjecture 1 follows from Bertrand's Postulate. - Charles R Greathouse IV, Feb 11 2022
Row products give A097448. - Omar E. Pol, Feb 17 2022

			Array begins:
  1;
  1,  2;
  1,  3;
  2;
  1,  5;
  2,  3;
  1,  7;
  2,  4;
  3;
  2,  5;
  1, 11;
  3,  4;
  1, 13;
...

Alois P. Heinz, Rows n = 1..5000
Omar E. Pol, Illustration of the divisors of the first 12 positive integers

Row n has length A169695(n).
Row sums give A207376.
Cf. A000005, A000290, A008578, A027750, A033676, A033677, A097448, A161901, A161904.

A207375row[n_] := ArrayPad[#, -Floor[(Length[#] - 1)/2]] & [Divisors[n]];
Array[A207375row, 50] (* Paolo Xausa, Apr 07 2025 *)

A007503 Number of subgroups of dihedral group: sigma(n) + d(n).

Original entry on oeis.org

2, 5, 6, 10, 8, 16, 10, 19, 16, 22, 14, 34, 16, 28, 28, 36, 20, 45, 22, 48, 36, 40, 26, 68, 34, 46, 44, 62, 32, 80, 34, 69, 52, 58, 52, 100, 40, 64, 60, 98, 44, 104, 46, 90, 84, 76, 50, 134, 60, 99, 76, 104, 56, 128, 76, 128, 84, 94, 62, 180, 64, 100, 110, 134

Offset: 1

N. J. A. Sloane, Robert G. Wilson v

nonn

Essentially first differences of A257644. - Franklin T. Adams-Watters, Nov 05 2015
Write D_{2n} as , then the subgroups are of the form  for d|n or  for d|n and 0 <= r < d. There are d(n) subgroups of the first type and sigma(n) subgroups of the second type. - Jianing Song, Jul 21 2022
Inverse Möbius transform of n+1 (see Arndt formula). - Wesley Ivan Hurt, Jul 05 2025

			a(4) = 10 since D_8 = <a, x | a^4 = x^2 = 1, x*a*x = a^(-1)> has 10 subgroups. The 6 subgroups {e}, {e,a^2}, {e,a,a^2,a^3}, {e,a^2,x,a^2*x}, {e,a^2,a*x,a^3*x} and D_8 are normal, and the 4 subgroups {e,x}, {e,a*x}, {e,a^2*x} and {e,a^3*x} are not. - _Jianing Song_, Jul 21 2022

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 1..1000
Keith Conrad, Dihedral Groups II
David W. Jensen and Eric R. Bussian, A Number-Theoretic Approach to Counting Subgroups of Dihedral Groups, Two-Year College Math. Jnl., 23 (1992), 150-152.
The Group Properties Wiki, Subgroup structure of dihedral groups
Index entries for sequences related to groups

Cf. A000005, A000203, A037852 (number of normal subgroups).

Cf. A027750, A257644 (cumulative sums, start=1).

a007503 = sum . map (+ 1) . a027750_row'
-- Reinhard Zumkeller, Nov 09 2015

with(numtheory): seq(sigma(n)+tau(n), n=1..56) ; # Zerinvary Lajos, Jun 04 2008

A007503[n_]:=DivisorSum[n,DivisorSigma[0,2^#]&]; Array[A007503,20] (* Enrique Pérez Herrero, Apr 14 2012 *)

a(n) = sumdiv(n,d, d+1 ); \\ Joerg Arndt, Apr 14 2013

G.f.: Sum_{k>=1} 1/(1-x^k)^2. - Benoit Cloitre, Apr 21 2003
G.f.: Sum_{i>=1} (1+i)*x^i/(1-x^i). - Jon Perry, Jul 03 2004
a(n) = Sum_{d|n} tau(p^d), where tau is A000005 and p any prime. - Enrique Pérez Herrero, Apr 14 2012
a(n) = Sum_{d divides n} d+1. - Joerg Arndt, Apr 14 2013
L.g.f.: -log(Product_{k>=1} (1 - x^k)^(1+1/k)) = Sum_{n>=1} a(n)*x^n/n. - Ilya Gutkovskiy, May 23 2018
a(n) = A000005(n) + A000203(n). - Omar E. Pol, Aug 19 2019
a(n) = A348219(n) + A386438(n). - Wesley Ivan Hurt, Jul 21 2025

A060640 If n = Product p_i^e_i then a(n) = Product (1 + 2p_i + 3p_i^2 + ... + (e_i+1)*p_i^e_i).

Original entry on oeis.org

1, 5, 7, 17, 11, 35, 15, 49, 34, 55, 23, 119, 27, 75, 77, 129, 35, 170, 39, 187, 105, 115, 47, 343, 86, 135, 142, 255, 59, 385, 63, 321, 161, 175, 165, 578, 75, 195, 189, 539, 83, 525, 87, 391, 374, 235, 95, 903, 162, 430, 245, 459, 107, 710, 253, 735, 273, 295, 119

Offset: 1

N. J. A. Sloane, Apr 17 2001

Equals row sums of triangle A143313. - Gary W. Adamson, Aug 06 2008
Equals row sums of triangle A127099. - Gary W. Adamson, Jul 27 2008
Sum of the divisors d2 from the ordered pairs of divisors of n, (d1,d2) with d1<=d2, such that d1|d2. - Wesley Ivan Hurt, Mar 22 2022

			a(4) = a(2^2) = 1 + (2)*(2) + (3)*(2^2) = 17;
a(6) = a(2)*a(3) = (1 + (2)*(2))*(1+(2)*(3)) = (5)*(7) = 35.
a(6) = tau(1) + 2*tau(2) + 3*tau(3) + 6*tau(6) = 1 + 2*2 + 3*2 + 6*4 = 35.

D. M. Burton, Elementary Number Theory, Allyn and Bacon Inc., Boston, MA, 1976, p. 120.

Seiichi Manyama, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Harry J. Smith)

Cf. A000005, A000203, A001001, A006171, A038040 (Mobius transform), A049060, A057660, A057723, A327960 (Dirichlet inverse).

Cf. also triangles A027750, A127099, A143313.

a060640 n = sum [d * a000005 d | d <- a027750_row n]
-- Reinhard Zumkeller, Feb 29 2012

A060640 := proc(n) local ans, i, j; ans := 1: for i from 1 to nops(ifactors(n)[2]) do ans := ans*(1+sum((j+1)*ifactors(n)[2][i][1]^j,j=1..ifactors(n)[2][i][2])): od: RETURN(ans) end:

a[n_] := Total[#*DivisorSigma[1, n/#] & /@ Divisors[n]];
a /@ Range[59] (* Jean-François Alcover, May 19 2011, after Vladeta Jovovic *)
f[p_, e_] := ((e + 1)*p^(e + 2) - (e + 2)*p^(e + 1) + 1)/(p - 1)^2; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Apr 10 2022 *)

j=[]; for(n=1,200,j=concat(j,sumdiv(n,d,n/d*sigma(d)))); j

a(n)=if(n<1,0,direuler(p=2,n,1/(1-X)/(1-p*X)^2)[n]) /* Ralf Stephan */

N=66; default(seriesprecision,N); x=z+O(z^(N+1))
c=sum(j=1,N,j*x^j); t=1/prod(j=1,N, eta(x^(j)));
t=log(t);t=serconvol(t,c);
Vec(t) /* Joerg Arndt, May 03 2008 */

{ for (n=1, 1000, write("b060640.txt", n, " ", direuler(p=2, n, 1/(1 - X)/(1 - p*X)^2)[n]); ) } /* Harry J. Smith, Jul 08 2009 */

def A060640(n) :
    sigma = sloane.A000203
    return add(sigma(k)*(n/k) for k in divisors(n))
[A060640(i) for i in (1..59)] # Peter Luschny, Sep 15 2012

a(n) = Sum_{d|n} d*tau(d), where tau(d) is the number of divisors of d, cf. A000005. a(n) = Sum_{d|n} d*sigma(n/d), where sigma(n)=sum of divisors of n, cf. A000203. - Vladeta Jovovic, Apr 23 2001
Multiplicative with a(p^e) = ((e+1)*p^{e+2} - (e+2)*p^{e+1} + 1) / (p-1)^2. Dirichlet g.f.: zeta(s)*zeta(s-1)^2. - Franklin T. Adams-Watters, Aug 03 2006
L.g.f.: Sum(A060640(n)*x^n/n) = -log( Product_{j>=1} P(x^j) ) where P(x) = Product_{k>=1} (1-x^k). - Joerg Arndt, May 03 2008
G.f.: Sum_{k>=1} k*tau(k)*x^k/(1 - x^k). - Ilya Gutkovskiy, Sep 06 2018
Sum_{k=1..n} a(k) ~ n^2/24 * ((4*gamma - 1)*Pi^2 + 2*Pi^2 * log(n) + 12*Zeta'(2)), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, Feb 01 2019

More terms from James Sellers, Vladeta Jovovic and Matthew Conroy, Apr 17 2001

cp's OEIS Frontend

A037278 Replace n with concatenation of its divisors.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

MATLAB

Magma

Mathematica

PARI

Python

Formula

Extensions

A054973 Number of numbers whose divisors sum to n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

PARI

PARI

Formula

Extensions

A336811 Irregular triangle read by rows T(n,k) in which the length of row n equals the partition number A000041(n-1) and every column k gives the positive integers A000027, with n >= 1 and k >= 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

A000013 Definition (1): Number of n-bead binary necklaces with beads of 2 colors where the colors may be swapped but turning over is not allowed.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

A161906 Triangle read by rows in which row n lists the divisors of n that are <= sqrt(n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Extensions

A063655 Smallest semiperimeter of integral rectangle with area n.

Views

Author

Keywords

Comments

Examples

Links

A060640 If n = Product p_i^e_i then a(n) = Product (1 + 2p_i + 3p_i^2 + ... + (e_i+1)*p_i^e_i).