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

A382204 Number of normal multiset partitions of weight n into constant blocks with a common sum.

Original entry on oeis.org

1, 1, 2, 3, 4, 4, 7, 5, 8, 8, 10, 8, 15, 9, 14, 15, 17, 13, 22, 14, 25, 21, 23, 19, 34, 24, 29, 28, 37, 27, 45, 29, 44, 38, 43, 43, 59, 40, 51, 48, 69, 48, 71, 52, 73, 69, 72, 61, 93, 72, 91, 77, 99, 78, 105, 95, 119, 95, 113, 96, 146, 107, 126, 123, 151, 130
Offset: 0

Views

Author

Gus Wiseman, Mar 26 2025

Keywords

Comments

We call a multiset or multiset partition normal iff it covers an initial interval of positive integers. The weight of a multiset partition is the sum of sizes of its blocks.

Examples

			The a(1) = 1 through a(6) = 7 multiset partitions:
  {1} {11}   {111}     {1111}       {11111}         {111111}
      {1}{1} {2}{11}   {11}{11}     {2}{11}{11}     {111}{111}
             {1}{1}{1} {2}{2}{11}   {2}{2}{2}{11}   {22}{1111}
                       {1}{1}{1}{1} {1}{1}{1}{1}{1} {11}{11}{11}
                                                    {2}{2}{11}{11}
                                                    {2}{2}{2}{2}{11}
                                                    {1}{1}{1}{1}{1}{1}
The a(1) = 1 through a(7) = 5 factorizations:
  2  4    8      16       32         64           128
     2*2  3*4    4*4      3*4*4      8*8          3*4*4*4
          2*2*2  3*3*4    3*3*3*4    9*16         3*3*3*4*4
                 2*2*2*2  2*2*2*2*2  4*4*4        3*3*3*3*3*4
                                     3*3*4*4      2*2*2*2*2*2*2
                                     3*3*3*3*4
                                     2*2*2*2*2*2
		

Crossrefs

Without a common sum we have A055887.
Twice-partitions of this type are counted by A279789.
Without constant blocks we have A326518.
For distinct block-sums and strict blocks we have A381718.
Factorizations of this type are counted by A381995.
For distinct instead of equal block-sums we have A382203.
For strict instead of constant blocks we have A382429.
A000670 counts patterns, ranked by A055932 and A333217, necklace A019536.
A001055 count multiset partitions of prime indices, strict A045778.
A089259 counts set multipartitions of integer partitions.
A255906 counts normal multiset partitions, row sums of A317532.
A321469 counts multiset partitions with distinct block-sums, ranks A326535.
Normal multiset partitions: A035310, A304969, A356945.
Set multipartitions: A116540, A270995, A296119, A318360.
Set multipartitions with distinct sums: A279785, A381806, A381870.
Constant blocks with distinct sums: A381635, A381636, A381716.

