纳西瑟斯数(水仙花数)是n位正整数,其各位数字的n次幂之和等于自身;0不计入,1~9均为1位纳西瑟斯数;判断时需先求位数(避免log10精度问题),再用整数幂运算防浮点误差与溢出,sum须用long long并及时剪枝。

直接判断:先算位数,再逐位求幂累加,最后比对原数——但要注意整型溢出和 0 的边界情况。
什么是纳西瑟斯数(水仙花数)?
纳西瑟斯数(Narcissistic number),也叫自幂数或水仙花数,指一个 n 位正整数,其各位数字的 n 次幂之和等于它本身。例如 153 是 3 位数,且 1^3 + 5^3 + 3^3 == 153,所以它是纳西瑟斯数。
注意:0 不算——虽然 0^1 == 0,但按惯例纳西瑟斯数定义为「正整数」;1 到 9 都是 1 位纳西瑟斯数(因为 1^1 == 1,…,9^1 == 9)。
如何用 C++ 计算位数并避免 pow 的精度陷阱
C++ 标准库的 pow(double, int) 返回 double,对大整数(如 7 位以上)容易因浮点舍入导致错误结果,比如 pow(9, 7) 可能返回 4782968.999999,floor 后变成 4782968,少 1。
组合式C++代码评审方案,融合静态分析、AI推理、多轮迭代评审和C++专项检查,适用于PR审查、增量代码审查、全项目评审和代码质量评分,触发词包括review cpp、cpp代码评审、C++review、代码审查。
立即学习“C++免费学习笔记(深入)”;
- 用整数循环自乘实现幂运算,例如写个
int ipow(int base, int exp),每次乘法后检查是否溢出(可用if (result > INT_MAX / base)提前截断) - 位数不要用
log10(n) + 1——对n == 1000这类边界,log10可能返回2.9999999,取整成 2 - 更稳妥的做法是转字符串:
to_string(n).length(),或用 while 循环除 10 计数
完整判断函数怎么写(含关键防错点)
以下是一个生产可用的 is_narcissistic 实现逻辑:
bool is_narcissistic(int n) {
if (n <= 0) return false; // 排除负数和 0
string s = to_string(n);
int len = s.length();
long long sum = 0; // 用 long long 防止中间求和溢出(如 9^9 * 10 ≈ 3.8e9,int 可能溢)
for (char c : s) {
int digit = c - '0';
long long power = 1;
for (int i = 0; i < len; ++i) {
if (power > INT_MAX / digit && digit != 0) return false; // 提前防乘法溢出
power *= digit;
}
sum += power;
if (sum > n) break; // 提前退出,优化性能
}
return sum == n;
}- 必须用
long long存sum:即使输入是int,像9474(4 位)的各位四次方和已达9^4 + 4^4 + 7^4 + 4^4 = 6561 + 256 + 2401 + 256 = 9474,但更大的 6 位数(如548834)求和过程可能超int范围 - 循环幂计算中,
digit == 0时跳过溢出检查(0^k恒为 0) - 一旦
sum > n就可提前返回false,尤其对高位数效率提升明显
实际使用时最容易忽略的三个点
很多人写完测试 153、371 就以为没问题,但上线后在边界上翻车:
-
INT_MAX以内最大的纳西瑟斯数是115132219018763992565095597973971522401(39 位),但 C++int或long long根本存不下——所以实际应用中,若需支持大数,得换string输入 + 自定义大数幂运算,不能依赖原生整型 - 输入是
0或负数时,多数人忘了加n 判断,导致 <code>to_string(0)得到"0",长度为 1,然后算0^1 == 0,误判为真 - 用
pow的代码在 x86_64 和 ARM 平台下结果可能不一致,因为pow底层调用不同数学库,浮点误差路径不同
真要稳定判断,别碰 pow,手写整数幂;别信 log10,老实用字符串或除法数位;别假设输入小,溢出检查得嵌在每一步乘法里。

















