cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 11-20 of 57 results. Next

A321832 a(n) = Sum_{d|n, n/d==1 (mod 4)} d^8 - Sum_{d|n, n/d==3 (mod 4)} d^8.

Original entry on oeis.org

1, 256, 6560, 65536, 390626, 1679360, 5764800, 16777216, 43040161, 100000256, 214358880, 429916160, 815730722, 1475788800, 2562506560, 4294967296, 6975757442, 11018281216, 16983563040, 25600065536, 37817088000, 54875873280, 78310985280, 110058536960, 152588281251, 208827064832
Offset: 1

Views

Author

N. J. A. Sloane, Nov 24 2018

Keywords

Links

Crossrefs

Cf. A101455, A258816.
Cf. A321543 - A321565, A321807 - A321836 for similar sequences.
Glaisher's E'_i (i=0..12): A002654, A050469, A050470, A050471, A050468, A321829, A321830, A321831, this sequence, A321833, A321834, A321835, A321836.

Programs

  • Mathematica
    s[n_,r_] := DivisorSum[n, #^8 &, Mod[n/#,4]==r &]; a[n_] := s[n,1] - s[n,3]; Array[a, 30] (* Amiram Eldar, Nov 26 2018 *)
s[n_] := If[OddQ[n], (-1)^((n-1)/2), 0]; (* A101455 *)
f[p_, e_] := (p^(8*e+8) - s[p]^(e+1))/(p^8 - s[p]); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 04 2023 *)
  • PARI
    apply( A321832(n)=factorback(apply(f->f[1]^(8*f[2]+8)\/(f[1]^8+f[1]%4-2),Col(factor(n)))), [1..50]) \\ M. F. Hasler, Nov 26 2018

Formula

G.f.: Sum_{k>=1} k^8*x^k/(1 + x^(2*k)). - Ilya Gutkovskiy, Nov 26 2018
Multiplicative with a(p^e) = round(p^(8e+8)/(p^8 + (p mod 4) - 2)). (Following R. Israel in A321833.) - M. F. Hasler, Nov 26 2018
Sum_{k=1..n} a(k) ~ c * n^9 / 9, where c = 277*Pi^9/8257536 (A258816). - Amiram Eldar, Nov 04 2023
a(n) = Sum_{d|n} (n/d)^8*sin(d*Pi/2). - Ridouane Oudra, Sep 27 2024

A321833 a(n) = Sum_{d|n, n/d==1 mod 4} d^9 - Sum_{d|n, n/d==3 mod 4} d^9.

Original entry on oeis.org

1, 512, 19682, 262144, 1953126, 10077184, 40353606, 134217728, 387400807, 1000000512, 2357947690, 5159518208, 10604499374, 20661046272, 38441425932, 68719476736, 118587876498, 198349213184, 322687697778, 512000262144, 794239673292
Offset: 1

Views

Author

N. J. A. Sloane, Nov 24 2018

Keywords

Links

Crossrefs

Cf. A101455.
Cf. A321807 - A321836 for similar sequences.
Glaisher's E'_i (i=0..12): A002654, A050469, A050470, A050471, A050468, A321829, A321830, A321831, A321832, this sequence, A321834, A321835, A321836.

Programs

  • Maple
    f:= n ->
mul(piecewise(t[1]=2,2^(9*t[2]), t[1] mod 4 = 1, (t[1]^(9*(t[2]+1))-1)/(t[1]^9-1), (t[1]^(9*(t[2]+1))+(-1)^t[2])/(t[1]^9+1)), t = ifactors(n)[2]):
map(f, [$1..100]); # Robert Israel, Nov 26 2018
  • Mathematica
    s[n_,r_] := DivisorSum[n, #^9 &, Mod[n/#,4]==r &]; a[n_] := s[n,1] - s[n,3]; Array[a, 30] (* Amiram Eldar, Nov 26 2018 *)
s[n_] := If[OddQ[n], (-1)^((n-1)/2), 0]; (* A101455 *)
f[p_, e_] := (p^(9*e+9) - s[p]^(e+1))/(p^9 - s[p]); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 04 2023 *)
  • PARI
    apply( a(n)=sumdiv(n,d,if(bittest(n\d,0),(2-n\d%4)*d^9)), [1..30]) \\ M. F. Hasler, Nov 26 2018

