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 31-40 of 73 results. Next

A228350 Triangle read by rows: T(j,k) is the k-th part in nonincreasing order of the j-th region of the set of compositions (ordered partitions) of n in colexicographic order, if 1<=j<=2^(n-1) and 1<=k<=A006519(j).

Original entry on oeis.org

1, 2, 1, 1, 3, 2, 1, 1, 1, 2, 1, 1, 4, 3, 2, 2, 1, 1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 1, 1, 2, 1, 1, 5, 4, 3, 3, 2, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 1, 1, 2, 1, 1, 4, 3, 2, 2, 1, 1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 1, 1, 2, 1, 1, 6, 5, 4, 4, 3, 3
Offset: 1

Views

Author

Omar E. Pol, Aug 20 2013

Keywords

Comments

Triangle read by rows in which row n lists the A006519(n) elements of the row A001511(n) of triangle A065120, n >= 1.
The equivalent sequence for integer partitions is A206437.

Examples

			---------------------------------------------------------
.              Diagram                Triangle
Compositions     of            of compositions (rows)
.   of 5       regions          and regions (columns)
----------------------------------------------------------
.             _ _ _ _ _
.         5  |_        |                                5
.       1+4  |_|_      |                              1 4
.       2+3  |_  |     |                            2   3
.     1+1+3  |_|_|_    |                          1 1   3
.       3+2  |_    |   |                        3       2
.     1+2+2  |_|_  |   |                      1 2       2
.     2+1+2  |_  | |   |                    2   1       2
.   1+1+1+2  |_|_|_|_  |                  1 1   1       2
.       4+1  |_      | |                4               1
.     1+3+1  |_|_    | |              1 3               1
.     2+2+1  |_  |   | |            2   2               1
.   1+1+2+1  |_|_|_  | |          1 1   2               1
.     3+1+1  |_    | | |        3       1               1
.   1+2+1+1  |_|_  | | |      1 2       1               1
.   2+1+1+1  |_  | | | |    2   1       1               1
. 1+1+1+1+1  |_|_|_|_|_|  1 1   1       1               1
.
Also the structure could be represented by an isosceles triangle in which the n-th diagonal gives the n-th region. For the composition of 4 see below:
.             _ _ _ _
.         4  |_      |                  4
.       1+3  |_|_    |                1   3
.       2+2  |_  |   |              2       2
.     1+1+2  |_|_|_  |            1   1       2
.       3+1  |_    | |          3               1
.     1+2+1  |_|_  | |        1   2               1
.     2+1+1  |_  | | |      2       1               1
.   1+1+1+1  |_|_|_|_|    1   1       1               1
.
Illustration of the four sections of the set of compositions of 4:
.                                      _ _ _ _
.                                     |_      |     4
.                                     |_|_    |   1+3
.                                     |_  |   |   2+2
.                       _ _ _         |_|_|_  | 1+1+2
.                      |_    |   3          | |     1
.             _ _      |_|_  | 1+2          | |     1
.     _      |_  | 2       | |   1          | |     1
.    |_| 1     |_| 1       |_|   1          |_|     1
.
.
Illustration of initial terms. The parts of the eight regions of the set of compositions of 4:
--------------------------------------------------------
\j:  1      2    3        4     5      6    7          8
k
--------------------------------------------------------
.  _    _ _    _    _ _ _     _    _ _    _    _ _ _ _
1 |_|1 |_  |2 |_|1 |_    |3  |_|1 |_  |2 |_|1 |_      |4
2        |_|1        |_  |2         |_|1        |_    |3
3                      | |1                       |   |2
4                      |_|1                       |_  |2
5                                                   | |1
6                                                   | |1
7                                                   | |1
8                                                   |_|1
.
Triangle begins:
1;
2,1;
1;
3,2,1,1;
1;
2,1;
1;
4,3,2,2,1,1,1,1;
1;
2,1;
1;
3,2,1,1;
1;
2,1;
1;
5,4,3,3,2,2,2,2,1,1,1,1,1,1,1,1;
...
.
Also triangle read by rows T(n,m) in which row n lists the parts of the n-th section of the set of compositions of the integers >= n, ordered by regions. Row lengths give A045623. Row sums give A001792 (see below):
[1];
[2,1];
[1],[3,2,1,1];
[1],[2,1],[1],[4,3,2,2,1,1,1,1];
[1],[2,1],[1],[3,2,1,1],[1],[2,1],[1],[5,4,3,3,2,2,2,2,1,1,1,1,1,1,1,1];
		

