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.

Showing 1-4 of 4 results.

A062011 a(n) = 2*tau(n) = 2*A000005(n).

Original entry on oeis.org

2, 4, 4, 6, 4, 8, 4, 8, 6, 8, 4, 12, 4, 8, 8, 10, 4, 12, 4, 12, 8, 8, 4, 16, 6, 8, 8, 12, 4, 16, 4, 12, 8, 8, 8, 18, 4, 8, 8, 16, 4, 16, 4, 12, 12, 8, 4, 20, 6, 12, 8, 12, 4, 16, 8, 16, 8, 8, 4, 24, 4, 8, 12, 14, 8, 16, 4, 12, 8, 16, 4, 24, 4, 8, 12, 12, 8, 16, 4, 20, 10, 8, 4, 24, 8, 8, 8, 16
Offset: 1

Views

Author

Ahmed Fares (ahmedfares(AT)my-deja.com), Jul 12 2001

Keywords

Comments

Old definition was "Number of cyclic subgroups of the group C_n X C_2 (where C_n is the cyclic group with n elements)."
More generally, the number of cyclic subgroups of the group C_n X C_m is Sum_{i|n, j|m} phi(i)*phi(j)/phi(lcm(i,j)), where phi=Euler totient function, cf. A000010. - Vladeta Jovovic, Jul 15 2001
Number of divisors of p*n, where p is any prime not dividing n. - Reinhard Zumkeller, May 17 2006
From Enrique Pérez Herrero, Jul 21 2011: (Start)
If p(x) is a polynomial with integer coefficients, and if r is an integer zero of p(x), then r is a divisor of the constant term c_0 of p(x). Under this theorem, p(x) can have a(c_0) possible integer roots.
a(n) is the number of integer divisors of n, while A000005(n) is the number of positive divisors. (End)
Number of solutions to the Diophantine equation i*j = n*i + j. - Robert G. Wilson v, Apr 10 2019
a(n) is also the number of times n appears in the triangle A333119, or equivalently, the number of positive integer solutions of the equation A333119(x, y) = n for y < x. - Stefano Spezia, Oct 05 2022

Crossrefs

Programs

Formula

a(n) = A000005(A087560(n)) = A000005(A119416(n)). - Reinhard Zumkeller, May 17 2006
L.g.f.: -log(Product_{k>=1} (1 - x^k)^(2/k)) = Sum_{n>=1} a(n)*x^n/n. - Ilya Gutkovskiy, Mar 18 2018

Extensions

More terms from Vladeta Jovovic, Jul 14 2001
Edited by N. J. A. Sloane, Sep 20 2018, replacing old definition (which was of course correct) with a simple formula.

A062249 a(n) = n + d(n), where d(n) = number of divisors of n, cf. A000005.

Original entry on oeis.org

2, 4, 5, 7, 7, 10, 9, 12, 12, 14, 13, 18, 15, 18, 19, 21, 19, 24, 21, 26, 25, 26, 25, 32, 28, 30, 31, 34, 31, 38, 33, 38, 37, 38, 39, 45, 39, 42, 43, 48, 43, 50, 45, 50, 51, 50, 49, 58, 52, 56, 55, 58, 55, 62, 59, 64, 61, 62, 61, 72, 63, 66, 69, 71, 69, 74, 69, 74, 73, 78, 73
Offset: 1

Views

Author

Ahmed Fares (ahmedfares(AT)my-deja.com), Jul 01 2001

Keywords

Comments

Number of cyclic subgroups of dihedral group with 2n elements.
a(n) is the n-th smallest number not a divisor of n. - J. Lowell, Apr 06 2008

Crossrefs

Cf. A064491 (iteration, start=1).

