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

A348733 a(n) = gcd(A003959(n), A034448(n)), where A003959 is multiplicative with a(p^e) = (p+1)^e and A034448 (usigma) is multiplicative with a(p^e) = (p^e)+1.

Original entry on oeis.org

1, 3, 4, 1, 6, 12, 8, 9, 2, 18, 12, 4, 14, 24, 24, 1, 18, 6, 20, 6, 32, 36, 24, 36, 2, 42, 4, 8, 30, 72, 32, 3, 48, 54, 48, 2, 38, 60, 56, 54, 42, 96, 44, 12, 12, 72, 48, 4, 2, 6, 72, 14, 54, 12, 72, 72, 80, 90, 60, 24, 62, 96, 16, 1, 84, 144, 68, 18, 96, 144, 72, 18, 74, 114, 8, 20, 96, 168, 80, 6, 2, 126, 84, 32
Offset: 1

Views

Author

Antti Karttunen, Nov 05 2021

Keywords

Comments

This is not multiplicative. The first point where a(m*n) = a(m)*a(n) does not hold for coprime m and n is 1444 = 2^2 * 19^2, where a(1444) = 10 != 1*2 = a(4)*a(361). See A348740 for the list of such positions.

Links

Crossrefs

Cf. A003959, A034448, A348732, A348734, A348735, A348740.
Cf. also A344695, A348047, A348503, A348946 for similar, almost multiplicative sequences.

Programs

  • Mathematica
    f1[p_, e_] := (p + 1)^e; f2[p_, e_] := p^e + 1; a[1] = 1; a[n_] := GCD[Times @@ f1 @@@ (f = FactorInteger[n]), Times @@ f2 @@@ f]; Array[a, 100] (* Amiram Eldar, Nov 05 2021 *)
  • PARI
    A003959(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1]++); factorback(f); };
