A069289 -id:A069289 - cp's OEIS frontend

A000593 Sum of odd divisors of n.

1, 1, 4, 1, 6, 4, 8, 1, 13, 6, 12, 4, 14, 8, 24, 1, 18, 13, 20, 6, 32, 12, 24, 4, 31, 14, 40, 8, 30, 24, 32, 1, 48, 18, 48, 13, 38, 20, 56, 6, 42, 32, 44, 12, 78, 24, 48, 4, 57, 31, 72, 14, 54, 40, 72, 8, 80, 30, 60, 24, 62, 32, 104, 1, 84, 48, 68, 18, 96, 48, 72, 13, 74, 38, 124

Offset: 1

N. J. A. Sloane

Denoted by Delta(n) or Delta_1(n) in Glaisher 1907. - Michael Somos, May 17 2013
A069289(n) <= a(n). - Reinhard Zumkeller, Apr 05 2015
A000203, A001227 and this sequence have the same parity: A053866. - Omar E. Pol, May 14 2016
For the g.f.s given below by Somos Oct 29 2005, Jovovic, Oct 11 2002 and Arndt, Nov 09 2010, see the Hardy-Wright reference, proof of Theorem 382, p. 312, with x^2 replaced by x. - Wolfdieter Lang, Dec 11 2016
a(n) is also the total number of parts in all partitions of n into an odd number of equal parts. - Omar E. Pol, Jun 04 2017
It seems that a(n) divides A000203(n) for every n. - Ivan N. Ianakiev, Nov 25 2017 [Yes, see the formula dated Dec 14 2017].
Also, alternating row sums of A126988. - Omar E. Pol, Feb 11 2018
Where a(n) shows the number of equivalence classes of Hurwitz quaternions with norm n (equivalence defined by right multiplication with one of the 24 Hurwitz units as in A055672), A046897(n) seems to give the number of equivalence classes of Lipschitz quaternions with norm n (equivalence defined by right multiplication with one of the 8 Lipschitz units). - R. J. Mathar, Aug 03 2025

			G.f. = x + x^2 + 4*x^3 + x^4 + 6*x^5 + 4*x^6 + 8*x^7 + x^8 + 13*x^9 + 6*x^10 + 12*x^11 + ...

Jean-Marie De Koninck and Armel Mercier, 1001 Problèmes en Théorie Classique des Nombres, Problème 496, pp. 69-246, Ellipses, Paris, 2004.
G. H. Hardy, Ramanujan: twelve lectures on subjects suggested by his life and work, Cambridge, University Press, 1940, p. 132.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, Fifth Edition, Clarendon Press, Oxford, 2003, p. 312.
Friedrich Hirzebruch, Thomas Berger, and Rainer Jung, Manifolds and Modular Forms, Vieweg, 1994, p. 133.
John Riordan, Combinatorial Identities, Wiley, 1968, p. 187.
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).

T. D. Noe, Table of n, a(n) for n = 1..10000
Francesca Aicardi, Matricial formulas for partitions, arXiv:0806.1273 [math.NT], 2008.
Michael Baake and Robert V. Moody, Similarity submodules and root systems in four dimensions, arXiv:math/9904028 [math.MG], 1999; Canad. J. Math. 51 (1999), 1258-1276.
John A. Ewell, On the sum-of-divisors function, Fib. Q., 45 (2007), 205-207.
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).
Heekyoung Hahn, Convolution sums of some functions on divisors, arXiv:1507.04426 [math.NT], 2015.
Kaya Lakein and Anne Larsen, A Proof of Merca's Conjectures on Sums of Odd Divisor Functions, arXiv:2107.07637 [math.NT], 2021.
Mircea Merca, The Lambert series factorization theorem, The Ramanujan Journal, January 2017, also here
Mircea Merca, Congruence identities involving sums of odd divisors function, Proceedings of the Romanian Academy, Series A, Volume 22, Number 2/2021, pp. 119-125.
Hossein Movasati and Younes Nikdelan, Gauss-Manin Connection in Disguise: Dwork Family, arXiv:1603.09411 [math.AG], 2016.
Yash Puri and Thomas Ward, Arithmetic and growth of periodic orbits, J. Integer Seqs., Vol. 4 (2001), Article 01.2.1.
N. J. A. Sloane, Transforms.
H. J. Stephen Smith, Report on the Theory of Numbers. — Part VI., Report of the 35 Meeting of the British Association for the Advancement of Science (1866). See p. 336.
Eric Weisstein's World of Mathematics, Odd Divisor Function.
Eric Weisstein's World of Mathematics, Partition Function Q.
Eric Weisstein's World of Mathematics, q-Pochhammer Symbol.
Index entries for "core" sequences.
Index entries for sequences mentioned by Glaisher.

