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

A098223 Integer quotients when sigma(sigma(x))/x is an integer.

Original entry on oeis.org

1, 2, 2, 3, 4, 2, 3, 7, 6, 8, 2, 6, 6, 9, 8, 6, 10, 10, 3, 8, 4, 6, 7, 8, 2, 9, 10, 8, 4, 10, 10, 7, 13, 8, 8, 8, 2, 6, 8, 14, 2, 9, 7, 8, 6, 9, 8, 13, 8, 15, 14, 6, 9, 9, 8, 10, 12, 14, 13, 8, 8, 11, 6, 14, 16, 12, 14, 12, 16, 15, 12, 18, 16, 11, 8, 22
Offset: 1

Views

Author

Labos Elemer, Oct 25 2004

Keywords

Comments

Below n=5x10^11, q=5 and 17 quotients do not appear; smallest numbers providing integer quotients = 1, 2, 3, 4,..., 16,... are as follows: 1, 2, 8, 15, ?, 42, 24, 60, 168, 480, 57669920, 2200380, 57120, 217278, 1058148, 7526400, ... - updated by Jud McCranie, Feb 08 2012
The above sequence is now A272930. - Franklin T. Adams-Watters, May 11 2016
See A019278 for the actual numbers x such that x | sigma(sigma(x)). - M. F. Hasler, Jul 03 2016

Crossrefs

