cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A352226 Consider a 2D sandpile model where each site with 2 or more grains, say at location (x, y), topples and transfers one grain of sand to the sites at locations (x+1, y) and (x, y+1). Let S(n) be the configuration after stabilization of a configuration with n grains at the origin. a(n) = Max_{ (x,y) in S(n) } (x+y).

This page as a plain text file.
%I A352226 #16 Mar 11 2022 07:42:51
%S A352226 0,1,1,3,3,3,3,5,5,5,5,7,7,7,7,9,9,9,9,9,9,9,9,11,11,11,11,13,13,13,
%T A352226 13,15,15,15,15,15,15,15,15,15,15,15,15,17,17,17,17,19,19,19,19,19,19,
%U A352226 19,19,19,19,19,19,21,21,21,21,23,23,23,23,23,23,23,23
%N A352226 Consider a 2D sandpile model where each site with 2 or more grains, say at location (x, y), topples and transfers one grain of sand to the sites at locations (x+1, y) and (x, y+1). Let S(n) be the configuration after stabilization of a configuration with n grains at the origin. a(n) = Max_{ (x,y) in S(n) } (x+y).
%C A352226 Sites containing 0 or 1 grain are stable. S(n) contains n elements.
%H A352226 Rémy Sigrist, <a href="/A352226/b352226.txt">Table of n, a(n) for n = 1..10000</a>
%H A352226 Rémy Sigrist, <a href="/A352226/a352226.png">Representation of the configuration for n = 100000</a>
%H A352226 Wikipedia, <a href="https://en.wikipedia.org/wiki/Abelian_sandpile_model#Sandpile_models_on_directed_graphs">Sandpile models on directed graphs</a>
%e A352226 For n = 15:
%e A352226 - S(15) corresponds to the following configuration:
%e A352226     4|  X X X
%e A352226     3|X   X   X
%e A352226     2|X     X X
%e A352226     1|X       X
%e A352226     0|X X X X
%e A352226      +---------
%e A352226       0 1 2 3 4
%e A352226 - x+y is maximized for (x,y) = (4,3) and (3,4),
%e A352226 - so a(15) = 3+4 = 7.
%o A352226 (PARI) a(n) = { my (s=[n]); for (k=-1, oo, if (vecmax(s)==0, return (k), s \= 2; s = concat(0, s) + concat(s, 0); if (#s>2 && s[1]==0, s = s[2..#s-1]))) }
%Y A352226 Cf. A180230, A349990, A350188, A351783.
%K A352226 nonn
%O A352226 1,4
%A A352226 _Rémy Sigrist_, Mar 08 2022