One of my favourite classical results using $p$-adic methods in elementary number theory is the theorem of Skolem-Mahler-Lech:
This is a theorem about linear recurrence sequences, which are sequences of integers where each term is a fixed linear combination of $n$ previous ones.
So fixing $n$ the sequence $s_i$ is defined by choosing the first $n$ terms
$$s_0,\ldots, s_{n-1}\in \mathbf Z$$
and a relation for all $k$
$$s_{k + n} = \sum_{i=0}^{n-1} a_i s_{k+i}$$
for fixed $a_i$.
Some examples are the Fibonacci sequence ($n = 2$,$s_0 = 0, s_1 = 1$, $a_0=a_1= 1$), and simpler things like any eventually periodic sequence, or the sequence $s_k = k$ (here $n=2$, $s_0 = 0, s_1=1$, $a_0 = -1, a_1= 2$).
We can make other such sequences easily by noting that that the sum of any two linear recurrence sequences is also a linear recurrence sequence.
An important fact about such sequences are that their generating functions
$$f_s = \sum_{k= 0}^\infty s_k x^k$$
are always rational functions of the variable $x$ (one polynomial divided by another), where the numerator defines the initial terms $s_0, \ldots, s_{n-1}$ and the denominator defines the recurrence relation.
Of the examples I mentioned above, the fibonacci sequence grows (exponentially), eventually periodic sequences are bounded, and the sequence $s_k=k$ also grows, just less quickly than fibonacci.
One question that one might then ask is:
What is the set of $k$ for which $s_k = 0$?
from these examples (and others) we might conjecture that this set is periodic, except for finitely many exceptions (after all we can always change finitely many terms of any linear recurrence sequence to make a sequence with the same behaviour eventually but with zeroes wherever we want at the start).
How might one go about proving this?
The first step of the proof is to use the rational generating function $f_s$ and write out its partial fraction decomposition over an algebraically closed field (like $\overline {\mathbf Q}$), this will be of the form
$$f_s = \sum_{i=1}^{\ell} \sum_{j = 1}^{r_i} \frac{\alpha_{ij}}{(x - \beta_{i})^j} $$
for some fixed roots $\beta_j$ of the original denominator of $f_s$.
Now using this decomposition
we have
$$f_s = \sum_{i=1}^{\ell} \sum_{j = 1}^{r_i} \alpha_{ij}{\left(\sum_{n=0}^\infty \beta_i^n x^n\right)^j} $$
what this gives is that
$$s_n = \text{some polynomial expression involving terms }\beta^n $$
For example for the fibonacci sequence this recovers Binet's formula
$$s_n = \frac{1}{\sqrt 5} \left(\frac{1+ \sqrt 5}2\right)^n-\frac{1}{\sqrt 5} \left(\frac{1- \sqrt 5}2\right)^n.$$
Or for the periodic sequence $0,1,0,1,0,1,\ldots$ this is
$$ s_n = 1^ n - (-1)^n$$
So we've written $s_n$ as a sum of exponential-type functions in $n$ with different bases, that we want to describe the zeroes of this function for $n \in \mathbf N$.
Now the magic part: the function $e^x$ is an analytic function, and on a bounded domain analytic functions have only finitely many zeroes (unless they are zero everywhere). This would give us a lot of control over the zeroes of $s_n$ if the naturals were bounded.
Which leads to the slightly strange question:
What if the natural numbers were bounded? And the functions $\beta^n$ were still analytic in some way?
Of course using the usual absolute value and metric on $\mathbf Q$ and $\mathbf C$ this is totally false.
But in the $p$-adic numbers this is true! The integers are all bounded ($p$-adically) by norm $\le 1$. So lets treat these functions as $p$-adic functions and control the zero sets in some way.
How does this prove the result? The functions $\beta^n$ are not $p$-adic analytic functions of $n$ on their own, but they are so on small enough $p$-adic disks though, what ends up happening is we get distinction between congruence classes of $n$ mod $p-1$ for some well-chosen $p$ such that in each congruence class there are only finitely many zeroes of $s_n$ or the function $s_n$ is identically zero on that congruence class.
This gives us the theorem mentioned above, that the zeroes of $s_n$ are periodic, except for finitely many exceptions.