Programs

  • Maple
    with(numtheory): A098223:=n->`if`(sigma(sigma(n)) mod n = 0, sigma(sigma(n))/n, NULL): seq(A098223(n), n=1..10^5); # Wesley Ivan Hurt, Oct 10 2014
  • Mathematica
    Select[DivisorSigma[1, DivisorSigma[1, #]]/# &@ Range[10^6], IntegerQ] (* Michael De Vlieger, May 11 2016 *)
  • PARI
    for(n=1,1e7, sigma(sigma(n))%n||print1(sigma(sigma(n))/n",")) \\ M. F. Hasler, Jul 03 2016

Formula

In order of appearance the sigma(sigma(A019278(n)))/A019278(n) quotients which are by definition integers.

A274397 Positive integers m such that sigma(m) is divisible by 5.

Original entry on oeis.org

8, 19, 24, 27, 29, 38, 40, 54, 56, 57, 58, 59, 72, 76, 79, 87, 88, 89, 95, 104, 108, 109, 114, 116, 118, 120, 128, 133, 135, 136, 139, 145, 149, 152, 158, 168, 171, 174, 177, 178, 179, 184, 189, 190, 199, 200, 203, 209, 216, 218, 228, 229, 232, 236, 237, 239, 247, 248, 261, 264, 266, 267, 269, 270, 278, 280, 285, 290, 295, 296, 297
Offset: 1

Views

Author

M. F. Hasler, Jul 02 2016

Keywords

Comments

See the subsequence A274685 of odd terms for a remark on frequent pairs of the form (30k-3, 30k-1).
If m is in the sequence and gcd(k,m)=1, then k*m is also in the sequence. One might call "primitive" those terms which are not of this form, i.e., not a "coprime" multiple of an earlier term. The primitive terms are the primes and powers of primes within the sequence, cf. below.
Integers m > 0 where an integer k exists such that A000203(m) = A008587(k). - Felix Fröhlich, Jul 02 2016
For any prime p <> 5 there is an exponent k in {1, 3, 4} (depending on whether p is in A030433, A003631 or A030430) such that p^k is in this sequence. Given these p^k, the sequence consists of all numbers of the form n*p^(q*(k+1)-1) where n is coprime to p and q >= 1. Otherwise said, all numbers m which have some prime factor p with multiplicity q*(k+1)-1, where k = k(p) in {1, 3, 4} as introduced before. - M. F. Hasler, Jul 10 2016

Examples

			Some values for a(2^k): We have a(2) = 19, a(4) = 27, a(8) = 54, a(16) = 87, a(32) = 145, a(64) = 270, a(128) = 488, a(256) = 919, a(512) = 1736, a(1024) = 3267, a(2048) = 6258, a(4096) = 12035, a(8192) = 23160, a(16384) = 44878, a(32768) = 87207, a(65536) = 169911, a(131072) = 332009, a(262144) = 650031, a(524288) = 1274569, a(1048576) = 2503510, a(2097152) = 4924370, a(4194304) = 9697475, a(8388608) = 19116191.
		

Crossrefs

Cf. A000203, A028983 (sigma even), A087943 (sigma = 3k), A248150 (sigma = 4k); A028982 (sigma is odd), A248151 (sigma is not divisible by 4); A272930 (sigma(sigma(k)) = nk).

Programs

  • Maple
    select(t -> numtheory:-sigma(t) mod 5 = 0, [$1..1000]); # Robert Israel, Jul 12 2016
  • Mathematica
    Select[Range[300], Divisible[DivisorSigma[1, #], 5]&] (* Jean-François Alcover, Apr 09 2019 *)
  • PARI
    is(n)=sigma(n)%5==0
    
  • PARI
    is(n)=for(i=1,#n=factor(n)~,n[1,i] != 5 && (n[2,i]+1) % [5,4,4,2][n[1,i]%5] == 0 && return(1))

Formula

lim_{n->oo} a(k)/k = 2 (conjectured; cf. Examples).

Extensions

Edited by M. F. Hasler, Jul 10 2016

A274685 Odd numbers n such that sigma(n) is divisible by 5.

Original entry on oeis.org

19, 27, 29, 57, 59, 79, 87, 89, 95, 109, 133, 135, 139, 145, 149, 171, 177, 179, 189, 199, 203, 209, 229, 237, 239, 247, 261, 267, 269, 285, 295, 297, 319, 323, 327, 343, 349, 351, 359, 377, 379, 389, 395, 399, 409, 413, 417, 419, 435, 437, 439, 445, 447, 449, 459, 475, 479, 493, 499
Offset: 1

Views

Author

M. F. Hasler, Jul 02 2016

Keywords

Comments

The subsequence of odd terms in A274397.
If n is in the sequence and gcd(n,m)=1 for some odd m, then n*m is also in the sequence. One might call "primitive" those terms which are not of this form, i.e., not a "coprime" multiple of an earlier term. The list of these primitive terms is (19, 27, 29, 59, 79, 89, 109, 139, 149, 179, 199, 229, 239, 269, 343, 349, 359, 379, 389, 409, 419, 439, 449, 479, 499, ...). The primitive terms are the primes and powers of primes within the sequence. If a prime power p^k (k >= 1) is in the sequence, then p^(m(k+1)-1) is in the sequence for any m >= 1, since 1+p+...+p^(m(k+1)-1) = (1+p+...+p^k)(1+p^(k+1)+...+p^((m-1)*(k+1))). For example, with the prime p=19 we also have all odd powers 19^3, 19^5, ..., and with 27 = 3^3, we also have 27^5, 27^9, ... in the sequence.
On the other hand, for any prime p <> 5 there is an exponent k in {1, 3, 4} such that p^k is in this sequence (and therewith all higher powers of the form given above).
One may notice that there are many pairs of the form (30k-3, 30k-1), e.g., 27,29; 57,59; 87,89; 177,179; 237,239; 295,299; ... Indeed, it is likely that 30k-1 is prime and in this case, if 10k-1 is also prime, then sigma(30k-3) = 40k is divisible by 5 and sigma(30k-1) = 30k is also divisible by 5.

Examples

			Some values of a(2^k): a(2) = 27, a(4) = 57, a(8) = 89, a(16) = 171, a(32) = 297, a(64) = 545, a(128) = 1029, a(256) = 1937, a(512) = 3625, a(1024) = 6939, a(2048) = 13257, a(4096) = 25483, a(8192) = 49319, a(16384) = 95695, a(32768) = 185991, a(65536) = 362725, a(131072) = 708887, a(262144) = 1388367, a(524288) = 2722639, a(1048576) = 5346681, a(2097152) = 10514679, a(4194304) = 20698531, a(8388608) = 40790203.
		

Crossrefs

Cf. A000203 (sigma), A028983 (sigma even), A087943 (sigma = 3k), A248150 (sigma = 4k); A028982 (sigma is odd), A248151 (sigma is not divisible by 4); A272930(sigma(sigma(k)) = nk).

Programs

  • Mathematica
    Select[Range[1, 500, 2], Divisible[DivisorSigma[1, #], 5] &] (* Michael De Vlieger, Jul 16 2016 *)
  • PARI
    is_A274685(n)=sigma(n)%5==0&&bittest(n,0)
    
  • PARI
    forstep(n=1,999,2,sigma(n)%5||print1(n","))

Formula

a(n) ~ 2n. - Charles R Greathouse IV, Jul 16 2016

A318242 a(n) is the least k such that A188999(A188999(k)) = n*k, where A188999 is the bi-unitary sigma function, or 0 if no such k exists.

Original entry on oeis.org

1, 2, 8, 15, 24, 42, 240, 648, 168, 480, 321408, 4320, 57120, 103680, 1827840, 23591520, 898128000, 374250240
Offset: 1

Views

Author

Michel Marcus, Aug 22 2018

Keywords

Comments

It is also known that a(20) = 11975040.
Then for higher indices n, we have:
a(19) <= 5235707393280;
a(21) <= 49110437376000;
a(22) <= 106780561395056640;
a(24) <= 1099525819392000;
a(25) <= 41252767395840;
a(26) <= 202768780032000.

Crossrefs

Cf. A272930 (analog for sigma), A318272 (analog for infinitary sigma).

A318272 a(n) is the least k such that A049417(A049417(k)) = n*k, where A049417 is the infinitary sigma function, or 0 if no such k exists.

Original entry on oeis.org

1, 2, 8, 15, 24, 60, 1321920, 17850, 285600, 763776
Offset: 1

Views

Author

Michel Marcus, Aug 23 2018

Keywords

Comments

It is also known that a(12) = 71400.
Then for higher indices n, we have:
a(11) <= 414230544000;
a(13) <= 2667897127526400;
a(14) <= 446464417259520;
a(15) <= 23613006957281280;
a(16) <= 22227004800;
a(17) <= 3134896756992000;
a(18) <= 15414516736819200;
a(20) <= 53129265653882880.
a(16) = 8420630400. - Giovanni Resta, Aug 25 2018

Crossrefs

Cf. A272930 (analog for sigma), A318242 (analog for bi-unitary sigma).

A327630 a(n) is the least integer k such that sigma(sigma(k)) >= n*k where sigma is A000203, the sum of divisors.

Original entry on oeis.org

1, 2, 6, 6, 24, 24, 24, 60, 120, 480, 540, 1560, 2520, 10920, 27720, 30240, 191520, 360360, 360360, 3243240, 5765760, 28828800, 109549440, 438197760, 766846080, 3834230400, 9081072000, 32974381440, 147516969600, 880887047040, 2802822422400
Offset: 1

Views

Author

Michel Marcus, Oct 07 2019

Keywords

Comments

From Daniel Suteu, Oct 10 2019: (Start)
a(32) <= 14814918518400,
a(33) <= 64464915715200,
a(34) <= 709114072867200,
a(35) <= 9881550651772800,
a(36) <= 76648784785372800,
a(37) <= 2376112328346556800. (End)

Crossrefs

Programs

  • PARI
    a(n) = my(k=1); while (sigma(sigma(k))/k < n, k++); k;

Extensions

a(23)-a(29) from Giovanni Resta, Oct 07 2019
a(30)-a(31) from Giovanni Resta, Feb 27 2020

A292413 a(n) is the smallest k such that sigma(sigma(k)) = sigma(sigma(k) - k) + n*k, or 0 if no such k exists.

Original entry on oeis.org

4, 3, 14, 54, 1496, 42520, 177263760, 11712677760
Offset: 1

Views

Author

Altug Alkan, Sep 16 2017

Keywords

Comments

For n <= 5, the values k <= 10^6 such that sigma(sigma(k)) = sigma(sigma(k) - k) + n*k are:
For n = 1; 4, 144, 16384.
For n = 2; 3, 7, 8, 31, 84, 127, 8191, 35743, 42048, 131071, 524287.
For n = 3; 14, 20, 372, 1081, 16246, 98292.
For n = 4; 54, 168, 308, 504, 994, 3870, 4994, 7800, 16488, 17360, 51995, 475664.
For n = 5; 1496, 1704, 14976, 21552, 379938, 854168.
a(10) = 16708685994, a(9) > 6*10^11. - Giovanni Resta, Sep 18 2017

Crossrefs

Programs

  • PARI
    a(n) = {my(k=2); while(sigma(sigma(k))-sigma(sigma(k)-k) != n*k, k++); k; }

Extensions

a(8) from Giovanni Resta, Sep 18 2017
Showing 1-7 of 7 results.