%I#11 2021年11月4日20:47:39
%S 15324526136942547763772580183384587390998110171025,
%电话:1233132513411377142115571573162917731805181320092057,
%电话:2061209721692225233132349242124252525252925972637264522853287329893141317733321335734253501350935736093681
%N奇数k,其中A064989(sigma(k))<A064988(k),形式为p^(1+4k)*r^2,其中p是形式1+4m的素数,r>1,gcd(p,r)=1。
%C显然,任何假设的奇完全数都不会出现在这个序列中,也不会出现在A348939中。
%C在1..2^20范围内的数字中,9644位于该序列中,3865位于A348939中。在<=2^25的数字中,229480按此顺序排列,88270在A348939中。
%H<a href=“/index/Pri#prime_indices”>根据素因式分解中的索引计算出的序列的索引项</a>
%H<a href=“/index/Si#SIGMAN”>与sigma(n)相关序列的索引条目</a>
%tq[n_]:=模块[{f=FactorInteger[n]},p=f[[;,1]];e=f[[;;,2]];odde=选择[e,OddQ];长度[e]>1&&Length[odde]==1&&Divisible[odde[[1]]-1,4]&&Divisible[p[[位置[e,odde[[1]]][[1,1]]]-1,4]];f[2,e_]:=1;f[p_,e_]:=下一素数[p,-1]^e;s[1]=1;s[n_]:=倍@@f@@FactorInteger[n];选择[范围[1,4000,2],q[#]&s[DivisorSigma[1,#]]<s[#]&](*_Amiram Eldar_,2021年11月4日*)
%o(PARI)
%o A064989(n)={my(f=因子(n));如果(n>1&&f[1,1]==2),f[1,2]=0);对于(i=1,#f~,f[i,1]=precprime(f[i、1]-1));因子回退(f)};
%o isA228058(n)=如果(!(n%2)||(ω(n)<2),0,my(f=因子(n),y=0);对于(i=1,#f~,if(1==(f[i,2]%4),if((1==y)||(1!=(f[i,1]%4)),return(0),y=1),if(f[i,2]%2,return(0)));(y) );
%o是A348748(n)=((n%2)&&(A064989(σ(n))<A064988(n;
%o是A348938(n)=(是A228058(n)&是A348748(n;
%A228058和A348748的Y交叉口。
%Y参见A064989、A326042、A348939。
%K nonn公司
%O 1,1号机组
%2021年11月4日,安提·卡图宁
|