OEIS哀悼西蒙斯感谢西蒙斯基金会支持包括OEIS在内的许多科学分支的研究。
登录
OEIS由OEIS基金会的许多慷慨捐赠者.

 

标志
提示
(来自的问候整数序列在线百科全书!)
A302042型 A032742美元基于埃拉托斯坦筛的非标准因式分解过程的模拟(A083221号). 20

%I#30 2018年4月5日20:34:17

%S 1,1,2,1,3,1,4,3,5,1,6,1,7,5,8,1,9,10,9,11,12,5,13,7,14,15,1,

%电话:16,15,17,7,18,1,19,11,20,12,21,23,1,24,7,25,26,1,27,25,28,

%U 27,29,1,30,1,31,13,32,1,33,1,34,33,35,1,36,1,37,17,38,11,39,1,40,39,41,1,42,35,44,1,45,49,46,45,47,13,48单位

%N A032742类似物,用于基于Eratosthenes筛的非标准因子分解过程(A083221)。

%类似于[A020639(n),A032742(n)],有序对[A020639(n),a(n)]对于每个n是唯一的。迭代n,a(n),a(a(n)),a(a(a(n))。。。,直到达到1,然后取每个项的最小素因子(A020639),得到一组按升序排列的素数,对于每个自然数n>=1是唯一的。置换对A250245/A250246在n的非标准素因式分解和n的普通素因式化之间映射。

%H Antti Karttunen,n的表,n=1..65537的a(n)</a>

%H<a href=“/index/Si#sieve”>筛子生成序列的索引条目</a>

%F对于n>1,a(n)=A250469^(r)(A078898(n)),其中r=A055396。

%F a(n)=A250245(A032742(A250246(n))。

%F a(n)=n-A302043(n)。

%e对于n=66,A020639(66)[其最小素因子]为2。因为A055396(66)=A000720(2)=1,所以a(66)就是A078898(66)=66/2=33。

%e对于n=33,A020639(33)=3和A055396(33)=2,则a(33)=A250469(A078898(33))=A250.469(6)=15。[15在数组A083221]中低于6。

%e对于n=15,A020639(15)=3和A055396(15)=2,则a(15)=A250469(A078898(15))=A250.469(3)=5。[5小于3是数组A083221]。

%e对于n=5,A020639(5)=5和A055396(5)=3,则a(5)=A250469。

%e收集A020639给出的素数,我们得到一组多因子:[2,3,3,5]。注意,2*3*3*5=90=A250246(66)。

%e如果我们从n=66开始,迭代映射n->A302044(n)[而不是n->A302042(n)],并对获得的每个项应用A020639,我们只得到每个素数的一个实例:[2,3,5]。然后将A302045应用于相同的项,我们得到了这些素数的相应指数(重数):[1,2,1]。

%o(PARI)

%o假设A250469及其逆运算A268674已预计算,则以下运算速度足够快:

%o A302042(n)=如果(1==n,n,my(k=0);而(n%2),n=A268674(n);k++);n=n/2;而(k>0,n=A250469(n);k——);(n) );

%o(PARI)

%o A020639(n)=如果(n>1,如果(n>n=系数(n,0)[1,1],n,系数(n)[1,1,1]),1);\\来自A020639

%o A078898(n)=如果(n<=1,n,my(spf=A020639(n),k=1,m=n/spf);而(m>1,如果(A020639(m)>=spf,k++);m--);(k) );

%如果我们将A078898作为A020639的序数变换进行预计算,则速度更快:

%o up_to=65537;

%o ordinal_transform(invec)={my(om=Map(),outvec=vector(length(invec)),pt);对于(i=1,length;

%o v078898=序号变换(向量(up_to,n,A020639(n)));

%o A078898(n)=v078898【n】;

%o A302042(n)=if((1==n)||isprime(n),1,my(c=A078898(n),p=素数(-1+素数pi(A020639(n))+素数pi(A020639(c))),d=A078898(c),k=0);而(d,k++;如果((1==k)||(A020639(k)>=p),d-=1);(k*p));

%Y参见A020639、A032742、A055396、A078898、A083221、A250245、A250269、A268674、A276151、A280496、A302032。

%Y参考以下类比:A302041(ω)、A253557(bigomega)、A302043、A302044、A302045(最小素数存在的指数)、A302046(素数特征滤波器)、A301050(Moebius mu)、A320251(τ)、02052(正方形的char.fun)、A30.239、A302055(阿里斯导数)。

%K nonn公司

%O 1,4型

%A _Antti Karttunen,2018年3月31日

查找|欢迎光临|维基|注册|音乐|地块2|演示|索引|浏览|更多|网络摄像头
贡献新序列。或评论|格式|样式表|变换|超级搜索|最近
OEIS社区|维护人OEIS基金会。

许可协议、使用条款、隐私政策。.

上次修改时间:美国东部夏令时2024年6月18日19:02。包含373486个序列。(在oeis4上运行。)