Crossrefs

Formula

T(j,k) = A065120(A001511(j)),k) = A001511(j) - A029837(k), 1<=k<=A006519(j), j>=1.

A088839 Numerator of sigma(4n)/sigma(n).

Original entry on oeis.org

7, 5, 7, 31, 7, 5, 7, 21, 7, 5, 7, 31, 7, 5, 7, 127, 7, 5, 7, 31, 7, 5, 7, 21, 7, 5, 7, 31, 7, 5, 7, 85, 7, 5, 7, 31, 7, 5, 7, 21, 7, 5, 7, 31, 7, 5, 7, 127, 7, 5, 7, 31, 7, 5, 7, 21, 7, 5, 7, 31, 7, 5, 7, 511, 7, 5, 7, 31, 7, 5, 7, 21, 7, 5, 7, 31, 7, 5, 7, 127, 7, 5, 7, 31, 7, 5, 7, 21, 7, 5, 7, 31
Offset: 1

Views

Author

Labos Elemer, Nov 04 2003

Keywords

Crossrefs

Programs

  • Maple
    f:= proc(n) local m;
      m:= padic:-ordp(n,2);
      if m::odd then (2^(m+3)-1)/3 else 2^(m+3)-1 fi
    end proc:
    map(f, [$1..200]); # Robert Israel, Nov 19 2017
  • Mathematica
    k=4; Table[Numerator[DivisorSigma[1, k*n]/DivisorSigma[1, n]], {n, 1, 128}]
  • PARI
    A088839(n) = numerator(sigma(4*n)/sigma(n)); \\ Antti Karttunen, Nov 18 2017

Formula

a(n) = (8*A006519(n)-1)/(1+2*A096268(n)). - Robert Israel, Nov 19 2017
From Amiram Eldar, Jan 06 2023: (Start)
a(n) = numerator(A193553(n)/A000203(n)).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k)/A088840(k) = 3*A065442 + 1 = 5.820085... . (End)

Extensions

Typo in definition corrected by Antti Karttunen, Nov 18 2017

A088840 Denominator of sigma(4n)/sigma(n).

Original entry on oeis.org

1, 1, 1, 7, 1, 1, 1, 5, 1, 1, 1, 7, 1, 1, 1, 31, 1, 1, 1, 7, 1, 1, 1, 5, 1, 1, 1, 7, 1, 1, 1, 21, 1, 1, 1, 7, 1, 1, 1, 5, 1, 1, 1, 7, 1, 1, 1, 31, 1, 1, 1, 7, 1, 1, 1, 5, 1, 1, 1, 7, 1, 1, 1, 127, 1, 1, 1, 7, 1, 1, 1, 5, 1, 1, 1, 7, 1, 1, 1, 31, 1, 1, 1, 7, 1, 1, 1, 5, 1, 1, 1, 7, 1, 1, 1, 21, 1, 1, 1, 7, 1, 1
Offset: 1

Views

Author

Labos Elemer, Nov 04 2003

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Denominator[DivisorSigma[1, 4*n]/DivisorSigma[1, n]], {n, 1, 128}]
    a[n_] := Module[{e = IntegerExponent[n, 2]}, (((-1)^e+2)*(2^(e+1)-1))/3]; Array[a, 100] (* Amiram Eldar, Oct 03 2023 *)
  • PARI
    A088840(n) = denominator(sigma(4*n)/sigma(n)); \\ Antti Karttunen, Nov 18 2017
    
  • PARI
    a(n) = {my(e = valuation(n, 2)); (((-1)^e+2) * (2^(e+1)-1))/3;} \\ Amiram Eldar, Oct 03 2023

Formula

From Amiram Eldar, Oct 03 2023: (Start)
Multiplicative with a(2^e) = (((-1)^e+2)*(2^(e+1)-1))/3 = A213243(e+1), and a(p^e) = 1 for an odd prime p.
a(n) = A213243(A007814(n+1)).
Dirichlet g.f.: ((8^s + 4^s + 2^(s+1))/(8^s + 4^s - 2^(s+2) - 4)) * zeta(s).
Sum_{k=1..n} a(k) = (2*n/(3*log(2))) * (log(n) + gamma - 1 + 7*log(2)/12), where gamma is Euler's constant (A001620). (End)

