登录
A071081号
n X n矩阵的行列式,其元素(i,j)等于第i个复合数,如果i=j,则为0。
1, 0, -16, 192, -1904, 16416, -134608, 1102920, -8971103, 69262338, -527129920, 4002967800, -30263030000, 218133853800, -1565386817920, 11130108480678, -75244171093875, 496516351214832, -3261752198331472, 21401161780748720, -140093238345715827, 914525302322457472
抵消
0,3
链接
罗伯特·伊斯雷尔,n=0..1000时的n,a(n)表
MAPLE公司
comps:=删除(isprime,[$4..11000]):
f: =proc(n)局部M;
M: =矩阵(n,n,(i,j)->`如果`(i=j,0,comps[abs(i-j)]);
线性代数:行列式(M)
结束进程:
f(0):=1:
地图(f,[0..25]美元); #罗伯特·伊斯雷尔,2024年12月2日
数学
复合[n_Integer]:=固定点[n+PrimePi[#]+1&,n+PricePi[n]+1];f[n_]:=Det[表[i==j,0,复合[Abs[i-j]],{i,1,n},{j,1,n}]];表[f[n],{n,1,20}]
黄体脂酮素
(Python)
从sympy导入矩阵,复合
定义A071081号(n) :如果j!=k,则返回矩阵(n,n,[组合(abs(j-k));如果j!=k,则为范围内j的0;如果范围内k的0,则为(n)]).det()#柴华武2024年7月1日
(PARI)a(n)=my(复合(n)=my(k=-1);而(-n+n+=-k+k=primepi(n),);n);matdet(矩阵(n,n,i,j,如果(i==j,0,复合(abs(i-j)))); \\路德·范托尔(Ruud H.G.van Tol)2024年7月14日
交叉参考
关键词
签名
作者
罗伯特·威尔逊v2002年5月26日
扩展
a(21)来自斯特凡诺·斯佩齐亚2024年6月27日
a(0)=1前面加阿洛伊斯·海因茨2024年7月1日
状态
经核准的