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

A232100 a(n) = the position of A232099(n) in A055926.

Original entry on oeis.org

41, 123, 205, 287, 369, 451, 533, 615, 697, 779, 861, 943, 1025, 1107, 1189, 1271, 1353, 1435, 1517, 1599, 1681, 1763, 1845, 1927, 2009, 2091, 2173, 2255, 2337, 2419, 2501, 2583, 2665, 2706, 2747, 2829, 2911, 2993, 3075, 3157, 3239, 3321, 3403, 3485, 3567, 3649
Offset: 1

Views

Author

Antti Karttunen, Nov 18 2013

Keywords

Comments

In range n=1..2215, each a(n) is a multiple of 41, with a(2215) = 177079 = 4319*41. In that range, all 2160 odd multiples 1, 3, 5, ..., 4319*41 occur, but only 55 even multiples, with a(34) = 66*41 = 2706 being the first one of them. The multipliers for even terms in that range seems to be given by 66*A047253(1..55) from 66*1 up to 66*A047253(55) = 66*65 = 4290. (Where A047253 gives the numbers not divisible by six).

Crossrefs

Formula

A055926(a(n)) = A232099(n) for all n.

A055874 a(n) = largest m such that 1, 2, ..., m divide n.

Original entry on oeis.org

1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2
Offset: 1

Views

Author

Leroy Quet, Jul 16 2000

Keywords

Comments

From Antti Karttunen, Nov 20 2013 & Jan 26 2014: (Start)
Differs from A232098 for the first time at n=840, where a(840)=8, while A232098(840)=7. A232099 gives all the differing positions. See also the comments at A055926 and A232099.
The positions where a(n) is an odd prime is given by A017593 up to A017593(34)=414 (so far all 3's), after which comes the first 7 at a(420). (A017593 gives the positions of 3's.)
(Continued on Jan 26 2014):
Only terms of A181062 occur as values.
A235921 gives such n where a(n^2) (= A235918(n)) differs from A071222(n-1) (= A053669(n)-1). (End)
a(n) is the largest m such that A003418(m) divides n. - David W. Wilson, Nov 20 2014
a(n) is the largest number of consecutive integers dividing n. - David W. Wilson, Nov 20 2014
A051451 gives indices where record values occur. - Gionata Neri, Oct 17 2015
Yuri Matiyasevich calls this the maximum inheritable divisor of n. - N. J. A. Sloane, Dec 14 2023

Examples

			a(12) = 4 because 1, 2, 3, 4 divide 12, but 5 does not.
		

Crossrefs

Programs

  • Haskell
    a055874 n = length $ takeWhile ((== 0) . (mod n)) [1..]
    -- Reinhard Zumkeller, Feb 21 2012, Dec 09 2010
    
  • Maple
    N:= 1000: # to get a(1) to a(N)
    A:= Vector(N,1);
    for m from 2 do
      Lm:= ilcm($1..m);
      if Lm > N then break fi;
      if Lm mod (m+1) = 0 then next fi;
      for k from 1 to floor(N/Lm) do
        A[k*Lm]:=m
      od
    od:
    convert(A,list); # Robert Israel, Nov 28 2014
  • Mathematica
    a[n_] := Module[{m = 1}, While[Divisible[n, m++]]; m - 2]; Array[a, 100] (* Jean-François Alcover, Mar 07 2016 *)
  • PARI
    a(n) = my(m = 1); while ((n % m) == 0, m++); m - 1; \\ Michel Marcus, Jan 17 2014
    
  • Python
    from itertools import count
    def A055874(n):
        for m in count(1):
            if n % m:
                return m-1 # Chai Wah Wu, Jan 02 2022
  • Scheme
    (define (A055874 n) (let loop ((m 1)) (if (not (zero? (modulo n m))) (- m 1) (loop (+ 1 m))))) ;; Antti Karttunen, Nov 18 2013
    

Formula

a(n) = A007978(n) - 1. - Antti Karttunen, Jan 26 2014
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = A064859 (Farhi, 2009). - Amiram Eldar, Jul 25 2022

