求一个数的真因数。

zophim

% t. G% J2 n+ `# ]" W
. u( F8 |4 h% h/ v8 E这个是一个函数，求一个数的真因数
% E9 H! H: s3 \  N5 ^function y=myfun(n)
   y1=[];
   for i=1:fix(n/2);
2 ?  D! z& B5 m: r( Z5 M       if mod(n,i)==0;
9 a! }) J" f1 e+ e+ W           y1=[y1,i,n/i];
       end
   end
   y=unique(y1);   
! k4 W% c/ w: V+ x# E: ^/ Hend
然后我用这个函数求1到1000的亲和数对
for  a=1:10000;
! M+ Y- t  ^- E8 x9 b- t
    b=sum(myfun(a)); %求a的真因数和  为b

$ F) i7 G. u0 v3 N* E/ r& x    c=sum(myfun(b)); %求b的真因数和
1 p) H3 k- X; s" a; M9 i
      if   c==a  
2 t- H, V% E3 w% \9 B7 J4 ]3 H
          disp(a)                 %如果a的真因数和b的真因数相等  则输出a

/ O& F2 ]9 c( l0 ?' j  z      end
# {+ V2 C; M1 @4 }* ?2 `1 L1 u
% C0 @* y' x* N" Z7 }( a$ Mend

, G" R1 [. f' v: k+ u

- S2 ~' l+ f, o8 y
2 P5 m9 {, T2 L1 C( i! ~    但是他不输出任何数，
& t1 _; W( T4 e3 c# Z4 p- V# C" d附上亲和数对
3 J, V9 a: `5 C- o& K0 q# R$ S7 i亲和数是这样一对正整数a和b，使得a的所有真因子的和等于b，而b的所有真因子的和等于a。
就是220的真因子是1、2、4、5、10、11、20、22、44、55、110，它们的和是284；284的真因子是1、2、4、71、142，其和恰好是220。这是一对亲和数。
请教大神，非常感谢！

Taio

你可以用sum函数，另外i增大时会和此前的n/i数值上发生重复，直接去掉后者
改成
! L  J) i! `4 E. K) W, y0 ey1 = [];
, ?7 w7 {9 b  L$ ?; c/ _/ ofor i = 1:fix(n/2)
   if mod(n,i) == 0
- }2 K' _8 [7 Z$ j# i8 H" H# Z       y1 = [y1,i];
   end
! D" T& @7 @* W9 q$ t* |end
y = sum( unique(y1) );   
9 v. U0 f% J0 }8 n# gend