Extensions

Typo in definition corrected by Antti Karttunen, Nov 18 2017

A118821 2-adic continued fraction of zero, where a(n) = 2 if n is odd, -A006519(n/2) otherwise.

Original entry on oeis.org

2, -1, 2, -2, 2, -1, 2, -4, 2, -1, 2, -2, 2, -1, 2, -8, 2, -1, 2, -2, 2, -1, 2, -4, 2, -1, 2, -2, 2, -1, 2, -16, 2, -1, 2, -2, 2, -1, 2, -4, 2, -1, 2, -2, 2, -1, 2, -8, 2, -1, 2, -2, 2, -1, 2, -4, 2, -1, 2, -2, 2, -1, 2, -32, 2, -1, 2, -2, 2, -1, 2, -4, 2, -1, 2, -2, 2, -1, 2, -8, 2, -1, 2, -2, 2, -1, 2, -4, 2, -1, 2, -2, 2, -1, 2, -16
Offset: 1

Views

Author

Paul D. Hanna, May 01 2006

Keywords

Comments

Limit of convergents equals zero; only the 6th convergent is indeterminate. Other 2-adic continued fractions of zero are A118824, A118827, A118830. A006519(n) is the highest power of 2 dividing n; A080277 = partial sums of A038712, where A038712(n) = 2*A006519(n) - 1.

Examples

			For n >= 1, convergents A118822(k)/A118823(k) are:
  at k = 4*n: -1/A080277(n);
  at k = 4*n+1: -2/(2*A080277(n)-1);
  at k = 4*n+2: -1/(A080277(n)-1);
  at k = 4*n-1: 0/(-1)^n.
Convergents begin:
  2/1, -1/-1, 0/-1, -1/1, -2/1, 1/0, 0/1, 1/-4,
  2/-7, -1/3, 0/-1, -1/5, -2/9, 1/-4, 0/1, 1/-12,
  2/-23, -1/11, 0/-1, -1/13, -2/25, 1/-12, 0/1, 1/-16,
  2/-31, -1/15, 0/-1, -1/17, -2/33, 1/-16, 0/1, 1/-32, ...
		

Crossrefs

Cf. A006519, A080277; convergents: A118822/A118823; variants: A118824, A118827, A118830; A100338.

