直接代数反解比循环快,关键在于用整数校验r*r==discriminant避免sqrt浮点误差;暴力循环最坏需√(2x)次,x=1e12时超百万次;解n=(−1+√(1+8x))/2后须验证r²是否等于判别式且(r−1)为正偶数。

直接用代数反解,比循环快得多;关键不是算不算得出来,而是怎么避免 sqrt 的浮点误差导致误判。
为什么不能用 for 循环暴力试?
三角数定义是形如 n * (n + 1) / 2 的正整数。暴力法从 n = 1 开始试,直到结果 ≥ 目标值 x。最坏情况要迭代约 √(2x) 次——x 是 1e12 时,就要跑上百万次。对单次判断尚可,但批量或嵌套调用时明显拖慢。
真正高效的做法是把问题倒过来:给定 x,是否存在正整数 n 满足 n² + n − 2x = 0?解这个方程,得到 n = (−1 + √(1 + 8x)) / 2。只要这个 n 是正整数,x 就是三角数。
如何安全计算 sqrt(1 + 8x)?
double 类型的 std::sqrt 在输入超过 2⁵³ ≈ 9e15 后会丢失整数精度,而 int 最大才约 2.1e9,所以对 int 输入本身是安全的;但若扩展到 long long 或 uint64_t,就必须换方法。
- 先用
long long discriminant = 1LL + 8LL * x防溢出(尤其当x是int但接近INT_MAX) - 用
std::sqrt算近似值,再取最接近的整数:long long r = std::llround(std::sqrt(discriminant)) - 必须用整数校验:
r * r == discriminant—— 这一步不能省,否则sqrt返回999999.999999或1000000.000001都会导致static_cast<long long></long>截断出错 - 最后检查
(r - 1) % 2 == 0 && (r - 1) / 2 > 0,确保解出的n是正整数
边界和类型适配要注意什么?
对 int 类型输入,1 + 8 * x 最大约 1.7e10,仍在 long long 范围内,double 的 sqrt 可靠。但一旦换成 uint64_t(最大约 1.8e19),1 + 8x 可达 1.4e20,远超 double 精确表示整数的上限 2⁵³。
此时必须改用整数开方,比如牛顿迭代法,或调用 std::sqrtl(long double 版本,部分平台提供更高精度)。否则哪怕只错 1,校验 r * r == discriminant 就会失败,整个判定就漏判。
另外,负数、零、非整数值不参与判定——函数开头加 if (x 是必要防护。
真正容易被忽略的,不是公式推导,而是那个整数平方校验——它不是“优化”,而是必须步骤;少这一行,大数下大概率出错,而且错误难以复现。
C++免费学习笔记(深入):立即使用
在学习笔记中,你将探索 C++ 的入门与实战技巧!