Programs

  • Mathematica
    allnorm[n_Integer]:=Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[mset_]:=Union[Sort[Sort/@(#/.x_Integer:>mset[[x]])]&/@sps[Range[Length[mset]]]];
    Table[Length[Join@@(Select[mps[#],SameQ@@Total/@#&&And@@SameQ@@@#&]&/@allnorm[n])],{n,0,5}]
  • PARI
    h(s,x)=my(t=0,p=1,k=1);while(s%k==0,p*=1/(1-x^(s/k))-1;t+=p;k+=1);t
    lista(n)=Vec(1+sum(s=1,n,h(s,x+O(x*x^n)))) \\ Christian Sievers, Apr 05 2025

Formula

G.f.: 1 + Sum_{s>=1} Sum_{k=1..A055874(s)} Product_{v=1..k} (1/(1-x^(s/v)) - 1). - Christian Sievers, Apr 05 2025

Extensions

Terms a(16) and beyond from Christian Sievers, Apr 04 2025

A328194 Maximum length of a divisibility chain of consecutive nontrivial divisors of n (greater than 1 and less than n).

Original entry on oeis.org

0, 0, 0, 1, 0, 1, 0, 2, 1, 1, 0, 1, 0, 1, 1, 3, 0, 2, 0, 2, 1, 1, 0, 1, 1, 1, 2, 2, 0, 1, 0, 4, 1, 1, 1, 1, 0, 1, 1, 2, 0, 2, 0, 2, 1, 1, 0, 1, 1, 2, 1, 2, 0, 2, 1, 2, 1, 1, 0, 1, 0, 1, 1, 5, 1, 2, 0, 2, 1, 1, 0, 1, 0, 1, 2, 2, 1, 2, 0, 2, 3, 1, 0, 1, 1, 1, 1, 3, 0, 1, 1, 2, 1, 1, 1, 1, 0, 2, 2, 3, 0, 2, 0, 3, 1
Offset: 1

Views

Author

Gus Wiseman, Oct 14 2019

Keywords

Comments

The nontrivial divisors of n are row n of A163870.

Examples

			The nontrivial divisors of 272 are {2, 4, 8, 16, 17, 34, 68, 136} with divisibility chains {{2, 4, 8, 16}, {17, 34, 68, 136}}, so a(272) = 4.
		

Crossrefs

Positions of 1's are A328028 without 1.
The version with all divisors allowed is A328162.
Allowing n as a divisor of n gives A328195.
Indices of terms greater than 1 are A328189.
The maximum run-length of divisors of n is A055874(n).

Programs

  • Mathematica
    Table[Switch[n,1,0,?PrimeQ,0,,Max@@Length/@Split[DeleteCases[Divisors[n],1|n],Divisible[#2,#1]&]],{n,100}]
  • PARI
    A328194(n) = if(1==n || isprime(n), 0, my(divs=divisors(n), rl=0,ml=1); for(i=2,#divs-1,if(!(divs[i]%divs[i-1]), rl++, ml = max(rl,ml); rl=1)); max(ml,rl)); \\ Antti Karttunen, Dec 07 2024

Extensions

Data section extended up to a(105) by Antti Karttunen, Dec 07 2024

A328457 Length of the longest run of divisors > 1 of n.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 5, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 1, 1, 1, 2, 1, 2, 1, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Oct 16 2019

Keywords

Crossrefs

Records occur at A328448.
Positions of 0's and 1's are A088725.
The version that looks at all divisors (including 1) is A055874.
The number of successive pairs of divisors > 1 of n is A088722(n).
The Heinz number of the multiset of run-lengths of divisors of n is A328166(n).
The longest run of nontrivial divisors of n is A328458(n).

Programs

  • Mathematica
    Table[If[n==1,0,Max@@Length/@Split[Rest[Divisors[n]],#2==#1+1&]],{n,100}]
  • PARI
    A328457(n) = { my(rl=0,pd=0,m=0); fordiv(n, d, if(d>1, if(d>(1+pd), m = max(m,rl); rl=0); pd=d; rl++)); max(m,rl); }; \\ Antti Karttunen, Feb 23 2023

Extensions

Data section extended up to a(105) by Antti Karttunen, Feb 23 2023

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

A071222 Smallest k such that gcd(n,k) = gcd(n+1,k+1).

Original entry on oeis.org

1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2, 1, 4, 1, 2, 1, 2
Offset: 0

Views

Author

Benoit Cloitre, Jun 10 2002

Keywords

Comments

a(n) = least m>0 such that gcd(n!+1+m,n-m) = 1. [Clark Kimberling, Jul 21 2012]
From Antti Karttunen, Jan 26 2014: (Start)
a(n-1)+1 = A053669(n) = Smallest k >= 2 coprime to n = Smallest prime not dividing n.
Note that a(n) is equal to A235918(n+1) for the first 209 values of n. The first difference occurs at n=210 and A235921 lists the integers n for which a(n) differs from A235918(n+1).
(End)

Crossrefs

One less than A053669(n+1).

Programs

  • Haskell
    a071222 n = head [k | k <- [1..], gcd (n + 1) (k + 1) == gcd n k]
    -- Reinhard Zumkeller, Oct 01 2014
  • Mathematica
    sgcd[n_]:=Module[{k=1},While[GCD[n,k]!=GCD[n+1,k+1],k++];k]; Array[sgcd,110] (* Harvey P. Dale, Jul 13 2012 *)
  • PARI
    for(n=1,140,s=1; while(gcd(s,n)
    				
  • Scheme
    (define (A071222 n) (let loop ((k 1)) (cond ((= (gcd n k) (gcd (+ n 1) (+ k 1))) k) (else (loop (+ 1 k)))))) ;; Antti Karttunen, Jan 26 2014
    

Formula

Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = A249270 - 1. - Amiram Eldar, Jul 26 2022

Extensions

Added a(0)=1. - N. J. A. Sloane, Jan 19 2014

A328162 Maximum length of a divisibility chain of consecutive divisors of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Oct 06 2019

Keywords

Examples

			The divisors of 968 split into consecutive divisibility chains {{1, 2, 4, 8}, {11, 22, 44, 88}, {121, 242, 484, 968}}, so a(968) = 4.
		

Crossrefs

Records occur at powers of 2 (A000079).
Taking only proper divisors gives A328194.
Taking only divisors > 1 gives A328195.
The maximum run-length among divisors of n is A055874.

Programs

  • Maple
    f:= proc(n) local F,L,d,i;
      F:= sort(convert(numtheory:-divisors(n),list));
      d:= nops(F);
      L:= Vector(d);
      L[1]:= 1;
      for i from 2 to d do
        if F[i] mod F[i-1] = 0 then L[i]:= L[i-1]+1
        else L[i]:= 1
        fi
      od;
      max(L)
    end proc:
    map(f, [$1..100]); # Robert Israel, Apr 20 2023
  • Mathematica
    Table[Max@@Length/@Split[Divisors[n],Divisible[#2,#1]&],{n,100}]

A328195 Maximum length of a divisibility chain of consecutive divisors of n greater than 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Oct 14 2019

Keywords

Comments

Also the maximum length of a divisibility chain of consecutive divisors of n less than n.
The divisors of n (except 1) are row n of A027749.

Examples

			The divisors of 272 greater than 1 are {2, 4, 8, 16, 17, 34, 68, 136, 272}, with divisibility chains {{2, 4, 8, 16}, {17, 34, 68, 136, 272}}, so a(272) = 5.
		

Crossrefs

Allowing 1 as a divisor gives A328162.
Forbidding n as a divisor of n gives A328194.
Positions of 1's are A000040 (primes).
Indices of terms greater than 1 are A002808 (composite numbers).
The maximum run-length of divisors of n is A055874(n).

Programs

  • Mathematica
    Table[If[n==1,0,Max@@Length/@Split[DeleteCases[Divisors[n],1],Divisible[#2,#1]&]],{n,100}]
  • PARI
    A328195(n) = if(1==n, 0, my(divs=divisors(n), rl=0,ml=1); for(i=2,#divs,if(!(divs[i]%divs[i-1]), rl++, ml = max(rl,ml); rl=1)); max(ml,rl)); \\ Antti Karttunen, Dec 07 2024

Extensions

Data section extended up to a(105) by Antti Karttunen, Dec 07 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

A232099 Numbers n such that {largest m such that 1, 2, ..., m divide n} is different from {largest m such that m! divides n^2}.

Original entry on oeis.org

840, 2520, 4200, 5880, 7560, 9240, 10920, 12600, 14280, 15960, 17640, 19320, 21000, 22680, 24360, 26040, 27720, 29400, 31080, 32760, 34440, 36120, 37800, 39480, 41160, 42840, 44520, 46200, 47880, 49560, 51240, 52920, 54600, 55440, 56280, 57960, 59640, 61320, 63000
Offset: 1

Views

Author

Antti Karttunen, Nov 18 2013

Keywords

Comments

Numbers n such that A055874(n) differs from A232098(n). (By the definition of the sequence).
This sequence is a subset of A055926. Please see there for a proof. From that follows that A055881(a(n))+1 is always composite (in range n=1..100000, only values 6, 8, 9 and 10 occur).
Also, incidentally, for the first five terms, n=1..5, a(n) = 70*A055926(n), then a(6)=77*A055926(6), and the next time the ratio A232099(n)/A055926(n) is integral is at n=21, where a(n) = 82*A055926(21), at n=41 (a(41) = 79*A055926(41) = 79*840 = 66360), at n=136, a(136) = 80*A055926(136) = 80*2772 = 221760 and at n=1489, where a(1489) = 80*A055926(1489) = 80 * 30492 = 2439360. The ratio seems to converge towards some value a little less than 80. Please see the plot generated by Plot2 in the links section.

Examples

			840 (= 3*5*7*8) is in the sequence as all natural numbers up to 8 divide 840, but the largest factorial that divides its square, 705600, is 7! (840^2 = 140 * 5040), and 7 differs from 8.
		

Crossrefs

Formula

For all n, a(n) = A055926(A232100(n)). [Follows from the definition of A232100, but cannot as such be used to compute the sequence. Use the given Scheme-program instead.]

A328448 Smallest number whose divisors > 1 have a longest run of length n, and 0 if none exists.

Original entry on oeis.org

2, 6, 12, 504, 60, 420, 840, 4084080, 2520, 21162960, 27720, 2059318800, 0, 360360, 720720, 8494326640800, 12252240, 281206918792800, 0, 0, 232792560, 409547311252279200, 5354228880, 619808900849199341280, 26771144400, 54749786241679275146400, 80313433200, 5663770990518545704800
Offset: 1

Views

Author

Gus Wiseman, Oct 16 2019

Keywords

Examples

			The runs of divisors of 504 (greater than 1) are {{2,3,4},{6,7,8,9},{12},{14},{18},{21},{24},{28},{36},{42},{56},{63},{72},{84},{126},{168},{252},{504}}, the longest of which has length 4, and 504 is the smallest number with this property, so a(4) = 504.
		

Crossrefs

The version that looks at all divisors (including 1) is A328449.
The longest run of divisors of n greater than 1 has length A328457.
Numbers whose divisors > 1 have no non-singleton runs are A088725.
The number of successive pairs of divisors of n is A129308(n).
The Heinz number of the multiset of run-lengths of divisors of n is A328166(n).

Extensions

Data corrected and extended by Giovanni Resta, Oct 18 2019
Previous Showing 21-30 of 43 results. Next