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.

A322541 Lesser of semi-unitary amicable numbers pair: numbers (m, n) such that susigma(m) = susigma(n) = m + n, where susigma(n) is the sum of the semi-unitary divisors of n (A322485).

Original entry on oeis.org

114, 366, 1140, 3660, 3864, 5016, 11040, 15210, 16104, 16536, 18480, 44772, 57960, 67158, 68640, 68880, 142290, 142310, 155760, 196248, 198990, 240312, 248040, 275520, 278160, 308220, 322080, 326424, 339822, 348840, 352632, 366792, 462330, 485760, 607920
Offset: 1

Views

Author

Amiram Eldar, Dec 14 2018

Keywords

Examples

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

Crossrefs

Cf. A002025, A002952, A322485, A322486, A322542.

Programs

  • Mathematica
    f[p_, e_] := (p^Floor[(e + 1)/2] - 1)/(p - 1) + p^e; s[n_] := If[n == 1, 1, Times @@ (f @@@ FactorInteger[n])] - n; seq = {}; Do[n = s[m]; If[n > m && s[n] == m, AppendTo[seq, m]], {m, 1, 1000000}]; seq
  • PARI
    susigma(n) = {my(f = factor(n)); for (k=1, #f~, my(p=f[k, 1], e=f[k, 2]); f[k, 1] = (p^((e+1)\2) - 1)/(p-1) + p^e; f[k, 2] = 1; ); factorback(f); } \\ A322485
isok(n) = my(m=susigma(n)-n); (m > n) && (susigma(m) == n + m); \\ Michel Marcus, Dec 15 2018

A322542 Larger of semi-unitary amicable numbers pair: numbers (m, n) such that susigma(m) = susigma(n) = m + n, where susigma(n) is the sum of the semi-unitary divisors of n (A322485).

Original entry on oeis.org

126, 378, 1260, 3780, 4584, 5544, 11424, 15390, 16632, 16728, 25296, 49308, 68760, 73962, 88608, 84336, 179118, 168730, 172560, 225096, 256338, 266568, 250920, 297024, 287280, 365700, 374304, 391656, 374418, 387720, 386568, 393528, 548550, 502656, 623280
Offset: 1

Views

Author

Amiram Eldar, Dec 14 2018

Keywords

Comments

The terms are ordered according to the order of their lesser counterparts (A322541).

Examples

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

Crossrefs

Cf. A002025, A002952, A322485, A322486, A322541.

Programs

  • Mathematica
    f[p_, e_] := (p^Floor[(e + 1)/2] - 1)/(p - 1) + p^e; s[n_] := If[n == 1, 1, Times @@ (f @@@ FactorInteger[n])] - n; seq = {}; Do[n = s[m]; If[n > m && s[n] == m, AppendTo[seq, n]], {m, 1, 1000000}]; seq
  • PARI
    susigma(n) = {my(f = factor(n)); for (k=1, #f~, my(p=f[k, 1], e=f[k, 2]); f[k, 1] = (p^((e+1)\2) - 1)/(p-1) + p^e; f[k, 2] = 1; ); factorback(f); } \\ A322485
lista(nn) = {for (n=1, nn, my(m=susigma(n)-n); if ((m > n) && (susigma(m) == n + m), print1(m, ", ")););} \\ Michel Marcus, Dec 15 2018

A322486 Semi-unitary perfect numbers: numbers k such that susigma(k) = 2k, where susigma(k) is the sum of the semi-unitary divisors of k (A322485).

Original entry on oeis.org

6, 60, 90, 264, 3960, 4560, 8736, 13770, 131040, 384384, 605880, 5765760, 20049120, 882161280, 23253135360
Offset: 1

Views

Author

Amiram Eldar, Dec 11 2018

Keywords

Comments

a(16) <= 1846273228800. - David A. Corneth, Dec 11 2018

Examples

			264 is in the sequence since its sum of semi-unitary divisors is susigma(264) = 528 = 2 * 264.

Crossrefs

Cf. A000396, A002827, A322485.

Programs

  • Mathematica
    f[p_, e_] := (p^Floor[(e+1)/2] - 1)/(p-1) + p^e; susigma[n_] := If[n==1, 1, Times @@ (f @@@ FactorInteger[n])]; aQ[n_] := susigma[n]==2n; Select[Range[10000], aQ]
  • PARI
    ssu(n) = {my(f = factor(n)); for (k=1, #f~, my(p=f[k,1], e=f[k,2]); f[k,1] = (p^((e+1)\2) - 1)/(p-1) + p^e; f[k,2] = 1;); factorback(f);} \\ A322485
isok(n) = ssu(n) == 2*n; \\ Michel Marcus, Dec 14 2018

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.

Links

Crossrefs

Cf. A000203, A034448, A049417, A072691, A188999, A322485, A324707, A324708, A324709.

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
Showing 1-4 of 4 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.