%I #93 Feb 17 2024 02:46:12
%S 1,2,2,4,2,8,2,10,5,11,2,45,2,14,14,36,2,81,2,92,18,20,2,458,7,23,23,
%T 156,2,742,2,202,26,29,26,2234,2,32,30,1370,2,1654,2,337,286,38,2,
%U 9676,9,407,38,454,2,3132,38,3065,42,47,2,73155,2,50,493,1828,44,5257,2,740,50,5066
%N Number of partitions of n into divisors of n.
%C From _Reinhard Zumkeller_, Dec 11 2009: (Start)
%C For odd primes p: a(p^2) = p + 2; for n > 1: a(A001248(n)) = A052147(n);
%C For odd primes p > 3, a(3*p) = 2*p + 4; for n > 2: a(A001748(n)) = A100484(n) + 4. (End)
%C From _Matthew Crawford_, Jan 19 2021: (Start)
%C For a prime p, a(p^3) = (p^3 + p^2 + 2*p + 4)/2;
%C For distinct primes p and q, a(p*q) = (p+1)*(q+1)/2 + 2. (End)
%H Alois P. Heinz, <a href="/A018818/b018818.txt">Table of n, a(n) for n = 1..10000</a> (first 1000 terms from T. D. Noe)
%H Douglas Bowman et al., <a href="https://www.jstor.org/stable/2325075">Problem 6640: Partitions of n into parts which are divisors of n</a>, American Mathematical Monthly, 99(3) (1992), 276-277.
%H Hansraj Gupta, <a href="http://www.dli.gov.in/rawdataupload/upload/insa/INSA_2/20005a84_1276.pdf">Partitions of n into divisors of m</a>, Indian J. Pure Appl. Math., 6(11) (1975), 1276-1286.
%H Hansraj Gupta, <a href="https://insa.nic.in/writereaddata/UpLoadedFiles/IJPAM/20005a84_1276.pdf">Partitions of n into divisors of m</a>, Indian J. Pure Appl. Math., 6(11) (1975), 1276-1286.
%H Martin Klazar, <a href="http://arxiv.org/abs/1808.08449">What is an answer? — remarks, results and problems on PIO formulas in combinatorial enumeration, part I</a>, arXiv:1808.08449 [math.CO], 2018.
%H Noah Lebowitz-Lockard and Joseph Vandehey, <a href="https://arxiv.org/abs/2402.08119">On the number of partitions of a number into distinct divisors</a>, arXiv:2402.08119 [math.NT], 2024. See p. 1.
%H Rémy Sigrist, <a href="/A018818/a018818.png">Colored logarithmic scatterplot of the first 10000 terms</a> (where the color is function of A000005(n)).
%F Coefficient of x^n in the expansion of 1/Product_{d|n} (1-x^d). - _Vladeta Jovovic_, Sep 28 2002
%F a(n) = 2 iff n is prime. - _Juhani Heino_, Aug 27 2009
%F a(n) = f(n,n,1), where f(n,m,k) = f(n,m,k+1) + f(n,m-k,k)*0^(n mod k) if k <= m, otherwise 0^m. - _Reinhard Zumkeller_, Dec 11 2009
%F Paul Erdős, Andrew M. Odlyzko, and the Editors of the AMM give bounds; see Bowman et al. - _Charles R Greathouse IV_, Dec 04 2012
%e The a(6) = 8 representations of 6 are 6 = 3 + 3 = 3 + 2 + 1 = 3 + 1 + 1 + 1 = 2 + 2 + 2 = 2 + 2 + 1 + 1 = 2 + 1 + 1 + 1 + 1 = 1 + 1 + 1 + 1 + 1 + 1.
%p A018818 := proc(n)
%p local a,p,w,el ;
%p a := 0 ;
%p for p in combinat[partition](n) do
%p w := true ;
%p for el in p do
%p if modp(n,el) <> 0 then
%p w := false;
%p break;
%p end if;
%p end do:
%p if w then
%p a := a+1 ;
%p end if;
%p end do:
%p a ;
%p end proc: # _R. J. Mathar_, Mar 30 2017
%t Table[d = Divisors[n]; Coefficient[Series[1/Product[1 - x^d[[i]], {i, Length[d]}], {x, 0, n}], x, n], {n, 100}] (* _T. D. Noe_, Jul 28 2011 *)
%o (Haskell)
%o a018818 n = p (init $ a027750_row n) n + 1 where
%o p _ 0 = 1
%o p [] _ = 0
%o p ks'@(k:ks) m | m < k = 0
%o | otherwise = p ks' (m - k) + p ks m
%o -- _Reinhard Zumkeller_, Apr 02 2012
%o (PARI) a(n)=numbpartUsing(n, divisors(n));
%o numbpartUsing(n, v, mx=#v)=if(n<1, return(n==0)); sum(i=1,mx, numbpartUsing(n-v[i],v,i)) \\ inefficient; _Charles R Greathouse IV_, Jun 21 2017
%o (Magma) [#RestrictedPartitions(n,{d:d in Divisors(n)}): n in [1..100]]; // _Marius A. Burtea_, Jan 02 2019
%Y Cf. A002577, A027750, A033630, A161148 (partitions in squared divisors), A171565, A210442, A211110, A225244,
%Y Cf. A000005, A001248, A001748, A052147, A100484.
%K nonn,nice,look
%O 1,2
%A _David W. Wilson_
|