6
$\begingroup$

I want to know if there exists a measure $\mu$ on the positive integers (equipped with the $\sigma-$algebra of all subsets) satisfying:

1) For all $n > 1$, $\mu(A_n) = 1/n$, where $A_n = \{n, 2n, \ldots \}$ is the multiples of $n$, and

2) The events $A_m, A_n$ are independent when $m$ and $n$ are relatively prime.

If such a measure existed, it would mimic the notion of drawing a "random integer" from all of $\mathbb{N}$, at least with respect to these divisibility properties.

1) and 2) hold approximately for the uniform distribution on $\mathbb{N} \cap [1,X]$ for large $X$; a similar weakening of the question is to use asymptotic density, i.e. define $\mu(A) = \lim_{n \to \infty} \frac{|A \cap \{1, 2, 3, \ldots, n\}|}{n}$, which satisfies 1) and 2) but is not countably additive.

The zeta distribution, given by (for any $s > 1$) the density $\nu_s(n) = \frac{1}{\zeta(s)} n^{-s}$, comes very close: it is an easy check that 2) holds, but $\nu_s(A_n) = n^{-s}$ for all $n$. (Interestingly, as $s \to 1^+$, for a fixed subset $A \subset \mathbb{N}, \nu_s(A)$ converges to the asymptotic density of $A$, if it exists.)

I suspect the answer is that no such measure exists: the fact that the zeta distribution fails makes me think 1 and 2 are somehow at odds.

I believe any $\mu$ satisfying 1) and 2) should have $\mu(\mathbb{N}) = \infty$, but I don't have a proof.

Thanks!

Edit: As people have pointed out, 1) clearly implies 2). (Whoops!) And @Zhoraster has given a nice proof that no measure can exist in this case. Here's what (I think) will be a harder question:

Can we find a measure $\mu$ satisfying

1a) $\mu(A_p) = 1/p$ for primes $p$, and

2a) The events $A_p$ and $A_q$ are independent when $p$ and $q$ are prime.

