年度呼吁:请向OEIS基金会捐款支持OEIS的持续开发和维护。现在是我们的第61年,我们有超过378000个序列,我们已经获得了11000条引文(通常说“多亏了OEIS才发现”)。
编辑
提出
Robert Israel,<a href=“/A369969型/b369969_1.txt“>n表,n=1..10000时为a(n)</a>
过滤器:=proc(n)局部F,np,t;
如果n mod 5=0,则返回假fi;
F: =系数(n)[2];
np:=加(n*t[2]/t[1],t=F);
np模块5=0
结束进程:
选择(过滤器,[1..1000]美元)#罗伯特·伊斯雷尔2024年12月5日
经核准的
Robert Israel,<a href=“/A161142号/b161142_1.txt“>n,a(n)表,n=1..10000</a>
M: =10^6:#对于<=M的术语
S: ={1}:
对于从1到ilog10(M)的d,do
x: =4/9*(10^d-1);T: ={}:
对于s do中的s
T: =T联合{seq(s*x^i,i=1..层(log[x](M/s))};
od;
S: =S由壬T;
日期:
排序(转换(S,列表))#罗伯特·伊斯雷尔2024年12月5日
Robert Israel,<a href=“/A161141号/b161141_1.txt“>n,a(n)表,n=1..10000</a>
M: =10^8:#对于<=M的术语
x: =3/9*(10^d-1);T: ={}:
Robert Israel,<a href=“/A161140型/b161140_1.txt“>n,a(n)的表,对于n=1..10000</a>
M: =10 ^8:
x: =2/9*(10^d-1);T: ={}:
Robert Israel,<a href=“/A367984型/b367984_1.txt“>n表,n=1..10000时为a(n)</a>
M: =10000:#对于<=M的术语
立方体:={seq(i^3,i=0..floor(M^(1/3))}:
四分之一:={seq(i^4,i=0..floor(M^(1/4)))}:
排序(转换(选择(`<=`,{seq(seq(a+b,a=立方),b=四)},M),list))#罗伯特·伊斯雷尔2024年12月5日