%I#17 2024年2月24日11:03:59
%序号1,2,8,20,40384102401267201013760486604807612661760,
%电话:47317450752016701625625331240364210021904040742624909465600,
%电话:23814785347489280093297677510746570752066941119654277247139849044593230639040844267520
%N a(N)=A203418(N)/A000178(N)。
%H G.C.Greubel,n的表格,n=1..140的a(n)</a>
%H R.Chapman,<a href=“https://www.maa.org/sites/default/files/Robin_Chapman27238.pdf“>取整数值的多项式,《数学杂志》,29(1996),121。
%t composite=选择[Range[100],CompositeQ];(*A002808*)
%tz=20;
%t f[j]:=复合[[j]];
%tv[n_]:=乘积[乘积[f[k]-f[j],{j,1,k-1}],{k,2,n}];
%t d[n_]:=乘积[(i-1)!,{i,1,n}];
%t表[v[n],{n,z}](*A203418*)
%t表格[v[n+1]/v[n],{n,z}](*A203419*)
%t表[v[n]/d[n],{n,z}](*此序列*)
%o(岩浆)
%o A002808:=[2..250]|非IsPrime(n)]中的n:n;
%o巴恩斯G:=func<n|(&*[阶乘(k):[0..n-2])中的k)>;
%o a:=函数1选择1其他(&*[(&*[A002808[k+2]-A002808[j+1]:[0..k]]中的j):[0..n-2]中的k)/BarnesG(n+1)>;
%o[1..40]]中的[a(n):n;//_G.C.Greubel,2024年2月24日
%o(SageMath)
%o A002808=[n代表(2..250)中的n,如果不是_prime(n)]
%o def BarnesG(n):返回乘积(范围(1,n-1)中j的阶乘(j))
%o定义a(n):返回产品(范围(k+1)中j的产品(A002808[k+1]-A002808[j])范围(n-1)中k的产品)/BarnesG(n+1)
%o[a(n)for n in range(1,41)]#_G.C.Greubel_,2024年2月24日
%Y参考A000040、A000178、A002808、A018252、A202808,A203417、A203418和A203419。
%K非n
%O 1,2号机组
%A _百灵鸟金伯利,2012年1月2日
|