求一个数的真因数。

zophim

. X2 T" H( \- U  h% n; X
) a+ R1 E% H$ [- k/ f这个是一个函数，求一个数的真因数
% k, S! s/ ]$ G1 qfunction y=myfun(n)
$ o* D- P- u/ Z% n+ V1 N   y1=[];
   for i=1:fix(n/2);
       if mod(n,i)==0;
% B) p% x0 A- a/ ]; y           y1=[y1,i,n/i];
* N; W* _; P' h& B; N. z; i& y       end
   end
   y=unique(y1);   
$ F3 y' `4 I+ ]end
% p3 f" Q+ Z& P, }3 H- z然后我用这个函数求1到1000的亲和数对
9 L' S1 b" t4 t& I$ _for  a=1:10000;

    b=sum(myfun(a)); %求a的真因数和  为b
8 v( }# ?+ V8 T
6 C+ x! i/ S# v    c=sum(myfun(b)); %求b的真因数和

: [1 I" c& ?3 E+ x' ]# g      if   c==a  
/ T2 \% h  J$ f) ?
4 w$ i4 P0 f0 I8 R+ A# s, v          disp(a)                 %如果a的真因数和b的真因数相等  则输出a

      end
- G5 P* H3 o5 @, Z* l
end




) ]" _' o! D2 x& p! j! |* t    但是他不输出任何数，
附上亲和数对
亲和数是这样一对正整数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数值上发生重复，直接去掉后者
6 J1 J# j% s( ?% U8 Y2 M, Y( G改成
y1 = [];
for i = 1:fix(n/2)
   if mod(n,i) == 0
       y1 = [y1,i];
   end
( s# l0 w# ^% {+ o1 n/ F# Y0 jend
y = sum( unique(y1) );   
end