%I#12 2024年2月24日11:03:52

%编号：1、2、1624011520132710402548039680001589212348416000，

%电话：51261633510506496000000892887435278044668870400000，

%U 50689719717698351557731837542400000000 1257651788315794213051651266651252326400000000

%前N个复合数的N Vandermonde行列式（A002808）。

%C如A203419所示，每个术语都可以划分其后继术语，如A203420所示，每一术语都可以被相应的超因子A000178（n）划分。

%H G.C.Greubel，n表，n=1..42的a（n）</a>

%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}]（*此序列*）

%t表格[v[n+1]/v[n]，{n，z}]（*A203419*）

%t表格[v[n]/d[n]，{n，z}]（*A203420*）

%o（岩浆）

%o A002808:=[2..250]|非IsPrime（n）]中的n：n；

%o a:=func<n|n eq 0选择1 else（&*[（&*[A002808[k+2]-A002808[j+1]：[0..k]]中的j）：[0..n-1]中的k）>；

%o[0..20]]中的[a（n）:n；//_G.C.Greubel，2024年2月24日

%o（SageMath）

%o A002808=[n代表（2..250）中的n，如果不是_prime（n）]

%o定义a（n）：返回产品（范围（k+1）中j的产品（A002808[k+1]-A002808[j]），范围（n）中k的产品）

%o[a（n）for n in range（15）]#_G.C.Greubel_，2024年2月24日

%Y参考A000040、A000178、A002808、A203419、A203420。

%K nonn公司

%O 1、2

%A _百灵鸟金伯利，2012年1月2日