У меня получается максимум

пары, перебором чисел от

до

в несколько проходов, на каждом проходе постоянна сумма показателей в разложении на простые множители:
-

простых чисел, больших

, в любом случае останутся без пары;
- числа, большие

, у которых сумма показателей в разложении на простые равна

ставим в пару с их наибольшим делителем, если он пока свободен (пример -

); иначе - с наименьшим, если он свободен (пример -

, поскольку

уже занято в паре

); если и наименьший занят - этому числу не повезло (это числа

)
- далее четыре таких же прохода для сумм показателей от

до

;
- в конце остается

чисел

, из которых складывается

пар, а одно число добавляется к числам-"неудачникам"

.
Доказать оптимальность не возьмусь.