Formula

G.f.: Sum_{k>=1} k^9*x^k/(1 + x^(2*k)). - Ilya Gutkovskiy, Nov 26 2018
From Robert Israel, Nov 26 2018: (Start) a(2^m) = 2^(9*m).
For prime p == 1 (mod 4), a(p^m) = (p^(9(m+1))-1)/(p^9-1).
For prime p == 3 (mod 4), a(p^m) = (p^(9(m+1))+(-1)^m)/(p^9+1). (End)
From Amiram Eldar, Nov 04 2023: (Start)
Multiplicative with a(p^e) = (p^(9*e+9) - A101455(p)^(e+1))/(p^9 - A101455(p)).
Sum_{k=1..n} a(k) ~ c * n^10 / 10, where c = beta(10) = 0.99998316402... and beta is the Dirichlet beta function. (End)
a(n) = Sum_{d|n} (n/d)^9*sin(d*Pi/2). - Ridouane Oudra, Sep 27 2024

A321834 a(n) = Sum_{d|n, n/d==1 mod 4} d^10 - Sum_{d|n, n/d==3 mod 4} d^10.

Original entry on oeis.org

1, 1024, 59048, 1048576, 9765626, 60465152, 282475248, 1073741824, 3486725353, 10000001024, 25937424600, 61916315648, 137858491850, 289254653952, 576640684048, 1099511627776, 2015993900450, 3570406761472, 6131066257800, 10240001048576
Offset: 1

Views

Author

N. J. A. Sloane, Nov 24 2018

Keywords

Links

Crossrefs

Cf. A101455.
Cf. A321807 - A321836 for similar sequences.
Glaisher's E'_i (i=0..12): A002654, A050469, A050470, A050471, A050468, A321829, A321830, A321831, A321832, A321833, this sequence, A321835, A321836.

Programs

  • Mathematica
    s[n_,r_] := DivisorSum[n, #^10 &, Mod[n/#,4]==r &]; a[n_] := s[n,1] - s[n,3]; Array[a, 30] (* Amiram Eldar, Nov 26 2018 *)
s[n_] := If[OddQ[n], (-1)^((n-1)/2), 0]; (* A101455 *)
f[p_, e_] := (p^(10*e+10) - s[p]^(e+1))/(p^10 - s[p]); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 04 2023 *)
  • PARI
    apply( a(n)=sumdiv(n, d, if(bittest(n\d,0),(2-n\d%4)*d^10)), [1..30]) \\ M. F. Hasler, Nov 26 2018

Formula

G.f.: Sum_{k>=1} k^10*x^k/(1 + x^(2*k)). - Ilya Gutkovskiy, Nov 26 2018
From Amiram Eldar, Nov 04 2023: (Start)
Multiplicative with a(p^e) = (p^(10*e+10) - A101455(p)^(e+1))/(p^10 - A101455(p)).
Sum_{k=1..n} a(k) ~ c * n^11 / 11, where c = beta(11) = 50521*Pi^11/14863564800 = 0.999994374973... and beta is the Dirichlet beta function. (End)
a(n) = Sum_{d|n} (n/d)^10*sin(d*Pi/2). - Ridouane Oudra, Sep 27 2024

A321835 a(n) = Sum_{d|n, n/d==1 mod 4} d^11 - Sum_{d|n, n/d==3 mod 4} d^11.

Original entry on oeis.org

1, 2048, 177146, 4194304, 48828126, 362795008, 1977326742, 8589934592, 31380882463, 100000002048, 285311670610, 743004176384, 1792160394038, 4049565167616, 8649707208396, 17592186044416, 34271896307634, 64268047284224, 116490258898218
Offset: 1

Views

Author

N. J. A. Sloane, Nov 24 2018

Keywords

Links

Crossrefs

Cf. A101455.
Cf. A321807 - A321836 for related sequences.
Glaisher's E'_i (i=0..12): A002654, A050469, A050470, A050471, A050468, A321829, A321830, A321831, A321832, A321833, A321834, this sequence, A321836.

