cp's OEIS Frontend

A362452 Gilbreath transform of {sigma(i)-i, i >= 1} (see sum of aliquot parts, A001065).

%I A362452 #47 Sep 27 2023 15:04:18
%S A362452 0,1,1,1,1,0,1,0,0,1,0,0,0,1,0,0,1,1,1,0,1,1,1,1,1,1,1,0,1,1,1,1,1,0,
%T A362452 1,0,1,0,1,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,0,0,1,0,1,0,1,0,1,1,0,1,1,
%U A362452 1,0,1,1,0,1,0,0,0,1,0,1,0,0,0,1,0,0,0,1,0,0,0,1,0,0,0,1,1,1,1,0,1,1,1,0,0,0,0,1,0,0,0,0,1,1,1,0,1,1,1,62,0,12,0,3,0,2,0,25,1
%N A362452 Gilbreath transform of {sigma(i)-i, i >= 1} (see sum of aliquot parts, A001065).
%C A362452 See A362451 for further information.
%C A362452 The first 50000 terms of the present sequence suggest that the terms are usually 0's and 1's, except for occasional "geysers". See A362458, A362459.
%C A362452 [It would be nice to have plots of larger numbers of initial terms.]
%H A362452 N. J. A. Sloane, <a href="/A362452/b362452.txt">Table of n, a(n) for n = 1..50000</a>
%H A362452 N. J. A. Sloane, New Gilbreath Conjectures, Sum and Erase, Dissecting Polygons, and Other New Sequences, Doron Zeilberger's Exper. Math. Seminar, Rutgers, Sep 14 2023: <a href="https://vimeo.com/866583736?share=copy">Video</a>, <a href="http://neilsloane.com/doc/EMSep2023.pdf">Slides</a>, <a href="http://neilsloane.com/doc/EMSep2023.Updates.txt">Updates</a>. (Mentions this sequence.)
%H A362452 Paolo Xausa, <a href="/A362452/a362452.txt">Table of n, a(n) for n = 1..1000000</a>
%H A362452 Paolo Xausa, <a href="/A362452/a362452.png">Logarithmic scatterplot for n = 1..1000000</a>
%H A362452 <a href="/index/Ge#Gilbreath">Index entries for sequences related to Gilbreath conjecture and transform</a>
%p A362452 # To get M terms of the Gilbreath transform of s:
%p A362452 GT := proc(s,M) local G,u,i;
%p A362452 u := [seq(s(i),i=1..M)];
%p A362452 G:=[s(1)];
%p A362452 for i from 1 to M-1 do
%p A362452 u:=[seq(abs(u[i+1]-u[i]),i=1..nops(u)-1)];
%p A362452 G:=[op(G),u[1]]; od:
%p A362452 G;
%p A362452 end;
%p A362452 # For the present sequence:
%p A362452 aliq := proc(n) numtheory[sigma](n) - n; end;
%p A362452 GT(aliq,150);
%t A362452 A362452[nmax_]:=Module[{d=DivisorSigma[1,Range[nmax]]-Range[nmax]},Join[{0},Table[First[d=Abs[Differences[d]]],nmax-1]]];A362452[200] (* _Paolo Xausa_, May 07 2023 *)
%o A362452 (PARI)
%o A362452 f(n) = sigma(n) - n
%o A362452 lista(nn) = my(v=apply(f, [1..nn]), list = List(), nb=nn); listput(list, v[1]); for (n=2, nn, nb--; my(w = vector(nb, k, abs(v[k+1]-v[k]))); listput(list, w[1]); v = w; ); Vec(list);
%o A362452 lista(200)
%Y A362452 Cf. A000005, A000040, A000203, A001065, A036262, A361897, A362450, A362451, A362458, A362459.
%K A362452 nonn
%O A362452 1,120
%A A362452 _N. J. A. Sloane_, May 03 2023
%E A362452 More than the usual number of terms are displayed in order to go out beyond the long initial 0,1 subsequence.