A154789 Where records occurs in A085543 for positive values of n.
1, 2, 7, 9, 15, 19, 21, 39, 59, 66, 75, 96, 124, 125, 180, 327, 358, 434, 699, 999, 1685, 1779, 3561, 4184, 4835
Offset: 1
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.
DeleteDuplicates[DivisorSigma[0,PartitionsP[Range[0,2000]]],GreaterEqual] (* The program generates the first 22 terms of the sequence. *) (* Harvey P. Dale, Dec 09 2024 *)
For n = 20 the divisors of 20 are 1, 2, 4, 5, 10, 20 and three of them are also partition numbers: 1, 2, 5, so a(20) = 3. For n = 42 the divisors of 42 are 1, 2, 3, 6, 7, 14, 21, 42 and five of them are also partition numbers: 1, 2, 3, 7, 42, so a(42) = 5.
p = {1}; Table[If[n >= Last@p, AppendTo[p, PartitionsP[1 + Length@p]]]; Length@Select[p, Mod[n, #] == 0 &], {n, 90}] (* Giovanni Resta, Jan 22 2014 *)
For n = 12 the divisors of 12 are 1, 2, 3, 4, 6, 12. There are 6 divisors of 12 and the number of partitions of 6 is A000041(6) = 11, so a(12) = 11.
Table[PartitionsP@ DivisorSigma[0, n], {n, 120}] (* Michael De Vlieger, Apr 25 2016 *)
a(n) = numbpart(numdiv(n)); \\ Michel Marcus, Apr 26 2016
DivisorSigma[0, PartitionsQ@ Range[0, 100]]
a(n) = numdiv(polcoeff(prod(k=1, n, 1 + x^k, 1 + x * O(x^n)), n)); \\ Michel Marcus, Jun 19 2015
21 belongs to this sequence as the partition number of 21 is 792, which has 22 nontrivial divisors (divisors different from 1 and itself), and 22 > 21.
Select[Range[1000], DivisorSigma[0, PartitionsP[#]] > # + 2 &] (* or *) First/@ Select[ Import[ "https://oeis.org/A085543/b085543.txt", "Table"], #[[2]] > #[[1]] + 2&] (* Giovanni Resta, Jul 16 2018 *)
select(n->numdiv(numbpart(n))-2 > n, [1..1000]) \\ Andrew Howroyd, Jul 15 2018
r="" for i in range(1000): l=Partitions(i+1).cardinality().divisors() if ((len(l)-2)>i+1): r=r+str(i+1)+"," print(r)
Comments