Programs

  • Mathematica
    s[n_,r_] := DivisorSum[n, #^11 &, Mod[n/#,4]==r &]; a[n_] := s[n,1] - s[n,3]; Array[a, 30] (* Amiram Eldar, Nov 26 2018 *)
s[n_] := If[OddQ[n], (-1)^((n-1)/2), 0]; (* A101455 *)
f[p_, e_] := (p^(11*e+11) - s[p]^(e+1))/(p^11 - s[p]); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 04 2023 *)
  • PARI
    apply( a(n)=sumdiv(n,d,if(bittest(n\d,0),(2-n\d%4)*d^11)), [1..30]) \\ M. F. Hasler, Nov 26 2018

Formula

G.f.: Sum_{k>=1} k^11*x^k/(1 + x^(2*k)). - Ilya Gutkovskiy, Nov 26 2018
From Amiram Eldar, Nov 04 2023: (Start)
Multiplicative with a(p^e) = (p^(11*e+11) - A101455(p)^(e+1))/(p^11 - A101455(p)).
Sum_{k=1..n} a(k) ~ c * n^12 / 12, where c = beta(12) = 0.99999812235..., and beta is the Dirichlet beta function. (End)
a(n) = Sum_{d|n} (n/d)^11*sin(d*Pi/2). - Ridouane Oudra, Sep 27 2024

A321810 Sum of 6th powers of odd divisors of n.

Original entry on oeis.org

1, 1, 730, 1, 15626, 730, 117650, 1, 532171, 15626, 1771562, 730, 4826810, 117650, 11406980, 1, 24137570, 532171, 47045882, 15626, 85884500, 1771562, 148035890, 730, 244156251, 4826810, 387952660, 117650, 594823322, 11406980, 887503682
Offset: 1

Views

Author

N. J. A. Sloane, Nov 24 2018

Keywords

Links

Crossrefs

Column k=6 of A285425.
Cf. A050999, A051000, A051001, A051002, A321811 - A321816 (analog for 2nd .. 12th powers).
Cf. A321543 - A321565, A321807 - A321836 for related sequences.
Cf. A000265, A013665, A013954.

Programs

  • Mathematica
    f[2, e_] := 1; f[p_, e_] := (p^(6*e + 6) - 1)/(p^6 - 1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 02 2022 *)
  • PARI
    apply( A321810(n)=sigma(n>>valuation(n,2),6), [1..30]) \\ M. F. Hasler, Nov 26 2018
  • Python
    from sympy import divisor_sigma
def A321810(n): return int(divisor_sigma(n>>(~n&n-1).bit_length(),6)) # Chai Wah Wu, Jul 16 2022

Formula

a(n) = A013954(A000265(n)) = sigma_6(odd part of n); in particular, a(2^k) = 1 for all k >= 0. - M. F. Hasler, Nov 26 2018
G.f.: Sum_{k>=1} (2*k - 1)^6*x^(2*k-1)/(1 - x^(2*k-1)). - Ilya Gutkovskiy, Dec 22 2018
From Amiram Eldar, Nov 02 2022: (Start)
Multiplicative with a(2^e) = 1 and a(p^e) = (p^(6*e+6)-1)/(p^6-1) for p > 2.
Sum_{k=1..n} a(k) ~ c * n^7, where c = zeta(7)/14 = 0.0720249... . (End)
a(n) + a(n/2)*2^6 = A013954(n) where a(.)=0 for non-integer arguments. - R. J. Mathar, Aug 15 2023

A321552 a(n) = Sum_{d|n} (-1)^(n/d+1)*d^7.

Original entry on oeis.org

1, 127, 2188, 16255, 78126, 277876, 823544, 2080639, 4785157, 9922002, 19487172, 35565940, 62748518, 104590088, 170939688, 266321791, 410338674, 607714939, 893871740, 1269938130, 1801914272, 2474870844, 3404825448, 4552438132, 6103593751, 7969061786, 10465138360, 13386707720, 17249876310
Offset: 1

Views

Author

N. J. A. Sloane, Nov 23 2018

Keywords

Links

Crossrefs

