login

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”).

A341150
Triangle read by rows, lower triangle of the all-zeros matrix modified step by step with T(0,0) = 1, then at each step, if T(n,k) has changed in previous step, increase T(n+k,k) and T(3n+1-2k,2n+1-k) by 1.
1
1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 1, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1
OFFSET
0
LINKS
Niccolò Castronuovo, On the number of fixed points of the map gamma, arXiv:2102.02739 [math.NT], 2021.
EXAMPLE
Triangle begins:
[1];
[0, 1];
[0, 1, 1];
[0, 1, 0, 1];
[0, 1, 1, 0, 1];
[0, 1, 0, 0, 1, 1];
[0, 1, 1, 1, 0, 0, 1].
PROG
(PARI) incr(m, nn, vincr) = if ((vincr[1]<nn) && (vincr[2]<nn), m[vincr[1]+1, vincr[2]+1]++); m;
tabl(nn) = {my(m = matrix(nn, nn), istep=0, vstep, nvstep); for (istep=0, nn, if (istep==0, vstep = [[0, 0]], vstep = nvstep); nvstep = List(); for (i=1, #vstep, my(step=vstep[i]); m = incr(m, nn, step); if (istep, listput(nvstep, [step[1]+step[2], step[2]])); listput(nvstep, [3*step[1]+1-2*step[2], 2*step[1]+1-step[2]]); ); ); m; }
out(nn) = my(m=tabl(nn)); for(i=1, nn, print(vector(i, k, m[i, k])); );
CROSSREFS
Row sums give A341151.
Sequence in context: A287722 A284588 A113704 * A244220 A283963 A131670
KEYWORD
nonn,tabl
AUTHOR
Michel Marcus, Feb 06 2021
STATUS
approved