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

A324709 Larger of tri-unitary amicable numbers pair: numbers (m, n) such that tsigma(m) = tsigma(n) = m + n, where tsigma(n) is the sum of the tri-unitary divisors of n (A324706).

Original entry on oeis.org

126, 846, 1260, 8460, 11760, 10856, 14595, 17700, 49308, 83142, 62700, 71145, 73962, 83904, 107550, 88730, 131100, 168730, 149952, 196650, 203432, 306612, 365700, 399592, 419256, 548550, 721962, 669688, 831420, 686072, 691256, 712216, 652664, 661824, 827700
Offset: 1

Views

Author

Amiram Eldar, Mar 11 2019

Keywords

Comments

The terms are ordered according to their lesser counterparts (A324708).

Examples

			126 is in the sequence since it is the larger of the amicable pair (114, 126): tsigma(114) = tsigma(126) = 114 + 126.
		

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := If[e == 3, (p^4-1)/(p-1), If[e==6, (p^8-1)/(p^2-1), p^e+1]]; tsigma[1]=1; tsigma[n_]:= Times @@ f @@@ FactorInteger[n]; s[n_] := tsigma[n] - n; seq={}; Do[m=s[n]; If[m>n && s[m]==n, AppendTo[seq, m]] ,{n,1,700000}]; seq

A324707 Tri-unitary perfect numbers: numbers k such that tsigma(k) = 2k, where tsigma(k) is the sum of the tri-unitary divisors of k (A324706).

Original entry on oeis.org

6, 60, 90, 36720, 47520, 8173440, 22276800, 126463680, 597542400, 4201148160, 287704872000, 1632485836800
Offset: 1

Views

Author

Amiram Eldar, Mar 11 2019

Keywords

Comments

Also in the sequence is 21623407345626345971712000.
a(13) > 5*10^12. - Giovanni Resta, Mar 14 2019

