%I#13 2020年1月25日05:56:36

%S 1,2,8,10,12,32112181361441521845121088214422722528，

%电话：2736324833124531245927936819298001180017176186321890422984，

%电话：3276832896336643468849024571527724885952131076400212400309168335376

%N双优缺完全数：双优缺数k，2*k-bsigma（k）是k的双优除数，其中bsigma是k的两优除数之和（A188999）。

%C A271816的双功能版本。

%C包括2的所有奇数幂（A004171）。

%H Amiram Eldar，n的表，n的a（n）=1.171</a>

%由于e112的双优除数之和是1+2+7+8+14+16+56+112=216，并且2*112-216=8是112的双劣除数，因此e112在序列中。

%t f[n_]：=选择[Divisors[n]，函数[d，互质Q[d，n/d]]]；bsigma[m_]：=除数和[m，#&，最后一个@交点[f@#，f[m/#]]==1&]；双单位除数Q[div_，n_]：=如果[Mod[#2，#1]==0，上次@应用[Intersection，Map[Select[Divisors[#]，Function[d，CoprimQ[d，#/d]]&，{#1，#2/#1}]==1，False]@@{div，n}；aQ[n_]：=模[{d=2n-bsigma[n]}，如果[d<=0，False，双单位除数Q[d，n]]]；s={}；做[If[aQ[n]，AppendTo[s，n]]，{n，1，10000}]；秒

%o（PARI）udivs（n）={my（d=除数（n））；选择（x->（gcd（x，n/x）==1），d）；}

%o gcud（n，m）=vecmax（setintersect（udivs（n），udives（m）））；

%o biudivs（n）=选择（x->（gcud（x，n/x）==1），除数（n））；

%o isok（n）=我的（divs=biudivs（n），sig=vecsum（divs））；（sig<2*n）&&vecsearch（divs，2*n-sig）；\\_米歇尔·马库斯，2018年4月27日

%Y参见A004171、A188999、A271816、A292982、A303359。

%K nonn公司

%O 1,2号机组

%A _Amiram Eldar_和_Michael De Vlieger_，2018年4月22日