Programs

  • Haskell
    a062249 n = a000005 n + n  -- Reinhard Zumkeller, Mar 29 2014
    
  • Maple
    with(numtheory):seq(n+tau(n), n=1..71) ; # Zerinvary Lajos, Jun 04 2008
  • Mathematica
    Table[n + DivisorSigma[0, n], {n, 100}] (* Indranil Ghosh, Apr 12 2017 *)
  • PARI
    a(n) = n + numdiv(n) \\ Harry J. Smith, Aug 03 2009
    
  • Python
    from sympy.ntheory import divisor_count
    [n + divisor_count(n) for n in range(101)] # Indranil Ghosh, Apr 12 2017

Formula

a(n) = n + A000005(n). - Omar E. Pol, Dec 12 2008
From Ilya Gutkovskiy, Apr 12 2017: (Start)
G.f.: x/(1 - x)^2 + Sum_{k>=1} x^k/(1 - x^k).
Dirichlet g.f.: zeta(s)^2 + zeta(s-1). (End)

Extensions

Formula and more terms from Vladeta Jovovic, Jul 03 2001

A216624 Square array read by antidiagonals, T(n,k) = sum_{c|n,d|k} gcd(c,d) for n>=1, k>=1.

Original entry on oeis.org

1, 2, 2, 2, 5, 2, 3, 4, 4, 3, 2, 8, 6, 8, 2, 4, 4, 6, 6, 4, 4, 2, 10, 4, 15, 4, 10, 2, 4, 4, 12, 6, 6, 12, 4, 4, 3, 11, 4, 16, 8, 16, 4, 11, 3, 4, 6, 8, 6, 8, 8, 6, 8, 6, 4, 2, 10, 10, 22, 4, 30, 4, 22, 10, 10, 2, 6, 4, 8, 9, 8, 8, 8, 8, 9, 8, 4, 6
Offset: 1

Views

Author

Peter Luschny, Sep 12 2012

Keywords

Comments

T(n,k) = number of subgroups of C_n X C_k. [Hampjes et al.] - N. J. A. Sloane, Feb 02 2013

Examples

			[----1---2---3---4---5---6---7---8---9--10--11--12]
[ 1] 1,  2,  2,  3,  2,  4,  2,  4,  3,  4,  2,  6
[ 2] 2,  5,  4,  8,  4, 10,  4, 11,  6, 10,  4, 16
[ 3] 2,  4,  6,  6,  4, 12,  4,  8, 10,  8,  4, 18
[ 4] 3,  8,  6, 15,  6, 16,  6, 22,  9, 16,  6, 30
[ 5] 2,  4,  4,  6,  8,  8,  4,  8,  6, 16,  4, 12
[ 6] 4, 10, 12, 16,  8, 30,  8, 22, 20, 20,  8, 48
[ 7] 2,  4,  4,  6,  4,  8, 10,  8,  6,  8,  4, 12
[ 8] 4, 11,  8, 22,  8, 22,  8, 37, 12, 22,  8, 44
[ 9] 3,  6, 10,  9,  6, 20,  6, 12, 23, 12,  6, 30
[10] 4, 10,  8, 16, 16, 20,  8, 22, 12, 40,  8, 32
[11] 2,  4,  4,  6,  4,  8,  4,  8,  6,  8, 14, 12
[12] 6, 16, 18, 30, 12, 48, 12, 44, 30, 32, 12, 90
.
Displayed as a triangular array:
1,
2,  2,
2,  5,  2,
3,  4,  4,  3,
2,  8,  6,  8, 2,
4,  4,  6,  6, 4,  4,
2, 10,  4, 15, 4, 10, 2,
4,  4, 12,  6, 6, 12, 4,  4,
3, 11,  4, 16, 8, 16, 4, 11, 3,
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    T:= (n, k)-> add(add(igcd(c,d), c=divisors(n)), d=divisors(k)):
    seq(seq(T(n, 1+d-n), n=1..d), d=1..14); # Alois P. Heinz, Sep 12 2012
    T:=proc(m,n) local d; add( d*tau(m*n/d^2), d in divisors(gcd(m,n))); end; # N. J. A. Sloane, Feb 02 2013
  • Mathematica
    t[n_, k_] := Sum[Sum[GCD[c, d], {c, Divisors[n]}], {d, Divisors[k]}]; Table[t[n-k+1, k], {n, 1, 12}, {k, 1, n}] // Flatten (* Jean-François Alcover, May 21 2013 *)
  • Sage
    def A216624(n, k) :
        cp = cartesian_product([divisors(n), divisors(k)])
        return reduce(lambda x,y: x+y, map(gcd, cp))
    for n in (1..12): [A216624(n,k) for k in (1..12)]

