%I#48 2022年2月18日23:08:49
%S 1,2,6,2,30,2,42,2,302,26,22730,2,6,2510,2798,2330,2138,22730,
%电话:2,6,2870,214322,2510,2,6,219190,2,6,13530,21806,2690,2282,
%U 246410,2,66,21590,2798,2870,2354,256786730,2,6510,264722,2,30,24686
%N 1后面是A027760,这是伯努利数分母的变体。
%C对于n>0,伯努利数的分母。B_n序列开始于1、-1/2、1/6、0/2、-1/30、0/2,1/42、0/2。。。这是克劳森定理建议的A027642的替代版本_Peter Luschny_,2009年4月29日
%C设f(n,k)=gcd{多项式(n;n1,…,nk)|n1+…+nk=n};则a(n)=f(n,n-n+1)/f(n,n-n),对于n>>n.-Mamuka Jibladze_,2017年3月7日
%H Antti Karttunen,n表,n=0..10080的a(n)</a>
%H托马斯·克劳森,<a href=“http://adsabs.harvard.edu/abs/1840AN…..17R.351“>《Lehrsatz aus einer Abhandlung ut ber die Bernoullischen Zahlen》,《美国国家科学院学报》17(22)(1840),351-352。
%H维基百科,<a href=“http://en.wikipedia.org/wiki/Bernoulli_number“>伯努利数</a>
%F a(n)是由cosh(x*z)*z/(1-exp(-z))计算得出的x=1生成的多项式的分母。分子见A176328_Peter Luschny_,2018年8月18日
%F a(n)=分母(和{j=0..n}(-1)^(n-j)*j*斯特林2(n,j)*B(j)),其中B是伯努利数A164555/A027642_Fabián Pereyra,2022年1月6日
%e公式部分中示例f.给出的合理值开始:1,1/2,7/6,3/2,59/30,5/2,127/42,7/2,119/30,…-_Peter Luschny_,2018年8月18日
%p克劳森:=proc(n)局部S,i;
%p S:=数值[除数](n);S:=映射(i->i+1,S);
%p S:=选择(isprime,S);mul(i,i=S)结束进程:
%p seq(克劳森(i),i=0..24);
%p#_Peter Luschny_,2009年4月29日
%p A141056:=程序(n)
%p如果n=0,则1,否则A027760(n)结束;
%p end程序:#_R.J.Mathar_,2013年10月28日
%t a[n_]:=和[Boole[PrimeQ[d+1]]/(d+1),{d,除数[n]}]//分母;表[a[n],{n,0,70}](*Jean-François Alcover_,2012年8月9日*)
%o(PARI)
%o A141056(n)=
%o个{
%o p=1;
%o如果(n>0,
%o代表(n,d,
%o r=d+1;
%o如果(isprime(r),p=p*r)
%o)
%o);
%o返回(p)
%o}(o)
%o表示(n=0,70,打印1(A141056(n),“,”);/*_Peter Luschny_,2012年5月7日*/
%Y参见A027760、A027642、A176328。
%Y参见A164555、A027642、A048993。
%K nonn公司
%0、2
%2008年8月1日,A Paul Curtz
%E由R.J.Mathar_延长,2009年11月22日
|