%I#24 2024年2月14日16:27:57
%S 7,16,0,16,25,0,7,16,10,70,0,0,52,70,016,34,0133160,0,0115,0,25,52,
%电话:0700232,0124160,0,34,70,0304,0160205,0,43,0,03221222,0,
%U 2401250,0,52106,0,07000,02324012,0,61124,01024520,01122340单位
%N数字和为7的N的最小倍数,如果不存在这样的数字,则为0,例如a(3k)=a(11k)=0。
%C 37的倍数没有数字和7,或者它们的第一个值大于10000000拉里·里维斯(larryr(AT)acm.org),2002年7月2日
%C没有3、11、37、101或271的倍数有数字和7_罗伯特·伊斯雷尔,2024年2月13日
%H Robert Israel,n的表,a(n)表示n=1..10000</a>
%F a(n)=n*A088396(n).-_R.J.Mathar,2019年8月6日
%p未完成:=真:V:=向量(300):V0:=选择(t->igcd(t,3*11*37*101*271)=1,{$1..300}):
%未完成时从0中为i1指定p
%p表示i2从0到i1,而未完成do
%p表示i3从0到i2,而未完成do
%i4从0到i3的p,而未完成do
%p表示i5从0到i4,而未完成do
%p表示i6从0到i5,而未完成do
%p表示i7从0到i6,而未完成do
%p v:=10^i1+10^i2+10^i 3+10^i 4+10^ i 5+10^ii 6+10^il 7;
%p dv:=数字理论:-除数(v);
%V0中s的p与dv do相交
%p V[s]:=V;
%p od;
%p V0:=V0减去dv;
%p未完成:=evalb(V0<>{});
%日期od od od od od:
%p转换(V,列表);#_罗伯特·伊斯雷尔,2024年2月13日
%Y参考A052221、A069521、A069522、A069523、A069524、A069525、A069527、A069528。
%K nonn,看,基础
%O 1,1号机组
%A _Amarnath Murthy,2002年4月1日
%E Larry Reeves提供的更多条款,2002年7月2日
|