Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%I #20 Dec 21 2021 21:11:55
%S 2,9,22,48,85,151,231,355,500,709,937,1267,1617,2069,2575,3193,3860,
%T 4686,5549,6593,7725,8985,10337,11961,13591,15464,17498,19714,22036,
%U 24690,27378,30382,33603,37023,40597,44733,48720,53152,57950,62978,68074,73898,79558
%N Convolution of sequence of primes with sequence sigma(n).
%C From _Omar E. Pol_, Dec 06 2021: (Start)
%C Antidiagonal sums of A272214.
%C Convolution of A000040 and A000203.
%C Convolution of A054541 and A024916.
%C Convolution of the nonzero terms of A007504 and A340793.
%C a(n) is also the volume of a tower or polycube in which the successive terraces are the symmetric representation of sigma(k), k = 1..n starting from the top, and the successive heights of the terraces are the prime numbers starting from the base. (End)
%H Robert Israel, <a href="/A086718/b086718.txt">Table of n, a(n) for n = 1..10000</a>
%p N:= 100: # to get a(1)..a(N)
%p P:= [seq(ithprime(i),i=1..N+1)]:
%p S:= [seq(numtheory:-sigma(i),i=1..N+1)]:
%p seq(add(P[i]*S[n-i],i=1..n-1),n=2..N+1); # _Robert Israel_, Sep 09 2020
%o (PARI) p=primes(30); s=vector(30,i, sigma(i)); conv(u,v)=local(w); w=vector(length(u),i,sum(j=1,i,u[j]*v[i+1-j])); w;
%o conv(p,s)
%Y Cf. A000040, A000203.
%Y Cf. A007504, A024916, A054541, A066186, A175254, A191831, A237593, A272214, A340793.
%K nonn
%O 1,1
%A _Jon Perry_, Jul 29 2003
%E More terms from _Robert Israel_, Sep 09 2020