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 21 results. Next

A057732 Numbers k such that 2^k + 3 is prime.

Original entry on oeis.org

1, 2, 3, 4, 6, 7, 12, 15, 16, 18, 28, 30, 55, 67, 84, 228, 390, 784, 1110, 1704, 2008, 2139, 2191, 2367, 2370, 4002, 4060, 4062, 4552, 5547, 8739, 17187, 17220, 17934, 20724, 22732, 25927, 31854, 33028, 35754, 38244, 39796, 40347, 55456, 58312, 122550, 205962, 235326, 363120, 479844, 685578, 742452, 1213815, 1434400, 1594947, 1875552, 1940812, 2205444
Offset: 1

Views

Author

G. L. Honaker, Jr., Oct 29 2000

Keywords

Comments

Some of the larger entries may only correspond to probable primes.
A number k is in this sequence iff A062709(k) is in A057733; this is the case iff A257273(k) is in A125246. - M. F. Hasler, Apr 27 2015

Examples

			For k = 6, 2^6 + 3 = 67 is prime.
For k = 28, 2^28 + 3 = 268435459 is prime.
		

References

  • Mike Oakes, posting to primenumbers(AT)yahoogroups.com on Jul 08 2001

Crossrefs

Cf. A019434 (primes 2^k+1), this sequence (2^k+3), A059242 (2^k+5), A057195 (2^k+7), A057196(2^k+9), A102633 (2^k+11), A102634 (2^k+13), A057197 (2^k+15), A057200 (2^k+17), A057221 (2^k+19), A057201 (2^k+21), A057203 (2^k+23).

