Does there exist $( x_n)_{n\in\mathbb{N}}\subset \mathbb{R}$ such that $\displaystyle\sum_{k=1}^{\infty} x_{kn} =1$ for every $n\in\mathbb{N}?$
I suspect so, but I am struggling to actually construct such a sequence. All I can say so far is that there must be infinitely many members of $( x_n)_{n\in\mathbb{N}}$ that are positive and infinitely many members of $( x_n)_{n\in\mathbb{N}}$ that are negative. Other than this, I have not made much progress.
Can someone sketch a construction or provide any hints of a disproof please?
A partial solution that may possibly get to a full solution with care (or of course may not, though I suspect the problem is solvable as per below)
One can construct an Abel summable sequence with the above properties (also in the Abel summability sense and even a bit more general) and maybe with care one can actually get a convergent sequence, st all the above series are convergent too by getting a sequence as below that satisfies Littlewood's Tauberian condition $n|x_n| \le C$; this, of course, implies $n|x_{nk}| \le C/k$ so all the subsequences $x_{nk}, n \ge 1, k$ fixed satisfy the Tauberian property so all the OP series converge (to $1$) if they are Abel summable (to $1$).
Let $S$ the set of all roots of unity which is countable, so there exists a non-negative integrable function on the unit circle $f \ge 0$ st if $f(re^{it})$ is its (nonnegative) harmonic extension to the unit disc (Poisson integral of $f$ or simply $f(re^{it})=\sum_{n \in \mathbb Z} a_nr^{|n|}e^{int}$ where $\sum a_ne^{int}$ is the Fourier series of $f$), then if $|w_k| \le 1, w_k \to w \in S$ we have $f(w_k) \to \infty$
This a general theorem valid for all sets of measure $0$ on the circle and is easily proven directly for $S$ by finding for every $m \ge 1$ an open circle set $G_m$ containing $S$ of circle measure at most $1/m^4$ and then define $f_m(w)=m^2, w \in G_m, f_m(w)=0$ otherwise and take $f=\sum f_m$.
Now considering $g$ the harmonic conjugate of $f$ in the unit disc and $F(z)=e^{-f-ig}$, we have that $F(0) \ne 0$ and $F(w_k) \to 0, w_k \to w \in S, |w_k|< 1$, while $F$ is analytic in the unit disc and bounded by $1$ in absolute value
Hence $G(z)=1-F(z)/F(0)$ satisfies $G(0)=0, G(w_k) \to 1, w_k \to w \in S, |w_k|< 1$ and if $G(z)= \sum_{n \ge 1} x_nz^n$ one has that if $k \ge 1$ is fixed $$\sum_{n \ge 1}x_{nk}z^{nk}=\frac{1}{k}\sum_{\omega^k=1}G(\omega z)$$
From the properties of $G$ it follows that if $z_m \to 1, |z_m| < 1$, then $$\sum_{n \ge 1}x_{nk}z_m^{nk} \to 1, m \to \infty$$ so the sequence $x_{n}$ satisfies the required properties in an Abelian sense (and a bit more as the convergence is allowed on any disc sequence going to $1$, not only radially or nontangentially)