Definition ijs (n: nat): list (nat * nat) :=
  concat (map (fun j => map (fun i => (i, j)) (nats 0 j)) (nats 0 n)).

Lemma sumOver_ijs_bound n: RsumOver (ijs n) (fun ij => 2 / INR (S (snd ij - fst ij))) <= 2 * INR n * INR (S (log2ceil n)).