Now it's not obvious that $1$ and $2$ hold for products of many primes: we do get $\mu(A_{pq}) = 1/pq$ for distinct primes $p$ and $q$, but that's it. I don't think the inclusion-exclusion proof will work as stated, but perhaps a similar idea can be used...

  • 0
    View $\mu$ as a bounded linear map $l^\infty \to \mathbb{R}$ where $A_n$ is the sequence indicating the multiples of $n$. Then $\mu$ is the [natural density](https://en.wikipedia.org/wiki/Natural_density). It exists by the Hahn-Banach theorem applied to $l^\infty$, but it may be very strange for weird subsets of $\mathbb{N}$. Also, you can restrict to the closed subspace $E = \{ (a_k) \in l^\infty, \lim_{N \to \infty} \frac{1}{N} \sum_{k=1}^N a_k \text{ exists } \}$ of $l^\infty$, in that case $\mu$ is well-defined.2017-02-14
  • 2
    As I said in the post, the 'natural density' or 'asymptotic density' is not a measure, since it isn't countably additive.2017-02-14
  • 0
    As a linear operator, it can be extended to any bounded sequence by Hahn-Banach. Not sure there is much more to say. And try to be more clear next time2017-02-14
  • 0
    The linear operator you're describing is perfectly valid, but it doesn't answer the question, since it doesn't represent a measure. For instance, the natural density of any finite subset of $\mathbb{N}$ (your operator $\mu$ applied to an element of $l^\infty$ with only finitely many non-zero terms) is zero, but the natural density of the infinite union of all even integers is 1/2. As you say, the 'measure' is finitely additive, (or 'linear' as an operator on $l^\infty$), but it isn't countably additive. What is exactly is unclear? I'm happy to clarify further.2017-02-15
  • 0
    Ok, you really want a probability measure on $\mathbb{N}$ such that $P(X \in A_n) = 1/n$, then use the Möbius function and independentness for isolating $P(X = n)$2017-02-15
  • 0
    This seems like a good idea: could you be a little more specific? I don't see how to use mobius inversion or the mobius function to solve for $P(X = n)$... do you have a particular formula or trick in mind?2017-02-15
  • 0
    Assuming $X \in A_n, X \in A_m$ are independent for every $n,m$ there are two formulas : $P(X=n) = \sum_{m=1}^\infty \mu(m) P(X \in A_{nm}) = \sum_{m=1}^\infty \frac{\mu(m)}{nm} = 0$ and $P(X=n) = \prod_p P(X \not \in A_{np}) = \prod_p (1-\frac{1}{np}) = 0$. I'd say you need to find a way to modify those using only $X \in A_n, X \in A_m$ independent when $gcd(n,m)=1$2017-02-15
  • 0
    @Richey do you want a probability or a measure ? I can show that such a probability cannot exist, but for a measure I don't know2017-02-18
  • 0
    @charMD Yes, I would be very interested a proof that there can't be a probability measure satisfying those properties, but I'm also curious to know if they are satisfied for some measure, possibly with infinite mass.2017-02-18
  • 1
    @Richey By the way, isn't $(2)$ a consequence of $(1)$ ?2017-02-18

2 Answers 2

4

Such measure doesn't exist. The idea is to show that the only measure satisfying 1 (which, by the way, obviously implies 2) is $$ \mu(\{n\}) = \frac{1}{n\zeta(1)} = 0, n\ge 1, $$ which is a contradiction.

So let $\mu$ satisfy 1. Denote by $\mathbb{P}$ the set of primes. Then, for any $n\in \mathbb N$, we have by the inclusion-exclusion formula (all series below are divergent, but one can proceed by finite approximations) $$ \mu(\{n\}) = \mu(A_n) - \mu\Big(\bigcup_{p\in \mathbb{P}} A_{np}\Big) \\ = \frac{1}{n} - \sum_{p\in \mathbb P}\mu(A_{np}) + \sum_{\substack{p_1,p_2\in \mathbb P\\p_1\neq p_2 }}\mu(A_{np_1}\cap A_{np_2}) - \sum_{\text{distinct } p_1,p_2,p_3\in \mathbb P}\mu(A_{np_1}\cap A_{np_2}\cap A_{np_3}) + \dots \\ = \frac1n + \sum_{k=1}^\infty (-1)^k \sum_{\text{distinct } p_1,p_2,\dots,p_k\in \mathbb P}\mu\Big(\bigcap_{j=1}^k A_{np_j}\Big)\\ = \frac1n + \sum_{k=1}^\infty (-1)^k \sum_{\text{distinct } p_1,p_2,\dots,p_k\in \mathbb P}\mu\Big( A_{np_1p_2\cdots p_k}\Big) \\ = \frac1n + \sum_{k=1}^\infty (-1)^k \sum_{\text{distinct } p_1,p_2,\dots,p_k\in \mathbb P}\frac1{np_1p_2\cdots p_k} = \frac1n\prod_{p\in \mathbb P}\Big(1-\frac{1}{p}\Big) = 0, $$ as claimed.


Similarly, it can be shown that the unique measure satisfying $\mu_s(A_n) = n^{-s}$ with $s>1$ is $$ \mu_s(\{n\})=\frac1n\prod_{p\in \mathbb P}\Big(1-\frac{1}{p^s}\Big) = \frac{1}{n^s \zeta(s)}, n\ge 1. $$

  • 0
    Great answer which is obviously more appropriate than mine ; minor detail : on your seventh line you wrote twice $\mu (A_n)-\mu \Big( \bigcup \limits_{p \in \mathbb{P}} A_{np} \Big)$2017-02-18
  • 0
    @charMD, thanks, fixed.2017-02-18
  • 0
    Nice solution! I had thought to use this inclusion-exclusion formula, but didn't realize it's (very simple) connection with the zeta function. I added an edit weakening the assumptions a little bit... do you think a similar proof will work in the case where we only assume 1) and 2) hold for distinct primes $p$ and $q$?2017-02-18
  • 0
    @JRichey If you only assume 1 and 2 for primes then, I think, such measure may be easily constructed.2017-02-19
  • 0
    @Zhoraster What do you have in mind? I don't have any ideas on how to construct one.2017-02-19
  • 0
    @JRichey, no quick solution, but I think one can achieve this by assigning some positive values to $pqr$ only, where $p,q,r$ are different primes. The equations one needs to solve are $\sum_{r:r\neq p, r\neq q} \mu(pqr) = \frac{1}{pq}$, $\sum_{q,r:q\neq p, r\neq p} \mu(pqr) = \frac{1}{p}$. Heuristically, for the first $N$ primes there are about $N^2$ equations with $N^3$ variables, so there is a plenty of space to make this possible.2017-02-19
  • 0
    OK thanks, yeah that sounds like it will work. I'll give it a try.2017-02-21