A055926 Numbers k such that {largest m such that 1, 2, ..., m divide k} is different from {largest m such that m! divides k}; numbers k which are either odd multiples of 12 or the largest m such that (m-1)! divides k is a composite number > 5.

Original entry on oeis.org

12, 36, 60, 84, 108, 120, 132, 156, 180, 204, 228, 240, 252, 276, 300, 324, 348, 360, 372, 396, 420, 444, 468, 480, 492, 516, 540, 564, 588, 600, 612, 636, 660, 684, 708, 732, 756, 780, 804, 828, 840, 852, 876, 900, 924, 948, 960, 972, 996, 1020, 1044, 1068
Offset: 1

Views

Author

Leroy Quet, Jul 16 2000

Keywords

Comments

From Antti Karttunen, Nov 20 - Dec 06 2013: (Start)
This sequence has several interpretations:
Numbers k such that A055874(k) differs from A055881(k). [Leroy Quet's original definition of the sequence. Note that A055874(k) >= A055881(k) for all k.]
Numbers k such that {largest m such that m! divides k^2} is different from {largest m such that m! divides k}, i.e., numbers k for which A232098(k) > A055881(k).
Numbers k which are either 12 times an odd number (A073762) or the largest m such that (m-1)! divides k is a composite number > 5 (A232743).
Please see my attached notes for the proof of the equivalence of these interpretations.
Additional implications based on that proof:
A232099 is a subset of this sequence.
A055881(a(n))+1 is always composite. In the range n = 1..17712, only values 4, 6, 8, 9 and 10 occur.
The new definition can be also rephrased by saying that the sequence contains all the positive integers k whose factorial base representation of (A007623(k)) either ends as '...200' (in which case k is an odd multiple of 12, 12 = '200', 36 = '1200', 60 = '2200', ...) or the number of trailing zeros + 2 in that representation is a composite number greater than or equal to 6, e.g. 120 = '10000' (in other words, A055881(k) is one of the terms of A072668 after the initial 3). Together these conditions also imply that all the terms are divisible by 12.
(End)

Examples

			12 is included because 3! is the largest factorial to divide 12, but 1, 2, 3 and 4 all divide 12. Equally, 12 is included because it is one of the terms of A073762, or equally, because its factorial base representation ends with digits '...200': A007623(12) = 200.
840 (= 3*5*7*8) is included because the largest factorial which divides 840 is 5! (840 = 7*120), but all positive integers up to 8 divide 840. Equally, 840 is included because it is one of the terms of A232743 as 5+1 = 6 is a composite number larger than 5. Note that A007623(840) = 110000.
		

Crossrefs

Union of A073762 and A232743. Equivalently, setwise difference of A232742 and A017593. Subset: A232099.

Extensions

More terms from Antti Karttunen, Dec 01 2013

A232098 a(n) is the largest m such that m! divides n^2; a(n) = A055881(n^2).

Original entry on oeis.org

1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 4, 1, 2, 1
Offset: 1

Views

Author

Antti Karttunen, Nov 18 2013

Keywords

Comments

For all n, A055881(n) <= a(n), and probably also a(n) <= A055874(n).
Moreover, a(n) > A055881(n) if and only if A055874(n) > A055881(n), thus A055926 gives (also) all the positions where this sequence differs from A055881. Please see Comments section in A055926 for the proof.
Differs from A055874 for the first time at n=840, where a(840)=7, while A055874(840)=8. A232099 gives all the positions where such differences occur.

Crossrefs

Programs

  • Mathematica
    Module[{nn=10,fct},fct=Table[{f,f!},{f,nn}];Table[Select[fct,Mod[n^2,#[[2]]]==0&][[-1,1]],{n,90}]] (* Harvey P. Dale, Aug 11 2024 *)
  • Scheme
    (define (A232098 n) (A055881 (A000290 n)))

Formula

a(n) = A055881(A000290(n)) = A055881(n^2).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{k>=1} 1/A065887(k) = 1.78672776922161809767... . - Amiram Eldar, Jan 01 2024

A235921 Numbers n such that smallest number not dividing n^2 (A236454) is different from smallest prime not dividing n (A053669).

Original entry on oeis.org

210, 630, 1050, 1470, 1890, 2310, 2730, 3150, 3570, 3990, 4410, 4830, 5250, 5670, 6090, 6510, 6930, 7350, 7770, 8190, 8610, 9030, 9450, 9870, 10290, 10710, 11130, 11550, 11970, 12390, 12810, 13230, 13650, 14070, 14490, 14910, 15330, 15750, 16170, 16590, 17010
Offset: 1

Views

Author

Antti Karttunen and Michel Marcus, Jan 17 2014

Keywords

Comments

Equivalent definition is: numbers n such that {the largest m such that 1, 2, ..., m divide n^2 = A055874(n^2) = A235918(n)} is different from {the smallest k such that gcd(n-1,k) = gcd(n,k+1) = A071222(n-1)}.
All terms are multiples of 210 = 2*3*5*7, the fourth primorial, A002110(4).
The first term which is an even multiple of 210 (i.e., 210 times an even number), is 446185740 = 2124694 * 210 = 2*223092870 = 2*A002110(9) = 2*A034386(23). Note that 23 is the 9th prime, and 223092870 is its primorial. Thus this sequence differs from its subsequence, A236432, "the odd multiples of 210" = (2n-1)*210, for the first time at n = 1062348, where a(n) = 446185740, while A236432(n) = 446185950.
Note that a more comprehensive description for which terms are included is still lacking. Compare for example to the third definition of A055926.
At least we know the following:
If a number is not divisible by 210, then it cannot be a member, as then it is "missing" (i.e., not divisible by) one of those primes, 2, 3, 5 or 7, and thus its square is also "missing" the same prime. In more detail, this follows because:
If the least nondividing prime is 2, then A053669(n) = A236454(n) = 2. If the least nondividing prime is 3, then A053669(n) = A236454(n) = 3.
If the least nondividing prime is 5 (so 2 and 3 are present), then as 2|n and 4|(n^2), we have A053669(n) = A236454(n) = 5.
If the least nondividing prime is 7, but 2, 3 and 5 are present, then we have A053669(n) = A236454(n) = 7.
On the other hand, when n is an odd multiple of 210 (= 2*3*5*7), i.e., (2k+1)*210, so that its prime factorization is of the form 2*3*5*7*{zero or more additional odd prime factors}, then A053669(n) must be at least 11, the next prime after 7, which is certainly different from A236454(n) = A007978(n^2) which must be 8, as then 4 is the highest power of 2 dividing n^2.
In contrast to that, when n is an even multiple of 210, so that its prime factorization is of the form 2*2*3*5*7*{zero or more additional prime factors}, then also all the composites 8, 9, 10, 12, 14, 15, 16, 18 and 20 divide n^2, thus if A053669(n) is any prime from 11 to 19, A236454(n) will return the same result.
However, if n is of the form k*446185740, where k is not a multiple of 3, so that the prime factorization of n is 2*2*3*5*7*11*13*17*19*23*{zero or more additional prime factors, all different from 3}, then A053669(n) must be at least 29 (next prime after 23), but A236454(n) = 27, because then 9 is the highest power of 3 dividing n^2.
The pattern continues indefinitely: If n is of the form (2k+1)*2*3*200560490130, where 200560490130 = A002110(11), so that n has a prime factorization of the form 2*2*3*3*5*7*11*13*17*19*23*29*31*{zero or more additional odd prime factors}, then A053669(n) must be at least 37, while A236454(n) = 32 = 2^5, because then 16 is the highest power of 2 dividing n^2.

Examples

			210 (= 2*3*5*7) is a member, because A053669(210)=11, while A236454(210) = A007978(210*210) = A007978(44100) = 8.
446185740 (= 2*2*3*5*7*11*13*17*19*23) is a member, because A053669(446185740) = 29, while A236454(446185740) = 27, as there is only one 3 present in 446185740, which means that its square is only divisible by 9, but not by 27 = 3^3.
		

Crossrefs

Showing 1-5 of 5 results.