First of all, you cannot conclude that an operator is compact if it is the strong limit of finite rank operators, for example in any Hilbert space the identity operator is the strong limit of a net of finite-rank projections, but of course it is compact if and only if the dimension of our space is finite.
As Timur Bakiev says, the multiplication operator $T:\ell^p\to\ell^p$ is compact if and only if the sequence of $\ell^\infty$ that acts on $\ell^p$ belongs to $c_0$.
You already showed the converse: In your notation, $$\|(T_n-T)x\|^p=\sum_{k=n+1}^\infty|a_k|^p\cdot|x_k|^p\leq\sup_{k\geq n+1}|a_k|^p\cdot\|x\|^p$$
so $\|T_n-T\|\leq\sup_{k\geq n+1}|a_k|\to0$.
Conversely, suppose that $T$ is compact. We will use the following lemma from real analysis: if $(t_n)$ is a sequence of real numbers and every subsequence of $(t_n)$ has a further subsequence converging to $t$, then $(t_n)$ converges to $t$. See this post for example: Every subsequence of $x_n$ has a further subsequence which converges to $x$. Then the sequence $x_n$ converges to $x$.
By definition of compactness, if $(x_n)\subset\ell^p$ is a bounded sequence, $(Tx_n)$ has a convergent subsequence in $\ell^p$. Fix indexes $n_1<n_2<n_3<\dots$ and consider the elements $e_k\in\ell^p$ having $0$ everywhere and $1$ in the $n_k$ position. Then $\|e_k\|=1$, so this sequence is bounded and therefore $(Te_k)$ has a convergent subsequence in $\ell^p$. To simplify my notation, I will assume that $n_1=1, n_2=2,\dots$ but the method below can be applied to arbitrary subsequences too, you can verify this yourself easily.
So as I said, let $e_n$ be the sequence having $0$ everywhere and $1$ in the $n$-th position, so $\|e_n\|=1$ for all $n$. Note that $Te_n$ is the sequence that has $0$ everywhere and $a_n$ in the $n$-th position, so $\|Te_n\|_p=|a_n|$. Now $(Te_n)$ has a convergent subsequence, say $(Te_{n_k})$ and say that $Te_{n_k}\to y\in\ell^p$, thus $\|Te_{n_k}-y\|_p^p\to0$. If $y=(y_i)$, then
$$\|Te_{n_k}-y\|_p^p=\sum_{i=1}^{n_k-1}|y_i|^p+|a_{n_k}-y_{n_k}|^p+\sum_{i=n_k+1}^\infty|y_i|^p$$
Fix an integer $m$. For $\varepsilon>0$ choose $k>m$ so large so that $\|Te_{n_k}-y\|_p^p<\varepsilon$.
Then $|y_m|^p\leq\sum_{i\neq n_k}|y_i|^p+|a_{n_k}-y_{n_k}|^p=\|Te_{n_k}-y\|_p^p<\varepsilon$ and as $\varepsilon>0$ was arbitrary, we conclude that $y_m=0$, so $y=0$ since $m$ was arbitrary before fixing. We conclude that $|a_{n_k}|\to0$.
Now as I said, this can be done for any subsequence (this is why I fixed $n_1<n_2<\dots$ in the beginning, nevermind that I eventually went wLOG to the case $n_1=1,n_2=2\dots$), so the lemma will yield that $|a_n|\to0$, i.e. $a_n\in c_0$.