Неправильно разложена на множители разность кубов:
Правильная формула:
x³−y³=(x−y)(x²+xy+y²)
Ничего не меняет, суть та же.
Вот уточнил- изменив плюс на минус:
Алгоритм Ферма:
Если для числа n мы найдем такие x и y, что:

(x и y сравнимы по модулю n)
и при этом
x ≠ y (mod n)
Затем

То есть разность кубов делится на n , которая раскладывается:


делится на n ,
Но

не делится на n , а с какой то вероятностью, и

, не делится на n,
(с какой то вероятностью, может делиться на n) ,
значит, если в случае, когда оба не делятся, то с какой то вероятностью получим,
1 < НОД(n, x - y) < n ,
Тогда, этот НОД , даст первый нетривиальный делитель n , а второй, если подставить трехчлен,
PS "
то с какой то вероятностью получим" - для разности квадратов, эта вероятность больше 1/2,
доказательство видел в книге Василенко. Если в случае с кубами (а скорее всего так и есть),
то же самое с приемлемой точностью, то принцип факторизации тот же.
Значит если найти случаи, что с кубами быстрее алгоритм где то будет работать,
(не знаю, может быть, метод просеивания, или более быстрое нахождение "B-гладких чисел",
или часть алгоритма в шаге линейной алгебры, и т.п.), то весь принцип с кубами
может оказаться быстрее чем с квадратами, и мы улучшили бы ECM, SIQS, GNFS ,
Просто интересное наблюдение, но нигде не видел каких то исследований по этому поводу