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

A319717 Filter sequence combining the largest proper divisor of n (A032742) with modulo 6 residue of the smallest prime factor, A010875(A020639(n)), and a single bit A319710(n) telling whether the smallest prime factor is unitary.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 5, 11, 7, 12, 13, 14, 5, 15, 7, 16, 17, 18, 5, 19, 20, 21, 22, 23, 5, 24, 7, 25, 26, 27, 28, 29, 7, 30, 31, 32, 5, 33, 7, 34, 35, 36, 5, 37, 38, 39, 40, 41, 5, 42, 43, 44, 45, 46, 5, 47, 7, 48, 49, 50, 51, 52, 7, 53, 54, 55, 5, 56, 7, 57, 58, 59, 60, 61, 7, 62, 63, 64, 5, 65, 66, 67, 68, 69, 5, 70, 71, 72, 73, 74, 75, 76, 7, 77, 78, 79, 5, 80, 7, 81, 82, 83, 5, 84, 7, 85, 86, 87, 5, 88, 89, 90, 91, 92, 93, 94, 95
Offset: 1

Views

Author

Antti Karttunen, Oct 04 2018

Keywords

Comments

Restricted growth sequence transform of triple [A010875(A020639(n)), A032742(n), A319710(n)] (with a separate value allotted for a(1)), or equally, of ordered pair [A319716(n), A319710(n)].
In addition to A319716, this filter sequence also records in the value of a(n) also the fact whether the smallest prime factor of n is unitary or not. This information is enough to determine the modulo 6 residues of all the divisors of n, thus sequences like A002324 are essentially functions of this sequence. Moreover, a lot of other information is immediately (and unavoidably) present, for example the exact prime signature of n, including also the relative order of exponents.
Any such filtering sequence can be perceived also in terms of what information it leaves out from a(n) that would be needed to reconstruct whole n from each a(n). If the whole n could be reconstructed from a(n) each time, then sequence a would be injective, and would be useless for filtering, because then it would match with any sequence. In this filter, what is left out is only the exact identity of the smallest prime factor, although its residue class mod 6 is retained. However, when the smallest prime factor is 2 or 3, this can be seen from that residue value, so for any number x in A047229, both A020639(x) and A032742(x) are known, and as x = A020639(x)*A032742(x), it means such numbers must occur in their own singleton equivalence classes.
Likewise, for any n in A283050, even if not divisible by 2 or 3, when we have A319710(n) stored in the triple as 1, this immediately gives away the exact identity of the smallest prime factor, which is equal to A014673(n) = A020639(A032742(n)) in these cases.
Thus there is a substantial subset of N (containing at least the union of A047229 and A283050) which is actually in the "blind sector" of this filter, "where anything goes", as this sequence obtains only unique values in that subdomain.
There is a related filter sequence A319996, which operates by "cleaving n from its high end" (by storing the residue class of the largest prime factor, A006530, instead of the smallest, together with n/A006530(n)), which has its own blind spots, but fortunately, they do not fully coincide with the blind spots of this filter. Naturally, any sequence like A002324 should match both to this sequence and A319996.
For all i, j:
a(i) = a(j) => A002324(i) = A002324(j),
a(i) = a(j) => A067029(i) = A067029(j),
a(i) = a(j) => A071178(i) = A071178(j),
a(i) = a(j) => A077462(i) = A077462(j) => A101296(i) = A101296(j),
a(i) = a(j) => A319716(i) = A319716(j) => A319690(i) = A319690(j).

Examples

			For n = 65 = 5*13 and 143 = 11*13, the smallest prime factor is of the form 6k+5,  doesn't occur more than once in the factorization, and the largest proper divisor is the same number (13) in both cases, thus a(65) = a(143) (= 51, a running count value allotted by rgs-transform for this equivalence class).
For n = 1805 (5*19^2), 3971 (11*19^2), 6137 (17*19^2), it's like above, but the largest proper divisor is in all three cases 361 = 19^2, thus a(1805) = a(3971) = a(6137) (= 1405).
Note that such nontrivial equivalence classes may only contain numbers that are 5-rough, A007310, with no prime factors 2 or 3, and also, they may not contain numbers from A283050. See the comments section.
		

