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.

Previous Showing 51-60 of 72 results. Next

A262908 a(n) = largest k such that A049820(k + A262509(n)) <= A262509(n).

Original entry on oeis.org

53, 49, 69, 55, 53, 31, 47, 39, 25, 35, 31, 39, 37, 51, 33, 43, 33, 69, 65, 57, 43, 41, 57, 49, 33, 33, 43, 41, 37, 33, 37, 39, 35, 27, 41, 27, 43, 75, 177, 171, 173, 155, 45, 133, 107, 121, 111, 139, 78, 119, 123, 47, 65, 79, 77, 97, 81, 151, 149, 145, 111, 197, 375, 71, 59, 81, 259, 257
Offset: 1

Views

Author

Antti Karttunen, Oct 08 2015

Keywords

Comments

For all nonzero terms a(n), A263083(n) = a(n) + A262509(n) and A155043(A263083(n)) < A155043(A262509(n)) because at each A262509(n) the "distance to zero", A155043 obtains a unique value A262508(n), thus no A049820-iteration trajectory starting from any k larger than A262509(n) and using a greater or equal number of steps to reach zero may bypass A262509(n) [i.e., without going through A262509(n)], because then A262508(n) would not be unique anymore. See also comments in A262909.

Crossrefs

Programs

Formula

Other identities. For all n >= 1:
a(n) <= A262909(n).

A263090 Primes p for which A049820(x) = p has a solution.

Original entry on oeis.org

2, 3, 5, 11, 17, 23, 29, 31, 41, 47, 53, 59, 61, 71, 73, 83, 89, 97, 101, 107, 127, 137, 139, 149, 151, 157, 173, 179, 181, 191, 197, 199, 211, 223, 227, 233, 239, 263, 269, 277, 281, 283, 311, 317, 331, 337, 347, 349, 367, 373, 389, 409, 419, 421, 431, 433, 443, 449, 457, 461, 467, 521, 523, 541, 547, 557, 569, 577, 587, 593, 599, 601, 607, 613, 617, 619, 631, 641, 643
Offset: 1

Views

Author

Antti Karttunen, Oct 11 2015

Keywords

Comments

Primes p that there is at least one such k for which k - d(k) = p, where d(k) is the number of divisors of k (A000005).

Examples

			2 is present, as we have 6 - d(6) = 6 - 4 = 2.
3 is present, as we have 5 - d(5) = 3. The same holds for all lesser twin primes (A001359).
		

Crossrefs

Complement among primes: A263091.
Intersection of A000040 and A236562.
Cf. A001359 (a subsequence).
Cf. also A263094.

Programs

  • Mathematica
    lim = 10000; s = Select[Sort@ DeleteDuplicates@ Table[n - DivisorSigma[0, n], {n, lim}], PrimeQ]; Take[s, 79] (* Michael De Vlieger, Oct 13 2015 *)
  • PARI
    allocatemem(123456789);
    uplim1 = 2162160 + 320; \\ = A002182(41) + A002183(41).
    v060990 = vector(uplim1);
    for(n=3, uplim1, v060990[n-numdiv(n)]++);
    A060990 = n -> if(!n,2,v060990[n]);
    n=0; forprime(p=2, 131071, if((A060990(p) > 0), n++; write("b263090.txt", n, " ", p)));
    
  • Scheme
    ;; With Antti Karttunen's IntSeq-library.
    (define A263090 (MATCHING-POS 1 1 (lambda (n) (and (= 1 (A010051 n)) (not (zero? (A060990 n)))))))

A345147 a(1)=1. If a(n) is a novel term, a(n+1) = A000005(a(n)). If a(n) is a repeat term, a(n+1) = a(n)+m, where m is the least prior term (which once used in this way is not used again).

Original entry on oeis.org

1, 1, 2, 2, 3, 2, 4, 3, 5, 2, 4, 6, 4, 7, 2, 5, 7, 11, 2, 6, 8, 4, 8, 12, 6, 11, 16, 5, 11, 16, 22, 4, 10, 4, 8, 12, 19, 2, 9, 3, 5, 8, 13, 2, 10, 12, 20, 6, 14, 4, 10, 14, 22, 31, 2, 12, 14, 24, 8, 18, 6, 14, 20, 31, 42, 8, 19, 27, 4, 16, 20, 32, 6, 18, 24, 36
Offset: 1

Views

Author

Keywords

Comments

