login
a(1) = 1; for n>1, a(n) = n/(largest Fermi-Dirac factor of n).
6

%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