Comment:

Here I consider a particular case .

$N=p_1p_2$ ⇒ $\sigma_N=p_1+p_2+p_1p_2+1$

Where $P_1$ and $p_2$ are primes.

Suppose $p_1p_2 +p_1+p_2=p$ is prime, then:

$\sigma_p=\sigma_N=p_1+p_2+p_1p_2+1=p+1 $

There can be infinitely many couples of this type, for example:

$p_1=2$, $p_2=13$ gives $N=26$ and $p=41$, or $p_1=3$ and $p_2=17$ gives $N=51$ and $p=71$

Now consider couple (33, 35), we try to find a third number with sum of divisors equal to that of $33$:

$N=33=3\times 11$ ⇒ $\sigma_N=3+11+33+1=48$

$p=3+11+33=47$ ⇒ $\sigma_p=47+1=48$

$\sigma_{q=35}=5+7+35+1=48$

so we have a triple $(33, 35, 47)$

There can be infinitely many triples of this type too.This experiment shows that the probability of existence of a triple with small difference is not zero if N is large enough.

For $\sigma(n)=\sigma(n+1)=\sigma(n+2)$, if we want to use this model, we have to look for a large number like $N=p_1p_2p_3 . . .$such that it gives numerous options for combination to construct two new numbers,the difference between three numbers $N$, $p$ and $q$ is as small as possible.This difference might never be equal especially equal to $1$. algorithm can be more efficient if one of these numbers is prime.

UPDATE:

I checked this by a simple proram which finds integer solutions for x and y in following Diophantine equation:

$xy+x+y=p$

Where p is a primes which is the known of equation. This is the results:

1- Some primes such as $9973, 9967, 9949, 9817, 9721 . . .$ give no solution.

2- Some primes such as $9923$ give only one prime solution for x and y:

$x=2, y=3307$ such that: $9923=2\times 3307+2+3307$ and we have:

$\sigma_{9923}=\sigma_{2\times 3307=6614}=9924$

3-Some primes give more than one prime solution for x and y:

$71$ gives: $(x, y, N)=(2, 23, 46), (3, 17, 51), (5, 11, 55)$

The differences are interesting:

$55-51=4$

$51-46=5$

$9719$ gives: $(x, y, N)=(89, 107, 9523),(53, 179, 9487), (11, 809, 8899), (5, 1619, 8095)$

The program is simple for Python:

for x in range (2, 5000):

```
for y in range (2, 5000):
if x*y+x+y==p:
print x, y
```

I gave p manually. I am not good enough to write a program that my computer do it automatically, also my computer is too weak for this calculations. But I am quite sure we can find triples or quadruple or even family of numbers with larger number of members which have equal sum of divisors, among these members some make an arithmetic progression.