Sum_{k>=1} k^b*x^k/(1 + x^k): A000593 (b=1), A078306 (b=2), A078307 (b=3), A284900 (b=4), A284926 (b=5), A284927 (b=6), this sequence (b=7), A321553 (b=8), A321554 (b=9), A321555 (b=10), A321556 (b=11), A321557 (b=12).
Cf. A321543 - A321565, A321807 - A321836 for similar sequences.
Cf. A013666.

Programs

  • Mathematica
    f[p_, e_] := (p^(7*e + 7) - 1)/(p^7 - 1); f[2, e_] := (63*2^(7*e + 1) + 1)/127; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 50] (* Amiram Eldar, Nov 11 2022 *)
  • PARI
    apply( A321552(n)=sumdiv(n, d, (-1)^(n\d-1)*d^7), [1..30]) \\ M. F. Hasler, Nov 26 2018

Formula

G.f.: Sum_{k>=1} k^7*x^k/(1 + x^k). - Seiichi Manyama, Nov 23 2018
From Amiram Eldar, Nov 11 2022: (Start)
Multiplicative with a(2^e) = (63*2^(7*e+1)+1)/127, and a(p^e) = (p^(7*e+7) - 1)/(p^7 - 1) if p > 2.
Sum_{k=1..n} a(k) ~ c * n^8, where c = 127*zeta(8)/1024 = 0.124529... . (End)

A321557 a(n) = Sum_{d|n} (-1)^(n/d+1)*d^12.

Original entry on oeis.org

1, 4095, 531442, 16773119, 244140626, 2176254990, 13841287202, 68702695423, 282430067923, 999755863470, 3138428376722, 8913939907598, 23298085122482, 56680071092190, 129746582562692, 281406240452607, 582622237229762, 1156551128144685, 2213314919066162, 4094999772632494
Offset: 1

Views

Author

N. J. A. Sloane, Nov 23 2018

Keywords

Links

Crossrefs

Cf. A321543 - A321565, A321807 - A321836 for similar sequences.
Cf. A013671.

Programs

  • Mathematica
    f[p_, e_] := (p^(12*e + 12) - 1)/(p^12 - 1); f[2, e_] := (2047*2^(12*e + 1) + 1)/4095; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 50] (* Amiram Eldar, Nov 11 2022 *)
  • PARI
    apply( A321557(n)=sumdiv(n, d, (-1)^(n\d-1)*d^12), [1..30]) \\ M. F. Hasler, Nov 26 2018

Formula

G.f.: Sum_{k>=1} k^12*x^k/(1 + x^k). - Seiichi Manyama, Nov 25 2018
From Amiram Eldar, Nov 11 2022: (Start)
Multiplicative with a(2^e) = (2047*2^(12*e+12)+1)/4095, and a(p^e) = (p^(12*e+12) - 1)/(p^12 - 1) if p > 2.
Sum_{k=1..n} a(k) ~ c * n^13, where c = 315*zeta(13)/4096 = 0.0769137... . (End)

A321816 Sum of 12th powers of odd divisors of n.

Original entry on oeis.org

1, 1, 531442, 1, 244140626, 531442, 13841287202, 1, 282430067923, 244140626, 3138428376722, 531442, 23298085122482, 13841287202, 129746582562692, 1, 582622237229762, 282430067923, 2213314919066162, 244140626, 7355841353205284, 3138428376722
Offset: 1

Views

Author

N. J. A. Sloane, Nov 24 2018

Keywords

Links

Crossrefs

Column k=12 of A285425.
Cf. A050999, A051000, A051001, A051002, A321810 - A321815 (analog for 2nd .. 11th powers).
Cf. A321543 - A321565, A321807 - A321836 for related sequences.
Cf. A000265, A013671, A013960.

