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-10 of 11 results. Next

A349169 Numbers k such that k * gcd(sigma(k), A003961(k)) is equal to the odd part of {sigma(k) * gcd(k, A003961(k))}, where A003961 shifts the prime factorization one step towards larger primes, and sigma is the sum of divisors function.

Original entry on oeis.org

1, 15, 105, 3003, 3465, 13923, 45045, 264537, 459459, 745875, 1541475, 5221125, 8729721, 10790325, 14171625, 29288025, 34563375, 57034575, 71430975, 99201375, 109643625, 144729585, 205016175, 255835125, 295708875, 356080725, 399242025, 419159475, 449323875, 928602675, 939495375, 1083656925, 1941623775, 1962350685, 2083228875
Offset: 1

Views

Author

Antti Karttunen, Nov 10 2021

Keywords

Comments

Numbers k such that A348990(k) [= k/gcd(k, A003961(k))] is equal to A348992(k), which is the odd part of A349162(k), thus all terms must be odd, as A348990 preserves the parity of its argument.
Equally, numbers k for which gcd(A064987(k), A191002(k)) is equal to A000265(gcd(A064987(k), A341529(k))).
Also odd numbers k for which A348993(k) = A319627(k).
Odd terms of A336702 are given by the intersection of this sequence and A349174.
Conjectures:
(1) After 1, all terms are multiples of 3. (Why?)
(2) After 1, all terms are in A104210, in other words, for all n > 1, gcd(a(n), A003961(a(n))) > 1. Note that if we encountered a term k with gcd(k, A003961(k)) = 1, then we would have discovered an odd multiperfect number.
(3) Apart from 1, 15, 105, 3003, 13923, 264537, all other terms are abundant. [These apparently are also the only terms that are not Zumkeller, A083207. Note added Dec 05 2024]
(4) After 1, all terms are in A248150. (Cf. also A386430).
(5) After 1, all terms are in A348748.
(6) Apart from 1, there are no common terms with A349753.
Note: If any of the last four conjectures could be proved, it would refute the existence of odd perfect numbers at once. Note that it seems that gcd(sigma(k), A003961(k)) < k, for all k except these four: 1, 2, 20, 160.
Questions:
(1) For any term x here, can 2*x be in A349745? (Partial answer: at least x should be in A191218 and should not be a multiple of 3). Would this then imply that x is an odd perfect number? (Which could explain the points (1) and (4) in above, assuming the nonexistence of opn's).

Crossrefs

Programs

  • Mathematica
    Select[Range[10^6], #1/GCD[#1, #3] == #2/(2^IntegerExponent[#2, 2]*GCD[#2, #3]) & @@ {#, DivisorSigma[1, #], Times @@ Map[NextPrime[#1]^#2 & @@ # &, FactorInteger[#]]} &] (* Michael De Vlieger, Nov 11 2021 *)
  • PARI
    A000265(n) = (n >> valuation(n, 2));
    A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
    isA349169(n) = { my(s=sigma(n),u=A003961(n)); (n*gcd(s,u) == A000265(s)*gcd(n,u)); }; \\ (Program simplified Nov 30 2021)

Formula

For all n >= 1, A007949(A000203(a(n))) = A007949(a(n)). [sigma preserves the 3-adic valuation of the terms of this sequence] - Antti Karttunen, Nov 29 2021

Extensions

Name changed and comment section rewritten by Antti Karttunen, Nov 29 2021

A068390 Numbers k such that sigma(k) = 4*phi(k).

Original entry on oeis.org

14, 105, 248, 418, 1485, 3135, 3596, 3956, 4064, 5396, 8636, 20026, 23374, 25714, 35074, 35343, 39105, 41656, 55154, 56134, 56536, 71145, 74613, 87087, 124605, 150195, 175305, 192855, 263055, 393104, 413655, 421005, 474548, 604012, 697851, 711988, 819772
Offset: 1

Views

Author

Benoit Cloitre, Mar 03 2002

Keywords

Comments

If 2^p-1 is a prime (Mersenne prime) greater than 3 then 2^(p-2)*(2^p-1) is in the sequence. So for n>1, 2^(A000043(n)-2)*(2^A000043(n)-1) is in the sequence. - Farideh Firoozbakht, Feb 23 2005
Theorem: If m>0, k is an integer and p=2^(m+2)+k-1 is a prime number then n=2^m*p is a solution to the equation sigma(x) = 4*phi(x)-k. The previous comment is the special case k=0. - Farideh Firoozbakht, Oct 01 2014

References

  • D. S. Mitrinovic et al., Handbook of Number Theory, Kluwer, p. 88.

Crossrefs

Subsequence of A248150 (sigma(k) is divisible by 4).

Programs

  • Magma
    [n: n in [1..10^6] | SumOfDivisors(n) eq 4*EulerPhi(n)]; // Vincenzo Librandi, Sep 25 2017
  • Mathematica
    Select[Range[900000],DivisorSigma[1,#]==4EulerPhi[#]&] (* Harvey P. Dale, Nov 29 2013 *)
  • PARI
    for(n=1,300000, if(sigma(n)==4*eulerphi(n),print1(n,",")))
    

Extensions

More terms from Carl Najafi, Aug 16 2011

A349745 Numbers k for which k * gcd(sigma(k), A003961(k)) is equal to sigma(k) * gcd(k, A003961(k)), where A003961 shifts the prime factorization one step towards larger primes, and sigma is the sum of divisors function.

Original entry on oeis.org

1, 120, 216, 672, 2464, 22176, 228480, 523776, 640640, 837760, 5581440, 5765760, 7539840, 12999168, 19603584, 33860736, 38342304, 71344000, 95472000, 102136320, 197308800, 220093440, 345080736, 459818240, 807009280, 975576960, 1476304896, 1510831360, 1773584640
Offset: 1

Views

Author

Antti Karttunen, Nov 29 2021

Keywords

Comments

Numbers k for which k * A342671(k) = A000203(k) * A322361(k).
Numbers k such that gcd(A064987(k), A191002(k)) = gcd(A064987(k), A341529(k)).
Obviously, all odd terms in this sequence must be squares.
All the terms k of A005820 that satisfy A007949(k) < A007814(k) [i.e., whose 3-adic valuation is strictly less than their 2-adic valuation] are also terms of this sequence. Incidentally, the first six known terms of A005820 satisfy this condition, while on the other hand, any hypothetical odd 3-perfect number would be excluded from this sequence. Also, as a corollary, any hypothetical 3-perfect numbers of the form 4u+2 must not be multiples of 3 if they are to appear here. Similarly for any k which occurs in A349169, for 2*k to occur in this sequence, it shouldn't be a multiple of 3 and k should also be a term of A191218. See question 2 and its partial answer in A349169.
From Antti Karttunen, Feb 13-20 2022: (Start)
Question: Are all terms/2 (A351548) abundant, from n > 1 onward?
Note that of the 65 known 5-multiperfect numbers, all others except these three 1245087725796543283200, 1940351499647188992000, 4010059765937523916800 are also included in this sequence. The three exceptions are distinguished by the fact that their 3 and 5-adic valuations are equal. In 62 others the former is larger.
If k satisfying the condition were of the form 4u+2, then it should be one of the terms of A191218 doubled as only then both k and sigma(k) are of the form 4u+2, with equal 2-adic valuations for both. More precisely, one of the terms of A351538.
(End)

Crossrefs

Cf. also A349169, A349746, A351458, A351549 for other variants.
Subsequence of A351554 and also of its subsequence A351551.
Cf. A351459 (subsequence, intersection with A351458), A351548 (terms halved).

Programs

  • Mathematica
    f1[p_, e_] := (p^(e + 1) - 1)/(p - 1); f2[p_, e_] := NextPrime[p]^e; q[1] = True; q[n_] := n * GCD[(s = Times @@ f1 @@@ (f = FactorInteger[n])), (r = Times @@ f2 @@@ f)] == s*GCD[n, r]; Select[Range[10^6], q] (* Amiram Eldar, Nov 29 2021 *)
  • PARI
    A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
    isA349745(n) = { my(s=sigma(n),u=A003961(n)); (n*gcd(s,u) == (s*gcd(n,u))); };

Formula

For all n >= 1, A007814(A000203(a(n))) = A007814(a(n)). [sigma preserves the 2-adic valuation of the terms of this sequence]

A105824 a(n) = sigma(n) mod 4.

Original entry on oeis.org

1, 3, 0, 3, 2, 0, 0, 3, 1, 2, 0, 0, 2, 0, 0, 3, 2, 3, 0, 2, 0, 0, 0, 0, 3, 2, 0, 0, 2, 0, 0, 3, 0, 2, 0, 3, 2, 0, 0, 2, 2, 0, 0, 0, 2, 0, 0, 0, 1, 1, 0, 2, 2, 0, 0, 0, 0, 2, 0, 0, 2, 0, 0, 3, 0, 0, 0, 2, 0, 0, 0, 3, 2, 2, 0, 0, 0, 0, 0, 2, 1, 2, 0, 0, 0, 0, 0, 0, 2, 2, 0, 0, 0, 0, 0, 0, 2, 3, 0, 1, 2, 0, 0, 2, 0
Offset: 1

Views

Author

Shyam Sunder Gupta, May 05 2005

Keywords

Crossrefs

Sequences sigma(n) mod k: A053866 (k=2), A074941 (k=3), A105824 (k=4), A105825 (k=5), A084301 (k=6), A105826 (k=7), A105827 (k=8).

Programs

  • Maple
    A105824:= n-> (numtheory[sigma](n) mod 4):
    seq (A105824(n), n=1..105); # Jani Melik, Jan 26 2011
  • Mathematica
    Table[Mod[DivisorSigma[1, n], 4], {n, 100}] (* Wesley Ivan Hurt, Nov 07 2017 *)
  • PARI
    a(n)=sigma(n)%4

Formula

a(n) = A010873(A000203(n)). - Antti Karttunen, Nov 07 2017

A115414 Odd abundant numbers not divisible by 3.

Original entry on oeis.org

5391411025, 26957055125, 28816162375, 33426748355, 34393484125, 37739877175, 40342627325, 48150877775, 50866790975, 53356378075, 55959128225, 59305521275, 60711976325, 61164628525, 63395557225, 64899009175, 67275433225, 68972878975, 70088343325, 74922022175, 75665665075
Offset: 1

Views

Author

Sergio Pimentel, Mar 08 2006

Keywords

Comments

An odd abundant number (A005231) not divisible by 3 must have at least 7 distinct prime factors (e.g., 5^4*7^2*11^2*13*17*19*23) and be >= 5*29#/3# = 5^2*7*11*13*17*19*23*29 = 5391411025 = A047802(2) = a(1). This is most easily seen by writing the relative abundancy A(N) = sigma(N)/2N = sigma[-1](N) as A(Product p_i^e_i) = (1/2)*Product (p_i-1/p_i^e_i)/(p_i-1) < (1/2)*Product p_i/(p_i-1). See A064001 for odd abundant numbers not divisible by 5. - M. F. Hasler, Jul 27 2016
This is not a subsequence of A248150. For example, 81324229811825 and 37182145^2 = 1382511906801025 are terms, with sigma(.) == 2 (mod 4) and sigma(.) == 3 (mod 4) respectively. - Amiram Eldar, Aug 24 2020

Examples

			a(1)=5391411025 because it is the smallest abundant number (sigma(n)/n =~ 2.003) that is not divisible by 2 or 3.
		

Crossrefs

Programs

Extensions

Added missing term 55959128225 and a(14)-a(16) from Donovan Johnson, Dec 29 2008
a(17)-a(20) from Donovan Johnson, Dec 01 2011
More terms from M. F. Hasler, Jul 28 2016

A248151 Numbers n such that the sum of the divisors of n is not divisible by 4.

Original entry on oeis.org

1, 2, 4, 5, 8, 9, 10, 13, 16, 17, 18, 20, 25, 26, 29, 32, 34, 36, 37, 40, 41, 45, 49, 50, 52, 53, 58, 61, 64, 68, 72, 73, 74, 80, 81, 82, 89, 90, 97, 98, 100, 101, 104, 106, 109, 113, 116, 117, 121, 122, 128, 136, 137, 144, 146, 148, 149, 153, 157, 160, 162, 164, 169, 173, 178, 180
Offset: 1

Views

Author

M. F. Hasler, Oct 02 2014

Keywords

Comments

The complement of A248150; contains in particular A028982, the squares and twice the squares, for which sigma is odd.
The terms cannot have more than one odd prime factor to an odd power. Therefore this set has asymptotic density zero. The sequence grows faster than .75 n log(n). In particular: a(10) = 17, a(100) = 298, a(10^3) = 4724, a(10^4) = 66533, a(10^5) = 872434, a(10^6) = 10819205. - M. F. Hasler, Apr 26 2017

Crossrefs

Cf. A000203 (sum of divisors), A028982, A028983, A248150.
Subsequence of A285801.

Programs

  • Mathematica
    Select[Range[200],Mod[DivisorSigma[1,#],4]!=0&] (* Harvey P. Dale, Apr 03 2025 *)
  • PARI
    for(n=1,999,sigma(n)%4&&print1(n","))

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

A285800 Numbers having more than one odd prime factor to an odd power in their prime factorization.

Original entry on oeis.org

15, 21, 30, 33, 35, 39, 42, 51, 55, 57, 60, 65, 66, 69, 70, 77, 78, 84, 85, 87, 91, 93, 95, 102, 105, 110, 111, 114, 115, 119, 120, 123, 129, 130, 132, 133, 135, 138, 140, 141, 143, 145, 154, 155, 156, 159, 161, 165, 168, 170, 174, 177, 182, 183, 185, 186
Offset: 1

Views

Author

M. F. Hasler, Apr 26 2017

Keywords

Comments

The sequence is of asymptotic density one, a(n) ~ n.

Examples

			15 = 3*5, 21 = 3*7, 30 = 2*15, 33 = 3*11 are the smallest positive integers having at least two prime factors to an odd power in their factorization.
a(10) = 57, a(100) = 287, a(10^3) = 1950, a(10^4) = 15701, a(10^5) = 138540, a(10^6) = 1284998.
		

Crossrefs

Subsequence of A248150.
Complement of A285801.

Programs

  • Maple
    s800:=[]; s801:=[1];
    for n from 2 to 1000 do
    c:=0;
    t2:=ifactors(n)[2];
    for t3 in t2 do if t3[1]>2 and (t3[2] mod 2 = 1) then c:=c+1; fi; od:
    if c <= 1 then s801:=[op(s801),n]; else s800:=[op(s800),n]; fi;
    od:
    s800; # A285800
    s801; # A285801 - N. J. A. Sloane, Sep 30 2017
  • PARI
    is(n)=1<#select(t->bittest(t,0),factor(n>>valuation(n,2))[,2])

A386430 Odd numbers k such that there are no prime factors p of sigma(k) such that p does not divide A003961(k) and the valuation(k, p) is different from valuation(sigma(k), p), where A003961 is fully multiplicative with a(p) = nextprime(p), and sigma is the sum of divisors function.

Original entry on oeis.org

1, 3, 7, 15, 21, 27, 31, 33, 57, 69, 87, 91, 93, 105, 127, 141, 177, 189, 195, 217, 231, 237, 273, 285, 301, 381, 399, 447, 465, 483, 495, 513, 567, 573, 597, 609, 627, 651, 717, 775, 819, 837, 861, 889, 903, 987, 1023, 1029, 1149, 1185, 1239, 1311, 1365, 1419, 1431, 1437, 1455, 1497, 1561, 1653, 1659, 1687, 1743
Offset: 1

Views

Author

Antti Karttunen, Aug 22 2025

Keywords

Comments

Conjecture: After the initial 1, and apart from any hypothetical odd perfect numbers, all other terms are in A248150, i.e., sigma(k) == 0 (mod 4). This would imply (with the same caveat), that this sequence has no common terms with A228058 and no squares larger than one. This is true at least for the first 709203 terms (terms in range [1..2^34]).
Terms k such that A162642(k) = 1 are rare: 3, 7, 27, 31, 127, 567, 775, 8191, 27783, 131071, 524287, 2147483647, ... (odd terms of A387160).

Examples

			a(386548) = 5919068925 = 3^4 * 5^2 * 7^2 * 11^2 * 17 * 29. sigma(5919068925) = 15355618740 = 2^2 * 3^4 * 5 * 7 * 11^2 * 19^2 * 31. The "don't care primes" is given by A003961(A007947(5919068925))) = 2947945 = 5*7*11*13*19*31, thus only odd prime factor that matters here is 3, which in case has the same exponent (4) in both n = 5919068925 and sigma(n). In a way, this number is very close to satisfying Euler's criterion for odd perfect numbers (A228058), except that it has two unitary prime factors of the form 4k+1, instead of just one, apart from the square factor. Both n/17 and n/29 are in A228058.
		

Crossrefs

Odd terms of A351554.
Cf. A349169 (subsequence).

Programs

  • PARI
    A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
    A351555(n) = { my(s=sigma(n),f=factor(s),u=A003961(n)); sum(k=1,#f~,if((f[k,1]%2) && 0!=(u%f[k,1]), (valuation(n,f[k,1])!=f[k,2]), 0)); };
    isA386430(n) = ((n%2) && (0==A351555(n)));

Formula

{k | k odd, A351555(k) = 0}.

A292762 Numbers of the form p^k or 2*p^k, where p is a prime == 3 mod 4 and k is odd.

Original entry on oeis.org

3, 6, 7, 11, 14, 19, 22, 23, 27, 31, 38, 43, 46, 47, 54, 59, 62, 67, 71, 79, 83, 86, 94, 103, 107, 118, 127, 131, 134, 139, 142, 151, 158, 163, 166, 167, 179, 191, 199, 206, 211, 214, 223, 227, 239, 243, 251, 254, 262, 263, 271, 278, 283, 302, 307, 311, 326, 331, 334, 343, 347, 358, 359, 367, 379, 382
Offset: 1

Views

Author

N. J. A. Sloane, Sep 26 2017

Keywords

Comments

Numbers m such that sigma(m) == 0 mod 4 and phi(m) == 2 mod 4.

Crossrefs

Intersection of A097987 and A248150.

Programs

  • Mathematica
    Do[If[Mod[DivisorSigma[1,n],4]==0 && Mod[EulerPhi[n],4]==2,Print[n]],{n,1,10^3}] (* Vincenzo Librandi, Oct 02 2017 *)
  • PARI
    isok(m) = ((sigma(m) % 4) == 0) && ((eulerphi(m) % 4) == 2); \\ Michel Marcus, Oct 02 2017
    
  • PARI
    upto(n) = {my(l=List()); forprime(p=3, n, if(p%4==3, forstep(e=1, logint(n,p), 2, listput(l,p^e); if(2*p^e <= n, listput(l,2*p^e))))); listsort(l); l} \\ David A. Corneth, Oct 02 2017

Formula

As 22 = 2 * 11^1, which is of the form 2 * p^k with p = 11 = 2 * 4 + 3 == 3 mod 4 and k = 1 which is odd, 22 is a term. - David A. Corneth, Oct 02 2017
Showing 1-10 of 11 results. Next