%I#12 2021年7月31日22:17:06
%S 3,18,33,48,83,981131631782432582732883353418513528593,
%电话:627642657707722768788289796211371251126612981313,
%电话:132813311378139314581506155315681633180818751922193720022177240324182433248324982546256325932682658
%N个数字是三个四分之一次幂的和
%C第60项与A003337和A047714不同,因为2673=2^4+4^4+7^4=3^4+6^4+6 ^4,见A309762。
%H David Consiglio,Jr.,<a href=“/A34188/b344188.txt”>n表,n=1..20000的a(n)</a>
%e 33是这个序列的成员,因为33=1^4+2^4+2 ^4
%o(Python)
%o从itertools导入combinations_with_replacement作为cwr
%o从集合导入defaultdict
%o keep=默认dict(λ:0)
%o power_terms=[x**4表示x在范围(1,50)内]
%o对于cwr中的pos(power_terms,3):
%o总计=总和(pos)
%o保持[tot]+=1
%o rets=已排序([k代表k,v在keep.items()中,如果v==1])
%o表示x在范围内(len(rets)):
%o打印(rets[x])
%Y参见A003337、A025395、A344187、A344189、A344 192、A34 4641。
%K nonn公司
%O 1,1号机组
%2021年5月11日,小阿维德·Consiglio
|