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.

A121572 Subprimorials: inverse binomial transform of primorials (A002110).

This page as a plain text file.
%I A121572 #23 Mar 13 2017 10:42:09
%S A121572 1,1,3,17,119,1509,18799,342397,6340263,151918421,4619754311,
%T A121572 140219120601,5396354613583,221721908976697,9431597787000999,
%U A121572 447473598316521449,24163152239530299719,1444153946379288324477,87200644323074509092943,5929294512595059362045041
%N A121572 Subprimorials: inverse binomial transform of primorials (A002110).
%C A121572 By analogy with subfactorials, which are the inverse binomial transform of the factorials.
%H A121572 Alois P. Heinz, <a href="/A121572/b121572.txt">Table of n, a(n) for n = 0..350</a>
%F A121572 a(n) = Sum_{k=0..n} (-1)^(n-k) C(n,k) Prime(k)#, where p# is p primorial and Prime(0)# = 1.
%F A121572 A007318^(-1) * A002110. - _Gary W. Adamson_, Dec 14 2007
%e A121572 a(3) = 30 - 3*6 + 3*2 - 1 = 17.
%p A121572 b:= proc(n) option remember; `if`(n=0, 1, ithprime(n)*b(n-1)) end:
%p A121572 a:= n-> add(binomial(n,k)*b(k)*(-1)^(n-k), k=0..n):
%p A121572 seq(a(n), n=0..20);  # _Alois P. Heinz_, Sep 19 2016
%t A121572 b[n_] := b[n] = If[n==0, 1, Prime[n]*b[n-1]]; a[n_] := Sum[Binomial[n, k]* b[k]*(-1)^(n-k), {k, 0, n}]; Table[a[n], {n, 0, 20}] (* _Jean-François Alcover_, Mar 13 2017, after _Alois P. Heinz_ *)
%Y A121572 Cf. A002110, A000166, A136104.
%Y A121572 See A079266 for a different definition of subprimorial.
%K A121572 nonn
%O A121572 0,3
%A A121572 _Franklin T. Adams-Watters_, Aug 08 2006
%E A121572 More terms from _R. J. Mathar_, Sep 18 2007
%E A121572 Edited by _N. J. A. Sloane_, May 15 2008 at the suggestion of _R. J. Mathar_