Let conditions 0 and 1 respectively pertain to novel a(n) and extant a(n) respectively. For n>3 condition 0 results in a fall from a(n) to a(n+1), whereas condition 1 results in an increase, which may be sustained over several consecutive terms as the least unused term is applied to the result of a sum which is the same as a prior term, therefore prompting a further application of condition 1, until a novel term is reached.
A distinguishing feature of the scatterplot of this sequence is that the data are highly concentrated around several quasi-radial, apparently straight line "spokes" emanating from the origin, which appear to be equally spaced. The lower spoke relates to applications of condition 0, producing the values d(k) of novel terms k, which being small by comparison to k, are distributed widely across the lower zone of the plot. Higher spokes relate to consecutive terms generated by condition 1, namely repeated addition of the least unused term m.
Since m appears in multiples and increments by small amounts, successive additions result in roughly equally spaced consecutive terms. Repeated application of condition 1 continues until a novel term k is reached, followed by a fall to d(k), from which the process repeats, carrying forward any remaining copies of m until all are used, at which point m increments, causing the spokes to diverge gradually as the sequence advances.
Records beyond a(1) arise consequent only to repetitions of condition 1, and not to the appearance of a novel term. With expansion of the data the spokes are gradually obscured from the bottom up by horizontal striations resulting from the fact that many different numbers have the same number of divisors, and that with greater values of a(n), greater values of d(a(n)) can appear. As the lower spokes are obscured, so new spokes appear above, corresponding to new records.
The greatest fall from a(n) to a(n+1) is when a(n) is a prime record, then a(n+1)=2. This sometimes creates a remarkable disturbance in the sequence, involving many chaotic reversals before recovery to the more stable pattern described above is achieved. The first conspicuous example of this is at a(243)=89. A repeat of condition 0 can only be consequent to a novel highly divisible term k for which d(k) is also novel. For n<=2^19 the only occurrence of this is at a(422)=144, followed by d(144)=15. For more details see links.

Examples

			a(1)=1 is a novel term so a(2)=d(1)=1. Since a(2) is a repeat term and 1 is the least unused prior term, a(3)=a(2)+1=2. Then since 2 is a novel term, a(4)=d(a(3))=d(2)=2; and so on.
		

Crossrefs

