使用SAS判断水仙花数

data mywork.sxh;

*int是取整

*mod的是取余
do num = 100 to 999;
b100 = int(num/100);
b10 = int(mod(num,100)/10);
b = mod(num,10);
i = b100**3+b10**3+b**3;
if num =  i then output;
end;
run;
proc print;
run;