求一个数的真因数。

zophim

2 Z7 H" g! j# T0 `) ?' W7 X这个是一个函数，求一个数的真因数
; u* O) O1 d3 ?# C7 x9 K* Efunction y=myfun(n)
   y1=[];
   for i=1:fix(n/2);
       if mod(n,i)==0;
           y1=[y1,i,n/i];
       end
   end
   y=unique(y1);   
3 Z. k7 _1 _6 p3 N3 J7 }7 W3 V+ j) Rend
9 i" `; D1 X* X( t% p/ P3 b( H9 W, H然后我用这个函数求1到1000的亲和数对
for  a=1:10000;

    b=sum(myfun(a)); %求a的真因数和  为b
6 ~# z/ I5 b- I% _$ r7 c
    c=sum(myfun(b)); %求b的真因数和

      if   c==a  

          disp(a)                 %如果a的真因数和b的真因数相等  则输出a
9 _, z$ n. `- Q; v( g$ t5 \/ d
      end
; E) S; \: f, O
end


9 c, U" @4 ?+ ]; P) n/ n, I% d, W

+ ?% d# z9 i0 G8 Q8 y$ r4 i    但是他不输出任何数，
附上亲和数对
% h5 d" A0 D3 z/ K2 W亲和数是这样一对正整数a和b，使得a的所有真因子的和等于b，而b的所有真因子的和等于a。
! _5 X! O. B8 I就是220的真因子是1、2、4、5、10、11、20、22、44、55、110，它们的和是284；284的真因子是1、2、4、71、142，其和恰好是220。这是一对亲和数。
; D) O( _/ ?: q1 ~) g2 p+ G请教大神，非常感谢！

Taio

3 ^  K" j- S$ _0 Q" i2 p你可以用sum函数，另外i增大时会和此前的n/i数值上发生重复，直接去掉后者
! }7 V, k& c1 p& I2 ~) ?改成
y1 = [];
+ M- z% t, }. l' f  I( ^4 V2 L# F/ yfor i = 1:fix(n/2)
4 f) ?  a5 d4 P# ^. x   if mod(n,i) == 0
2 h0 t* r9 [8 W" s3 M3 [* T9 A       y1 = [y1,i];
   end
end
' k! b0 L& F- \/ J3 [. S, e) ?y = sum( unique(y1) );   
end