登录
OEIS由OEIS基金会的许多慷慨捐赠者.

 

标志
提示
(来自的问候整数序列在线百科全书!)
A068571号 快乐数字的数量<=10^n。 5
1, 3, 20, 143, 1442, 14377, 143071, 1418854, 14255667, 145674808, 1492609148, 15091199357, 149121303586, 1443278000870, 13770853279685, 130660965862333, 1245219117260664, 12024696404768025, 118226055080025491, 1183229962059381238, 12005034444292997294 (列表;图表;参考;;历史;文本;内部格式)
抵消
0,2
链接
布莱恩·沃尔夫,n=0..1000时的n,a(n)表(前122个术语来自Lambert Klasen)
贾斯汀·吉尔默,关于快乐数的密度,arXiv:1110.3836[math.NT],2011-2015。
Lambert Klasen,Xeen3d的happy.html.
沃尔特·施耐德,快乐数字.
埃里克·魏斯坦的数学世界,快乐数字.
配方奶粉
发件人布莱恩·沃尔夫2012年11月6日:(开始)
h(n,x)=h(n-1,x-0^2)+h(n-1,x-1^2)=…+h(n-1,x-9^2),其中h()是长度为n的数字的个数,其数字的平方和为x。
a(n)是所有h(n,0<x<=81*n)的总和,其中x是一个快乐数,加上1代表10^n。(结束)
例子
对于n=0,h(0,0)=1,h(0,x>0)=0。
黄体脂酮素
(PARI)ssd(n)=n=数字(n);总和(i=1,#n,n[i]^2)
快乐(n)=while(n>6,n=ssd(n));n==1
a(n)=我的(f=n!,s,d);forvec(v=向量(9,i,[0,n]),d=向量(9,i,如果(i>8,n,v[i+1])-v[i]);如果(快乐(总和(i=1,9,d[i]*i^2)),s+=f/prod(i=1,9,d[i]!)/v[1]!),1); 秒+1\\查尔斯·格里特豪斯四世2016年11月1日
交叉参考
囊性纤维变性。A007770号,A124095型,A219667型.
关键词
非n,基础
作者
萨沙·库尔兹2002年3月26日
扩展
更多术语来自Lambert Klasen(Lambert.Klasen,AT)gmx.de),2004年11月3日
状态
经核准的

查找|欢迎光临|维基|注册|音乐|地块2|演示|索引|浏览|更多|网络摄像头
贡献新序列。或评论|格式|样式表|变换|超级搜索|最近
OEIS社区|维护人OEIS基金会。

许可协议、使用条款、隐私政策。.

最后修改时间:美国东部时间2024年6月22日01:41。包含373561个序列。(在oeis4上运行。)