The well known scaling property for the Dirac $\delta$-function (distribution) states (see the Wikpedia entry)
$$\int_{-\infty}^\infty dx\:\delta(a x) =\frac{1}{|a|}$$
Now we can generally consider $\delta(x)$ as the limiting case of a probability density function
$$p(x,\epsilon)=\frac{1}{\epsilon \sqrt{\pi}} e^{-x^2/\epsilon^2}$$
with
$$\ \int_{-\infty}^\infty dx\:p(x,\epsilon) = 1$$
and
$$\delta(x)=\lim_{\epsilon\to 0}\:p(x,\epsilon)$$
Assume now we spread the events over a different $x$-range by doing the linear scale transformation $x\rightarrow ax$, so $p(x,\epsilon)$ transforms to
$$p(ax,\epsilon)=\frac{1}{\epsilon \sqrt{\pi}} e^{-a^2x^2/\epsilon^2}$$
However, this function is not normalized anymore, as
$$\ \int_{-\infty}^\infty dx\:p(a x,\epsilon) = \frac{1}{|a|}$$
Imagine $p(x)$ describes the way you spend a given amount of money over time, then here, by increasing the length of time ($|a|<1$), one would have miraculously gained money! Obviously, this is not possible. In order to qualify as a probability density function $p(a x)$ must be properly normalized (this is expressed in the 'Location-Scale Family' theorem of probability theory, see also this link) i.e.
$$p(ax,\epsilon)=\frac{|a|}{\epsilon \sqrt{\pi}} e^{-a^2x^2/\epsilon^2}$$
With this normalized function, the increased period of time the money is spent is compensated by the smaller amounts you have to spend, so you correctly end up with the same total spent, i.e.
$$\ \int_{-\infty}^\infty dx\:p(a x,\epsilon) = 1$$
However, if one generates the $\delta$-function from this, one finds
$$\int_{-\infty}^\infty dx\: \delta(a x)=\lim_{\epsilon\to 0}\: \int_{-\infty}^\infty dx\: p(a x,\epsilon) =1 $$
rather than $\int_{-\infty}^\infty dx\: \delta(a x)= \frac{1}{|a|}$ as stated throughout in the literature.
Can anybody see anything wrong with my argument?
EDIT: My argument above used the (by definition) normal property of the generating Gaussian distribution function if one generally interprets the Dirac $\delta$-function as the limit of a probability density function. In this sense it was an additional constraint that enables the use of a continuous mathematical function (in this case the Gaussian function) to approximate the actually discrete probability distribution. However, this can actually also be derived without resorting to probability theory by considering $\delta$ as a distribution that acts as a 'unit impulse' function when applied to a test function $f(x)$ i.e.
$$\int_{-\infty}^\infty dx\: \delta(x)f(x) = f(0) $$
Also with this interpretation, $\delta(x)$ can be interpreted as the limit of $\delta$-sequences, for instance
$$\delta_n(x)=\frac{n}{\sqrt{\pi}} e^{-n^2x^2}$$
with the notation as defined in the link. If we set $n=\frac{a}{\epsilon}$ this becomes instead
$$\delta_{a,\epsilon}(x)=\frac{a}{\epsilon\sqrt{\pi}} e^{-a^2x^2/\epsilon^2}$$
which means
$$\ \int_{-\infty}^\infty dx\:\delta_{a,\epsilon}(x) = 1 \neq \frac{1}{|a|}$$