求一个数的真因数。

zophim

- {3 e4 v% ?4 L* [6 H
5 l& o6 _% _3 c4 J- g. C* N这个是一个函数，求一个数的真因数
function y=myfun(n)
$ u5 z8 _: C) `8 N$ g  `6 r) `   y1=[];
2 d; Y( J7 E# f( k3 b7 \& p   for i=1:fix(n/2);
       if mod(n,i)==0;
2 k; P( s1 X# {2 d; s3 J. e# O/ V           y1=[y1,i,n/i];
       end
; q( L9 E9 Q/ M& i   end
8 D  ~) G( |9 a* D3 Q6 m$ x   y=unique(y1);   
end
然后我用这个函数求1到1000的亲和数对
+ ?# ~! k# }6 B4 O' L! C( _4 W! efor  a=1:10000;

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

6 ?9 t5 O. U9 `2 I  }9 i    c=sum(myfun(b)); %求b的真因数和

+ U: {! r+ V" p6 v- k) h      if   c==a  
9 u( L5 D) U& ^$ s
          disp(a)                 %如果a的真因数和b的真因数相等  则输出a

# i2 O- Z8 |! p" \      end
5 G* l0 X7 ~8 M0 X, I
3 t5 B+ A% l- O5 N7 B/ N$ d8 J3 D# ^9 N$ Wend
$ ]2 o; M. P5 x5 r' M
5 n  z2 W- }+ {: j& ?
" S3 m( d) y; t+ P

" ^; J- c8 G9 v0 O  t    但是他不输出任何数，
) J! J+ _3 C4 k7 s2 l附上亲和数对
' A, U$ D8 Z& v( [/ 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

) P3 {1 ^* m% d% r! E你可以用sum函数，另外i增大时会和此前的n/i数值上发生重复，直接去掉后者
( P. }7 h! Q  a改成
$ P" u/ l  u. C3 @. Jy1 = [];
7 ]  \8 t$ W1 m+ V6 Z( Yfor i = 1:fix(n/2)
  d9 H  j; x- X7 t   if mod(n,i) == 0
, y" N) J6 K" m, \& i1 B7 Y* B       y1 = [y1,i];
   end
end
* g+ Q5 l8 _1 z7 K( y& ~y = sum( unique(y1) );   
end