Crossrefs

Cf. also A320004 (analogous sequence for modulo 4 residues).
Differs from A319707 for the first time at n=143, where a(143) = 51, differs from A319716 for the first time at n=121, where a(121) = 95.

Programs

  • PARI
    up_to = 100000;
    rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om,invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om,invec[i],i); outvec[i] = u; u++ )); outvec; };
    A032742(n) = if(1==n,n,n/vecmin(factor(n)[,1]));
    A286476(n) = if(1==n,n,(6*A032742(n) + (n % 6)));
    A319710(n) = ((n>1)&&(factor(n)[1,2]>1));
    v319717 = rgs_transform(vector(up_to,n,[A286476(n),A319710(n)]));
    A319717(n) = v319717[n];

A253569 Composite numbers n = p_i * p_j * p_k * ... * p_u, p_i <= p_j <= p_k <= ... <= p_u, where each successive prime factor (when sorted into a nondecreasing order) is greater than the square of the previous: (p_i)^2 < p_j, (p_j)^2 < p_k, etc.

Original entry on oeis.org

10, 14, 22, 26, 33, 34, 38, 39, 46, 51, 57, 58, 62, 69, 74, 82, 86, 87, 93, 94, 106, 111, 118, 122, 123, 129, 134, 141, 142, 145, 146, 155, 158, 159, 166, 177, 178, 183, 185, 194, 201, 202, 205, 206, 213, 214, 215, 218, 219, 226, 235, 237, 249, 254, 262, 265, 267, 274, 278, 290
Offset: 1

Views

Author

Keywords

Comments

Numbers n = A020639(n) * A014673(n) * A054576(n), for which A020639(n)^2 < A014673(n) and either A054576(n) = 1 or A032742(n) satisfies the same condition (is the term of this sequence).

Examples

			290 = 2*5*29 is a member, because 2^2 < 5 and 5^2 < 29.
		

Crossrefs

Complement: A253567.
Subsequence of A002808, A005117, A088381, A251727, A245729 and A253785.
A138511 is a subsequence, from which this sequence differs for the first time at n=60, where A138511(60) = 291, while here a(60) = 290.

