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

A330901 Numbers k such that k and k+2 have the same deficiency (A033879).

Original entry on oeis.org

2, 6497, 12317, 91610, 133787, 181427, 404471, 439097, 485237, 1410119, 2696807, 6220607, 6827369, 6954767, 9770027, 10302419, 10449347, 10887977, 11014007, 16745387, 18959111, 25883519, 27334469, 39508037, 40311149, 40551617, 42561437, 44592209, 47717471, 48912107
Offset: 1

Views

Author

Amiram Eldar, May 01 2020

Keywords

Comments

Are 2 and 91610 the only even terms?
Are there any abundant numbers (A005101) in this sequence?
Numbers k such that k and k+1 have the same deficiency are 1, 145215, and no more below 10^13 (they are a subset of A112645).
Up to a(2214) = 2001876242879 there are no further even terms nor abundant terms. - Giovanni Resta, May 01 2020

Examples

			2 is a term since 2 and 4 have the same deficiency: A033879(2) = 2*2 - sigma(2) = 4 - 3 = 1, and A033879(4) = 2*4 - sigma(4) = 8 - 7 = 1.
		

Crossrefs

Programs

  • Mathematica
    def[n_] := 2*n - DivisorSigma[1, n]; Select[Range[10^5], def[#] == def[# + 2] &]
    SequencePosition[Table[2n-DivisorSigma[1,n],{n,48920000}],{x_,,x}][[;;,1]] (* Harvey P. Dale, Apr 26 2025 *)
  • PARI
    j1=1;j2=1;for(k=3,50000000,j=k+k-sigma(k);if(j==j1,print1(k-2,", "));j1=j2;j2=j) \\ Hugo Pfoertner, May 01 2020

A112646 Special solutions to abs(sigma(x+1) - sigma(x)) = 2 where x + 1 = 3^y.

Original entry on oeis.org

8, 26, 2186, 1594322, 7509466514979724803946715958257546, 13915193059764305937984450503671774362956903094026
Offset: 1

Views

Author

Labos Elemer, Sep 28 2005

Keywords

Comments

The relevant computed exponents providing -1+3^z solutions for A112645 are 2, 3, 7, 13, 71, 103.
The next term, if it exists, is too large to include here. - N. J. A. Sloane, Jul 12 2009
Contribution from David Wasserman, Jul 06 2009: No more terms < 10^86.

Crossrefs

Programs

  • Mathematica
    ta={{0}};Do[s=Abs[DivisorSigma[1, 3^n]-DivisorSigma[1, -1+3^n]]; If[Equal[s, 2], ta=Append[ta, n];Print[ -1+3^n]], {n, 1, 100000000}]; ta=Delete[ta, 1]

A112647 a(n)=x is the smallest solution to abs(sigma(x+1)-sigma(x))=n, or 0 if no solution exists.

Original entry on oeis.org

14, 2, 1, 3, 6, 9, 5, 7, 62
Offset: 0

Views

Author

Labos Elemer, Sep 28 2005

Keywords

Comments

The known values plus conjectured zero values are: 14, 2, 1, 3, 6, 9, 5, 7, 62, 0, 13, 25, 22, 16, 12, 32, 11, 0, 104, 18, 837, 17, 19, 63, 46, 0, 28, 0, 116, 24, 58, 31, 2222, 0, 39, 242, 23, 0, 147, 0, 30, 675, 29, 35, 52, 0, 777, 0, 40, 0, 435, 0, 42, 36, 41, 0, 91, 0, 67, 0, 65, 99, 0, 195, 110, 80, 53, 48, 124, 0, 243, 0, 70, 97.
The first unknown value is a(9).
The zero values are based on a search up to 10000000.
While it is known that not all m values satisfy sigma(x) = m (see A007369), it is more difficult to determine those numbers which cannot be a difference of sigma(u)-sigma(w) for some u and w.
No solutions to abs(sigma(x+1)-sigma(x)) = n with x < 2*10^8 for n = 9, 17, 25, 27, 33, 37, 39, 45, 47, 49, 51, 55, 57, 59, 62, 69, 71. - Robert Israel, May 24 2016
Except for a(62) = 1159742042, all the terms a(n)>0 with n <= 100 are either smaller than 10^6 or greater than 2*10^12. - Giovanni Resta, Oct 29 2019

Examples

			n=5: least solution is 9 because sigma for 9 and 9+1=10 are 13 and 13+5=18.
		

Crossrefs

Programs

  • MATLAB
    N = 2*10^8; % to search sigma(n) for n <= N
    M = 100;    % to get a(1) to a(M)
    Sigma = ones(1,N);
    for n=2:N
      inds = [n:n:N];
      Sigma(inds) = Sigma(inds) + n;
    end
    DSigma = abs(Sigma(2:end) - Sigma(1:end-1));
    A = zeros(1,M);
    for v = 1:M
      r = find(DSigma == v,1,'first');
      if numel(r) > 0
        A(v) = r;
      end
    end
    A   % Robert Israel, May 24 2016
  • Mathematica
    f[x_] :=Abs[DivisorSigma[1, n+1] - DivisorSigma[1, n]]; t=Table[0, {258}]; Do[s=f[n]; If[s<258 && t[[s+1]]==0, t[[s+1]]=n], {n, 10^7}]; t (* edited by Giovanni Resta, Oct 29 2019 *)

Extensions

Entry revised by N. J. A. Sloane, May 25 2016
a(0) prepended by Giovanni Resta, Oct 29 2019
Showing 1-3 of 3 results.