|
%I #6 Jul 21 2013 03:19:59
%S 3,4,5,7,12,13,16,18,19,22,23,30,31,38,39,40,42,43,46,48,49,50,51,52,
%T 53,56,57,58,60,61,68,69,70,72,73,80,81,82,84,85,86,87,88,89,92,93,94,
%U 95,96,98,99,100,102,103,106,108,110,111,112,113,121,122,123,124,125,126
%N "Difference" of the sequences of triangular numbers (A000217) and the primes (cf. A092858).
%C Here the complement of the sequence of primes (1 and the composites) is "added" to the sequence of triangulars, according to the definition outlined in A092858.
%H Ferenc Adorjan, <a href="http://web.axelero.hu/fadorjan/aronsf.pdf">Binary mapping of monotonic sequences and the Aronson function</a>
%o (PARI) {sdif(a,b)= /*Returns the "difference" of monotonic sequences a and b */ return(mtinv(mt(a)+mt(compl(b)))) /* the functions mt(a) and mtinv(r) are defined in A051006 and A092855, respectively */ } {compl(v)=/* Returns the complement of v monotonic positive sequence */ local(n,p=0,vv=[]);n=matsize(v)[2];for(i=1,n, for(j=p+1,v[i]-1,vv=concat (vv,j));p=v[i]);return(vv)}
%Y Cf. A092855, A051006, A092857, A092858, A092860, A092861, A092862, A092863, A092874.
%K easy,nonn
%O 1,1
%A Ferenc Adorjan (fadorjan(AT)freemail.hu)
|
