login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A362451 Gilbreath transform of {sigma(i), i >= 1} (cf. A000203). 14
1, 2, 1, 1, 1, 2, 1, 0, 0, 1, 0, 4, 0, 3, 0, 2, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 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, 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, 68, 0, 14, 0, 7, 0, 2, 0, 21, 1, 8, 1, 9, 1, 0, 1, 18, 0, 7, 0, 2, 0, 1, 0, 13, 1, 1, 1, 2, 1, 1 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
Given a sequence {u(i), i >= o} with offset o, its absolute difference sequence is the sequence {v(i) = |u(i+1)-u(i)|, i >= o}.
The Gilbreath transform of a sequence s = {s(i), i >= o} is constructed as follows.
Form an array A in which the initial row is s and each subsequence row is the absolute difference sequence of the previous row. The sequence of leading terms of the rows of A is the Gilbreath transform of s.
If "absolute difference sequence" is changed to the familiar "first differences", instead of the Gilbreath transform we get the usual inverse binomial transform.
It appears that the terms are mostly 0's and 1's, with occasional eruptions of "geysers". See A362456, A362457.
LINKS
N. J. A. Sloane, Transforms
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: Video, Slides, Updates. (Mentions this sequence.)
EXAMPLE
We give two examples. (1) For the Gilbreath transform of the sequence of primes (cf. A000040), the array A is given in A036262. The Gilbreath transform begins {2, 1, 1, 1, 1, ...}, and the famous Gilbreath conjecture is that every term after the initial 2 is equal to 1.
(2) For the Gilbreath transform of {tau(i), i >= 1} (cf. A000005), the array A is given in A362450, and the Gilbreath transform is given in A361897. The authors of the latter sequence conjecture that its terms are just 0's and 1's.
See A362452 for a further example.
MAPLE
# To get M terms of the Gilbreath transform of s, assuming offset is 1:
GT := proc(s, M) local G, u, i;
u := [seq(s[i], i=1..M)];
G:=[s[1]];
for i from 1 to M-1 do
u:=[seq(abs(u[i+1]-u[i]), i=1..nops(u)-1)];
G:=[op(G), u[1]]; od:
G;
end;
# For the present sequence:
GT(numtheory[sigma], 150);
# See link for a more comprehensive Maple program
MATHEMATICA
A362451[nmax_]:=Module[{d=DivisorSigma[1, Range[nmax]]}, Join[{1}, Table[First[d=Abs[Differences[d]]], nmax-1]]]; A362451[200] (* Paolo Xausa, May 07 2023 *)
PROG
(PARI)
lista(nn) = my(v=apply(sigma, [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);
lista(200) \\ (based on PARI program in A361897
CROSSREFS
Sequence in context: A281185 A260683 A337683 * A092673 A243842 A367405
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, May 03 2023
EXTENSIONS
More than the usual number of terms are displayed in order to go out beyond the long initial 0,1 subsequence.
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 05:02 EDT 2024. Contains 371235 sequences. (Running on oeis4.)