Programs

  • Mathematica
    Array[-2^(IntegerExponent[#, 2] - 1) /. -1/2 -> 2 &, 96] (* Michael De Vlieger, Nov 02 2018 *)
  • PARI
    a(n)=local(p=+2,q=-1);if(n%2==1,p,q*2^valuation(n/2,2))

A118824 2-adic continued fraction of zero, where a(n) = -2 if n is odd, A006519(n/2) otherwise.

Original entry on oeis.org

-2, 1, -2, 2, -2, 1, -2, 4, -2, 1, -2, 2, -2, 1, -2, 8, -2, 1, -2, 2, -2, 1, -2, 4, -2, 1, -2, 2, -2, 1, -2, 16, -2, 1, -2, 2, -2, 1, -2, 4, -2, 1, -2, 2, -2, 1, -2, 8, -2, 1, -2, 2, -2, 1, -2, 4, -2, 1, -2, 2, -2, 1, -2, 32, -2, 1, -2, 2, -2, 1, -2, 4, -2, 1, -2, 2, -2, 1, -2, 8, -2, 1, -2, 2, -2, 1, -2, 4, -2, 1, -2, 2, -2, 1, -2, 16, -2, 1, -2, 2, -2, 1
Offset: 1

Views

Author

Paul D. Hanna, May 01 2006

Keywords

Comments

Limit of convergents equals zero; only the 6th convergent is indeterminate. Other 2-adic continued fractions of zero are: A118821, A118827, A118830. A006519(n) is the highest power of 2 dividing n; A080277 = partial sums of A038712, where A038712(n) = 2*A006519(n) - 1.

Examples

			For n >= 1, convergents A118825(k)/A118826(k):
  at k = 4*n: 1/A080277(n);
  at k = 4*n+1: 2/(2*A080277(n)-1);
  at k = 4*n+2: 1/(A080277(n)-1);
  at k = 4*n-1: 0.
Convergents begin:
  -2/1, -1/1, 0/-1, -1/-1, 2/1, 1/0, 0/1, 1/4,
  -2/-7, -1/-3, 0/-1, -1/-5, 2/9, 1/4, 0/1, 1/12,
  -2/-23, -1/-11, 0/-1, -1/-13, 2/25, 1/12, 0/1, 1/16,
  -2/-31, -1/-15, 0/-1, -1/-17, 2/33, 1/16, 0/1, 1/32, ...
		

Crossrefs

Cf. A006519, A080277; convergents: A118825/A118826; variants: A118821, A118827, A118830; A100338.

Programs

  • Mathematica
    Array[If[OddQ@ #, -2, 2^(IntegerExponent[#, 2] - 1)] &, 102] (* Michael De Vlieger, Nov 06 2018 *)
  • PARI
    a(n)=local(p=-2,q=+1);if(n%2==1,p,q*2^valuation(n/2,2))

A118827 2-adic continued fraction of zero, where a(n) = 1 if n is odd, otherwise -2*A006519(n/2).

Original entry on oeis.org

1, -2, 1, -4, 1, -2, 1, -8, 1, -2, 1, -4, 1, -2, 1, -16, 1, -2, 1, -4, 1, -2, 1, -8, 1, -2, 1, -4, 1, -2, 1, -32, 1, -2, 1, -4, 1, -2, 1, -8, 1, -2, 1, -4, 1, -2, 1, -16, 1, -2, 1, -4, 1, -2, 1, -8, 1, -2, 1, -4, 1, -2, 1, -64, 1, -2, 1, -4, 1, -2, 1, -8, 1, -2, 1, -4, 1, -2, 1, -16, 1, -2, 1, -4, 1, -2, 1, -8, 1, -2, 1, -4, 1, -2, 1, -32, 1, -2, 1
Offset: 1

Views

Author

Paul D. Hanna, May 01 2006

Keywords

Comments

Limit of convergents equals zero; only the 6th convergent is indeterminate. Other 2-adic continued fractions of zero are: A118821, A118824, A118830. A006519(n) is the highest power of 2 dividing n; A080277 = partial sums of A038712, where A038712(n) = 2*A006519(n) - 1.
Multiplicative because both A006519 and A165326 are. - Andrew Howroyd, Aug 01 2018

Examples

			For n >= 1, convergents A118828(k)/A118829(k):
  at k = 4*n: -1/(2*A080277(n));
  at k = 4*n+1: -1/(2*A080277(n)-1);
  at k = 4*n+2: -1/(2*A080277(n)-2);
  at k = 4*n-1: 0.
Convergents begin:
  1/1, -1/-2, 0/-1, -1/2, -1/1, 1/0, 0/1, 1/-8,
  1/-7, -1/6, 0/-1, -1/10, -1/9, 1/-8, 0/1, 1/-24,
  1/-23, -1/22, 0/-1, -1/26, -1/25, 1/-24, 0/1, 1/-32,
  1/-31, -1/30, 0/-1, -1/34, -1/33, 1/-32, 0/1, 1/-64, ...
		

Crossrefs

Programs

  • Mathematica
    Array[If[OddQ@ #, 1, -2*2^(IntegerExponent[#, 2] - 1)] &, 99] (* Michael De Vlieger, Nov 06 2018 *)
  • PARI
    a(n)=local(p=+1,q=-2);if(n%2==1,p,q*2^valuation(n/2,2))

Formula

a(n) = A165326(n) * A006519(n). - Andrew Howroyd, Aug 01 2018
From Amiram Eldar, Oct 28 2023: (Start)
Multiplicative with a(2^e) = -2^e, and a(p^e) = 1 for an odd prime p.
Dirichlet g.f.: zeta(s) * (1 - 2^(1-s) + 1/(2-2^s)).
Sum_{k=1..n} a(k) ~ (-1/(2*log(2))) * n *(log(n) + gamma - log(2)/2 - 1), where gamma is Euler's constant (A001620). (End)

A118830 2-adic continued fraction of zero, where a(n) = -1 if n is odd, 2*A006519(n/2) otherwise.

Original entry on oeis.org

-1, 2, -1, 4, -1, 2, -1, 8, -1, 2, -1, 4, -1, 2, -1, 16, -1, 2, -1, 4, -1, 2, -1, 8, -1, 2, -1, 4, -1, 2, -1, 32, -1, 2, -1, 4, -1, 2, -1, 8, -1, 2, -1, 4, -1, 2, -1, 16, -1, 2, -1, 4, -1, 2, -1, 8, -1, 2, -1, 4, -1, 2, -1, 64, -1, 2, -1, 4, -1, 2, -1, 8, -1, 2, -1, 4, -1, 2, -1, 16, -1, 2, -1, 4, -1, 2, -1, 8, -1, 2, -1, 4, -1, 2, -1, 32, -1, 2, -1
Offset: 1

Views

Author

Paul D. Hanna, May 01 2006

Keywords

Comments

Limit of convergents equals zero; only the 6th convergent is indeterminate. Other 2-adic continued fractions of zero are: A118821, A118824, A118827. A006519(n) is the highest power of 2 dividing n; A080277 = partial sums of A038712, where A038712(n) = 2*A006519(n) - 1.

Examples

			For n >= 1, convergents A118831(k)/A118832(k):
  at k = 4*n: 1/(2*A080277(n));
  at k = 4*n+1: 1/(2*A080277(n)-1);
  at k = 4*n+2: 1/(2*A080277(n)-2);
  at k = 4*n-1: 0.
Convergents begin:
  -1/1, -1/2, 0/-1, -1/-2, 1/1, 1/0, 0/1, 1/8,
  -1/-7, -1/-6, 0/-1, -1/-10, 1/9, 1/8, 0/1, 1/24,
  -1/-23, -1/-22, 0/-1, -1/-26, 1/25, 1/24, 0/1, 1/32,
  -1/-31, -1/-30, 0/-1, -1/-34, 1/33, 1/32, 0/1, 1/64, ...
		

Crossrefs

Cf. A006519, A080277; convergents: A118831/A118832; variants: A118821, A118824, A118827; A100338.

Programs

  • Mathematica
    Array[If[OddQ@ #, -1, 2^IntegerExponent[#, 2]] &, 99] (* Michael De Vlieger, Nov 06 2018 *)
  • PARI
    a(n)=local(p=-1,q=+2);if(n%2==1,p,q*2^valuation(n/2,2))

A385006 The sum of the biquadratefree divisors of n.

Original entry on oeis.org

1, 3, 4, 7, 6, 12, 8, 15, 13, 18, 12, 28, 14, 24, 24, 15, 18, 39, 20, 42, 32, 36, 24, 60, 31, 42, 40, 56, 30, 72, 32, 15, 48, 54, 48, 91, 38, 60, 56, 90, 42, 96, 44, 84, 78, 72, 48, 60, 57, 93, 72, 98, 54, 120, 72, 120, 80, 90, 60, 168, 62, 96, 104, 15, 84, 144
Offset: 1

Views

Author

Amiram Eldar, Jun 15 2025

Keywords

Comments

First differs from A365682 and A366992 at n = 32.
The number of these divisors is A252505(n), and the largest of them is A058035(n).

Crossrefs

The sum of divisors d of n such that d is: A000593 (odd), A033634 (exponentially odd), A035316 (square), A038712 (power of 2), A048250 (squarefree), A072079 (3-smooth), A073185 (cubefree), A113061 (cube), A162296 (nonsquarefree), A183097 (powerful), A186099 (5-rough), A353900 (exponentially 2^n), A385005 (cubefull), this sequence (biquadratefree).

Programs

  • Mathematica
    f[p_, e_] := (p^Min[e+1, 4] - 1)/(p - 1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    a(n) = {my(f = factor(n), p, e); prod(i = 1, #f~, p = f[i,1]; e = f[i,2]; (p^min(e+1, 4) - 1)/(p - 1));}

Formula

Multiplicative with a(p^e) = (p^min(e+1, 4) - 1)/(p - 1).
In general, the sum of the k-free (numbers that are not divisible by a k-th power larger than 1) divisors of n is multiplicative with a(p^e) = (p^min(e+1, k) - 1)/(p - 1).
Dirichlet g.f.: zeta(s) * zeta(s-1) /zeta(4*s-4).
In general, the sum of the k-free divisors of n has Dirichlet g.f.: zeta(s)*zeta(s-1)/zeta(k*s-k).
Sum_{k=1..n} a(k) ~ (15/(2*Pi^2)) * n^2.
In general, the sum of the k-free divisors of n has an average order (Pi^2/(12*zeta(k))) * n^2.

A385047 The sum of the unitary divisors of n that are powers of 2.

Original entry on oeis.org

1, 3, 1, 5, 1, 3, 1, 9, 1, 3, 1, 5, 1, 3, 1, 17, 1, 3, 1, 5, 1, 3, 1, 9, 1, 3, 1, 5, 1, 3, 1, 33, 1, 3, 1, 5, 1, 3, 1, 9, 1, 3, 1, 5, 1, 3, 1, 17, 1, 3, 1, 5, 1, 3, 1, 9, 1, 3, 1, 5, 1, 3, 1, 65, 1, 3, 1, 5, 1, 3, 1, 9, 1, 3, 1, 5, 1, 3, 1, 17, 1, 3, 1, 5, 1, 3
Offset: 1

Views

Author

Amiram Eldar, Jun 16 2025

Keywords

Crossrefs

The unitary analog of A038712.
The sum of unitary divisors of n that are: A092261 (squarefree), A192066 (odd), A358346 (exponentially odd), A358347 (square), A360720 (powerful), A371242 (cubefree), A380396 (cube), A383763 (exponentially squarefree), A385043 (exponentially 2^n), A385045 (5-rough), A385046 (3-smooth), this sequence (power of 2), A385048 (cubefull), A385049 (biquadratefree).

Programs

  • Mathematica
    a[n_] := If[OddQ[n], 1, 2^IntegerExponent[n, 2] + 1]; Array[a, 100]
  • PARI
    a(n) = if(n%2, 1, 2^valuation(n, 2)+1);

Formula

Multiplicative with a(2^e) = 2^e + 1, and a(p^e) = 1 for an odd prime p.
a(n) = A034448(n) / A192066(n).
a(n) = A059841(n) + A006519(n), i.e., a(n) = A006519(n) + 1 if n is even, and 1 is n is odd.
Dirichlet g.f.: zeta(s) * ((1-1/2^(2*s-1))/(1-1/2^(s-1))).
Sum_{k=1..n} a(k) ~ (n/(2*log(2))) * (log(n) + gamma - 1 + 5*log(2)/2), where gamma is Euler's constant (A001620).

A088842 Denominator of the quotient sigma(7n)/sigma(n).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 57, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 57, 1, 1, 1, 1, 1, 1, 8
Offset: 1

Views

Author

Labos Elemer, Nov 04 2003

Keywords

Comments

Sum of powers of 7 dividing n. - Amiram Eldar, Nov 27 2022

Crossrefs

Cf. A000203 (sigma), A001620, A088841 (numerators), A283078 (sigma(7n)).

Programs

  • Mathematica
    Table[Denominator[DivisorSigma[1, 7*n]/DivisorSigma[1, n]], {n, 1, 128}] (* corrected by Ilya Gutkovskiy, Dec 15 2020 *)
    a[n_] := (7^(IntegerExponent[n, 7] + 1) - 1)/6; Array[a, 100] (* Amiram Eldar, Nov 27 2022 *)
  • PARI
    a(n) = denominator(sigma(7*n)/sigma(n)); \\ Michel Marcus, Dec 15 2020
    
  • PARI
    a(n) = (7^(valuation(n, 7) + 1) - 1)/6; \\ Amiram Eldar, Nov 27 2022

Formula

G.f.: Sum_{k>=0} 7^k * x^(7^k) / (1 - x^(7^k)). - Ilya Gutkovskiy, Dec 15 2020
From Amiram Eldar, Nov 27 2022: (Start)
Multiplicative with a(7^e) = (7^(e+1)-1)/6, and a(p^e) = 1 for p != 7.
Dirichlet g.f.: zeta(s) / (1 - 7^(1 - s)).
Sum_{k=1..n} a(k) ~ n*log_7(n) + (1/2 + (gamma - 1)/log(7))*n, where gamma is Euler's constant (A001620). (End)
Previous Showing 31-40 of 73 results. Next