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-19 of 19 results.

A363891 Sum of divisors of 3*n-2 of form 3*k+2.

Original entry on oeis.org

0, 2, 0, 7, 0, 10, 0, 13, 5, 16, 0, 19, 0, 35, 0, 25, 0, 28, 16, 31, 0, 42, 0, 56, 0, 40, 0, 43, 22, 65, 0, 49, 0, 77, 0, 55, 0, 80, 28, 61, 11, 64, 0, 98, 0, 95, 0, 73, 34, 76, 0, 104, 0, 147, 0, 85, 0, 88, 40, 91, 0, 125, 28, 140, 0, 114, 0, 103, 46, 140, 0, 109, 0, 192, 0, 115, 0
Offset: 1

Views

Author

Seiichi Manyama, Jun 26 2023

Keywords

Crossrefs

Programs

  • Mathematica
    a[n_] := DivisorSum[3*n - 2, # &, Mod[#, 3] == 2 &]; Array[a, 100] (* Amiram Eldar, Jun 26 2023 *)
  • PARI
    a(n) = sumdiv(3*n-2, d, (d%3==2)*d);

Formula

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

A082050 Sum of divisors of n that are not of the form 3k+1.

Original entry on oeis.org

0, 2, 3, 2, 5, 11, 0, 10, 12, 7, 11, 23, 0, 16, 23, 10, 17, 38, 0, 27, 24, 13, 23, 55, 5, 28, 39, 16, 29, 61, 0, 42, 47, 19, 40, 86, 0, 40, 42, 35, 41, 88, 0, 57, 77, 25, 47, 103, 0, 57, 71, 28, 53, 119, 16, 80, 60, 31, 59, 153, 0, 64, 96, 42, 70, 121, 0, 87, 95, 56, 71, 190
Offset: 1

Views

Author

Ralf Stephan, Apr 02 2003

Keywords

Crossrefs

Programs

  • Mathematica
    sd[n_]:=Total[Select[Divisors[n],!IntegerQ[(#-1)/3]&]]; Array[sd,80] (* Harvey P. Dale, May 04 2011 *)
  • PARI
    for(n=1,100,print1(sumdiv(n,d,if(d%3!=1,d))","))
    
  • PARI
    N = 66;  x = 'x + O('x^N);
    gf = sum(n=1,N, (3*n-1)*x^(3*n-1)/(1-x^(3*n-1)) + (3*n)*x^(3*n)/(1-x^(3*n)) );
    v = Vec(gf);  concat([0],v)
    \\ Joerg Arndt, May 17 2013

Formula

a(A004611(n)) = 0.
G.f.: Sum_{k>=1} x^(2*k)*(2+3*x^k+x^(3*k))/(1-x^(3*k))^2. - Vladeta Jovovic, Apr 11 2006
Sum_{k=1..n} a(k) ~ c * n^2, where c = Pi^2/18 = 0.548311... (A086463). - Amiram Eldar, Jan 06 2024

A082051 Sum of divisors of n that are not of the form 3k+2.

Original entry on oeis.org

1, 1, 4, 5, 1, 10, 8, 5, 13, 11, 1, 26, 14, 8, 19, 21, 1, 37, 20, 15, 32, 23, 1, 50, 26, 14, 40, 40, 1, 65, 32, 21, 37, 35, 8, 89, 38, 20, 56, 55, 1, 80, 44, 27, 73, 47, 1, 114, 57, 36, 55, 70, 1, 118, 56, 40, 80, 59, 1, 141, 62, 32, 104, 85, 14, 131, 68, 39, 73, 88, 1, 185
Offset: 1

Views

Author

Ralf Stephan, Apr 02 2003

Keywords

Crossrefs

Programs

  • Mathematica
    sd[n_]:= Total[Select[Divisors[n], !IntegerQ[(# - 2) / 3]&]]; Array[sd, 100] (* Vincenzo Librandi, May 17 2013 *)
  • PARI
    for(n=1,100,print1(sumdiv(n,d,if(d%3!=2,d))","))
    
  • PARI
    N = 66;  x = 'x + O('x^N);
    gf = sum(n=1,N, (3*n-2)*x^(3*n-2)/(1-x^(3*n-2)) + (3*n)*x^(3*n)/(1-x^(3*n)) );
    v = Vec(gf)
    \\ Joerg Arndt, May 17 2013

Formula

a(A003627(n)) = 1.
G.f.: Sum_{k>=1} x^k*(1 + 3*x^(2*k) + 2*x^(3*k))/(1 - x^(3*k))^2. - Ilya Gutkovskiy, Sep 12 2019
Sum_{k=1..n} a(k) ~ c * n^2, where c = Pi^2/18 = 0.548311... (A086463). - Amiram Eldar, Jan 06 2024

A082052 Sum of divisors of n that are not of the form 4k+1.

Original entry on oeis.org

0, 2, 3, 6, 0, 11, 7, 14, 3, 12, 11, 27, 0, 23, 18, 30, 0, 29, 19, 36, 10, 35, 23, 59, 0, 28, 30, 55, 0, 66, 31, 62, 14, 36, 42, 81, 0, 59, 42, 84, 0, 74, 43, 83, 18, 71, 47, 123, 7, 62, 54, 84, 0, 110, 66, 119, 22, 60, 59, 162, 0, 95, 73, 126, 0, 110, 67, 108, 26, 138, 71
Offset: 1

Views

Author

Ralf Stephan, Apr 02 2003

Keywords

Comments

a(A004613(n))=0.

Crossrefs

Programs

  • Mathematica
    sd[n_]:= Total[Select[Divisors[n], !IntegerQ[(# - 1) / 4]&]]; Array[sd, 100] (* Vincenzo Librandi, May 17 2013 *)
    Table[DivisorSum[n,#&,(!IntegerQ[(#-1)/4]&)],{n,80}] (* Harvey P. Dale, Nov 30 2019 *)
  • PARI
    for(n=1,100,print1(sumdiv(n,d,if(d%4!=1,d))","))

Formula

G.f.: Sum_{k>=1} x^(2*k)*(2 + 3*x^k + 4*x^(2*k) + 2*x^(4*k) + x^(5*k))/(1 - x^(4*k))^2. - Ilya Gutkovskiy, Sep 12 2019

A363902 Expansion of Sum_{k>0} x^(2*k) / (1 - x^(3*k))^2.

Original entry on oeis.org

0, 1, 0, 1, 2, 1, 0, 4, 0, 3, 4, 1, 0, 6, 2, 4, 6, 1, 0, 10, 0, 5, 8, 4, 2, 10, 0, 6, 10, 3, 0, 15, 4, 7, 14, 1, 0, 14, 0, 13, 14, 6, 0, 20, 2, 9, 16, 4, 0, 20, 6, 10, 18, 1, 6, 28, 0, 11, 20, 10, 0, 22, 0, 15, 24, 5, 0, 30, 8, 20, 24, 4, 0, 26, 2, 14, 30, 10, 0, 40, 0, 15, 28, 6, 8
Offset: 1

Views

Author

Seiichi Manyama, Jun 27 2023

Keywords

Crossrefs

Programs

  • Mathematica
    a[n_] := DivisorSum[n, # + 1 &, Mod[#, 3] == 2 &]/3; Array[a, 100] (* Amiram Eldar, Jun 27 2023 *)
  • PARI
    a(n) = sumdiv(n, d, (d%3==2)*(d+1))/3;

Formula

a(n) = (1/3) * Sum_{d|n, d==2 mod 3} (d+1) = (A001822(n) + A078182(n))/3.
G.f.: Sum_{k>0} k * x^(3*k-1) / (1 - x^(3*k-1)).

A082053 Sum of divisors of n that are not of the form 4k+3.

Original entry on oeis.org

1, 3, 1, 7, 6, 9, 1, 15, 10, 18, 1, 25, 14, 17, 6, 31, 18, 36, 1, 42, 22, 25, 1, 57, 31, 42, 10, 49, 30, 54, 1, 63, 34, 54, 6, 88, 38, 41, 14, 90, 42, 86, 1, 73, 60, 49, 1, 121, 50, 93, 18, 98, 54, 90, 6, 113, 58, 90, 1, 150, 62, 65, 31, 127, 84, 130, 1, 126, 70, 102, 1, 192
Offset: 1

Views

Author

Ralf Stephan, Apr 02 2003

Keywords

Comments

a(A002145(n))=1.

Crossrefs

Programs

  • Mathematica
    sd[n_]:= Total[Select[Divisors[n], !IntegerQ[(# - 3) / 4]&]]; Array[sd, 100] (* Vincenzo Librandi, May 17 2013 *)
  • PARI
    for(n=1,100,print1(sumdiv(n,d,if(d%4!=3,d))","))

Formula

G.f.: Sum_{k>=1} x^k*(1 + 2*x^k + 4*x^(3*k) + 3*x^(4*k) + 2*x^(5*k))/(1 - x^(4*k))^2. - Ilya Gutkovskiy, Sep 12 2019

A272715 Numbers equal to the sum of their proper divisors d such that d mod 3 = 2.

Original entry on oeis.org

280, 34960, 4960000, 15887872, 20527600, 105056320, 44781248512
Offset: 1

Views

Author

Giovanni Resta, May 05 2016

Keywords

Comments

The numbers equal to the sum of their proper divisors which are a multiple of 3 are the perfect numbers (A000396) multiplied by 3.

Examples

			The proper divisors of 280 which are congruent to 2 mod 3 are 2, 5, 8, 14, 20, 35, 56, and 140. Since their sum is 280, 280 is a term.
		

Crossrefs

Cf. A272716, A000396, A078182 (sum of proper and improver divisors).

Programs

  • Mathematica
    Select[Range[40000], # == Plus @@ Select[ Most@ Divisors[#], Mod[#, 3] == 2 &] &]
  • PARI
    is(n)=sumdiv(n, d, if(d%3==2, d, 0))==if(n%3==2, 2*n, n) \\ Charles R Greathouse IV, May 09 2016

A326503 Expansion of Sum_{k>=1} x^k * (1 - x^(2*k)) / (1 + x^k + x^(2*k))^2.

Original entry on oeis.org

1, -1, 1, 3, -4, -1, 8, -5, 1, 4, -10, 3, 14, -8, -4, 11, -16, -1, 20, -12, 8, 10, -22, -5, 21, -14, 1, 24, -28, 4, 32, -21, -10, 16, -32, 3, 38, -20, 14, 20, -40, -8, 44, -30, -4, 22, -46, 11, 57, -21, -16, 42, -52, -1, 40, -40, 20, 28, -58, -12, 62, -32, 8, 43, -56, 10
Offset: 1

Views

Author

Ilya Gutkovskiy, Sep 12 2019

Keywords

Crossrefs

Cf. A002129, A050457, A078181, A078182, A078708, A162397 (Moebius transform), A326401.

Programs

  • Mathematica
    nmax = 66; CoefficientList[Series[Sum[x^k (1 - x^(2 k))/(1 + x^k + x^(2 k))^2, {k, 1, nmax}], {x, 0, nmax}], x] // Rest
    Table[DivisorSum[n, # &, MemberQ[{1}, Mod[#, 3]] &] - DivisorSum[n, # &, MemberQ[{2}, Mod[#, 3]] &], {n, 1, 66}]
    f[p_, e_] := If[Mod[p, 3] == 1, (p^(e + 1) - 1)/(p - 1), ((-p)^(e + 1) - 1)/(-p - 1)]; f[3, e_] := 1; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 28 2023 *)
  • PARI
    a(n)={sumdiv(n, d, d*((d+1)%3-1))} \\ Andrew Howroyd, Sep 12 2019

Formula

a(n) = Sum_{d|n, d==1 (mod 3)} d - Sum_{d|n, d==2 (mod 3)} d.
a(n) = A078181(n) - A078182(n).
Multiplicative with a(3^e) = 1, a(p^e) = (p^(e+1)-1)/(p-1) if p == 1 (mod 3) and a(p^e) = ((-p)^(e+1)-1)/(-p-1) if p == 2 (mod 3). - Amiram Eldar, Nov 28 2023

A374018 Expansion of Product_{k>=1} 1 / (1 - x^(3*k-1))^2.

Original entry on oeis.org

1, 0, 2, 0, 3, 2, 4, 4, 7, 6, 13, 10, 19, 18, 27, 30, 42, 44, 63, 66, 91, 100, 130, 144, 187, 206, 263, 294, 364, 412, 506, 568, 696, 782, 943, 1070, 1273, 1444, 1713, 1936, 2285, 2586, 3027, 3428, 3996, 4516, 5243, 5924, 6841, 7730, 8895, 10030, 11512, 12966, 14825, 16696
Offset: 0

Views

Author

Ilya Gutkovskiy, Jun 25 2024

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 55; CoefficientList[Series[Product[1/(1 - x^(3 k - 1))^2, {k, 1, nmax}], {x, 0, nmax}], x]

Formula

a(0) = 1; a(n) = (2/n) * Sum_{k=1..n} A078182(k) * a(n-k).
a(n) = Sum_{k=0..n} A035386(k) * A035386(n-k).
a(n) ~ exp(2*Pi*sqrt(n)/3) * Pi^(4/3) / (3^(3/2) * Gamma(1/3)^2 * n^(11/12)). - Vaclav Kotesovec, Jun 25 2024
Previous Showing 11-19 of 19 results.