A105441 -id:A105441 - cp's OEIS frontend

A001227 Number of odd divisors of n.

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

Offset: 1

N. J. A. Sloane

Also (1) number of ways to write n as difference of two triangular numbers (A000217), see A136107; (2) number of ways to arrange n identical objects in a trapezoid. - Tom Verhoeff
Also number of partitions of n into consecutive positive integers including the trivial partition of length 1 (e.g., 9 = 2+3+4 or 4+5 or 9 so a(9)=3). (Useful for cribbage players.) See A069283. - Henry Bottomley, Apr 13 2000
This has been described as Sylvester's theorem, but to reduce ambiguity I suggest calling it Sylvester's enumeration. - Gus Wiseman, Oct 04 2022
a(n) is also the number of factors in the factorization of the Chebyshev polynomial of the first kind T_n(x). - Yuval Dekel (dekelyuval(AT)hotmail.com), Aug 28 2003
Number of factors in the factorization of the polynomial x^n+1 over the integers. See also A000005. - T. D. Noe, Apr 16 2003
a(n) = 1 if and only if n is a power of 2 (see A000079). - Lekraj Beedassy, Apr 12 2005
Number of occurrences of n in A049777. - Philippe Deléham, Jun 19 2005
For n odd, n is prime if and only if a(n) = 2. - George J. Schaeffer (gschaeff(AT)andrew.cmu.edu), Sep 10 2005
Also number of partitions of n such that if k is the largest part, then each of the parts 1,2,...,k-1 occurs exactly once. Example: a(9)=3 because we have [3,3,2,1],[2,2,2,2,1] and [1,1,1,1,1,1,1,1,1]. - Emeric Deutsch, Mar 07 2006
Also the number of factors of the n-th Lucas polynomial. - T. D. Noe, Mar 09 2006
Lengths of rows of triangle A182469;
Denoted by Delta_0(n) in Glaisher 1907. - Michael Somos, May 17 2013
Also the number of partitions p of n into distinct parts such that max(p) - min(p) < length(p). - Clark Kimberling, Apr 18 2014
Row sums of triangle A247795. - Reinhard Zumkeller, Sep 28 2014
Row sums of triangle A237048. - Omar E. Pol, Oct 24 2014
A069288(n) <= a(n). - Reinhard Zumkeller, Apr 05 2015
A000203, A000593 and this sequence have the same parity: A053866. - Omar E. Pol, May 14 2016
a(n) is equal to the number of ways to write 2*n-1 as (4*x + 2)*y + 4*x + 1 where x and y are nonnegative integers. Also a(n) is equal to the number of distinct values of k such that k/(2*n-1) + k divides (k/(2*n-1))^(k/(2*n-1)) + k, (k/(2*n-1))^k + k/(2*n-1) and k^(k/(2*n-1)) + k/(2*n-1). - Juri-Stepan Gerasimov, May 23 2016, Jul 15 2016
Also the number of odd divisors of n*2^m for m >= 0. - Juri-Stepan Gerasimov, Jul 15 2016
a(n) is odd if and only if n is a square or twice a square. - Juri-Stepan Gerasimov, Jul 17 2016
a(n) is also the number of subparts in the symmetric representation of sigma(n). For more information see A279387 and A237593. - Omar E. Pol, Nov 05 2016
a(n) is also the number of partitions of n into an odd number of equal parts. - Omar E. Pol, May 14 2017 [This follows from the g.f. Sum_{k >= 1} x^k/(1-x^(2*k)). - N. J. A. Sloane, Dec 03 2020]

			G.f. = q + q^2 + 2*q^3 + q^4 + 2*q^5 + 2*q^6 + 2*q^7 + q^8 + 3*q^9 + 2*q^10 + ...
From _Omar E. Pol_, Nov 30 2020: (Start)
For n = 9 there are three odd divisors of 9; they are [1, 3, 9]. On the other hand there are three partitions of 9 into consecutive parts: they are [9], [5, 4] and [4, 3, 2], so a(9) = 3.
Illustration of initial terms:
                              Diagram
   n   a(n)                         _
   1     1                        _|1|
   2     1                      _|1 _|
   3     2                    _|1  |1|
   4     1                  _|1   _| |
   5     2                _|1    |1 _|
   6     2              _|1     _| |1|
   7     2            _|1      |1  | |
   8     1          _|1       _|  _| |
   9     3        _|1        |1  |1 _|
  10     2      _|1         _|   | |1|
  11     2    _|1          |1   _| | |
  12     2   |1            |   |1  | |
...
a(n) is the number of horizontal line segments in the n-th level of the diagram. For more information see A286001. (End)