Programs

  • Magma
    [n: n in [0..1000] | IsPrime(2^n+3)]; // Vincenzo Librandi, Apr 27 2015
  • Mathematica
    Select[Range[10000], PrimeQ[2^# + 3] &] (* Vincenzo Librandi, Apr 27 2015 *)
  • PARI
    for(n=1, 2200, if(isprime(2^n+3), print1(n, ", ")));
    
  • PARI
    for (n=1, 2, if (isprime(2^n+3), print1(n, ", "))); for(n=3, 100000, N=2^n+3 ; S=(N-5)/2 ; x=S ; for(j=1, n-1, x=Mod(x^2-2, N)) ; if(x==S , print1(n, ", "))) \\ produces terms corresponding to probable primes, see formula; Tony Reix, Aug 27 2015
    

Formula

Here is an LLT-like algorithm, using a cycle of the digraph x^2-2 modulo N, that finds terms of this sequence generating a PRP (PRobable Prime) of A057733 numbers: N=2^k+3; S0=(N-5)/2; s(0)=S0; s(i+1)=s(i)^2-2 modulo N; if s(k-1) == S0 then N is prime. - Tony Reix, Aug 27 2015

Extensions

More terms from Jason Earls, Jul 18 2001 and Mike Oakes, Jul 28 2001
a(47)-a(50) from Donovan Johnson 2006, verified by Paul Bourdelais, Mar 22 2012
a(51) is a probable prime based on trial factoring to 1E9 and PRP testing base 3,5,7 (PFGW v3.3.1). Discovered by Paul Bourdelais, Apr 09 2012
a(52)-a(54) from Paul Bourdelais, Jun 18 2019
a(55) from Paul Bourdelais, Jul 16 2019
a(56) from Paul Bourdelais, Apr 22 2020
a(57) from Paul Bourdelais, Jun 12 2020
a(58) from Paul Bourdelais, Aug 04 2020

A077374 Odd numbers m whose abundance by absolute value is at most 10, that is, -10 <= sigma(m) - 2m <= 10.

Original entry on oeis.org

1, 3, 5, 7, 9, 11, 15, 21, 315, 1155, 8925, 32445, 442365, 815634435
Offset: 1

Views

Author

Jason Earls, Nov 30 2002

Keywords

Comments

Apart from {1, 3, 5, 7, 9, 11, 15, 21, 315}, subset of A088012. Probably finite. - Charles R Greathouse IV, Mar 28 2011
a(15) > 10^13. - Giovanni Resta, Mar 29 2013
The abundance of the given terms a(1..14) is: (-1, -2, -4, -6, -5, -10, -6, -10, -6, -6, 6, 6, 6, -6). See also A171929, A188263 and A188597 for numbers with abundancy sigma(n)/n close to 2. - M. F. Hasler, Feb 21 2017
a(15) > 10^22. - Wenjie Fang, Jul 13 2017

Examples

			sigma(32445) = 64896 and 32445*2 = 64890, which makes the odd number 32445 six away from perfection: A(32445) = 6 and hence in this sequence.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[1, 10^6, 2], -10 <= DivisorSigma[1, #] - 2 # <= 10 &] (* Michael De Vlieger, Feb 22 2017 *)
  • PARI
    forstep(n=1,442365,2,if(abs(sigma(n)-2*n)<=10,print1(n,",")))

Extensions

a(14) from Farideh Firoozbakht, Jan 12 2004

A141548 Numbers n whose deficiency is 6.

Original entry on oeis.org

7, 15, 52, 315, 592, 1155, 2102272, 815634435
Offset: 1

Views

Author

Keywords

Comments

a(9) > 10^12. - Donovan Johnson, Dec 08 2011
a(9) > 10^13. - Giovanni Resta, Mar 29 2013
a(9) > 10^18. - Hiroaki Yamanouchi, Aug 21 2018
For all k in A059242, the number m = 2^(k-1)*(2^k+5) is in this sequence. This yields further terms 2^46*(2^47+5), 2^52*(2^53+5), 2^140*(2^141+5), ... All even terms known so far and the initial 7 = 2^0*(2^1+5) are of this form. All odd terms beyond a(2) are of the form a(n) = a(k)*p*q, k < n. We have proved that there is no further term of this form with the a(k) given so far. - M. F. Hasler, Apr 23 2015
A term n of this sequence multiplied by a prime p not dividing it is abundant if and only if p < sigma(n)/6 = n/3-1. For the even terms 592 and 2102272, there is such a prime near this limit (191 resp. 693571) such that n*p is a primitive weird number, cf. A002975. For a(3)=52, the largest such prime, 11, is already too small. Odd weird numbers do not exist within these limits. - M. F. Hasler, Jul 19 2016
Any term x of this sequence can be combined with any term y of A087167 to satisfy the property (sigma(x)+sigma(y))/(x+y) = 2, which is a necessary (but not sufficient) condition for two numbers to be amicable. - Timothy L. Tiffin, Sep 13 2016

Examples

			a(1) = 7, since 2*7 - sigma(7) = 14 - 8 = 6. - _Timothy L. Tiffin_, Sep 13 2016
		

Crossrefs

Cf. A087485 (odd terms).
Cf. A000203, A033880, A005100; A191363 (deficiency 2), A125246 (deficiency 4), A141548 (deficiency 6), A125247 (deficiency 8), A101223 (deficiency 10), A141549 (deficiency 12), A141550 (deficiency 14), A125248 (deficiency 16), A223608 (deficiency 18), A223607 (deficiency 20).
Cf. A087167 (abundance 6).

Programs

  • Magma
    [n: n in [1..9*10^6] | (SumOfDivisors(n)-2*n) eq -6]; // Vincenzo Librandi, Sep 14 2016
  • Mathematica
    lst={};Do[If[n==Plus@@Divisors[n]-n+6,AppendTo[lst,n]],{n,10^4}];Print[lst];
    Select[Range[1, 10^8], DivisorSigma[1, #] - 2 # == - 6 &] (* Vincenzo Librandi, Sep 14 2016 *)
  • PARI
    is(n)=sigma(n)==2*n-6 \\ Charles R Greathouse IV, Apr 23 2015, corrected by M. F. Hasler, Jul 18 2016
    

Extensions

a(8) from Donovan Johnson, Dec 08 2011

A191363 Numbers m such that sigma(m) = 2*m - 2.

Original entry on oeis.org

3, 10, 136, 32896, 2147516416
Offset: 1

Views

Author

Luis H. Gallardo, May 31 2011

Keywords

Comments

Let k be a nonnegative integer such that F(k) = 2^(2^k) + 1 is prime (a Fermat prime A019434), then m = (F(k)-1)*F(k)/2 appears in the sequence.
Conjecture: a(1)=3 is the only odd term of the sequence.
Conjecture: All terms of the sequence are of the above form derived from Fermat primes.
The sequence has 5 (known) terms in common with sequences A055708 (k-1 | sigma(k)) and A056006 (k | sigma(k)+2) since {a(n)} is a subsequence of both.
The first five terms of the sequence are respectively congruent to 3, 4, 4, 4, 4 modulo 6.
After a(5) there are no further terms < 8*10^9.
Up to m = 1312*10^8 there are no further terms in the class congruent to 4 modulo 6.
a(6) > 10^12. - Donovan Johnson, Dec 08 2011
a(6) > 10^13. - Giovanni Resta, Mar 29 2013
a(6) > 10^18. - Hiroaki Yamanouchi, Aug 21 2018
See A125246 for numbers with deficiency 4, i.e., sigma(m) = 2*m - 4, and A141548 for numbers with deficiency 6. - M. F. Hasler, Jun 29 2016 and Jul 17 2016
A term m of this sequence multiplied by a prime p not dividing it is abundant if and only if p < m-1. For each of a(2..5) there is such a prime near this limit (here: 7, 127, 30197, 2147483647) such that a(k)*p is a primitive weird number, cf. A002975. - M. F. Hasler, Jul 19 2016
Any term m of this sequence can be combined with any term j of A088831 to satisfy the property (sigma(m) + sigma(j))/(m+j) = 2, which is a necessary (but not sufficient) condition for two numbers to be amicable. [Proof: If m = a(n) and j = A088831(k), then sigma(m) = 2m-2 and sigma(j) = 2j+2. Thus, sigma(m) + sigma(j) = (2m-2) + (2j+2) = 2m + 2j = 2(m+j), which implies that (sigma(m) + sigma(j))/(m+j) = 2(m+j)/(m+j) = 2.] - Timothy L. Tiffin, Sep 13 2016
At least the first five terms are a subsequence of A295296 and of A295298. - David A. Corneth, Antti Karttunen, Nov 26 2017
Conjectures: all terms are second hexagonal numbers (A014105). There are no terms with middle divisors. - Omar E. Pol, Oct 31 2018
The symmetric representation of sigma(m) of each of the 5 numbers in the sequence consists of 2 parts of width 1 that meet at the diagonal (subsequence of A246955). - Hartmut F. W. Hoft, Mar 04 2022
The first five terms coincide with the sum of two successive terms of A058891. The same is not true for a(6), if such exists. - Omar E. Pol, Mar 03 2023

Examples

			For n=1, a(1) = 3 since sigma(3) = 4 = 2*3 - 2.
		

Crossrefs

Cf. A000203, A002975, A056006, A055708, A088831 (abundance 2).
Cf. A033880, A125246 (deficiency 4), A141548 (deficiency 6), A125247 (deficiency 8), A125248 (deficiency 16).
Cf. A058891.

Programs

  • Magma
    [n: n in [1..9*10^6] | (SumOfDivisors(n)-2*n) eq -2]; // Vincenzo Librandi, Sep 15 2016
  • Mathematica
    ok[n_] := DivisorSigma[1,n] == 2*n-2; Select[ Table[ 2^(2^k-1) * (2^(2^k)+1), {k, 0, 5}], ok] (* Jean-François Alcover, Sep 14 2011, after conjecture *)
    Select[Range[10^6], DivisorSigma[1, #] == 2 # - 2 &] (* Michael De Vlieger, Sep 14 2016 *)
  • PARI
    zp(a,b) = {my(c,c1,s); c = a; c1 = 2*c-2;
    while(c
    				
  • PARI
    a(k)=(2^2^k+1)<<(2^k-1) \\ For k<6. - M. F. Hasler, Jul 27 2016
    

Formula

a(n) = (A019434(n)-1)*A019434(n)/2 for all terms known so far. - M. F. Hasler, Jun 29 2016

A088832 Numbers k whose abundance is 4: sigma(k) - 2*k = 4.

Original entry on oeis.org

12, 70, 88, 1888, 4030, 5830, 32128, 521728, 1848964, 8378368, 34359083008, 66072609790, 549753192448, 259708613909470, 2251799645913088, 9223372026117357568
Offset: 1

Views

Author

Labos Elemer, Oct 28 2003

Keywords

Comments

If 2^m-5 is prime (A059608) then n=2^(m-1)*(2^m-5) is in the sequence. 12, 88, 1888, 32128, 521728, 8378368 & 34359083008 are such terms. See comments in A088831. - Farideh Firoozbakht, Feb 15 2008
Any term x of this sequence can be combined with any term y of A125246 to satisfy the property (sigma(x)+sigma(y))/(x+y) = 2, which is a necessary (but not sufficient) condition for two numbers to be amicable. - Timothy L. Tiffin, Sep 13 2016
Also contains 865268370658615254581248 = 2^23 * 16823249 * 6131278669. - Max Alekseyev, May 29 2025

Examples

			Abundances of terms in A045769: {-5,4,4,4,4,4,4,4,4,4} so A045769(1)=9 is not here.
		

Crossrefs

Subsequence of A045769.
Cf. A033880, A045768, A088830, A059608, A125246 (deficiency 4).

Programs

Formula

Solutions to sigma(x)-2*x=4.

Extensions

One more terms from Farideh Firoozbakht, Feb 15 2008
a(11)-a(12) from Donovan Johnson, Dec 23 2008
a(13) from Donovan Johnson, Dec 08 2011
a(14)-a(15) from Hiroaki Yamanouchi, Aug 23 2018
a(16) from Max Alekseyev, May 29 2025

A101223 Numbers m whose deficiency is 10, or: sigma(m) = 2m - 10.

Original entry on oeis.org

11, 21, 26, 68, 656, 2336, 8768, 133376, 528896, 34360918016, 35184409837568, 576460757135261696
Offset: 1

Views

Author

Vassil K. Tintschev (tinchev(AT)sunhe.jinr.ru), Dec 15 2004

Keywords

Comments

a(13) > 10^18. - Hiroaki Yamanouchi, Aug 21 2018
A subsequence of A274556. a(11) <= b(23) = 35184409837568 ~ 3.5*10^13, since b(k) := 2^(k-1)*(2^k+9) is in this sequence for all k in A057196 (2^k+9 is prime). All known terms except a(2) = 21 are of that form. - M. F. Hasler, Jul 18 2016
Any term x of this sequence can be combined with any term y of A223609 to satisfy the property (sigma(x)+sigma(y))/(x+y) = 2, which is a necessary (but not sufficient) condition for two numbers to be amicable. - Timothy L. Tiffin, Sep 13 2016

Examples

			The divisors of 68 are {1, 2, 4, 17, 34, 68} and so sigma(68) = 1 + 2 + 4 + 17+ 24 + 68 = 126 = 2*68 - 10; thus, the deficiency of 68 is 10 so 68 is a term of the sequence.
		

Crossrefs

Cf. A033879, A033880, A125246 (deficiency 4), A141548 (deficiency 6), A125247 (deficiency 8), A125248 (deficiency 16).
Cf. also A274556.
Cf. A223609 (abundance 10).

Programs

  • Magma
    [n: n in [1..9*10^6] | (SumOfDivisors(n)) eq 2*n-10]; // Vincenzo Librandi, Sep 15 2016
  • Mathematica
    Select[ Range[ 85000000], DivisorSigma[1, # ] + 10 == 2# &]

Extensions

Edited and extended by Robert G. Wilson v, Dec 15 2004
a(10) from Donovan Johnson, Dec 23 2008
Edited by M. F. Hasler, Jul 18 2016
a(11)-a(12) from Hiroaki Yamanouchi, Aug 21 2018

A125248 Numbers n whose abundance sigma(n)-2n = -16. Numbers n whose deficiency is 16.

Original entry on oeis.org

17, 38, 92, 170, 248, 752, 988, 2528, 8648, 12008, 34688, 63248, 117808, 526688, 531968, 820808, 1292768, 1495688, 2095208, 2112512, 3477608, 4495808, 8419328, 12026888, 13192768, 16102808, 26347688, 29322008, 33653888, 169371008
Offset: 1

Views

Author

Jason G. Wurtzel, Nov 25 2006

Keywords

Comments

When p=2^k+15 is prime (cf. A057197), then 2^(k-1)*p is in this sequence. The terms { 17, 38, 92, 248, 752, 2528, 34688, 531968, 2112512, 8419328, 537116672, 2147975168, ...} are of this from, with k in {1, 2, 3, 4, 5, 6, 8, 10, 11, 12, 15, 16, ...} = A057197. - M. F. Hasler, Jul 18 2016
Any term x of this sequence can be combined with any term y of A141547 to satisfy the property (sigma(x)+sigma(y))/(x+y) = 2, which is a necessary (but not sufficient) condition for two numbers to be amicable. - Timothy L. Tiffin, Sep 13 2016

Examples

			The abundance of 38 = (1+2+19+38)-76 = -16
		

Crossrefs

Cf. A000203, A033880, A005100; A191363 (deficiency 2), A125246 (deficiency 4), A141548 (deficiency 6), A125247 (deficiency 8), A101223 (deficiency 10), A141549 (deficiency 12), A141550 (deficiency 14), A125248 (this), A223608 (deficiency 18), A223607 (deficiency 20); A141547 (abundance 16).

Programs

  • Magma
    [n: n in [1..9*10^6] | (SumOfDivisors(n)-2*n) eq -16]; // Vincenzo Librandi, Sep 14 2016
  • Mathematica
    Select[Range[1, 10^6], DivisorSigma[1, #] - 2 # == - 16 &] (* Vincenzo Librandi, Sep 14 2016 *)
  • PARI
    for(n=1,1000000,if(((sigma(n)-2*n)==-16),print1(n,",")))
    

Extensions

a(17) to a(30) from Klaus Brockhaus, Nov 29 2006

A141549 Numbers k whose deficiency is 12: 2k - sigma(k) = 12.

Original entry on oeis.org

13, 45, 76, 688, 8896, 133888, 537051136, 35184418226176, 144115191028645888, 2305843021024854016
Offset: 1

Views

Author

Keywords

Comments

Numbers n whose abundance is -12. No other terms up to n=100,000,000. - Jason G. Wurtzel, Aug 24 2010
For all k in A102633, the number 2^(k-1)*(2^k+11) is in this sequence. So far all terms except a(2) are of this form. For k = 55, 71, this yields terms 649037107316853651724695645454336, 2787593149816327892704951291908936712585216. - M. F. Hasler, Apr 23 2015; edited by Max Alekseyev, May 27 2025
Any term x = a(m) can be combined with any term y = A141545(n) to satisfy the property (sigma(x)+sigma(y))/(x+y) = 2. Although this property is a necessary condition for two numbers to be amicable, it is not a sufficient one. So far, these two sequences have not produced an amicable pair. However, if one is ever found, then it will exhibit x-y = 12. - Timothy L. Tiffin, Sep 13 2016
a(11) > 10^20. - Max Alekseyev, May 27 2025

Examples

			a(1) = 13, since 2*13 - sigma(13) = 26 - 14 = 12. - _Timothy L. Tiffin_, Sep 13 2016
		

Crossrefs

Cf. A000203, A033880, A005100; A191363 (deficiency 2), A125246 (deficiency 4), A141548 (deficiency 6), A125247 (deficiency 8), A101223 (deficiency 10), A141550 (deficiency 14), A125248 (deficiency 16), A223608 (deficiency 18), A223607 (deficiency 20); A141545 (abundance 12).

Programs

  • Magma
    [n: n in [1..9*10^6] | (SumOfDivisors(n)-2*n) eq -12]; // Vincenzo Librandi, Sep 14 2016
  • Mathematica
    lst={};Do[If[n==Plus@@Divisors[n]-n+12,AppendTo[lst,n]],{n,10^4}];Print[lst];
    Select[Range[1, 10^8], DivisorSigma[1, #] - 2 # == - 12 &] (* Vincenzo Librandi, Sep 14 2016 *)
  • PARI
    for(n=1, 10^8, if(((sigma(n)-2*n)==-12), print1(n, ", "))) \\ Jason G. Wurtzel, Aug 24 2010
    

Extensions

a(7) from Donovan Johnson, Dec 08 2011
a(8)-a(9) from Hiroaki Yamanouchi, Aug 21 2018
a(10) from Max Alekseyev, May 27 2025

A141550 Numbers n whose deficiency is 14.

Original entry on oeis.org

27, 34, 232, 34432, 549762629632
Offset: 1

Views

Author

Keywords

Comments

a(6) > 10^12. - Donovan Johnson, Dec 08 2011
a(6) > 10^13. - Giovanni Resta, Mar 29 2013
a(6) > 10^18. - Hiroaki Yamanouchi, Aug 21 2018
a(6) <= b(38) = 37778931864743868104704 = 3.77789*10^22, since b(k) = 2^(k-1)*(2^k+13) is in this sequence for all k in A102634, i.e., 2^k+13 is prime. All known terms except a(1) = 27 are of this form: a(2..5) = b(k) with k = 2, 4, 8, 20, and k = 38 yields the next larger term of this form. - M. F. Hasler, Jul 18 2016
Any term x of this sequence can be combined with any term y of A141546 to satisfy the property (sigma(x)+sigma(y))/(x+y) = 2, which is a necessary (but not sufficient) condition for two numbers to be amicable. - Timothy L. Tiffin, Sep 13 2016

Examples

			a(1) = 27, since 2*27 - sigma(27) = 54 - 40 = 14. - _Timothy L. Tiffin_, Sep 13 2016
		

Crossrefs

Cf. A000203, A033880, A005100; A191363 (deficiency 2), A125246 (deficiency 4), A141548 (deficiency 6), A125247 (deficiency 8), A101223 (deficiency 10), A141549 (deficiency 12), A125248 (deficiency 16); A141546 (abundance 14).

Programs

  • Magma
    [n: n in [1..9*10^6] | (SumOfDivisors(n)-2*n) eq -14]; // Vincenzo Librandi, Sep 14 2016
  • Mathematica
    lst={};Do[If[n==Plus@@Divisors[n]-n+14,AppendTo[lst,n]],{n,10^4}];Print[lst];
    Select[Range[1, 10^8], DivisorSigma[1, #] - 2 # == - 14 &] (* Vincenzo Librandi, Sep 14 2016 *)

Extensions

a(5) from Donovan Johnson, Dec 08 2011

A257273 a(n) = 2^(n-1)*(2^n+3).

Original entry on oeis.org

2, 5, 14, 44, 152, 560, 2144, 8384, 33152, 131840, 525824, 2100224, 8394752, 33566720, 134242304, 536920064, 2147581952, 8590131200, 34360131584, 137439739904, 549757386752, 2199026401280, 8796099313664, 35184384671744, 140737513521152, 562950003752960, 2251799914348544, 9007199456067584
Offset: 0

Views

Author

M. F. Hasler, Apr 27 2015

Keywords

Comments

a(n) is in A125246 <=> n is in A057732 <=> A062709(n) is in A057733.
These are also the row sum of the triangle A146769: For n>=1, a(n-1) is the sum of row n of A146769.

Crossrefs

Programs

  • Magma
    [2^(n-1)*(2^n+3): n in [0..35]]; // Vincenzo Librandi, Apr 27 2015
    
  • Mathematica
    Table[2^(n - 1) (2^n + 3), {n, 0, 30}] (* Bruno Berselli, Apr 27 2015 *)
    CoefficientList[Series[(2 - 7 x)/((1 - 4 x) (1 - 2 x)), {x, 0, 30}], x] (* Vincenzo Librandi, Apr 27 2015 *)
    LinearRecurrence[{6,-8},{2,5},30] (* Harvey P. Dale, Dec 21 2024 *)
  • PARI
    a(n)=2^(n-1)*(2^n+3)
    
  • PARI
    Vec((2-7*x)/((1-4*x)*(1-2*x)) + O(x^100)) \\ Colin Barker, Apr 27 2015

Formula

G.f.: (2-7*x)/((1-4*x)*(1-2*x)). - Vincenzo Librandi, Apr 27 2015
a(n) = 6*a(n-1)-8*a(n-2). - Colin Barker, Apr 27 2015
Showing 1-10 of 21 results. Next