整数N的所有因子之和等于N的倍数,则N称为多因子完备数,求[1,1000]之间所有多因子完备数.

vfp:::::::
set talk off
clear
k=0
for n=1 to 1000
s=0
for i=1 to n
if mod(n,i)=0
s=s+i
endif
endfor
if mod(s,n)=0
k=k+1
endif
endfor
? k
set talk on
return