Programs

  • Mathematica
    a[n_] := DivisorSum[n, #^12&, OddQ[#]&]; Array[a, 20] (* Amiram Eldar, Dec 07 2018 *)
  • PARI
    apply( A321816(n)=sigma(n>>valuation(n,2),12), [1..30]) \\ M. F. Hasler, Nov 26 2018
  • Python
    from sympy import divisor_sigma
def A321816(n): return int(divisor_sigma(n>>(~n&n-1).bit_length(),12)) # Chai Wah Wu, Jul 16 2022

Formula

a(n) = A013960(A000265(n)) = sigma_12(odd part of n); in particular, a(2^k) = 1 for all k >= 0. - M. F. Hasler, Nov 26 2018
G.f.: Sum_{k>=1} (2*k - 1)^12*x^(2*k-1)/(1 - x^(2*k-1)). - Ilya Gutkovskiy, Dec 22 2018
From Amiram Eldar, Nov 02 2022: (Start)
Multiplicative with a(2^e) = 1 and a(p^e) = (p^(12*e+12)-1)/(p^12-1) for p > 2.
Sum_{k=1..n} a(k) ~ c * n^13, where c = zeta(13)/26 = 0.0384662... . (End)

A321553 a(n) = Sum_{d|n} (-1)^(n/d+1)*d^8.

Original entry on oeis.org

1, 255, 6562, 65279, 390626, 1673310, 5764802, 16711423, 43053283, 99609630, 214358882, 428360798, 815730722, 1470024510, 2563287812, 4278124287, 6975757442, 10978587165, 16983563042, 25499674654, 37828630724, 54661514910, 78310985282, 109660357726, 152588281251, 208011334110, 282472589764
Offset: 1

Views

Author

N. J. A. Sloane, Nov 23 2018

Keywords

Links

Crossrefs

Cf. A321543 - A321565, A321807 - A321836 for similar sequences.
Cf. A013667.

Programs

  • Mathematica
    Table[Total[(-1)^(n/#+1) #^8&/@Divisors[n]],{n,30}] (* Harvey P. Dale, May 05 2021 *)
f[p_, e_] := (p^(8*e + 8) - 1)/(p^8 - 1); f[2, e_] := (127*2^(8*e + 1) + 1)/255; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 50] (* Amiram Eldar, Nov 11 2022 *)
  • PARI
    apply( A321553(n)=sumdiv(n, d, (-1)^(n\d-1)*d^8), [1..30]) \\ M. F. Hasler, Nov 26 2018

Formula

G.f.: Sum_{k>=1} k^8*x^k/(1 + x^k). - Seiichi Manyama, Nov 23 2018
From Amiram Eldar, Nov 11 2022: (Start)
Multiplicative with a(2^e) = (127*2^(8*e+1)+1)/255, and a(p^e) = (p^(8*e+8) - 1)/(p^8 - 1) if p > 2.
Sum_{k=1..n} a(k) ~ c * n^9, where c = 85*zeta(9)/768 = 0.110899... . (End)

A321554 a(n) = Sum_{d|n} (-1)^(n/d+1)*d^9.

Original entry on oeis.org

1, 511, 19684, 261631, 1953126, 10058524, 40353608, 133955071, 387440173, 998047386, 2357947692, 5149944604, 10604499374, 20620693688, 38445332184, 68584996351, 118587876498, 197981928403, 322687697780, 510998308506, 794320419872, 1204911270612, 1801152661464, 2636771617564, 3814699218751
Offset: 1

Views

Author

N. J. A. Sloane, Nov 23 2018

Keywords

Links

Crossrefs

Cf. A321543 - A321565, A321807 - A321836 for similar sequences.
Cf. A013668.

Programs

  • Mathematica
    f[p_, e_] := (p^(9*e + 9) - 1)/(p^9 - 1); f[2, e_] := (255*2^(9*e + 1) + 1)/511; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 50] (* Amiram Eldar, Nov 11 2022 *)
  • PARI
    apply( A321554(n)=sumdiv(n, d, (-1)^(n\d-1)*d^9), [1..30]) \\ M. F. Hasler, Nov 26 2018

Formula

G.f.: Sum_{k>=1} k^9*x^k/(1 + x^k). - Seiichi Manyama, Nov 24 2018
From Amiram Eldar, Nov 11 2022: (Start)
Multiplicative with a(2^e) = (255*2^(9*e+1)+1)/511, and a(p^e) = (p^(9*e+9) - 1)/(p^9 - 1) if p > 2.
Sum_{k=1..n} a(k) ~ c * n^10, where c = 511*zeta(10)/5120 = 0.0999039... . (End)
Previous Showing 11-20 of 57 results. Next

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.