B. C. Berndt, Ramanujan's Notebooks Part V, Springer-Verlag, see p. 487 Entry 47.
L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 1, p. 306.
J. W. L. Glaisher, On the representations of a number as the sum of two, four, six, eight, ten, and twelve squares, Quart. J. Math. 38 (1907), 1-62 (see p. 4).
Ronald. L. Graham, Donald E. Knuth, and Oren Patashnik, Concrete Mathematics, 2nd ed. (Addison-Wesley, 1994), see exercise 2.30 on p. 65.
P. A. MacMahon, Combinatory Analysis, Cambridge Univ. Press, London and New York, Vol. 1, 1915 and Vol. 2, 1916; see vol. 2, p 28.

N. J. A. Sloane, Table of n, a(n) for n = 1..10000
K. S. Brown's Mathpages, Partitions into Consecutive Integers.
Atli Fannar Franklín, Pattern avoidance enumerated by inversions, arXiv:2410.07467 [math.CO], 2024. See pp. 2, 18.
J. W. L. Glaisher, On the representations of a number as the sum of two, four, six, eight, ten, and twelve squares, Quart. J. Math. 38 (1907), 1-62 (see p. 4 and p. 8).
A. Heiligenbrunner, Sum of adjacent numbers (in German).
Christian Kassel and Christophe Reutenauer, Pairs of intertwined integer sequences, arXiv:2507.15780 [math.NT], 2025. See p. 11.
Gerzson Kéri, The factorization of compressed Chebyshev polynomials and other polynomials related to multiple-angle formulas, Annales Univ. Sci. Budapest (Hungary, 2022) Sect. Comp., Vol. 53, 93-108.
Mircea Merca, Combinatorial interpretations of a recent convolution for the number of divisors of a positive integer, Journal of Number Theory, Volume 160, March 2016, Pages 60-75, function tau_o(n).
M. A. Nyblom, On the representation of the integers as a difference of nonconsecutive triangular numbers, Fibonacci Quarterly 39:3 (2001), pp. 256-263.
R. C. Read, Letter to N. J. A. Sloane, Oct. 29, 1976.
N. J. A. Sloane, Transforms.
T. Verhoeff, Rectangular and Trapezoidal Arrangements, J. Integer Sequences, Vol. 2 (1999), Article 99.1.6.
Eric Weisstein's World of Mathematics, Binomial Number and Odd Divisor Function.
Eric Weisstein's World of Mathematics, q-Polygamma Function.
Index entries for "core" sequences.
Index entries for sequences mentioned by Glaisher.

