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.

A241646 Numbers m such that the GCD of the x's that satisfy sigma(x) = m is 1.

Original entry on oeis.org

1, 12, 18, 24, 31, 32, 42, 48, 54, 56, 60, 72, 80, 84, 90, 96, 98, 104, 108, 114, 120, 128, 132, 140, 144, 152, 156, 168, 180, 182, 192, 216, 224, 228, 234, 240, 248, 252, 264, 270, 272, 280, 288, 294, 308, 312, 324, 336, 342, 360, 372, 384, 390, 408, 420
Offset: 1

Views

Author

Michel Marcus, Apr 26 2014

Keywords

Examples

			We have sigma(6) = sigma(11) = 12, and gcd(6, 11) = 1, hence 12 is in the sequence.
For x in [20, 26, 41], sigma(x) = 42, and gcd(20, 26, 41) = 1, hence 42 is here.
		

Crossrefs

Programs

  • Maple
    N:= 10^4: # for terms <= N
    V:= Vector(N):
    for x from 1 to N do
      s:= numtheory:-sigma(x);
      if s <= N then
        if V[s] = 0 then V[s]:= x
        else V[s]:= igcd(V[s], x)
        fi
      fi
    od: select(t -> V[t]=1, [$1..N]);  # Robert Israel, Aug 18 2019
  • PARI
    is(k) = gcd(invsigma(k)) == 1; \\ Amiram Eldar, Dec 19 2024, using Max Alekseyev's invphi.gp

A241647 Numbers m such that the GCD of the x's that satisfy sigma(x) = m is 2.

Original entry on oeis.org

3, 126, 186, 399, 924, 1350, 1386, 1530, 1806, 2106, 2646, 2652, 2814, 2916, 3066, 3150, 3654, 3870, 4662, 4914, 6162, 6426, 6846, 6882, 6930, 7098, 7566, 7620, 8190, 8910, 9270, 10842, 11076, 12222, 12870, 14586, 14910, 15210, 15246, 15930, 16506, 17010
Offset: 1

Views

Author

Michel Marcus, Apr 26 2014

Keywords

Examples

			We have sigma(68) = sigma(82) = 126, and gcd(68, 82) = 2, hence 126 is in the sequence.
On the other hand, for x in [20, 26, 41], sigma(x) = 42, and gcd(20, 26, 41) = 1, hence 42 is not here, although gcd(20, 26) is 2.
		

Crossrefs

Programs

  • Maple
    N:= 10^5: # for terms <= N
    V:= Vector(N):
    for x from 1 to N do
      s:= numtheory:-sigma(x);
      if s <= N then
        if V[s] = 0 then V[s]:= x
        else V[s]:= igcd(V[s], x)
        fi
      fi
    od: select(t -> V[t]=2, [$1..N]); # Robert Israel, Aug 18 2019
  • PARI
    is(k) = gcd(invsigma(k)) == 2; \\ Amiram Eldar, Dec 19 2024, using Max Alekseyev's invphi.gp

A241648 Numbers m such that the GCD of the x's that satisfy sigma(x) = m is 3.

Original entry on oeis.org

4, 124, 320, 392, 416, 800, 1352, 1520, 2912, 2960, 3536, 3872, 5720, 5936, 6320, 7112, 8216, 9176, 9912, 10472, 11816, 12152, 12896, 13280, 14960, 15176, 16080, 16400, 16536, 18032, 18392, 18560, 19136, 19880, 20000, 21632, 21680, 21920, 22736, 23120, 23816
Offset: 1

Views

Author

Michel Marcus, Apr 26 2014

Keywords

Examples

			We have sigma(48) = sigma(75) = 124, and gcd(48, 75) = 3, hence 124 is in the sequence.
Likewise, we have sigma(x) = 2912 for x = [1116, 1236, 1701, 2007, 2181], with gcd 3.
		

Crossrefs

Programs

  • Maple
    N:= 10^5: # for terms <= N
    V:= Vector(N):
    for x from 1 to N do
      s:= numtheory:-sigma(x);
      if s <= N then
        if V[s] = 0 then V[s]:= x
        else V[s]:= igcd(V[s], x)
        fi
      fi
    od: select(t -> V[t]=3, [$1..N]); # Robert Israel, Aug 18 2019
  • PARI
    is(k) = gcd(invsigma(k)) == 3; \\ Amiram Eldar, Dec 19 2024, using Max Alekseyev's invphi.gp

A241649 Numbers m such that the GCD of the x's that satisfy sigma(x) = m is 4.

Original entry on oeis.org

7, 210, 378, 630, 1904, 3570, 6188, 6510, 7154, 9296, 9800, 10220, 12446, 13664, 14378, 17654, 17780, 18536, 19110, 19376, 19530, 20034, 20580, 21266, 23240, 23310, 24150, 24584, 25298, 26754, 27930, 28938, 29106, 29610, 30380, 31640, 34146, 34230, 34664
Offset: 1

Views

Author

Michel Marcus, Apr 26 2014

Keywords

Examples

			sigma(104) = sigma(116) = 210, and gcd(104, 116) = 4, hence 210 is in the sequence.
Likewise 6510 is obtained with sigma of [2600, 2900, 3464, 3716], with gcd 4.
		

Crossrefs

Programs

  • Maple
    N:= 10^5: # for terms <= N
    V:= Vector(N):
    for x from 1 to N do
      s:= numtheory:-sigma(x);
      if s <= N then
        if V[s] = 0 then V[s]:= x
        else V[s]:= igcd(V[s],x)
        fi
      fi
    od:
    select(t -> V[t]=4, [$1..N]); # Robert Israel, Aug 18 2019
  • PARI
    is(k) = gcd(invsigma(k)) == 4; \\ Amiram Eldar, Dec 19 2024, using Max Alekseyev's invphi.gp
Showing 1-4 of 4 results.