1

Such a measure can not be a probability (i.e. we can't have $\mu(\mathbb{N})=1$).

Ad absurdum, suppose we have a probability $P$ satisfying $(1)$, $(2)$.

Let $(p_n)_{n \ge 1}$ be an ordered enumeration of prime numbers. There is no positive integer which has an infinity of prime divisors, and thus $$\bigcap \limits_{n \in \mathbb{N}} \bigcup \limits_{k \ge n} A_{p_k} = \varnothing.$$

Hence $P \Big(\bigcap \limits_{n \in \mathbb{N}} \bigcup \limits_{k \ge n} A_{p_k} \Big) = 0$, so $$P \Big( \bigcup \limits_{n \in \mathbb{N}} \bigcap \limits_{k \ge n} \overline{A_{p_k}} \Big) > 0.$$

However, $P$ is continuous below so \begin{align*} P \Big( \bigcup \limits_{n \in \mathbb{N}} \bigcap \limits_{k \ge n} \overline{A_{p_k}} \Big) & = \lim \limits_{n \to \infty} P \Big( \bigcap \limits_{k \ge n} \overline{A_{p_k}} \Big) \\ & = \lim \limits_{n \to \infty} \prod \limits_{k \ge n} P (\overline{A_{p_k}}) \end{align*}

because $P$ is a probability, the $A_{p_k}$ are mutually independant, and thus so are the $\overline{A_{p_k}}$. Hence \begin{align*} P \Big( \bigcup \limits_{n \in \mathbb{N}} \bigcap \limits_{k \ge n} \overline{A_{p_k}} \Big) & = \lim \limits_{n \to \infty} \prod \limits_{k \ge n} \Big( 1 - P (A_{p_k}) \Big) \\ & \le \lim \limits_{n \to \infty} \prod \limits_{k \ge n} e^{-P(A_{p_k})} \end{align*}

because $1-x \le e^{-x}$ (convexity inequality). We deduce finally :

$$P \Big( \bigcup \limits_{n \in \mathbb{N}} \bigcap \limits_{k \ge n} \overline{A_{p_k}} \Big) \le \lim \limits_{n \to \infty} e^{-\sum \limits_{k=1}^n P(A_{p_k})}$$

But $$\sum \limits_{k=1}^n P(A_k) = \sum \limits_{k=1}^n \frac{1}{p_k} \underset{n \to \infty}{\longrightarrow} +\infty$$

(see Divergence sum of reciprocal of primes), so we conclude

$$P \Big( \bigcup \limits_{n \in \mathbb{N}} \bigcap \limits_{k \ge n} \overline{A_{p_k}} \Big) = 0$$

which is absurd. Hence, there is no probability satisfying your conditions.


See Zhoraster's post which provides a better and clearer answer to your question.

  • 1
    There are some typos, but essentially it is correct. In fact, you show that for each $n$ the measure of numbers not dividing $p_n,p_{n+1}\dots$ must be zero; as a result, by letting $n\to\infty$, the measure of whole $\mathbb N$ must be zero. The assumption that $\mu$ is a probability measure is not significant, as one can modify the argument for the set of even numbers, not the whole $n$ (my proof is essentially in the very same spirit).2017-02-18