登录
A008935号
如果2n=总和2^e(k),则a(n)=总和e(k”^2。
6
1, 4, 5, 9, 10, 13, 14, 16, 17, 20, 21, 25, 26, 29, 30, 25, 26, 29, 30, 34, 35, 38, 39, 41, 42, 45, 46, 50, 51, 54, 55, 36, 37, 40, 41, 45, 46, 49, 50, 52, 53, 56, 57, 61, 62, 65, 66, 61, 62, 65, 66, 70, 71, 74, 75, 77, 78, 81, 82, 86, 87, 90, 91, 49, 50, 53, 54, 58, 59, 62
抵消
1,2
链接
配方奶粉
一般公式:1/(1-x)*Sum_{k>=0}(k+1)^2*x^2^k/(1+x^2*k)。 -拉尔夫·斯蒂芬2003年6月23日
例子
为了得到a(5),我们写10=2+8=2^1+2^3,所以a(5”)=1^2+3^2=10。
MAPLE公司
a: =n->(l->加(l[i]*i^2,i=1..nops(l)))(转换(n,基数,2)):
seq(a(n),n=1..80); #阿洛伊斯·海因茨2019年11月20日
数学
a[n_]:=总计[Flatten[Position[Reverse[Integer Digits[n,2]],1]]^2];表[a[n],{n,1,70}](*Jean-François Alcover 2011年3月21日*)
黄体脂酮素
(C)
#包括<stdio.h>
#包括<stdlib.h>
#define USAGE“用法:”A008935号num’,其中num是序列中所需结束值的索引。\n“
#定义MAX 1023
#定义SHIFT_MAX 9
int main(int argc,char*argv[]){无符号短掩码,i,j,end;无符号长和;if(argc<2){fprintf(stderr,USAGE);return EXIT_FAILURE;}end=atoi(argv[1]);end=(end>=最大值)?MAX:结束;
fprintf(标准输出,“值:”);对于(i=1;i<=end;i++){sum=0;mask=1;对于(j=0;j<SHIFT_MAX;j++,mask*=2)if(i&mask)sum+=(j+1)*(j+1;返回EXIT_SUCCESS; }
(哈斯克尔)
a008935=f 1,其中
f k x | x==0=0
|r==0=f(k+1)x'
|否则=k^2+f(k+1)x',其中(x',r)=divMod x 2
(Python)
a=λn:总和(((k+1)**2)*((n>>k)&1)范围内的k(0,n.bit_length()))
打印([a(n)代表范围(1,68)中的n])#达里奥·克拉维乔2024年12月27日
交叉参考
给予A003995号如果排序并删除重复项。
关键词
非n,美好的,容易的
作者
扩展
由Larry Reeves(larryr(AT)acm.org)修订和扩展,2000年3月22日
状态
经核准的