This question is related to Gromov's theorem in that it is about the growth rate with respect to some set of generators of $\mathbb Z$. Gromov's Theorem, along with the Bass-Guivarc'h formula, tells us that for any finite symmetric set of generators $S$ of $\mathbb Z$, if we let $$B_n = \{x \in \mathbb Z : \exists k \leq n, s_1, \ldots, s_k \in S \text{ such that }n = s_1 + \cdots + s_k\}$$ and $c(n) = \left|B_n\right|$, then $c(n) \in \Theta(n)$ since $\operatorname{rank} \mathbb Z = 1$.
Computational evidence leads me to expect that there is a stronger condition on $c(n)$: I conjecture that $c(n)$ is necessarily linear for large enough $n$, that is, that there are some integers $N, a, b$ such that for all $n \geq N$, we have $c(n) = an + b$.
Is this conjecture true or not? Can anyone think of a proof either way?
EDIT: Computing the growth rates induced from various generating sets makes me think that in fact we have $c(n)$ linear for all $n$ such that $\{-\max S, \ldots, -1, 0, 1, \ldots, \max S\} \subseteq B_{n-1}$, or equivalently for all $n > \max d^{-1}[\{-\max S, \ldots, 0, \ldots, \max S\}]$, where $d : \mathbb Z \to \mathbb N$ is given by $d(k) = \min\{n \in \mathbb N : k \in B_n\}$; this gives a very precise bound on where the linearity must start, and in most cases this bound is actually sharp.
EDIT 2: Allow me to describe some of my observations. The Mathematica code I used is provided below.
$S = \pm \{4, 7\}$. Here are some plots which tell about the behavior of the $B_n$ and of $c(n)$:
The upper plot, with all the dots, shows the nonnegative elements of each $B_n$ for $n = 0, 1, \ldots, 10$. ($B_n$ is always symmetric, so we do not need to show both sides.) Note that $\max S = 7$, and observe from this plot that $B_n$ contains $\{-7, -6, \ldots, 6, 7\}$ for $n \geq 5$. Now, the lower-left plot is just $c(n)$ vs. $n$, and the lower-right plot is of the differences $c(n) - c(n - 1)$. Gromov's Theorem says that $c(n) \in \Theta(n)$, and indeed the lower-left plot shows that $c(n)$ looks linear. The lower-right plot strongly suggests that $c(n)$ is exactly linear for $n > 5$, with formula determined to be $c(n) = 14n - 9$. Looking back at the various $B_n$, one may observe that the right "edge" of each $B_n$ is exactly the same for $n \geq 5$, which is itself close to a proof that $c(n)$ is linear past this point.$S = \pm \{1, 2, 9\}$. Here are the plots:
Here again my conjectures seem to be confirmed: the lower-right plot suggests that $c(n)$ is eventually linear with slope 18, and furthermore $\max S = 9$, and we see from the upper plot that for all $n \geq 3$, $\{-9, -8, \ldots, 8, 9\} \subseteq B_n$, and from the lower-right plot again that $c(n)$ is linear on $n \geq 3$. (It in fact seems that the bound on where linearity begins is sharp exactly when $1 \not\in S$, and needs to be shifted left by one when $1 \in S$.)
Here is the code. gens specifies the positive portion of the generators you wish to use (e.g. if gens = {2, 3}, then $S = \{-3, -2, 2, 3\}$), while iters and countsiters just say how much data to show: iters is the number of balls $B_n$ to display in the first plot, and countsiters is the number of ball-sizes $c(n)$ to display in the lower plots.
(* input parameters *)
gens := {6, 7}
iters := 10
countsiters := 20
(* computation *)
MinkowskiSum[lst__] := DeleteDuplicates[Total /@ Tuples[{lst}]]
MinkowskiProduct[lst__] :=
DeleteDuplicates[Apply[#1 #2 &, Tuples[{lst}], {1}]]
Clear[b]
allgens := MinkowskiProduct[gens, {-1, 0, 1}]
b[-1] := {}
b[0] := {0}
b[n_] := b[n] = MinkowskiSum[b[n - 1], allgens]
c[n_] := Length[b[n]]
(* plots *)
ListPlot[Table[{#, -n} & /@ Select[b[n], # >= 0 &], {n, 0, iters}],
AspectRatio -> 3*iters/(Max[b[iters]] - Min[b[iters]]),
GridLines -> {Range[0, iters Max[allgens]], Range[-iters, 0]},
Ticks -> {Range[0, iters Max[allgens], 5], Range[-iters, 0]},
ImageSize -> Full]
{DiscretePlot[c[n], {n, 0, countsiters}, ImageSize -> Large],
DiscretePlot[c[n] - c[n - 1], {n, 0, countsiters},
ImageSize -> Large]}