cp's OEIS Frontend

A255243 Number x such that x | A255242(x).

%I A255243 #38 Apr 29 2019 20:50:01
%S A255243 1,4,10,16,18,64,80,96,195,256,462,576,768,880,1024,2560,3120,3136,
%T A255243 4096,6656,16384,40704,53248,57344,64000,65536,67896,78864,80640,
%U A255243 101376,103680,120320,120336,125440,126208,139264,147968,195840,217600,225280,250624,262144
%N A255243 Number x such that x | A255242(x).
%C A255243 For 4, 10, 195 we have x = A255242(x).
%H A255243 Amiram Eldar, <a href="/A255243/b255243.txt">Table of n, a(n) for n = 1..50</a>
%e A255243 For a(1) = 1 we have sigma(1) - 1 = 0 and 0 / 1 = 0.
%e A255243 Aliquot parts of a(2) = 4 are 1, 2 and their sum is 3.
%e A255243 Let us repeat the calculation with 1 and 2: 1 => 0; 2 => 1.
%e A255243 Their sum is 1. Finally, 3 + 1 = 4 and 4 / 4 = 1.
%e A255243 Aliquot parts of a(3) = 10 are 1, 2, 5. Their sum is 8.
%e A255243 Let us repeat the calculation with 1, 2 and 5: 1 => 0; 2 => 1; 5 => 1.
%e A255243 Their sum is 2. Finally, 8 + 2 = 10 and 10 / 10 = 1.
%e A255243 Aliquot parts of a(4) = 16 are 1, 2, 4, 8. Their sum is 15.
%e A255243 Let us repeat the calculation with 1, 2, 4 and 8: 1 => 0; 2 => 1; 4 => 1, 2; 8 => 1, 2, 4.
%e A255243 Their sum is 1 + 1 + 2 + 1 + 2 + 4 = 11.
%e A255243 Repeat the calculation with 1, 1, 2, 1, 2, 4: 1 => 0; 1 => 0; 2 => 1; 1 => 0; 2 => 1;  4 => 1, 2.
%e A255243 Their sum is 1 + 1 + 1 + 2 = 5.
%e A255243 Repeat the calculation with 1, 1, 1, 2: 1 => 0; 1 => 0; 1 => 0; 2 => 1; Their sum is 1.
%e A255243 Finally, 15 + 11 + 5 + 1 = 32 and 32 / 16 = 2.
%p A255243 with(numtheory): P:=proc(q) local a,b,c,k,n,t,v;
%p A255243 for n from 1 to q do b:=0; a:=sort([op(divisors(n))]); t:=nops(a)-1;
%p A255243 while add(a[k],k=1..t)>0 do b:=b+add(a[k],k=1..t); v:=[];
%p A255243 for k from 2 to t do c:=sort([op(divisors(a[k]))]); v:=[op(v),op(c[1..nops(c)-1])]; od;
%p A255243 a:=v; t:=nops(a); od; if type(b/n,integer) then print(n); fi; od; end: P(10^9);
%t A255243 f[s_] := Flatten[Most[Divisors[#]] & /@ s]; a[n_] := Total@Flatten[FixedPointList[ f, {n}]] - n; Select[Range[10000], Divisible[a[#], #] &] (* _Amiram Eldar_, Apr 06 2019 *)
%Y A255243 Cf. A001065, A255242.
%K A255243 nonn
%O A255243 1,2
%A A255243 _Paolo P. Lava_, Feb 19 2015
%E A255243 a(32)-a(42) from _Amiram Eldar_, Apr 06 2019