A034448(n) = { my(f = factor(n)); prod(k=1, #f~, 1+(f[k, 1]^f[k, 2])); };
A348733(n) = gcd(A003959(n), A034448(n));

Formula

a(n) = gcd(A003959(n), A034448(n)).
a(n) = gcd(A003959(n), A348732(n)) = gcd(A034448(n), A348732(n)).
a(n) = A003959(n) / A348734(n) = A034448(n) / A348735(n).

A348503 a(n) = gcd(sigma(n), usigma(n)), where sigma is the sum of divisors function, A000203, and usigma is the unitary sigma, A034448.

Original entry on oeis.org

1, 3, 4, 1, 6, 12, 8, 3, 1, 18, 12, 4, 14, 24, 24, 1, 18, 3, 20, 6, 32, 36, 24, 12, 1, 42, 4, 8, 30, 72, 32, 3, 48, 54, 48, 1, 38, 60, 56, 18, 42, 96, 44, 12, 6, 72, 48, 4, 1, 3, 72, 14, 54, 12, 72, 24, 80, 90, 60, 24, 62, 96, 8, 1, 84, 144, 68, 18, 96, 144, 72, 15, 74, 114, 4, 20, 96, 168, 80, 6, 1, 126, 84, 32
Offset: 1

Views

Author

Antti Karttunen, Oct 29 2021

Keywords

Comments

This is not multiplicative. The first point where a(m*n) = a(m)*a(n) does not hold for coprime m and n is 72 = 8*9, where a(72) = 15 != 3*1 = a(8)*a(9).

Links

Crossrefs

Cf. A000203, A034448, A048146, A348504, A348505.
Differs from A344695 for the first time at n=72, where a(72) = 15, while A344695(72) = 3.
Differs from A348047 for the first time at n=27, where a(27) = 4, while A348047(27) = 8.

Programs

  • Mathematica
    f1[p_, e_] := p^e + 1; f2[p_, e_] := (p^(e + 1) - 1)/(p - 1); a[1] = 1; a[n_] := GCD[Times @@ f1 @@@ (fct = FactorInteger[n]), Times @@ f2 @@@ fct]; Array[a, 100] (* Amiram Eldar, Oct 29 2021 *)
  • PARI
    A034448(n) = { my(f=factorint(n)); prod(k=1, #f~, 1+(f[k, 1]^f[k, 2])); }; \\ After code in A034448
A348503(n) = gcd(sigma(n), A034448(n));

Formula

a(n) = gcd(A000203(n), A034448(n)).
a(n) = gcd(A000203(n), A048146(n)) = gcd(A034448(n), A048146(n)).
a(n) = A000203(n) / A348504(n) = A034448(n) / A348505(n).

A348946 a(n) = gcd(sigma(n), A348944(n)), where A348944 is the arithmetic mean of A003959 and A034448, and sigma is the sum of divisors function.

Original entry on oeis.org

1, 3, 4, 7, 6, 12, 8, 3, 13, 18, 12, 28, 14, 24, 24, 1, 18, 39, 20, 42, 32, 36, 24, 12, 31, 42, 2, 56, 30, 72, 32, 3, 48, 54, 48, 1, 38, 60, 56, 18, 42, 96, 44, 84, 78, 72, 48, 4, 57, 93, 72, 98, 54, 6, 72, 24, 80, 90, 60, 168, 62, 96, 104, 1, 84, 144, 68, 126, 96, 144, 72, 3, 74, 114, 124, 140, 96, 168, 80, 6, 1, 126
Offset: 1

Views

Author

Antti Karttunen, Nov 05 2021

Keywords

Comments

This is not multiplicative. The first point where a(m*n) = a(m)*a(n) does not hold for coprime m and n is 36 = 2^2 * 3^2, where a(36) = 1 <> 91 = 7*13 = a(4)*a(9).

Links

Crossrefs

Cf. A000203, A003959, A034448, A348944, A348945, A348947, A348948.
Cf. also A344695, A348047, A348503, A348733.

Programs

  • Mathematica
    f1[p_, e_] := (p^(e + 1) - 1)/(p - 1); f2[p_, e_] := (p + 1)^e; f3[p_, e_] := p^e + 1; a[1] = 1; a[n_] := GCD[Times @@ f1 @@@ (f = FactorInteger[n]), (Times @@ f2 @@@ f + Times @@ f3 @@@ f)/2]; Array[a, 100] (* Amiram Eldar, Nov 05 2021 *)
  • PARI
    A003959(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1]++); factorback(f); };
A034448(n) = { my(f = factor(n)); prod(k=1, #f~, 1+(f[k, 1]^f[k, 2])); };
A348946(n) = gcd(sigma(n), ((1/2)*(A003959(n)+A034448(n))));

Formula

a(n) = gcd(A000203(n), A348944(n)).
a(n) = gcd(A000203(n), A348945(n)) = gcd(A348944(n), A348945(n));
a(n) = A348944(n) / A348947(n) = A000203(n) / A348948(n).

A348048 a(n) = sigma(n) / gcd(sigma(n), A003959(n)), where A003959 is multiplicative with a(p^e) = (p+1)^e and sigma is the sum of divisors function.

Original entry on oeis.org

1, 1, 1, 7, 1, 1, 1, 5, 13, 1, 1, 7, 1, 1, 1, 31, 1, 13, 1, 7, 1, 1, 1, 5, 31, 1, 5, 7, 1, 1, 1, 7, 1, 1, 1, 91, 1, 1, 1, 5, 1, 1, 1, 7, 13, 1, 1, 31, 57, 31, 1, 7, 1, 5, 1, 5, 1, 1, 1, 7, 1, 1, 13, 127, 1, 1, 1, 7, 1, 1, 1, 65, 1, 1, 31, 7, 1, 1, 1, 31, 121, 1, 1, 7, 1, 1, 1, 5, 1, 13, 1, 7, 1, 1, 1, 7, 1, 57, 13
Offset: 1

Views

Author

Antti Karttunen, Oct 21 2021

Keywords

Comments

Not multiplicative. For example, a(196) = 133 != a(4) * a(49).

Links

Crossrefs

Cf. A000203, A003959, A348029, A348047, A348049.
Cf. also A344696.

Programs

  • Mathematica
    f[p_, e_] := (p + 1)^e; a[1] = 1; a[n_] := (s = DivisorSigma[1, n]) / GCD[s, Times @@ f @@@ FactorInteger[n]]; Array[a, 100] (* Amiram Eldar, Oct 21 2021 *)
  • PARI
    A003959(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1]++); factorback(f); };
A348048(n) = { my(u=sigma(n)); (u/gcd(u, A003959(n))); };

Formula

a(n) = A000203(n) / A348047(n) = A000203(n) / gcd(A000203(n), A003959(n)).

A348049 a(n) = A003959(n) / gcd(sigma(n), A003959(n)), where A003959 is multiplicative with a(p^e) = (p+1)^e and sigma is the sum of divisors function.

Original entry on oeis.org

1, 1, 1, 9, 1, 1, 1, 9, 16, 1, 1, 9, 1, 1, 1, 81, 1, 16, 1, 9, 1, 1, 1, 9, 36, 1, 8, 9, 1, 1, 1, 27, 1, 1, 1, 144, 1, 1, 1, 9, 1, 1, 1, 9, 16, 1, 1, 81, 64, 36, 1, 9, 1, 8, 1, 9, 1, 1, 1, 9, 1, 1, 16, 729, 1, 1, 1, 9, 1, 1, 1, 144, 1, 1, 36, 9, 1, 1, 1, 81, 256, 1, 1, 9, 1, 1, 1, 9, 1, 16, 1, 9, 1, 1, 1, 27, 1, 64
Offset: 1

Views

Author

Antti Karttunen, Oct 21 2021

Keywords

Comments

Not multiplicative. For example, a(196) = 192 != a(4) * a(49).

Links

Crossrefs

Cf. A000203, A003959, A005117 (positions of 1's), A348029, A348047, A348048.
Cf. also A344697.

Programs

  • Mathematica
    f[p_, e_] := (p + 1)^e; a[1] = 1; a[n_] := (m = Times @@ f @@@ FactorInteger[n]) / GCD[m, DivisorSigma[1, n]]; Array[a, 100] (* Amiram Eldar, Oct 21 2021 *)
  • PARI
    A003959(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1]++); factorback(f); };
A348049(n) = { my(u=A003959(n)); (u/gcd(u, sigma(n))); };

Formula

a(n) = A003959(n) / A348047(n) = A003959(n) / gcd(A000203(n), A003959(n)).

A348984 a(n) = gcd(sigma(n), A325973(n)), where A325973 is the arithmetic mean of {sum of squarefree divisors} and {sum of unitary divisors}.

Original entry on oeis.org

1, 3, 4, 1, 6, 12, 8, 3, 1, 18, 12, 4, 14, 24, 24, 1, 18, 3, 20, 6, 32, 36, 24, 12, 1, 42, 8, 8, 30, 72, 32, 9, 48, 54, 48, 1, 38, 60, 56, 18, 42, 96, 44, 12, 6, 72, 48, 4, 1, 3, 72, 14, 54, 24, 72, 24, 80, 90, 60, 24, 62, 96, 8, 1, 84, 144, 68, 18, 96, 144, 72, 3, 74, 114, 4, 20, 96, 168, 80, 6, 1, 126, 84, 32, 108
Offset: 1

Views

Author

Antti Karttunen, Nov 06 2021

Keywords

Comments

This is not multiplicative. The first point where a(m*n) = a(m)*a(n) does not hold for coprime m and n is 108 = 4*27, where a(108) = 4, although a(4) = 1 and a(27) = 8.

Links

Crossrefs

Cf. A000203, A048250, A034448, A325973, A325974, A348985, A348986.
Differs from A348047 for the first time at n=108, where a(108) = 4, while A348047(108) = 8.
Cf. also A348733, A348946.

Programs

  • Mathematica
    f1[p_, e_] := p + 1; f2[p_, e_] := p^e + 1; s[1] = 1; s[n_] := (Times @@ f1 @@@ (f = FactorInteger[n]) + Times @@ f2 @@@ f)/2; a[n_] := GCD[s[n], DivisorSigma[1, n]]; Array[a, 100] (* Amiram Eldar, Nov 06 2021 *)
  • PARI
    A325973(n) = (1/2)*sumdiv(n, d, d*(issquarefree(d) + (1==gcd(d, n/d))));
A348984(n) = gcd(sigma(n), A325973(n));

Formula

a(n) = gcd(A000203(n), A325973(n)).
a(n) = gcd(A000203(n), A325974(n)) = gcd(A325973(n), A325974(n)).
a(n) = A000203(n) / A348985(n) = A325973(n) / A348986(n).
Showing 1-6 of 6 results.

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

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