Examples

			36720 is in the sequence since its sum of tri-unitary divisors is A324706(36720) = 73440 = 2 * 36720.
		

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := If[e == 3, (p^4-1)/(p-1), If[e==6, (p^8-1)/(p^2-1), p^e+1]]; tsigma[1]=1; tsigma[n_]:= Times @@ f @@@ FactorInteger[n]; Select[Range[50000], tsigma[#]==2# &]

Extensions

a(11)-a(12) from Giovanni Resta, Mar 14 2019

A324706 The sum of the tri-unitary divisors of n.

Original entry on oeis.org

1, 3, 4, 5, 6, 12, 8, 15, 10, 18, 12, 20, 14, 24, 24, 17, 18, 30, 20, 30, 32, 36, 24, 60, 26, 42, 40, 40, 30, 72, 32, 33, 48, 54, 48, 50, 38, 60, 56, 90, 42, 96, 44, 60, 60, 72, 48, 68, 50, 78, 72, 70, 54, 120, 72, 120, 80, 90, 60, 120, 62, 96, 80, 85, 84, 144
Offset: 1

Views

Author

Amiram Eldar, Mar 11 2019

Keywords

Comments

A divisor d of n is tri-unitary if the greatest common bi-unitary divisor of d and n/d is 1.

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := If[e == 3, (p^4-1)/(p-1), If[e==6, (p^8-1)/(p^2-1), p^e+1]]; a[1]=1; a[n_]:= Times @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    A324706(n) = { my(f = factor(n)); prod(i=1, #f~, if(3==f[i,2], sigma(f[i,1]^f[i,2]), if(6==f[i,2], ((f[i,1]^8)-1)/((f[i,1]^2)-1), 1+(f[i,1]^f[i,2])))); }; \\ Antti Karttunen, Mar 12 2019

Formula

Multiplicative with a(p^3) = 1 + p + p^2 + p^3, a(p^6) = 1 + p^2 + p^4 + p^6, and a(p^e) = 1 + p^e otherwise.
Sum_{k=1..n} a(k) ~ c * n^2, where c = (Pi^2/12) * Product_{p prime} (1 - 1/p^3 + 1/p^4 - 2/p^6 + 2/p^8 - 1/p^9 - 1/p^12 + 1/p^13) = 0.72189237802... . - Amiram Eldar, Nov 24 2022

A348343 Smaller member of a noninfinitary amicable pair: numbers (k, m) such that nisigma(k) = m and nisigma(m) = k, where nisigma(k) is the sum of the noninfinitary divisors of k (A348271).

Original entry on oeis.org

336, 1792, 5376, 6096, 21504, 32004, 97536, 34062336, 64512000, 118008576, 30064771072
Offset: 1

Views

Author

Amiram Eldar, Oct 13 2021

Keywords

Comments

The larger counterparts are in A348344.

Examples

			336 is a term since A348271(336) = 448 and A348271(448) = 336.
		

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := Module[{b = IntegerDigits[e, 2], m}, m = Length[b]; Product[If[b[[j]] > 0, 1 + p^(2^(m - j)), 1], {j, 1, m}]]; isigma[1] = 1; isigma[n_] := Times @@ f @@@ FactorInteger[n]; s[n_] := DivisorSigma[1,n] - isigma[n]; seq={}; Do[m=s[n]; If[m>n && s[m]==n, AppendTo[seq, n]], {n,1,10^4}]; seq

A335385 The number of tri-unitary divisors of n.

Original entry on oeis.org

1, 2, 2, 2, 2, 4, 2, 4, 2, 4, 2, 4, 2, 4, 4, 2, 2, 4, 2, 4, 4, 4, 2, 8, 2, 4, 4, 4, 2, 8, 2, 2, 4, 4, 4, 4, 2, 4, 4, 8, 2, 8, 2, 4, 4, 4, 2, 4, 2, 4, 4, 4, 2, 8, 4, 8, 4, 4, 2, 8, 2, 4, 4, 4, 4, 8, 2, 4, 4, 8, 2, 8, 2, 4, 4, 4, 4, 8, 2, 4, 2, 4, 2, 8, 4, 4, 4
Offset: 1

Views

Author

Amiram Eldar, Jun 04 2020

Keywords

Comments

A divisor d of k is tri-unitary if the greatest common bi-unitary divisor of d and k/d is 1.
Differs from A037445 at n = 32, 96, 128, 160, 224, ...

Examples

			a(4) = 2 since 4 has 2 tri-unitary divisors, 1 and 4. 2 is not a tri-unitary divisor of 4 since the greatest common bi-unitary divisor of 2 and 4/2 = 2 is 2 and not 1.
		

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := If[e == 3 || e == 6, 4, 2]; a[1] = 1; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100]
  • PARI
    a(n) = vecprod(apply(x -> if(x == 3 || x == 6, 4, 2), factor(n)[, 2])); \\ Amiram Eldar, Dec 18 2023

Formula

Multiplicative with a(p^e) = 4 if e = 3 or 6, and a(p^e) = 2 otherwise.

A357495 Lesser of a pair of amicable numbers k < m such that s(k) = m and s(m) = k, where s(k) = A162296(k) - k is the sum of aliquot divisors of k that have a square factor.

Original entry on oeis.org

880, 10480, 20080, 24928, 42976, 69184, 110565, 252080, 267712, 489472, 566656, 569240, 603855, 626535, 631708, 687424, 705088, 741472, 786896, 904365, 1100385, 1234480, 1280790, 1425632, 1749824, 1993750, 2012224, 2401568, 2439712, 2496736, 2542496, 2573344, 2671856
Offset: 1

Views

Author

Amiram Eldar, Oct 01 2022

Keywords

Comments

Analogous to amicable numbers (A002025 and A002046) with nonsquarefree divisors.
The larger counterparts are in A357496.
Both members of each pair are necessarily nonsquarefree numbers.

Examples

			880 is a term since s(880) = 1136 and s(1136) = 880.
		

Crossrefs

Programs

  • Mathematica
    s[n_] := Module[{f = FactorInteger[n], p, e}, p = f[[;; , 1]]; e = f[[;; , 2]]; Times @@ ((p^(e + 1) - 1)/(p - 1)) - Times @@ (p + 1) - n]; seq = {}; Do[m = s[n]; If[m > n && s[m] == n, AppendTo[seq, n]], {n, 2, 3*10^6}]; seq

A371419 Lesser member of Carmichael's variant of amicable pair: numbers k < m such that s(k) = m and s(m) = k, where s(k) = A371418(k).

Original entry on oeis.org

12, 48, 112, 160, 192, 448, 1984, 12288, 28672, 126976, 196608, 458752, 520192, 786432, 1835008, 2031616, 8126464, 8323072, 33292288, 536805376, 2147221504, 3221225472, 7516192768, 33285996544, 34359476224, 136365211648
Offset: 1

Views

Author

Amiram Eldar, Mar 23 2024

Keywords

Comments

Analogous to amicable numbers (A002025 and A002046) with the largest aliquot divisor of the sum of divisors (A371418) instead of the sum of aliquot divisors (A001065).
Carmichael (1921) proposed this function (A371418) for the purpose of studying periodic chains that are formed by repeatedly applying the mapping x -> A371418(x). The chains of cycle 2 are analogous to amicable numbers.
Carmichael noted that if q < p are two different Mersenne exponents (A000043), then 2^(p-1)*(2^q-1) and 2^(q-1)*(2^p-1) are an amicable pair. With the 51 Mersenne exponents that are currently known it is possible to calculate 51 * 50 / 2 = 1275 amicable pairs. (160, 189) is a pair that is not of this "Mersenne form". Are there any other pairs like it? There are no other such pairs with lesser member below a(26).
a(27) <= 8795019280384.
The greater counterparts are in A371420.

Examples

			12 is a term since A371418(12) = 14 > 12, and A371418(14) = 12.
		

Crossrefs

Programs

  • Mathematica
    r[n_] := n/FactorInteger[n][[1, 1]]; s[n_] := r[DivisorSigma[1, n]]; seq = {}; Do[m = s[n]; If[m > n && s[m] == n, AppendTo[seq, n]], {n, 1, 10^6}]; seq
  • PARI
    f(n) = {my(s = sigma(n)); if(s == 1, 1, s/factor(s)[1, 1]);}
    lista(nmax) = {my(m); for(n = 1, nmax, m = f(n); if(m > n && f(m) == n, print1(n, ", ")));}

A348602 Smaller member of a nonexponential amicable pair: numbers (k, m) such that nesigma(k) = m and nesigma(m) = k, where nesigma(k) is the sum of the nonexponential divisors of k (A160135).

Original entry on oeis.org

198, 18180, 142310, 1077890, 1156870, 1511930, 1669910, 2236570, 2728726, 3776580, 4246130, 4532710, 5123090, 5385310, 6993610, 7288930, 8619765, 8754130, 8826070, 9478910, 10254970, 14426230, 17041010, 17257695, 21448630, 30724694, 34256222, 35361326, 37784810
Offset: 1

Views

Author

Amiram Eldar, Oct 25 2021

Keywords

Comments

The larger counterparts are in A348603.

Examples

			198 is a term since A160135(198) = 204 and A160135(204) = 198.
		

Crossrefs

Programs

  • Mathematica
    esigma[n_] := Times @@ (Sum[First[#]^d, {d, Divisors[Last[#]]}] &) /@ FactorInteger[n]; s[n_] := DivisorSigma[1, n] - esigma[n]; seq = {}; Do[m = s[n]; If[m > n && s[m] == n, AppendTo[seq, n]], {n, 1, 1.7*10^6}]; seq
Showing 1-8 of 8 results.