login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A334370 Expansion of e.g.f. Product_{k>=1} (1 + x^prime(k) / prime(k)!). 5

%I #26 Jul 20 2023 10:38:23

%S 1,0,1,1,0,11,0,22,56,36,2640,1,8712,79,72436,360465,48608,49008961,

%T 794376,4232764,7753140,942565890,18198334,14799637777,10577976,

%U 366619314900,2785137222400,1475339135400,1065920156634060,3765722000041,5869315258699050

%N Expansion of e.g.f. Product_{k>=1} (1 + x^prime(k) / prime(k)!).

%C a(n) is the number of functions f:[n]-> [n] such that the number of elements that are mapped to i is either 0 or the i-th prime. a(5) = 11: (33333), (11222), (12122), (12212), (12221), (21122), (21212), (21221), (22112), (22121), (22211). - _Alois P. Heinz_, Jul 18 2023

%H Alois P. Heinz, <a href="/A334370/b334370.txt">Table of n, a(n) for n = 0..836</a>

%p b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0, b(n, i-1)+

%p (p-> `if`(p>n, 0, b(n-p, i-1)*binomial(n, p)))(ithprime(i))))

%p end:

%p a:= n-> b(n, numtheory[pi](n)):

%p seq(a(n), n=0..30); # _Alois P. Heinz_, Jul 18 2023

%t nmax = 30; CoefficientList[Series[Product[(1 + x^Prime[k]/Prime[k]!), {k, 1, nmax}], {x, 0, nmax}], x] Range[0, nmax]!

%t a[n_] := a[n] = If[n == 0, 1, (n - 1)! Sum[DivisorSum[k, -#/(-#!)^(k/#) &, PrimeQ[#] &] a[n - k]/(n - k)!, {k, 1, n}]]; Table[a[n], {n, 0, 30}]

%o (PARI) my(N=40, x='x+O('x^N)); Vec(serlaplace(prod(k=1, N, 1+isprime(k)*x^k/k!))) \\ _Seiichi Manyama_, Feb 27 2022

%Y Cf. A000040, A000720, A007837, A032310, A115278, A178682, A190476, A319113, A364344.

%K nonn

%O 0,6

%A _Ilya Gutkovskiy_, May 11 2020

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 01:36 EDT 2024. Contains 371264 sequences. (Running on oeis4.)