這題 1+√(N+1) 是 N 的因數, 顯然 N + 1 會是個完全平方數.
而 N 是個 4 位數的數字, 介於 1000, 9999只要檢測 32 ~ 99 之間的數值
for i in range(32, 100):
n = i * i - 1
u = int((n + 1) ** 0.5) + 1
# 檢查 u 是否為 n 的第 4 個因數