Programs

  • Mathematica
    Block[{a = {1}, s = {}}, Do[If[FreeQ[#2, #1], AppendTo[a, DivisorSigma[0, a[[-1]]] ], AppendTo[a, a[[-1]] + First[s] ]; Set[s, Rest@ s]] & @@ {First[#1], #2} & @@ TakeDrop[a, -1]; Set[s, Insert[s, a[[-2]], LengthWhile[s, # < a[[-2]] &] + 1]], 105]; a] (* Michael De Vlieger, Jun 15 2021 *)

A091800 Largest n-digit number with maximal number of distinct prime divisors.

Original entry on oeis.org

6, 90, 990, 9870, 99330, 930930, 9699690, 99981420, 999068070, 9592993410, 99978788910, 999890501610, 9814524629910, 99999887777790, 999192361827660, 9999999768941490, 99992911041433410, 997799870344687410, 9999847102571786460, 99987077573596883670, 999999011467253427630, 9999928946485603635510
Offset: 1

Views

Author

Amarnath Murthy, Feb 21 2004

Keywords

Examples

			a(4) = 9870 as the largest number of distinct prime factors any 4-digit number can have and any number 9871 <= k <= 9999 has fewer than 5 prime factors. - _David A. Corneth_, Aug 19 2025
		

Crossrefs

Programs

  • Mathematica
    a[n_] := Module[{k=0, p=1, r=1, t=10^n}, While[r < t, p = NextPrime[p]; r *= p; k++]; k--; m = t-1; While[PrimeNu[m] != k, m--]; m]; Array[a, 8] (* Amiram Eldar, Mar 03 2020 *)
  • Python
    from sympy import nextprime, factorint
    def A091800(n: int) -> int:
        k, p, r, t = 0, 1, 1, 10**n
        while r < t:
           p = nextprime(p)
           r *= p
           k += 1
        m = t - 1
        while len(factorint(m)) != k - 1: m -= 1
        return m # John Reimer Morales, Aug 18 2025
    
  • Python
    # see linked program

Extensions

Edited, corrected and extended by Ray Chandler, Feb 23 2004
a(10)-a(12) from Amiram Eldar, Mar 03 2020
a(13) from Giovanni Resta, Mar 04 2020
a(14) onwards from John Reimer Morales and David A. Corneth, Aug 19 2025

A140645 Digit sums of highly composite numbers.

Original entry on oeis.org

1, 2, 4, 6, 3, 6, 9, 12, 6, 3, 9, 6, 9, 9, 12, 9, 15, 9, 9, 18, 9, 9, 9, 9, 18, 18, 9, 18, 18, 18, 18, 18, 18, 18, 36, 18, 27, 18, 18, 18, 18
Offset: 1

Views

Author

Parthasarathy Nambi, Jul 08 2008

Keywords

Comments

Except for the first three terms, all others are multiples of 3.

Examples

			2162160 is a highly composite number whose digit sum is 18.
		

Crossrefs

Programs

  • Mathematica
    HCN=NestList[Function[last,Module[{d = DivisorSigma[0, last]},NestWhile[# + 1 &, last, DivisorSigma[0, #] <= d &]]], 1,40]; DigitSum/@HCN (* James C. McMahon, Jul 12 2025 *)

Formula

a(n) = A007953(A002182(n)). - Amiram Eldar, Nov 13 2019

A172438 Numbers k such that tau(k^2+1) - tau(k^2) = 1 where the function tau(k) is the number of positive divisors of k.

Original entry on oeis.org

1, 3, 5, 11, 19, 27, 29, 59, 61, 71, 79, 101, 125, 131, 139, 181, 199, 242, 243, 271, 333, 349, 379, 387, 409, 423, 449, 461, 477, 521, 569, 571, 603, 631, 641, 661, 739, 747, 751, 772, 788, 821, 881, 929, 991, 1017, 1031, 1039, 1051, 1058, 1069, 1075, 1083
Offset: 1

Views

Author

Michel Lagneau, Feb 02 2010

Keywords

Comments

Square roots of perfect squares in A055927. [Juri-Stepan Gerasimov, Apr 06 2011]

Examples

			k=1, tau(2) - tau(1) = 2 - 1 = 1.
k=3, tau(10) - tau(9) = 4 - 3 = 1.
k=5, tau(26) - tau(25) = 4 - 3 = 1.
k=387, tau(149770)- tau(149769) = 16 - 15 = 1.
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 840.
  • T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 38.
  • D. S. Mitrinovic et al., Handbook of Number Theory, Kluwer, Chap. II. (For inequalities, etc.)

Crossrefs

Programs

  • Magma
    [m:m in [1..1100]| #Divisors(m^2+1) - #Divisors(m^2) eq 1]; // Marius A. Burtea, Jul 12 2019
  • Maple
    with(numtheory): for n from 1 to 100000 do; if tau(n^2+1)-tau(n^2)= 1 then print(n); else fi ; od;
  • Mathematica
    dsQ[n_]:=Module[{n2=n^2},DivisorSigma[0,n2+1]-DivisorSigma[0,n2]==1]; Select[Range[1200],dsQ] (* Harvey P. Dale, May 05 2011 *)
    Select[Sqrt[#]&/@Flatten[Position[Partition[DivisorSigma[0,Range[1200000]],2,1],?(#[[2]]-#[[1]]==1&),1,Heads->False]],IntegerQ] (* _Harvey P. Dale, Apr 09 2022 *)

A262504 First differences of A262502.

Original entry on oeis.org

1, 2, 4, 4, 6, 6, 8, 9, 10, 12, 12, 12, 12, 12, 16, 16, 16, 16, 18, 18, 18, 20, 20, 20, 20, 20, 20, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 30, 30, 30, 30, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 36, 36, 36, 36, 36, 36, 36, 36, 36, 36, 36, 36, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 48
Offset: 0

Views

Author

Antti Karttunen, Sep 24 2015

Keywords

Comments

Which of the terms of A002183 occur here, and for how many times each?

Crossrefs

Programs

Formula

a(n) = A262502(n+1) - A262502(n).

A273353 Number of divisors of A067128(n).

Original entry on oeis.org

1, 2, 2, 3, 4, 4, 4, 6, 6, 6, 8, 8, 9, 10, 12, 12, 12, 12, 12, 12, 16, 16, 18, 20, 20, 24, 24, 24, 24, 24, 24, 24, 24, 24, 30, 32, 32, 36, 36, 40, 40, 48, 48, 48, 48, 48, 48, 48, 48, 60, 64, 64, 72, 72, 72, 80, 80, 84, 90, 96, 96, 96, 96, 96, 96, 96, 96, 96, 100, 108, 120, 120, 120, 128, 128, 144, 144, 144, 144, 144, 160
Offset: 1

Views

Author

David A. Corneth, May 20 2016

Keywords

Comments

Is a(n + 1) / a(n) ~ 1 for large n?
Every term in this sequence also appears in A002183, where every element of this sequence occurs exactly once.
In A067128 it is asked if A034287 = A067128. If that is the case then this sequence is also the number of divisors of A034287.

Crossrefs

Programs

  • Mathematica
    s = {}; dmax = 0; Do[d = DivisorSigma[0, n]; If[d >= dmax, AppendTo[s, d]; dmax = d], {n, 1, 10^6}]; s (* Amiram Eldar, Jun 07 2019 *)
  • PARI
    is_a067128(n) = my(nd=numdiv(n)); for(k=1, n-1, if(numdiv(k) > nd, return(0))); return(1)
    for(n=1, 50000, if(is_a067128(n), print1(numdiv(n), ", "))) \\ Felix Fröhlich, May 24 2016

Formula

a(n) = A000005(A067128(n)).

A323393 a(n) is the number of divisors of A323392(n) in Eisenstein integers.

Original entry on oeis.org

1, 2, 3, 6, 9, 10, 12, 15, 24, 36, 40, 48, 60, 72, 80, 96, 100, 144, 160, 192, 240, 288, 320, 324, 336, 384, 400, 432, 480, 576, 640, 648, 768, 960, 1152, 1280, 1296, 1344, 1536, 1600, 1728, 1920, 2160, 2560, 2592, 2880, 3200, 3456, 3600, 3840, 4320, 4608, 5120, 5760, 6144, 6400, 7200, 7680
Offset: 1

Views

Author

Jianing Song, Jan 13 2019

Keywords

Comments

Records in A319442.
Analog of A002183 and A302249, which list the records of number of divisors in rational integers and Gaussian integers respectively.
It seems that 15 is the largest odd term.

Examples

			252 has 60 divisors up to association in Eisenstein integers, more than any previous positive integers, so 60 is a term.
		

Crossrefs

For the numbers whose number of divisors set new records see A323392.

Programs

  • Mathematica
    f[p_, e_] := Switch[Mod[p, 3], 0, 2*e + 1, 1, (e + 1)^2, 2, e + 1]; eisNumDiv[1] = 1; eisNumDiv[n_] := Times @@ f @@@ FactorInteger[n]; seq = {}; emax = 0; Do[eis = eisNumDiv[n]; If[eis > emax, emax = eis; AppendTo[seq, eis]], {n, 1, 10^6}]; seq (* Amiram Eldar, Mar 02 2020 *)
  • PARI
    my(r=0, t); for(n=1, 10^6, t=A319442(n); if(t>r, r=t; print1(r, ", ")));

Formula

a(n) = A319442(A323392(n)).

A349607 a(n) is the number of divisors of the n-th superabundant number (A004394).

Original entry on oeis.org

1, 2, 3, 4, 6, 8, 9, 10, 12, 16, 18, 20, 24, 30, 32, 36, 40, 48, 60, 72, 80, 90, 96, 120, 144, 160, 180, 192, 216, 224, 240, 288, 320, 360, 384, 432, 448, 480, 576, 640, 720, 768, 864, 896, 960, 1152, 1280, 1344, 1440, 1440, 1536, 1728, 1792, 1920, 2304, 2688
Offset: 1

Views

Author

Amiram Eldar, Nov 23 2021

Keywords

Comments

First differs from A002183 at n = 20.
The first position n where a(n+1) = a(n) is n = 49: a(49) = a(50) = 1440.
The first position n where a(n+1) < a(n) is n = 173: a(173) = 5308416 and a(174) = 5160960.

Examples

			a(1) = A000005(A004394(1)) = A000005(1) = 1.
a(10) = A000005(A004394(10)) = A000005(120) = 16.
		

Crossrefs

Programs

  • Mathematica
    s = {}; rm = 0; Do[r = DivisorSigma[1, n]/n; If[r > rm, rm = r; AppendTo[s, DivisorSigma[0, n]]], {n, 1, 10^5}]; s

Formula

a(n) = A000005(A004394(n)).
Previous Showing 51-60 of 72 results. Next