%I #15 Nov 26 2020 03:19:09
%S 1,1,1,1,1,2,1,2,1,2,1,3,1,2,3,1,1,2,1,4,3,2,1,6,1,2,3,4,1,6,1,2,3,2,
%T 5,4,1,2,3,8,1,6,1,4,5,2,1,3,1,2,3,4,1,6,5,8,3,2,1,12,1,2,7,4,5,6,1,4,
%U 3,10,1,8,1,2,3,4,7,6,1,5,1,2,1,12,5,2,3,8,1,10,7,4,3,2,5,6,1,2,9,4,1,6,1,8,15
%N a(1) = 1; for n>1, a(n) = n/(largest Fermi-Dirac factor of n).
%C For n > 1, a(n) = the smallest positive number d such that n/d is a "Fermi-Dirac prime", a term of A050376.
%H Antti Karttunen, <a href="/A302776/b302776.txt">Table of n, a(n) for n = 1..65537</a>
%F a(n) = n / A223491(n).
%F a(n) = A302023(A052126(A302024(n))).
%t f[p_, e_] := p^(2^Floor[Log2[e]]); a[n_] := n / Max @@ f @@@ FactorInteger[n]; Array[a, 100] (* _Amiram Eldar_, Nov 26 2020 *)
%o (PARI)
%o A209229(n) = (n && !bitand(n,n-1));
%o A302777(n) = A209229(isprimepower(n));
%o A302776(n) = if(1==n,n,fordiv(n, d, if(A302777(n/d), return(d))));
%Y Cf. A050376, A052331, A223491, A302023, A302024, A302777, A302792.
%Y Cf. A084400 (gives the positions of 1's).
%Y Cf. also A052126, A284600.
%K nonn
%O 1,6
%A _Antti Karttunen_, Apr 13 2018