|
%I M0065 #82 Jan 09 2023 01:31:13
%S 2,1,1,-1,3,-9,23,-53,115,-237,457,-801,1213,-1389,445,3667,-15081,
%T 41335,-95059,195769,-370803,652463,-1063359,1570205,-1961755,1560269,
%U 1401991,-11023119,36000427,-93408425,214275735,-450374071
%N Inverse binomial transform of primes.
%C a(n) is the (n-1)-st difference of the first n primes. Although the magnitude of the terms appears to grow exponentially, a plot shows that the sequence a(n)/2^n has quite a bit of structure. See A082594 for an interesting application. - _T. D. Noe_, May 09 2003
%C Graph this divided by A122803 using plot2! - _Franklin T. Adams-Watters_
%C From _Robert G. Wilson v_, Jan 28 2020: (Start)
%C a(n) is odd for all n>1 and a(2n) is negative for all n>1.
%C As opposed to A331573, there are terms where abs(a(n)) >= abs(a(n+1)). (End)
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%H Robert G. Wilson v, <a href="/A007442/b007442.txt">Table of n, a(n) for n = 1..3321</a> (first 1000 from Franklin T. Adams-Watters)
%H Vaclav Kotesovec, <a href="/A007442/a007442.jpg">Plot of |a(n)|^(1/n) for n = 1..10000</a>
%H T. D. Noe, <a href="http://www.sspectra.com/math/A007442.gif">Plot of A007442</a>
%H N. J. A. Sloane, <a href="/transforms.txt">Transforms</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/BinomialTransform.html">Binomial Transform</a>
%F a(n) = Sum_{k=0..n-1} (-1)^(n-k-1) binomial(n-1, k) prime(k+1).
%F a(n) = A095195(n,n-1). - _Alois P. Heinz_, Sep 25 2013
%F G.f.: Sum_{k>=1} prime(k)*x^k/(1 + x)^k. - _Ilya Gutkovskiy_, Apr 23 2019
%e a(4) = 7 - 3*5 + 3*3 - 2 = -1.
%t Diff[lst_List] := Table[lst[[i + 1]] - lst[[i]], {i, Length[lst] - 1}]; n=1000; dt = Prime[Range[n]]; a = Range[n]; a[[1]] = 2; Do[dt = Diff[dt]; a[[i]] = dt[[1]], {i, 2, n}]; a
%t u = Table[Prime[Range[k]], {k, 1, 100}];Flatten[Table[Differences[u[[k]], k - 1], {k, 1, 100}]] (* _Clark Kimberling_, May 15 2015 *)
%t t = Array[Prime, 30]; f[x_] := Rest[x] - Most[x];
%t Flatten[Last /@ (NestList[f, t[[1 ;; #]], (# - 1)] & /@ Range[1, 29])] (* _Horst H. Manninger_, Mar 22, 2021 *)
%o (PARI) vector(50, n, sum(k=0, n-1,(-1)^(n-k-1)*binomial(n-1, k)*prime(k+1))) \\ _Altug Alkan_, Oct 17 2015
%Y Cf. A007443, A082594, A095195.
%K sign,easy
%O 1,1
%A _N. J. A. Sloane_
|
