Given a Noetherian integral domain $A$ and a finitely generated torsion $A$-module $M$, we can define the divisor, or content of $M$ to be $\mathrm{div}(M)= \sum_{\mathrm{ht}(P)=1} \ell(M_P) [P]$, where the sum ranges over all height one prime ideals of $A$, and $M_P$ is the localization of $M$ at $P$, which is an $A_P$-module of finite length $\ell(M_P)$. The above sum makes sense because the only nonzero terms correspond to associated prime ideals of $M$ with height $1$.
Now for the question, as in the title. Let $a,b\in A$ be such that $\mathrm{div}(A/aA)=\mathrm{div}(A/bA)$. Can we deduce that $a,b$ differ only by a unit of $A$?
This is trivially true if $A$ is a UFD (factorial domain) and maybe also if $A$ is integrally closed, but is it true for every Noetherian integral domain? The question should be easy to settle once one has studied (which I have not yet done) Chapter 7 of Bourbaki's Commutative Algebra, for example.