求一个数的真因数。

zophim

$ r2 N, B' w' T0 l* q# {: b
这个是一个函数，求一个数的真因数
& m, a+ J$ {% T, M) }6 hfunction y=myfun(n)
; o! y; ^# Z* W8 V( l; y% v   y1=[];
6 f& |. N: y1 ?: V  u0 j   for i=1:fix(n/2);
4 A$ m  z& \1 T+ m: {       if mod(n,i)==0;
           y1=[y1,i,n/i];
       end
   end
   y=unique(y1);   
% x0 }( A) J# l0 Wend
然后我用这个函数求1到1000的亲和数对
for  a=1:10000;
2 q7 T4 N+ P: }5 s
$ L4 L; R/ s- h* a! J9 I. y8 |    b=sum(myfun(a)); %求a的真因数和  为b

    c=sum(myfun(b)); %求b的真因数和

- v/ {( H3 g7 i% A      if   c==a  
  R' g4 v) W$ G* |# n5 c, `! j4 i! M
8 u( C' v/ Z6 Y, Y) I( W4 }4 y          disp(a)                 %如果a的真因数和b的真因数相等  则输出a
3 s7 z5 w+ d/ K3 J& D/ P, q
      end
+ g! z% x. m1 |9 {* q3 s
6 Y# Z6 G! D$ r& P+ R8 |" p8 Lend

4 d9 s$ a4 ]# ]/ K7 N$ V
% M# @1 E6 P+ O+ p% @+ D/ Y" K% w
; z! e1 m# G- D+ D" T# A! a
  {9 L# H- Y$ B  l4 m    但是他不输出任何数，
% ~' ^9 q1 b$ Q% o附上亲和数对
亲和数是这样一对正整数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

3 ]  L" M1 _- P- N" r你可以用sum函数，另外i增大时会和此前的n/i数值上发生重复，直接去掉后者
改成
" m; N, x( r" c( j# L& Zy1 = [];
for i = 1:fix(n/2)
   if mod(n,i) == 0
- n3 x0 [) H  B9 g. w       y1 = [y1,i];
   end
( p4 m4 n( \) Q# X$ `  w6 ^: i5 Zend
y = sum( unique(y1) );   
end