%I#20 2023年7月23日13:42:39

%S 1,1,1,1,1,1,1,1,1,1,1,1,,1,1,11,1,1,1，1,1，1,1,1，1,1，1,1,1,2,1,1，

%T 2,1,1,1,1,2,1,1,1,1,1,2,1,1,2,1,1,1,1,1,1,1,1,1,1,11,1,1,2,1,1,1,1，

%U 1,2,1,1,1,1,2,1,1,2,1,2,1,1,2,1,1,2,1,2,1,1,1,3,1,1,1,1,2,1,1,1,1,1,1,1,1,1,2,2,2,1,1,2,2,2,2,2

%N形式为4*k+1的N的除数最多为sqrt（N）。

%F G.F.：和{k>=0}x ^（4*k+1）^2/（1-x^（4]k+1））。

%t表格[Count[Divisors[n]，_？（#<=Sqrt[n]&MemberQ[{1}，Mod[#，4]]&）]，{n，100}]

%t nmax=100；系数列表[级数[Sum[x^（4k+1）^2/（1-x^

%Y参考A001826、A038548、A364387。

%K nonn公司

%O 1,25号机组

%A _Ilya Gutkovskiy，2023年7月21日