Cf. A000005, A000203, A000265, A001227, A006128, A050999, A051000, A051001, A051002, A065442, A078471 (partial sums), A069289, A247837 (subset of the primes).

Cf. A301799, A301800.

a000593 = sum . a182469_row  -- Reinhard Zumkeller, May 01 2012, Jul 25 2011

m:=50; R:=PowerSeriesRing(Integers(), m); Coefficients(R!( (&+[j*x^j/(1+x^j): j in [1..2*m]])  )); // G. C. Greubel, Nov 07 2018

[&+[d:d in Divisors(n)|IsOdd(d)]:n in [1..75]]; // Marius A. Burtea, Aug 12 2019

A000593 := proc(n) local d,s; s := 0; for d from 1 by 2 to n do if n mod d = 0 then s := s+d; fi; od; RETURN(s); end;

Table[a := Select[Divisors[n], OddQ[ # ]&]; Sum[a[[i]], {i, 1, Length[a]}], {n, 1, 60}] (* Stefan Steinerberger, Apr 01 2006 *)
f[n_] := Plus @@ Select[ Divisors@ n, OddQ]; Array[f, 75] (* Robert G. Wilson v, Jun 19 2011 *)
a[ n_] := If[ n < 1, 0, Sum[ -(-1)^d n / d, {d, Divisors[ n]}]]; (* Michael Somos, May 17 2013 *)
a[ n_] := If[ n < 1, 0, DivisorSum[ n, -(-1)^# n / # &]]; (* Michael Somos, May 17 2013 *)
a[ n_] := If[ n < 1, 0, Sum[ Mod[ d, 2] d, {d, Divisors[ n]}]]; (* Michael Somos, May 17 2013 *)
a[ n_] := If[ n < 1, 0, Times @@ (If[ # < 3, 1, (#^(#2 + 1) - 1) / (# - 1)] & @@@ FactorInteger @ n)]; (* Michael Somos, Aug 15 2015 *)
Array[Total[Divisors@ # /. d_ /; EvenQ@ d -> Nothing] &, {75}] (* Michael De Vlieger, Apr 07 2016 *)
Table[SeriesCoefficient[n Log[QPochhammer[-1, x]], {x, 0, n}], {n, 1, 75}] (* Vladimir Reshetnikov, Nov 21 2016 *)
Table[DivisorSum[n,#&,OddQ[#]&],{n,80}] (* Harvey P. Dale, Jun 19 2021 *)

{a(n) = if( n<1, 0, sumdiv( n, d, (-1)^(d+1) * n/d))}; /* Michael Somos, May 29 2005 */

N=66; x='x+O('x^N); Vec( serconvol( log(prod(j=1,N,1+x^j)), sum(j=1,N,j*x^j)))  /* Joerg Arndt, May 03 2008, edited by M. F. Hasler, Jun 19 2011 */

s=vector(100);for(n=1,100,s[n]=sumdiv(n,d,d*(d%2)));s /* Zak Seidov, Sep 24 2011*/

a(n)=sigma(n>>valuation(n,2)) \\ Charles R Greathouse IV, Sep 09 2014

from math import prod
from sympy import factorint
def A000593(n): return prod((p**(e+1)-1)//(p-1) for p, e in factorint(n).items() if p > 2) # Chai Wah Wu, Sep 09 2021

[sum(k for k in divisors(n) if k % 2) for n in (1..75)] # Giuseppe Coppoletta, Nov 02 2016

Inverse Moebius Transform of [0, 1, 0, 3, 0, 5, ...].
Dirichlet g.f.: zeta(s)*zeta(s-1)*(1-2^(1-s)).
a(2*n) = A000203(2*n)-2*A000203(n), a(2*n+1) = A000203(2*n+1). - Henry Bottomley, May 16 2000
a(2*n) = A054785(2*n) - A000203(2*n). - Reinhard Zumkeller, Apr 23 2008
Multiplicative with a(p^e) = 1 if p = 2, (p^(e+1)-1)/(p-1) if p > 2. - David W. Wilson, Aug 01 2001
a(n) = Sum_{d divides n} (-1)^(d+1)*n/d, Dirichlet convolution of A062157 with A000027. - Vladeta Jovovic, Sep 06 2002
Sum_{k=1..n} a(k) is asymptotic to c*n^2 where c=Pi^2/24. - Benoit Cloitre, Dec 29 2002
G.f.: Sum_{n>0} n*x^n/(1+x^n). - Vladeta Jovovic, Oct 11 2002
G.f.: (theta_3(q)^4 + theta_2(q)^4 -1)/24.
G.f.: Sum_{k>0} -(-x)^k / (1 - x^k)^2. - Michael Somos, Oct 29 2005
a(n) = A050449(n)+A050452(n); a(A000079(n))=1; a(A005408(n))=A000203(A005408(n)). - Reinhard Zumkeller, Apr 18 2006
From Joerg Arndt, Nov 09 2010: (Start)
G.f.: Sum_{n>=1} (2*n-1) * q^(2*n-1) / (1-q^(2*n-1)).
G.f.: deriv(log(P)) = deriv(P)/P where P = Product_{n>=1} (1 + q^n). (End)
Dirichlet convolution of A000203 with [1,-2,0,0,0,...]. - R. J. Mathar, Jun 28 2011
a(n) = Sum_{k = 1..A001227(n)} A182469(n,k). - Reinhard Zumkeller, May 01 2012
G.f.: -1/Q(0), where Q(k) = (x-1)*(1-x^(2*k+1)) + x*(-1 +x^(k+1))^4/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Apr 30 2013
a(n) = Sum_{k=1..n} k*A000009(k)*A081362(n-k). - Mircea Merca, Feb 26 2014
a(n) = A000203(n) - A146076(n). - Omar E. Pol, Apr 05 2016
a(2*n) = a(n). - Giuseppe Coppoletta, Nov 02 2016
a(n) = n * [x^n] log((-1; x)inf), where (a; q)_inf is the q-Pochhammer symbol. - _Vladimir Reshetnikov, Nov 21 2016
From Wolfdieter Lang, Dec 11 2016: (Start)
G.f.: Sum_{n>=1} x^n*(1+x^(2*n))/(1-x^(2*n))^2, from the second to last equation of the proof to Theorem 382 (with x^2 -> x) of the Hardy-Wright reference, p. 312.
a(n) = Sum_{d|n} (-d)*(-1)^(n/d), commutating factors of the D.g.f. given above by Jovovic, Oct 11 2002. See also the a(n) version given by Jovovic, Sep 06 2002. (End)
a(n) = A000203(n)/A038712(n). - Omar E. Pol, Dec 14 2017
a(n) = A000203(n)/(2^(1 + (A183063(n)/A001227(n))) - 1). - Omar E. Pol, Nov 06 2018
a(n) = A000203(2n) - 2*A000203(n). - Ridouane Oudra, Aug 28 2019
From Peter Bala, Jan 04 2021: (Start)
a(n) = (2/3)*A002131(n) + (1/3)*A002129(n) = (2/3)*A002131(n) + (-1)^(n+1)*(1/3)*A113184(n).
a(n) = A002131(n) - (1/2)*A146076; a(n) = 2*A002131(n) - A000203(n). (End)
a(n) = A000203(A000265(n)) - John Keith, Aug 30 2021
Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k)/A000203(k) = A065442 - 1 = 0.60669... . - Amiram Eldar, Dec 14 2024

A069288 Number of odd divisors of n <= sqrt(n).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 2, 1, 2, 1, 1, 3, 1, 1, 2, 1, 2, 2, 1, 1, 2, 2, 1, 2, 1, 1, 3, 1, 1, 2, 2, 2, 2, 1, 1, 2, 2, 2, 2, 1, 1, 3, 1, 1, 3, 1, 2, 2, 1, 1, 2, 3, 1, 2, 1, 1, 3, 1, 2, 2, 1, 2, 3, 1, 1, 3, 2, 1, 2, 1, 1, 4

Offset: 1

Reinhard Zumkeller, Mar 14 2002

nonn

a(n) = #{d : d = A182469(n,k), d <= A000196(n), k=1..A001227(n)}. - Reinhard Zumkeller, Apr 05 2015

			From _Gus Wiseman_, Feb 11 2021: (Start)
The inferior odd divisors for selected n are the columns below:
n: 1    9   30   90  225  315  630  945 1575 2835 4410 3465 8190 6930
  --------------------------------------------------------------------
   1    3    5    9   15   15   21   27   35   45   63   55   65   77
        1    3    5    9    9   15   21   25   35   49   45   63   63
             1    3    5    7    9   15   21   27   45   35   45   55
                  1    3    5    7    9   15   21   35   33   39   45
                       1    3    5    7    9   15   21   21   35   35
                            1    3    5    7    9   15   15   21   33
                                 1    3    5    7    9   11   15   21
                                      1    3    5    7    9   13   15
                                           1    3    5    7    9   11
                                                1    3    5    7    9
                                                     1    3    5    7
                                                          1    3    5
                                                               1    3
                                                                    1
(End)

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

Cf. A000005, A000196, A001227, A069289, A182469.
Positions of first appearances are A334853.
A055396 selects the least prime index.
A061395 selects the greatest prime index.
- Odd -
A000009 counts partitions into odd parts (A066208).
A026424 lists numbers with odd Omega.
A027193 counts odd-length partitions.
A067659 counts strict partitions of odd length (A030059).
- Inferior divisors -
A033676 selects the greatest inferior divisor.
A033677 selects the least superior divisor.
A038548 counts inferior divisors.
A060775 selects the greatest strictly inferior divisor.
A063538 lists numbers with a superior prime divisor.
A063539 lists numbers without a superior prime divisor.
A063962 counts inferior prime divisors.
A064052 lists numbers with a properly superior prime divisor.
A140271 selects the least properly superior divisor.
A217581 selects the greatest inferior divisor.
A333806 counts strictly inferior prime divisors.
Cf. A001055, A244991, A300272, A340101, A340607, A340832, A340854/A340855.

a069288 n = length $ takeWhile (<= a000196 n) $ a182469_row n
-- Reinhard Zumkeller, Apr 05 2015

odn[n_]:=Count[Divisors[n],?(OddQ[#]&&#<=Sqrt[n ]&)]; Array[odn,100] (* _Harvey P. Dale, Nov 04 2017 *)

a(n) = my(ir = sqrtint(n)); sumdiv(n, d, (d % 2) * (d <= ir)); \\ Michel Marcus, Jan 14 2014

G.f.: Sum_{n>=1} 1/(1-q^(2*n-1)) * q^((2*n-1)^2). [Joerg Arndt, Mar 04 2010]

A333807 Sum of odd divisors of n that are < sqrt(n).

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 9, 1, 1, 4, 1, 6, 4, 1, 1, 4, 6, 1, 4, 1, 1, 9, 1, 1, 4, 1, 6, 4, 1, 1, 4, 6, 8, 4, 1, 1, 9, 1, 1, 11, 1, 6, 4, 1, 1, 4, 13, 1, 4, 1, 1, 9, 1, 8, 4, 1, 6, 4, 1, 1, 11, 6, 1, 4, 1, 1, 18

Offset: 1

Ilya Gutkovskiy, Apr 05 2020

nonn

Cf. A000593, A069289, A070039, A333805, A333808, A333810.

Table[DivisorSum[n, # &, # < Sqrt[n] && OddQ[#] &], {n, 1, 90}]
nmax = 90; CoefficientList[Series[Sum[(2 k - 1) x^(2 k (2 k - 1))/(1 - x^(2 k - 1)), {k, 1, nmax}], {x, 0, nmax}], x] // Rest

G.f.: Sum_{k>=1} (2*k - 1) * x^(2*k*(2*k - 1)) / (1 - x^(2*k - 1)).

A333782 G.f.: Sum_{k>=1} (-1)^(k + 1) * k * x^(k^2) / (1 - x^k).

Original entry on oeis.org

1, 1, 1, -1, 1, -1, 1, -1, 4, -1, 1, 2, 1, -1, 4, -5, 1, 2, 1, -5, 4, -1, 1, -2, 6, -1, 4, -5, 1, 7, 1, -5, 4, -1, 6, -8, 1, -1, 4, 0, 1, -4, 1, -5, 9, -1, 1, -8, 8, 4, 4, -5, 1, -4, 6, 2, 4, -1, 1, -3, 1, -1, 11, -13, 6, -4, 1, -5, 4, 11, 1, -16, 1, -1, 9, -5, 8, -4, 1, -8

Offset: 1

Ilya Gutkovskiy, Apr 05 2020

sign

Excess of sum of odd divisors of n that are <= sqrt(n) over sum of even divisors of n that are <= sqrt(n).

Cf. A002129, A066839, A069289, A333781.

nmax = 80; CoefficientList[Series[Sum[(-1)^(k + 1) k x^(k^2)/(1 - x^k), {k, 1, nmax}], {x, 0, nmax}], x] // Rest

a(n) = sumdiv(n, d, if (d^2<=n, if (d%2, d, -d))); \\ Michel Marcus, Apr 05 2020

A333751 Sum of nonprime divisors of n that are <= sqrt(n).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 5, 1, 1, 1, 5, 1, 1, 1, 5, 1, 1, 1, 5, 1, 1, 1, 11, 1, 1, 1, 5, 1, 7, 1, 5, 1, 1, 1, 11, 1, 1, 1, 5, 1, 7, 1, 5, 1, 1, 1, 11, 1, 1, 1, 13, 1, 7, 1, 5, 1, 1, 1, 19, 1, 1, 1, 5, 1, 7, 1, 13, 10, 1, 1, 11, 1, 1, 1, 13, 1, 16

Offset: 1

Ilya Gutkovskiy, Apr 03 2020

nonn

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A018252, A023890, A066839, A069289, A069293, A097974, A333748, A333752, A333753.

f:= proc(n) convert(select(t -> not isprime(t) and t^2 <= n, numtheory:-divisors(n)),`+`) end proc:
map(f, [$1..100]); # Robert Israel, Sep 12 2024

Table[DivisorSum[n, # &, # <= Sqrt[n] && !PrimeQ[#] &], {n, 1, 90}]
nmax = 90; CoefficientList[Series[Sum[Boole[!PrimeQ[k]] k x^(k^2)/(1 - x^k), {k, 1, nmax}], {x, 0, nmax}], x] // Rest

a(n) = sumdiv(n, d, if ((d^2<=n) && !isprime(d), d)); \\ Michel Marcus, Apr 03 2020

G.f.: Sum_{k>=1} A018252(k) * x^(A018252(k)^2) / (1 - x^A018252(k)).

A333753 Sum of prime power divisors of n that are <= sqrt(n).

Original entry on oeis.org

0, 0, 0, 2, 0, 2, 0, 2, 3, 2, 0, 5, 0, 2, 3, 6, 0, 5, 0, 6, 3, 2, 0, 9, 5, 2, 3, 6, 0, 10, 0, 6, 3, 2, 5, 9, 0, 2, 3, 11, 0, 5, 0, 6, 8, 2, 0, 9, 7, 7, 3, 6, 0, 5, 5, 13, 3, 2, 0, 14, 0, 2, 10, 14, 5, 5, 0, 6, 3, 14, 0, 17, 0, 2, 8, 6, 7, 5, 0, 19, 12, 2, 0, 16, 5, 2, 3, 14, 0, 19

Offset: 1

Ilya Gutkovskiy, Apr 03 2020

nonn

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A023889, A066839, A069289, A069293, A097974, A246655, A333750, A333751, A333752.

f:= proc(n) local F,i,j,t;
  F:= ifactors(n)[2];
  t:= 0;
  for i from 1 to nops(F) do
    j:= min(F[i,2],ilog[F[i,1]^2](n));
    t:= t + (F[i,1]^j-1)*F[i,1]/(F[i,1]-1)
  od;
  t
end proc:
map(f, [$1..100]); # Robert Israel, Feb 15 2023

Table[DivisorSum[n, # &, # <= Sqrt[n] && PrimePowerQ[#] &], {n, 1, 90}]
nmax = 90; CoefficientList[Series[Sum[Boole[PrimePowerQ[k]] k x^(k^2)/(1 - x^k), {k, 1, nmax}], {x, 0, nmax}], x] // Rest

a(n) = sumdiv(n, d, if ((d^2<=n) && isprimepower(d), d)); \\ Michel Marcus, Apr 03 2020

G.f.: Sum_{p prime, k>=1} p^k * x^(p^(2*k)) / (1 - x^(p^k)).

A333752 Sum of squarefree divisors of n that are <= sqrt(n).

Original entry on oeis.org

1, 1, 1, 3, 1, 3, 1, 3, 4, 3, 1, 6, 1, 3, 4, 3, 1, 6, 1, 3, 4, 3, 1, 6, 6, 3, 4, 3, 1, 11, 1, 3, 4, 3, 6, 12, 1, 3, 4, 8, 1, 12, 1, 3, 9, 3, 1, 12, 8, 8, 4, 3, 1, 12, 6, 10, 4, 3, 1, 17, 1, 3, 11, 3, 6, 12, 1, 3, 4, 15, 1, 12, 1, 3, 9, 3, 8, 12, 1, 8, 4, 3, 1, 19, 6, 3, 4, 3, 1, 17

Offset: 1

Ilya Gutkovskiy, Apr 03 2020

nonn

Cf. A005117, A048250, A066839, A069289, A069293, A097974, A333749, A333751, A333753.

Table[DivisorSum[n, # &, # <= Sqrt[n] && SquareFreeQ[#] &], {n, 1, 90}]
nmax = 90; CoefficientList[Series[Sum[MoebiusMu[k]^2 k x^(k^2)/(1 - x^k), {k, 1, nmax}], {x, 0, nmax}], x] // Rest

a(n) = sumdiv(n, d, if ((d^2<=n) && issquarefree(d), d)); \\ Michel Marcus, Apr 03 2020

G.f.: Sum_{k>=1} mu(k)^2 * k * x^(k^2) / (1 - x^k).

A347173 Sum of squares of odd divisors of n that are <= sqrt(n).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 10, 1, 1, 10, 1, 1, 10, 1, 1, 10, 1, 1, 10, 1, 1, 10, 26, 1, 10, 1, 1, 35, 1, 1, 10, 1, 26, 10, 1, 1, 10, 26, 1, 10, 1, 1, 35, 1, 1, 10, 50, 26, 10, 1, 1, 10, 26, 50, 10, 1, 1, 35, 1, 1, 59, 1, 26, 10, 1, 1, 10, 75, 1, 10, 1, 1, 35, 1, 50, 10, 1, 26

Offset: 1

Ilya Gutkovskiy, Aug 21 2021

			a(18) = 10 as the odd divisors of 18 are the divisors of 9 which are 1, 3 and 9. Of those, 1 and 3 are <= sqrt(18) so we find the squares of 1 and 3 then add them i.e., a(18) = 1^2 + 3^2 = 10. - _David A. Corneth_, Feb 24 2024

David A. Corneth, Table of n, a(n) for n = 1..10000

Cf. A001157, A050999, A069288, A069289, A095118, A347161, A347174, A347175.

Table[DivisorSum[n, #^2 &, # <= Sqrt[n] && OddQ[#] &], {n, 1, 80}]
nmax = 80; CoefficientList[Series[Sum[(2 k - 1)^2 x^((2 k - 1)^2)/(1 - x^(2 k - 1)), {k, 1, nmax}], {x, 0, nmax}], x] // Rest

a(n) = sum(k=0, sqrtint(n), if ((k%2) && !(n%k), k^2)); \\ Michel Marcus, Aug 22 2021

a(n) = {
	my(s = sqrtint(n), res);
	n>>=valuation(n, 2);
	d = divisors(n);
	for(i = 1, #d,
		if(d[i] <= s,
			res += d[i]^2
		,
			return(res)
		)
	); res
} \\ David A. Corneth, Feb 24 2024

G.f.: Sum_{k>=1} (2*k - 1)^2 * x^((2*k - 1)^2) / (1 - x^(2*k - 1)).

A347174 Sum of cubes of odd divisors of n that are <= sqrt(n).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 28, 1, 1, 28, 1, 1, 28, 1, 1, 28, 1, 1, 28, 1, 1, 28, 126, 1, 28, 1, 1, 153, 1, 1, 28, 1, 126, 28, 1, 1, 28, 126, 1, 28, 1, 1, 153, 1, 1, 28, 344, 126, 28, 1, 1, 28, 126, 344, 28, 1, 1, 153, 1, 1, 371, 1, 126, 28, 1, 1, 28, 469, 1, 28, 1, 1, 153

Offset: 1

Ilya Gutkovskiy, Aug 21 2021

			a(18) = 28 as the odd divisors of 18 are the divisors of 9 which are 1, 3 and 9. Of those, 1 and 3 are <= sqrt(18) so we find the cubes of 1 and 3 then add them i.e., a(18) = 1^3 + 3^3 = 28. - _David A. Corneth_, Feb 24 2024

David A. Corneth, Table of n, a(n) for n = 1..10000

Cf. A001158, A051000, A069288, A069289, A280375, A347162, A347173, A347175.

Table[DivisorSum[n, #^3 &, # <= Sqrt[n] && OddQ[#] &], {n, 1, 75}]
nmax = 75; CoefficientList[Series[Sum[(2 k - 1)^3 x^((2 k - 1)^2)/(1 - x^(2 k - 1)), {k, 1, nmax}], {x, 0, nmax}], x] // Rest

a(n) = sum(k=0, sqrtint(n), if ((k%2) && !(n%k), k^3)); \\ Michel Marcus, Aug 22 2021

a(n) = {
	my(s = sqrtint(n), res);
	n>>=valuation(n, 2);
	d = divisors(n);
	for(i = 1, #d,
		if(d[i] <= s,
			res += d[i]^3
		,
			return(res)
		)
	); res
} \\ David A. Corneth, Feb 24 2024

G.f.: Sum_{k>=1} (2*k - 1)^3 * x^((2*k - 1)^2) / (1 - x^(2*k - 1)).

A347175 Sum of 4th powers of odd divisors of n that are <= sqrt(n).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 82, 1, 1, 82, 1, 1, 82, 1, 1, 82, 1, 1, 82, 1, 1, 82, 626, 1, 82, 1, 1, 707, 1, 1, 82, 1, 626, 82, 1, 1, 82, 626, 1, 82, 1, 1, 707, 1, 1, 82, 2402, 626, 82, 1, 1, 82, 626, 2402, 82, 1, 1, 707, 1, 1, 2483, 1, 626, 82, 1, 1, 82, 3027, 1, 82, 1, 1, 707

Offset: 1

Ilya Gutkovskiy, Aug 21 2021

			a(18) = 82 as the odd divisors of 18 are the divisors of 9 which are 1, 3 and 9. Of those, 1 and 3 are <= sqrt(18) so we find the sum of fourth powers of 1 and 3 then add them i.e., a(18) = 1^4 + 3^4 = 82. - _David A. Corneth_, Feb 24 2024

David A. Corneth, Table of n, a(n) for n = 1..10000

Cf. A001159, A051001, A069288, A069289, A347143, A347172, A347173, A347174.

Table[DivisorSum[n, #^4 &, # <= Sqrt[n] && OddQ[#] &], {n, 1, 75}]
nmax = 75; CoefficientList[Series[Sum[(2 k - 1)^4 x^((2 k - 1)^2)/(1 - x^(2 k - 1)), {k, 1, nmax}], {x, 0, nmax}], x] // Rest

a(n) = {
	my(s = sqrtint(n), res);
	n>>=valuation(n, 2);
	d = divisors(n);
	for(i = 1, #d,
		if(d[i] <= s,
			res += d[i]^4
		,
			return(res)
		)
	); res
} \\ David A. Corneth, Feb 24 2024

G.f.: Sum_{k>=1} (2*k - 1)^4 * x^((2*k - 1)^2) / (1 - x^(2*k - 1)).

cp's OEIS Frontend

A000593 Sum of odd divisors of n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

Python

Sage

Formula

A069288 Number of odd divisors of n <= sqrt(n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Formula

A333807 Sum of odd divisors of n that are < sqrt(n).

Views

Author

Keywords

Crossrefs

Programs

Mathematica

Formula

A333782 G.f.: Sum_{k>=1} (-1)^(k + 1) * k * x^(k^2) / (1 - x^k).

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

PARI

A333751 Sum of nonprime divisors of n that are <= sqrt(n).

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A333753 Sum of prime power divisors of n that are <= sqrt(n).

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A333752 Sum of squarefree divisors of n that are <= sqrt(n).

Views

Author

Keywords

Crossrefs