|
%I #12 Sep 23 2018 20:56:26
%S 1,2,2,3,2,3,2,4,3,4,2,4,2,4,4,5,2,5,2,4,4,4,2,5,3,4,4,4,2,7,2,6,4,4,
%T 3,6,2,4,4,5,2,5,2,6,6,4,2,6,3,6,4,6,2,7,4,4,4,4,2,9,2,4,5,7,4,5,2,6,
%U 4,6,2,8,2,4,6,6,2,6,2,6,5,4,2,8,4,4,4,6,2,9,2,6,4,4,4,7,2,6,6,6,2,6,2,6,8
%N Number of divisors of 2n which are also differences between consecutive divisors of 2n (ordered by size).
%C For odd numbers the intersection is empty.
%H Antti Karttunen, <a href="/A060741/b060741.txt">Table of n, a(n) for n = 1..65537</a>
%F a(n) = A000005(2n) - A060764(2n).
%e For n=35, 2n=70; divisors={1,2,5,7,10,14,35,70}; differences={1,3,2,3,4,21,35}; intersection={1,2,35}, so a(35)=3.
%t a[n_ ] := Length[Intersection[Drop[d=Divisors[2n], 1]-Drop[d, -1], d]]
%t Table[Length[Intersection[Divisors[2n],Differences[Divisors[2n]]]],{n,110}] (* _Harvey P. Dale_, Nov 22 2015 *)
%o (PARI)
%o A060764(n) = { my(divs=divisors(n), diffs=vecsort(vector(#divs-1,i,divs[i+1]-divs[i]), ,8), c=#divs); for(i=1,#diffs,if(!(n%diffs[i]),c--)); (c); };
%o A060741(n) = (numdiv(2*n) - A060764(2*n)); \\ _Antti Karttunen_, Sep 21 2018
%Y Cf. A000005, A060742, A060764.
%K nonn
%O 1,2
%A _Labos Elemer_, Apr 23 2001
%E Edited by _Dean Hickerson_, Jan 22 2002
|