Formula

T(n,n) = A060724(n) = sum_{d|n} d*tau((n/d)^2).
T(n,1) = T(1,n) = A000005(n) = tau(n).
T(n,2) = T(2,n) = A060710(n) = sum_{d|n} (3-[d is odd]) (Iverson bracket).
T(n+1,n) = A092517(n) = tau(n+1)*tau(n).
T(prime(n),1) = A007395(n) = 2.
T(prime(n),prime(n)) = A113935(n) = prime(n)+3.

A054584 Number of subgroups of the group generated by a^n=1, b^3=1 and ab=ba.

Original entry on oeis.org

2, 4, 6, 6, 4, 12, 4, 8, 10, 8, 4, 18, 4, 8, 12, 10, 4, 20, 4, 12, 12, 8, 4, 24, 6, 8, 14, 12, 4, 24, 4, 12, 12, 8, 8, 30, 4, 8, 12, 16, 4, 24, 4, 12, 20, 8, 4, 30, 6, 12, 12, 12, 4, 28, 8, 16, 12, 8, 4, 36, 4, 8, 20, 14, 8, 24, 4, 12, 12, 16, 4, 40, 4, 8, 18, 12, 8, 24, 4, 20, 18, 8, 4
Offset: 1

Views

Author

John W. Layman, Apr 12 2000

Keywords

Comments

Also the number of subgroups of the group C_n X C_3 (where C_n is the cyclic group of order n). Number of subgroups of the group C_n X C_m is Sum_{i|n,j|m} gcd(i,j).

Crossrefs

Programs

  • Haskell
    a054584 n = a000005 n + 3 * a079978 n * a000005 (a051176 n) + a035191 n
    -- Reinhard Zumkeller, Aug 27 2012
  • Maple
    for n from 1 to 500 do a := ifactors(n):s := 1:for k from 1 to nops(a[2]) do p := a[2][k][1]:e := a[2][k][2]: if p=3 then b := 2*e+1:else b := e+1:fi:s := s*b:od:printf(`%d,`,2*s); od:
  • Mathematica
    f[d_ /; Mod[d, 3] == 0] = 4; f[] = 2; a[n] := Total[f /@ Divisors[n]]; Table[a[n], {n, 1, 100}](* Jean-François Alcover, Nov 21 2011, after Michael Somos *)
    f[p_, e_] := e + 1; f[3, e_] := 2*e + 1; a[1] = 2; a[n_] := 2*Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 29 2022 *)
  • PARI
    a(n)=if(n<1, 0, sumdiv(n,d, (d%3==0)*2+2)) /* Michael Somos, Sep 20 2005 */
    

Formula

a(n) = tau(n)+3*tau(n/3)+A035191(n) if n is congruent to 0 mod 3 else tau(n)+A035191(n), where A035191(n) is the number of divisors of n that are not congruent to 0 mod 3.
a(n)/2 is multiplicative with a(3^e)=2e+1 and a(p^e)=e+1 for p<>3.
Moebius transform is period 3 sequence [2, 2, 4, ...]. - Michael Somos, Sep 20 2005
G.f.: Sum_{k>0} x^k(2+2*x^k+4*x^(2k))/(1-x^(3k)).
From Amiram Eldar, Nov 29 2022: (Start)
Dirichlet g.f.: 2 * zeta(s)^2 * (1 + 1/3^s).
Sum_{k=1..n} a(k) ~ 2*(4*n*log(n) + (8*gamma - 4 - log(3))*n)/3, where gamma is Euler's constant (A001620). (End)

Extensions

Additional comments from Vladeta Jovovic, Oct 25 2001
Showing 1-4 of 4 results.