Cf. A000005, A000079, A000593, A010054 (char. func.), A038547 (positions of first appearances), A050999, A051000, A051001, A051002, A051731, A054844, A069283, A069288, A109814, A115369, A118235, A118236, A125911, A136655, A183063, A183064, A237593, A247795, A272887, A273401, A279387, A286001.
Cf. A000203, A001620, A002162, A053866, A060831.
If this sequence counts gapless sets by sum (by Sylvester's enumeration), these sets are ranked by A073485 and A356956. See also A055932, A066311, A073491, A107428, A137921, A333217, A356224, A356841, A356845.
Dirichlet inverse is A327276.

a001227 = sum . a247795_row
-- Reinhard Zumkeller, Sep 28 2014, May 01 2012, Jul 25 2011

[NumberOfDivisors(n)/Valuation(2*n, 2): n in [1..100]]; // Vincenzo Librandi, Jun 02 2019

for n from 1 by 1 to 100 do s := 0: for d from 1 by 2 to n do if n mod d = 0 then s := s+1: fi: od: print(s); od:
A001227 := proc(n) local a,d;
    a := 1 ;
    for d in ifactors(n)[2] do
        if op(1,d) > 2 then
            a := a*(op(2,d)+1) ;
        end if;
    end do:
    a ;
end proc: # R. J. Mathar, Jun 18 2015

f[n_] := Block[{d = Divisors[n]}, Count[ OddQ[d], True]]; Table[ f[n], {n, 105}] (* Robert G. Wilson v, Aug 27 2004 *)
Table[Total[Mod[Divisors[n], 2]],{n,105}] (* Zak Seidov, Apr 16 2010 *)
f[n_] := Block[{d = DivisorSigma[0, n]}, If[ OddQ@ n, d, d - DivisorSigma[0, n/2]]]; Array[f, 105] (* Robert G. Wilson v *)
a[ n_] := Sum[  Mod[ d, 2], { d, Divisors[ n]}]; (* Michael Somos, May 17 2013 *)
a[ n_] := DivisorSum[ n, Mod[ #, 2] &]; (* Michael Somos, May 17 2013 *)
Count[Divisors[#],?OddQ]&/@Range[110] (* _Harvey P. Dale, Feb 15 2015 *)
(* using a262045 from A262045 to compute a(n) = number of subparts in the symmetric representation of sigma(n) *)
(* cl = current level, cs = current subparts count *)
a001227[n_] := Module[{cs=0, cl=0, i, wL, k}, wL=a262045[n]; k=Length[wL]; For[i=1, i<=k, i++, If[wL[[i]]>cl, cs++; cl++]; If[wL[[i]]Hartmut F. W. Hoft, Dec 16 2016 *)
a[n_] := DivisorSigma[0, n / 2^IntegerExponent[n, 2]]; Array[a, 100] (* Amiram Eldar, Jun 12 2022 *)

{a(n) = sumdiv(n, d, d%2)}; /* Michael Somos, Oct 06 2007 */

{a(n) = direuler( p=2, n, 1 / (1 - X) / (1 - kronecker( 4, p) * X))[n]}; /* Michael Somos, Oct 06 2007 */

a(n)=numdiv(n>>valuation(n,2)) \\ Charles R Greathouse IV, Mar 16 2011

a(n)=sum(k=1,round(solve(x=1,n,x*(x+1)/2-n)),(k^2-k+2*n)%(2*k)==0) \\ Charles R Greathouse IV, May 31 2013

a(n)=sumdivmult(n,d,d%2) \\ Charles R Greathouse IV, Aug 29 2013

from functools import reduce
from operator import mul
from sympy import factorint
def A001227(n): return reduce(mul,(q+1 for p, q in factorint(n).items() if p > 2),1) # Chai Wah Wu, Mar 08 2021

def A001227(n): return len([1 for d in divisors(n) if is_odd(d)])
[A001227(n) for n in (1..80)]  # Peter Luschny, Feb 01 2012

Dirichlet g.f.: zeta(s)^2*(1-1/2^s).
Comment from N. J. A. Sloane, Dec 02 2020: (Start)
By counting the odd divisors f n in different ways, we get three different ways of writing the ordinary generating function.  It is:
A(x) =  x + x^2 + 2*x^3 + x^4 + 2*x^5 + 2*x^6 + 2*x^7 + x^8 + 3*x^9 + 2*x^10 + ...
  = Sum_{k >= 1} x^(2*k-1)/(1-x^(2*k-1))
  = Sum_{k >= 1} x^k/(1-x^(2*k))
  = Sum_{k >= 1} x^(k*(k+1)/2)/(1-x^k)  [Ramanujan, 2nd notebook, p. 355.].
(This incorporates comments from Vladeta Jovovic, Oct 16 2002 and  Michael Somos, Oct 30 2005.) (End)
G.f.: x/(1-x) + Sum_{n>=1} x^(3*n)/(1-x^(2*n)), also L(x)-L(x^2) where L(x) = Sum_{n>=1} x^n/(1-x^n). - Joerg Arndt, Nov 06 2010
a(n) = A000005(n)/(A007814(n)+1) = A000005(n)/A001511(n).
Multiplicative with a(p^e) = 1 if p = 2; e+1 if p > 2. - David W. Wilson, Aug 01 2001
a(n) = A000005(A000265(n)). - Lekraj Beedassy, Jan 07 2005
Moebius transform is period 2 sequence [1, 0, ...] = A000035, which means a(n) is the Dirichlet convolution of A000035 and A057427.
a(n) = A113414(2*n). - N. J. A. Sloane, Jan 24 2006 (corrected Nov 10 2007)
a(n) = A001826(n) + A001842(n). - Reinhard Zumkeller, Apr 18 2006
Sequence = M*V = A115369 * A000005, where M = an infinite lower triangular matrix and V = A000005, d(n); as a vector: [1, 2, 2, 3, 2, 4, ...]. - Gary W. Adamson, Apr 15 2007
Equals A051731 * [1,0,1,0,1,...]; where A051731 is the inverse Mobius transform. - Gary W. Adamson, Nov 06 2007
a(n) = A000005(n) - A183063(n).
a(n) = d(n) if n is odd, or d(n) - d(n/2) if n is even, where d(n) is the number of divisors of n (A000005). (See the Weisstein page.) - Gary W. Adamson, Mar 15 2011
Dirichlet convolution of A000005 and A154955 (interpreted as a flat sequence). - R. J. Mathar, Jun 28 2011
a(A000079(n)) = 1; a(A057716(n)) > 1; a(A093641(n)) <= 2; a(A038550(n)) = 2; a(A105441(n)) > 2; a(A072502(n)) = 3. - Reinhard Zumkeller, May 01 2012
a(n) = 1 + A069283(n). - R. J. Mathar, Jun 18 2015
a(A002110(n)/2) = n, n >= 1. - Altug Alkan, Sep 29 2015
a(n*2^m) = a(n*2^i), a((2*j+1)^n) = n+1 for m >= 0, i >= 0 and j >= 0. a((2*x+1)^n) = a((2*y+1)^n) for positive x and y. - Juri-Stepan Gerasimov, Jul 17 2016
Conjectures: a(n) = A067742(n) + 2*A131576(n) = A082647(n) + A131576(n). - Omar E. Pol, Feb 15 2017
a(n) = A000005(2n) - A000005(n) = A099777(n)-A000005(n). - Danny Rorabaugh, Oct 03 2017
L.g.f.: -log(Product_{k>=1} (1 - x^(2*k-1))^(1/(2*k-1))) = Sum_{n>=1} a(n)*x^n/n. - Ilya Gutkovskiy, Jul 30 2018
G.f.: (psi_{q^2}(1/2) + log(1-q^2))/log(q), where psi_q(z) is the q-digamma function. - Michael Somos, Jun 01 2019
a(n) = A003056(n) - A238005(n). - Omar E. Pol, Sep 12 2021
Sum_{k=1..n} a(k) ~ n*log(n)/2 + (gamma + log(2)/2 - 1/2)*n, where gamma is Euler's constant (A001620). - Amiram Eldar, Nov 27 2022
Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k)/A000005(k) = log(2) (A002162). - Amiram Eldar, Mar 01 2023
a(n) = Sum_{i=1..n} (-1)^(i+1)*A135539(n,i). - Ridouane Oudra, Apr 13 2023

A093641 Numbers of form 2^i * prime(j), i>=0, j>0, together with 1.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 10, 11, 12, 13, 14, 16, 17, 19, 20, 22, 23, 24, 26, 28, 29, 31, 32, 34, 37, 38, 40, 41, 43, 44, 46, 47, 48, 52, 53, 56, 58, 59, 61, 62, 64, 67, 68, 71, 73, 74, 76, 79, 80, 82, 83, 86, 88, 89, 92, 94, 96, 97, 101, 103, 104, 106, 107, 109, 112

Offset: 1

Reinhard Zumkeller, Apr 07 2004

a(n) is either 1, prime, or of form 2a(m), m

1 and Heinz numbers of hook integer partitions. The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k). A hook is a partition of the form (n,1,1,...,1). - Gus Wiseman, Sep 15 2018

Numbers whose odd part is noncomposite. - Peter Munn, Aug 06 2020

			55 is not a member, as 5*11 is not of the form 2^i * prime.

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000

A093640(a(n)) = A000005(a(n)); A000040 and A000079 are subsequences.
A105440 is a subsequence, see also A105442. - Reinhard Zumkeller, Apr 09 2005
Cf. A078822, A007088.
Complement of A105441; A001221(a(n))<=2; A005087(a(n))<=1; A087436(a(n))<=1.
See also A105442.
Union of A038550 and A000079, see also A008578.
Cf. A082733, A153452, A296188, A296561, A300121, A304438, A305940, A317554.
Cf. A000265 (odd part), A008578 (noncomposite).

a093641 n = a093641_list !! (n-1)
a093641_list = filter ((<= 2) . a001227) [1..]
-- Reinhard Zumkeller, May 01 2012

hookQ[n_]:=MatchQ[DeleteCases[FactorInteger[n],{2,}],{}|{{,1}}];
Select[Range[100],hookQ] (* Gus Wiseman, Sep 15 2018 *)

upTo(lim)=my(v=List([1])); for(e=0, log(lim)\log(2), forprime(p=2, lim>>e, listput(v,p<Charles R Greathouse IV, Aug 21 2011

isok(m) = my(k=m/2^valuation(m,2)); (k == 1) || isprime(k); \\ Michel Marcus, Mar 16 2023

from sympy import primepi
def A093641(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n-1+x-sum(primepi(x>>i) for i in range(x.bit_length()))
    return bisection(f,n,n) # Chai Wah Wu, Feb 02 2025

A001227(a(n)) <= 2. - Reinhard Zumkeller, May 01 2012

Number A(x) of a(n) not exceeding x equals 1 + pi(x) + pi(x/2) + pi(x/4) + ..., where pi(x) is the number of primes <= x. If x goes to infinity, A(x)~2*x/log(x) and a(n)~n*log(n)/2 (n-->infinity). - Vladimir Shevelev, Feb 06 2014

A325164 Heinz numbers of integer partitions with Durfee square of length 2.

Original entry on oeis.org

9, 15, 18, 21, 25, 27, 30, 33, 35, 36, 39, 42, 45, 49, 50, 51, 54, 55, 57, 60, 63, 65, 66, 69, 70, 72, 75, 77, 78, 81, 84, 85, 87, 90, 91, 93, 95, 98, 99, 100, 102, 105, 108, 110, 111, 114, 115, 117, 119, 120, 121, 123, 126, 129, 130, 132, 133, 135, 138, 140

Offset: 1

Gus Wiseman, Apr 05 2019

nonn

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
Also positions of 2 in A257990.
First differs from A105441 in lacking 125.
The Durfee length 1 case is A093641. The enumeration of Durfee length 2 partitions by sum is given by A006918, while that of Durfee length 3 partitions is given by A117485.

			The sequence of terms together with their prime indices begins:
   9: {2,2}
  15: {2,3}
  18: {1,2,2}
  21: {2,4}
  25: {3,3}
  27: {2,2,2}
  30: {1,2,3}
  33: {2,5}
  35: {3,4}
  36: {1,1,2,2}
  39: {2,6}
  42: {1,2,4}
  45: {2,2,3}
  49: {4,4}
  50: {1,3,3}
  51: {2,7}
  54: {1,2,2,2}
  55: {3,5}
  57: {2,8}
  60: {1,1,2,3}

Gus Wiseman, Young diagrams corresponding to the first 36 terms.

Cf. A006918, A056239, A093641, A112798, A115994, A117485, A252464, A257990, A325163, A325170.

durf[n_]:=Length[Select[Range[PrimeOmega[n]],Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]][[#]]>=#&]];
Select[Range[100],durf[#]==2&]

A000094 Number of trees of diameter 4.

Original entry on oeis.org

0, 0, 0, 0, 1, 2, 5, 8, 14, 21, 32, 45, 65, 88, 121, 161, 215, 280, 367, 471, 607, 771, 980, 1232, 1551, 1933, 2410, 2983, 3690, 4536, 5574, 6811, 8317, 10110, 12276, 14848, 17941, 21600, 25977, 31146, 37298, 44542, 53132, 63218, 75131, 89089

Offset: 1

N. J. A. Sloane

nonn

Number of partitions of n-1 with at least two parts of size 2 or larger. - Franklin T. Adams-Watters, Jan 13 2006
Also equal to the number of partitions p of n-1 such that max(p)-min(p) > 1. Example: a(7)=5 because we have [5,1],[4,2],[4,1,1],[3,2,1] and [3,1,1,1]. - Giovanni Resta, Feb 06 2006
Also number of partitions of n-1 with at least two parts that are smaller than the largest part. Example: a(7)=5 because we have [4,1,1],[3,2,1],[3,1,1,1],[2,2,1,1,1] and [2,1,1,1,1]. - Emeric Deutsch, May 01 2006
Also number of regions of n-1 that do not contain 1 as a part, n >= 2 (cf. A186114, A206437). - Omar E. Pol, Dec 01 2011
Also rank of the last region of n-1 multiplied by -1, n >= 2 (cf. A194447). - Omar E. Pol, Feb 11 2012
Also sum of ranks of the regions of n-1 that contain emergent parts, n >= 2 (cf. A182699). For the definition of "regions of n" see A206437. - Omar E. Pol, Feb 21 2012

			From _Gus Wiseman_, Apr 12 2019: (Start)
The a(5) = 1 through a(9) = 14 partitions of n-1 with at least two parts of size 2 or larger, or non-hooks, are the following. The Heinz numbers of these partitions are given by A105441.
  (22)  (32)   (33)    (43)     (44)
        (221)  (42)    (52)     (53)
               (222)   (322)    (62)
               (321)   (331)    (332)
               (2211)  (421)    (422)
                       (2221)   (431)
                       (3211)   (521)
                       (22111)  (2222)
                                (3221)
                                (3311)
                                (4211)
                                (22211)
                                (32111)
                                (221111)
The a(5) = 1 through a(9) = 14 partitions of n-1 whose maximum part minus minimum part is at least 2 are the following. The Heinz numbers of these partitions are given by A307516.
  (31)  (41)   (42)    (52)     (53)
        (311)  (51)    (61)     (62)
               (321)   (331)    (71)
               (411)   (421)    (422)
               (3111)  (511)    (431)
                       (3211)   (521)
                       (4111)   (611)
                       (31111)  (3221)
                                (3311)
                                (4211)
                                (5111)
                                (32111)
                                (41111)
                                (311111)
The a(5) = 1 through a(9) = 14 partitions of n-1 with at least two parts that are smaller than the largest part are the following. The Heinz numbers of these partitions are given by A307517.
  (211)  (311)   (321)    (322)     (422)
         (2111)  (411)    (421)     (431)
                 (2211)   (511)     (521)
                 (3111)   (3211)    (611)
                 (21111)  (4111)    (3221)
                          (22111)   (3311)
                          (31111)   (4211)
                          (211111)  (5111)
                                    (22211)
                                    (32111)
                                    (41111)
                                    (221111)
                                    (311111)
                                    (2111111)
(End)

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).

Christian G. Bower, Table of n, a(n) for n = 1..500
J. Riordan, Enumeration of trees by height and diameter, IBM J. Res. Dev. 4 (1960), 473-478.
J. Riordan, The enumeration of trees by height and diameter, IBM Journal 4 (1960), 473-478. (Annotated scanned copy)
Miloslav Znojil, Perturbation theory near degenerate exceptional points, arXiv:2008.00479 [math-ph], 2020.
Index entries for sequences related to trees

Cf. A000041, A206437, A034853, A000147 (diameter 5).

Cf. A006918, A083751, A084835, A105441, A115720, A257990, A307516, A307517, A325164.

g:=x/product(1-x^j,j=1..70)-x-x^2/(1-x)^2: gser:=series(g,x=0,48): seq(coeff(gser,x,n),n=1..46); # Emeric Deutsch, May 01 2006
A000094 := proc(n)
    combinat[numbpart](n-1)-n+1 ;
end proc: # R. J. Mathar, May 17 2016

t=Table[PartitionsP[n]-n,{n,0,45}];
ReplacePart[t,0,1]
(* Clark Kimberling, Mar 05 2012 *)
CoefficientList[1/QPochhammer[x]-x/(1-x)^2-1+O[x]^50, x] (* Jean-François Alcover, Feb 04 2016 *)

a(n+1) = A000041(n)-n for n>0. - John W. Layman
G.f.: x/product(1-x^j,j=1..infinity)-x-x^2/(1-x)^2. - Emeric Deutsch, May 01 2006
G.f.: sum(sum(x^(i+j+1)/product(1-x^k, k=i..j), i=1..j-2), j=3..infinity). - Emeric Deutsch, May 01 2006
a(n+1) = Sum_{m=1..n} A083751(m). - Gregory Gerard Wojnar, Oct 13 2020

More terms from Franklin T. Adams-Watters, Jan 13 2006

A093642 Numbers not containing all divisors in their binary representation.

Original entry on oeis.org

9, 15, 18, 21, 25, 27, 30, 33, 35, 36, 39, 42, 45, 49, 50, 51, 54, 57, 60, 63, 65, 66, 69, 70, 72, 75, 77, 78, 81, 84, 85, 87, 90, 91, 93, 95, 98, 99, 100, 102, 105, 108, 110, 111, 114, 115, 117, 119, 120, 121, 123, 125, 126, 129, 130, 132, 133, 135, 138, 140, 141

Offset: 1

Reinhard Zumkeller, Apr 07 2004

			55 is not a member, as the binary representations of 5 ("101") and 11 ("1011") both appear in the binary representation of 55 ("110111").

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

Complement of A123345.
Subsequence of A105441. - Reinhard Zumkeller, Apr 09 2005
Cf. A000005, A078822, A007088, A027750, A093640.

import Data.List (unfoldr, isInfixOf)
a093642 n = a093642_list !! (n-1)
a093642_list = filter
  (\x -> not $ all (`isInfixOf` b x) $ map b $ a027750_row x) [1..] where
  b = unfoldr (\x -> if x == 0 then Nothing else Just $ swap $ divMod x 2)
-- Reinhard Zumkeller, Oct 27 2012

q[n_] := !AllTrue[Divisors[n], StringContainsQ[IntegerString[n, 2], IntegerString[#, 2]] &]; Select[Range[150], q] (* Amiram Eldar, Jun 05 2022 *)

from sympy import divisors
def ok(n):
    b = bin(n)[2:]
    return not all(bin(d)[2:] in b for d in divisors(n, generator=True))
print([k for k in range(119) if ok(k)]) # Michael S. Branicky, Jun 05 2022

A093640(a(n)) < A000005(a(n)).

A307517 Numbers with at least two not necessarily distinct prime factors less than the largest prime factor.

Original entry on oeis.org

12, 20, 24, 28, 30, 36, 40, 42, 44, 45, 48, 52, 56, 60, 63, 66, 68, 70, 72, 76, 78, 80, 84, 88, 90, 92, 96, 99, 100, 102, 104, 105, 108, 110, 112, 114, 116, 117, 120, 124, 126, 130, 132, 135, 136, 138, 140, 144, 148, 150, 152, 153, 154, 156, 160, 164, 165, 168

Offset: 1

Gus Wiseman, Apr 12 2019

nonn

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so these are Heinz numbers of integer partitions with at least two not necessarily distinct parts less than the largest part. The enumeration of these partitions by sum is given by A000094.

			The sequence of terms together with their prime indices begins:
   12: {1,1,2}
   20: {1,1,3}
   24: {1,1,1,2}
   28: {1,1,4}
   30: {1,2,3}
   36: {1,1,2,2}
   40: {1,1,1,3}
   42: {1,2,4}
   44: {1,1,5}
   45: {2,2,3}
   48: {1,1,1,1,2}
   52: {1,1,6}
   56: {1,1,1,4}
   60: {1,1,2,3}
   63: {2,2,4}
   66: {1,2,5}
   68: {1,1,7}
   70: {1,3,4}
   72: {1,1,1,2,2}
   76: {1,1,8}

Cf. A000094, A000245, A000961, A001222, A052126, A056239, A061395, A064989, A105441, A112798, A243055, A256617, A307516, A325196.

q:= n-> (l-> add(l[i][2], i=1..nops(l)-1)>1)(sort(ifactors(n)[2])):
select(q, [$1..200])[];  # Alois P. Heinz, Apr 12 2019

Select[Range[100],PrimeOmega[#/Power@@FactorInteger[#][[-1]]]>1&]

A325226 Number of prime factors of n that are less than the largest, counted with multiplicity.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Apr 12 2019

nonn

			The prime factors of 300 are {2,2,3,5,5} of which {2,2,3} are less than the largest, so a(300) = 3.

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

Positions of 0's are A000961. Positions of 1's are A325230. Positions of terms > 1 are A307517.

Cf. A000094, A000245, A001222, A052126, A053585, A061395, A062977, A064989, A071178, A105441, A108951, A256617.

Table[PrimeOmega[n/Power@@FactorInteger[n][[-1]]],{n,100}]

A071178(n) = if(1==n, 0, factor(n)[omega(n), 2]);
A325226(n) = (bigomega(n) - A071178(n)); \\ Antti Karttunen, Nov 17 2019

a(n) = A001222(n/A053585(n)).

a(n) = A001222(n) - A071178(n) = A062977(A108951(n)). - Antti Karttunen, Nov 17 2019

Data section extended up to term a(105) by Antti Karttunen, Nov 17 2019

A069562 Numbers, m, whose odd part (largest odd divisor, A000265(m)) is a nontrivial square.

Original entry on oeis.org

9, 18, 25, 36, 49, 50, 72, 81, 98, 100, 121, 144, 162, 169, 196, 200, 225, 242, 288, 289, 324, 338, 361, 392, 400, 441, 450, 484, 529, 576, 578, 625, 648, 676, 722, 729, 784, 800, 841, 882, 900, 961, 968, 1058, 1089, 1152, 1156, 1225, 1250, 1296, 1352, 1369

Offset: 1

Benoit Cloitre, Apr 18 2002

Previous name: sum(d|n,6d/(2+mu(d))) is odd, where mu(.) is the Moebius function, A008683.
From Peter Munn, Jul 06 2020: (Start)
Numbers that have an odd number of odd nonsquarefree divisors.
[Proof of equivalence to the name, where m denotes a positive integer:
(1) These properties are equivalent: (a) m has an even number of odd squarefree divisors; (b) m has a nontrivial odd part.
(2) These properties are equivalent: (a) m has an odd number of odd divisors; (b) the odd part of m is square.
(3) m satisfies the condition at the start of this comment if and only if (1)(a) and (2)(a) are both true or both false.
(4) The trivial odd part, 1, is a square, so (1)(b) and (2)(b) cannot both be false, which (from (1), (2)) means (1)(a) and (2)(a) cannot both be false.
(5) From (3), (4), m satisfies the condition at the start of this comment if and only if (1)(a) and (2)(a) are true.
(6) m satisfies the condition in the name if and only if (1)(b) and (2)(b) are true, which (from (1), (2)) is equivalent to (1)(a) and (2)(a) being true, and hence from (5), to m satisfying the condition at the start of this comment.]
(End)
Numbers whose sum of non-unitary divisors (A048146) is odd. - Amiram Eldar, Sep 16 2024

			To determine the odd part of 18, remove all factors of 2, leaving 9. 9 is a nontrivial square, so 18 is in the sequence. - _Peter Munn_, Jul 06 2020

David A. Corneth, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Odd part.

A000265, A008683 are used in definitions of this sequence.
Lists of numbers whose odd part satisfies other conditions: A028982 (square), A028983 (nonsquare), A029747 (less than 6), A029750 (less than 8), A036349 (even number of prime factors), A038550 (prime), A070776 U {1} (power of a prime), A072502 (square of a prime), A091067 (has form 4k+3), A091072 (has form 4k+1), A093641 (noncomposite), A105441 (composite), A116451 (greater than 4), A116882 (less than or equal to even part), A116883 (greater than or equal to even part), A122132 (squarefree), A229829 (7-rough), A236206 (11-rough), A260488\{0} (has form 6k+1), A325359 (proper prime power), A335657 (odd number of prime factors), A336101 (prime power).
Cf. A048146, A091476.

Select[Range[1000], (odd = #/2^IntegerExponent[#, 2]) > 1 && IntegerQ @ Sqrt[odd] &] (* Amiram Eldar, Sep 29 2020 *)

upto(n) = { my(res = List()); forstep(i = 3, sqrtint(n), 2, for(j = 0, logint(n\i^2, 2), listput(res, i^2<David A. Corneth, Sep 28 2020

Sum_{n>=1} 1/a(n) = 2 * Sum_{k>=1} 1/(2*k+1)^2 = Pi^2/4 - 2 = A091476 - 2 = 0.467401... - Amiram Eldar, Feb 18 2021

New name from Peter Munn, Jul 06 2020

A106405 Number of odd semiprimes dividing n.

Original entry on oeis.org

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

Offset: 1

Reinhard Zumkeller, May 02 2005

nonn

			a(105) = #{15, 21, 35} = #{3*5, 3*7, 5*7} = 3.

Harvey P. Dale, Table of n, a(n) for n = 1..1000

Cf. A000005, A001227, A046315, A086971, A093641, A105441, A106404.

Table[Count[Divisors[n],?(OddQ[#]&&PrimeOmega[#]==2&)],{n,120}] (* _Harvey P. Dale, May 05 2015 *)
a[n_] := Count[e = FactorInteger[n][[;; , 2]], ?(# > 1 &)] + (o = Length[e])*(o - 1)/2 - If[EvenQ[n], If[e[[1]] > 1, o, o - 1], 0]; Array[a, 100] (* _Amiram Eldar, Jun 30 2022 *)

a(n) = A086971(n) - A106404(n);

a(A046315(n)) = 1; a(A093641(n)) = 0; a(A105441(n)) > 0.

A307516 Numbers whose maximum prime index and minimum prime index differ by more than 1.

Original entry on oeis.org

10, 14, 20, 21, 22, 26, 28, 30, 33, 34, 38, 39, 40, 42, 44, 46, 50, 51, 52, 55, 56, 57, 58, 60, 62, 63, 65, 66, 68, 69, 70, 74, 76, 78, 80, 82, 84, 85, 86, 87, 88, 90, 91, 92, 93, 94, 95, 98, 99, 100, 102, 104, 105, 106, 110, 111, 112, 114, 115, 116, 117, 118

Offset: 1

Gus Wiseman, Apr 12 2019

nonn

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so these are Heinz numbers of integer partitions whose maximum and minimum parts differ by more than 1. The enumeration of these partitions by sum is given by A000094.

Differs from A069900 first at n = 43.

			The sequence of terms together with their prime indices begins:
   10: {1,3}
   14: {1,4}
   20: {1,1,3}
   21: {2,4}
   22: {1,5}
   26: {1,6}
   28: {1,1,4}
   30: {1,2,3}
   33: {2,5}
   34: {1,7}
   38: {1,8}
   39: {2,6}
   40: {1,1,1,3}
   42: {1,2,4}
   44: {1,1,5}
   46: {1,9}
   50: {1,3,3}
   51: {2,7}
   52: {1,1,6}
   55: {3,5}

Positions of numbers > 1 in A243055. Complement of A000961 and A256617.

Cf. A000094, A000245, A001222, A052126, A056239, A061395, A064989, A069900, A105441, A112798, A307517, A325196.

with(numtheory):
q:= n-> (l-> pi(l[-1])-pi(l[1])>1)(sort([factorset(n)[]])):
select(q, [$2..200])[];  # Alois P. Heinz, Apr 12 2019

Select[Range[100],PrimePi[FactorInteger[#][[-1,1]]]-PrimePi[FactorInteger[#][[1,1]]]>1&]

Showing 1-10 of 14 results. Next

cp's OEIS Frontend

A001227 Number of odd divisors of n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

PARI

Python

SageMath

Formula

A093641 Numbers of form 2^i * prime(j), i>=0, j>0, together with 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

PARI

Python

Formula

A325164 Heinz numbers of integer partitions with Durfee square of length 2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A000094 Number of trees of diameter 4.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A093642 Numbers not containing all divisors in their binary representation.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

Python

Formula

A307517 Numbers with at least two not necessarily distinct prime factors less than the largest prime factor.

Views

Author

Keywords

Comments

Examples