登录
满足0<a<b<=n的1到有理部分a/b的分区数。
2

%I#5 2012年3月30日18:58:15

%S 0,1,3,6,12,19,33,50,771111662293294486115824112014671956,

%电话253833134258551269938944112781425817838224022778834629,

%电话:42747528326484479706972723211886814439217547212168256750

%N将1划分为满足0<a<b<=N的有理部分a/b的数量。

%F a(n)=a(n-1)-1+(n划分为相对素部分的次数)。

%e a(1)=0。

%e a(2)=1表示1/2+1/2。

%e a(3)=3计数1/2+1/2、1/3+1/3+1/3和1/3+2/3。

%t p[n_]:=p[n]=整数分区[n];

%t l[n_]:=长度[p[n]];

%t p[n_,j_]:=p[n,j]=零件[p[n],j]

%t g[n_,j_]:=g[n,j]=应用[GCD,p[n,j]]

%t h[n_]:=h[n]=表[g[n,j],{j,1,l[n]}]

%t c[n_]:=c[n]=计数[h[n],1]

%t表[c[n],{n,0,45}](*A000837*)

%t s[n]:=总和[c[k],{k,1,n}]

%t表[s[n]-1,{n,1,45}](*A209489*)

%Y参考A000837。

%K nonn公司

%氧1,3

%A_Clark Kimberling_,2012年3月9日