Programs

  • Haskell
    a253569 n = a253569_list !! (n-1)
    a253569_list = filter f [1..] where
                        f x = (p ^ 2 < a020639 q) && (a010051' q == 1 || f q)
                              where q = div x p; p = a020639 x
    -- Antti Karttunen after Reinhard Zumkeller's code for A138511, Jan 09 2015
    a253569 n = a253569_list !! (n-1)
    a253569_list = filter (not . f''') a002808_list where
       f''' x = p ^ 2 > a020639 q || (a010051 q == 0 && f''' q)
                where q = div x p; p = a020639 x
    -- Reinhard Zumkeller, Jan 12 2015
    (Scheme, with Antti Karttunen's IntSeq-library)
    (define A253569 (MATCHING-POS 1 1 (lambda (n) (and (> (A001222 n) 1) (numbers-sparsely-distributed? (ifactor n))))))
    (define (numbers-sparsely-distributed? lista) (cond ((null? lista) #t) ((null? (cdr lista)) #t) ((> (A000290 (car lista)) (cadr lista)) #f) (else (numbers-sparsely-distributed? (cdr lista)))))
    ;; Antti Karttunen, Jan 16 2015
  • Mathematica
    cnQ[n_]:=CompositeQ[n]&&Union[Boole[#[[2]]>#[[1]]^2&/@Partition[Flatten[Table[ #[[1]], #[[2]]]&/@FactorInteger[n]],2,1]]]=={1}; Select[Range[300],cnQ] (* Harvey P. Dale, Jul 10 2023 *)

A281891 Square array A(n,k): number of integers having k or more factors less than prime(n+1) in their prime factorization, within any interval of primorial(n)^k positive integers.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 4, 1, 0, 1, 14, 22, 1, 0, 1, 46, 412, 162, 1, 0, 1, 146, 7072, 22164, 1830, 1, 0, 1, 454, 115432, 2744088, 2822340, 24270, 1, 0, 1, 1394, 1827592, 319881696, 3913037880, 496348740, 418350, 1, 0, 1, 4246, 28390552, 35924741232, 5079363328560, 9082206410040, 147569907780, 8040810, 1
Offset: 0

Views

Author

Peter Munn, Feb 08 2017

Keywords

Comments

Square array read by descending antidiagonals; A(n,k) with rows n >= 0, columns k >= 0. Prime factors are counted with multiplicity. Primorial(n) = A002110(n): product of first n primes.
Visualize the prime factorizations of the positive integers as a table with row headings giving each successive integer, and the primes of which the row heading is the product listed across the columns in nondecreasing order, repeated when necessary. Except for 1, which lacks prime factors, column 1 has the row heading's least prime factor, column 2 has a value for composite numbers but is blank for primes, and so on. This sequence measures precisely how frequently values up to and including the various primes occur in each column. This is possible because any given prime occurs cyclically in any given column, for the reason following.
The occurrence pattern of up to k factors of prime(n) in such prime factorizations has a fundamental period over the positive integers of prime(n)^k. The least common period for up to k factors of each of the first n primes is primorial(n)^k, and this covers everything that can affect the occurrence of prime(n) in the least k factors. Thus prime(n) is k-th least prime factor of integer m if and only if it is k-th least prime factor of m + primorial(n)^k.
Intermediate values in the calculation of this sequence appear in A281890.
If n > 0, A(n,1) = A053144(n) in accordance with the comment on A053144 dated Apr 08 2010.
A(2,k) = A027649(k) = 2*(3^k) - 2^k.

Examples

			The table starts:
   1     0         0             0             0           0        0 ...
   1     1         1             1             1           1        1 ...
   1     4        14            46           146         454     1394 ...
   1    22       412          7072        115432     1827592 28390552 ...
   1   162     22164       2744088     319881696 35924741232    ...
   1  1830   2822340    3913037880 5079363328560      ...
   1 24270 496348740 9082206410040       ...
   ...
Primes less than prime(2+1)=5 occur as second least factor 14 times in the prime factorizations of every interval of 36 = primorial(2)^2 positive integers (cf. A014673). Therefore, A(2,2) = 14.
		

Crossrefs

A079474 re-read as a square array gives values of primorial(n)^k = A002110(n)^k.
The values in the body of the factorization table described in the author's comments are in the irregular array A027746.
A096294 gives the equivalent array for integers expressed as a product of prime powers.

Formula

A(n,0) = 1 for n >= 0, A(0,k) = 0 for k >= 1.
A(n,k) = prime(n)^k * A(n-1,k) + A281890(n,k) for n >= 1, k >= 1.

Extensions

Edited by M. F. Hasler, Apr 14 2017

A085393 Difference between the largest and the smallest prime factor of the greatest proper divisor of n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 1, 0, 0, 0, 5, 0, 2, 0, 0, 0, 0, 0, 1, 0, 0, 0, 3, 0, 4, 0, 9, 2, 0, 0, 1, 0, 0, 0, 11, 0, 0, 0, 5, 0, 0, 0, 3, 0, 0, 4, 0, 0, 8, 0, 15, 0, 2, 0, 1, 0, 0, 0, 17, 0, 10, 0, 3, 0, 0, 0, 5, 0, 0, 0, 9, 0, 2, 0, 21, 0, 0, 0, 1, 0, 0, 8, 3, 0, 14
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Mathematica
    PrimeFactors[n_] := Flatten[ Table[ # [[1]], {1}] & /@ FactorInteger[n]]; f[n_] := Block[{gpd = Divisors[n][[ -2]]}, If[gpd == 1, 0, PrimeFactors[gpd][[ -1]] - PrimeFactors[gpd][[1]] ]]; Table[ If[n == 1, 0, f[n]], {n, 1, 102}]
    {1}~Join~Array[#[[-1, 1]] - #[[1, 1]] &@ FactorInteger@ Last@ Most@ Divisors@ # &, 101, 2] (* Michael De Vlieger, Dec 03 2017 *)

Formula

a(n) = A085392(n) - A014673(n).

A285323 a(n) = A065642(A065642(A019565(n))) / A019565(n).

Original entry on oeis.org

1, 4, 9, 3, 25, 4, 5, 3, 49, 4, 7, 3, 7, 4, 5, 3, 121, 4, 9, 3, 11, 4, 5, 3, 11, 4, 7, 3, 7, 4, 5, 3, 169, 4, 9, 3, 13, 4, 5, 3, 13, 4, 7, 3, 7, 4, 5, 3, 13, 4, 9, 3, 11, 4, 5, 3, 11, 4, 7, 3, 7, 4, 5, 3, 289, 4, 9, 3, 17, 4, 5, 3, 17, 4, 7, 3, 7, 4, 5, 3, 17, 4, 9, 3, 11, 4, 5, 3, 11, 4, 7, 3, 7, 4, 5, 3, 17, 4, 9, 3, 13, 4, 5, 3, 13, 4, 7, 3, 7, 4, 5, 3
Offset: 0

Views

Author

Antti Karttunen, Apr 19 2017

Keywords

Comments

After the initial a(0)=1, the third row of array A285321 divided by its first row. After 1, all terms are either primes or squares of primes. See A285110.
The sequence is completely determined by the positions of two least significant 1-bits of n: After initial zero, if n is a power of two (only one 1-bit present) or if prime(1+A285099(n)) > prime(1+A007814(n))^2, a(n) = prime(1+A007814(n))^2 = A020639(A019565(n))^2, otherwise a(n) = prime(1+A285099(n)) = A014673(A019565(n)).

Crossrefs

Programs

Formula

a(n) = A065642(A065642(A019565(n))) / A019565(n).

A343597 Numbers divisible by a 7-smooth composite number.

Original entry on oeis.org

4, 6, 8, 9, 10, 12, 14, 15, 16, 18, 20, 21, 24, 25, 27, 28, 30, 32, 35, 36, 40, 42, 44, 45, 48, 49, 50, 52, 54, 56, 60, 63, 64, 66, 68, 70, 72, 75, 76, 78, 80, 81, 84, 88, 90, 92, 96, 98, 99, 100, 102, 104, 105, 108, 110, 112, 114, 116, 117, 120, 124, 125, 126, 128, 130
Offset: 1

Views

Author

Peter Munn, Apr 21 2021

Keywords

Comments

Numbers divisible by at least one of 4, 6, 9, 10, 14, 15, 21, 25, 35, 49.
Exactly half of the first 10, first 100 and first 600 positive integers are divisible by a 7-smooth composite number; the largest 7-smooth divisor of the remaining numbers is 1, 2, 3, 5 or 7.
Intervals extending to hundreds of integers with exactly 50% membership of this sequence are far from rare, some notable examples being [3000, 3999], [8000, 8999], [20000, 20999], [21000, 21999] and [23000, 23999]. This reflects the asymptotic density of the corresponding set being close to 0.5, precisely 1847 / 3675 = 0.50258503... (and membership of the set has a periodic pattern). See A343598 for further information.

Examples

			33 = 11 * 3 has divisors 1, 3, 11, 33, of which only 33 is composite. 33 is not 7-smooth, as its prime factors include 11, which is greater than 7. So 33 is not in the sequence.
52 = 13 * 2 * 2 is divisible by 4, which is composite and 7-smooth, so 52 is in the sequence.
		

Crossrefs

Cf. A002473 (7-smooth numbers), A014673, A020639, A210615 (smallest smoothest composite divisor), A343598.
Subsequence of A080672.

Programs

  • Mathematica
    Select[Range[130], Plus @@ IntegerExponent[#, {2, 3, 5, 7}] > 1 &] (* Amiram Eldar, May 04 2021 *)

Formula

{a(n)} = {k : k >= 1, 2 <= A014673(k) <= 7}, where A014673(k) = lpf(k/lpf(k)), where lpf(m) = A020639(m), the least prime factor of m.
For n >= 1, a(22164 + n) = 44100 + a(n).
For n < 22164, a(22164 - n) = 44100 - a(n).
Previous Showing 11-16 of 16 results.