想定解法は について全探索して を二分探索です。だいたい計算回数は最大でも くらいで間に合います。
または,分母 を全探索し,分母 を固定したうえで, 以下で出来る限り近くなるように を で求め(切り捨て除算でできます), を ずつ足していく方法も可能です。この方法の方が実装上は簡単です。 の扱いに注意してください。
※追記(2022/11/06) ファレイ数列で二分探索を行うことで,分母の最大値を として, 回ほどの計算回数で答えが求まり,とても高速であるそうです。(startcppさんの提出を参考にしてください。)