OFFSET
1,3
COMMENTS
Let f(n) = A050376(n) be the n-th Fermi-Dirac prime. The FDH number of a strict integer partition (y_1, ..., y_k) is f(y_1) * ... * f(y_k).
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..16384
FORMULA
For all n >= 1, a(A050376(n)) = n. - Antti Karttunen, Feb 18 2023
EXAMPLE
45 is the FDH number of (6,4), which has GCD 2, so a(45) = 2.
MATHEMATICA
nn=200;
FDfactor[n_]:=If[n==1, {}, Sort[Join@@Cases[FactorInteger[n], {p_, k_}:>Power[p, Cases[Position[IntegerDigits[k, 2]//Reverse, 1], {m_}->2^(m-1)]]]]];
FDprimeList=Array[FDfactor, nn, 1, Union]; FDrules=MapIndexed[(#1->#2[[1]])&, FDprimeList];
GCD@@@Table[Reverse[FDfactor[n]/.FDrules], {n, nn}]
PROG
(PARI) A319826(n) = { my(i=1, g=0, x=A052331(n)); while(x, if(x%2, g = gcd(g, i)); x>>=1; i++); (g); }; \\ (Uses the program given in A052331) - Antti Karttunen, Feb 18 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Sep 28 2018
EXTENSIONS
Secondary definition added by Antti Karttunen, Feb 18 2023
STATUS
approved