%I#34 2015年8月28日17:48:15
%S 5,8,9,18,21,24,26,27,36,44,45,50,54,60,62,63,72,80,81,86,90108116,
%电话:117126132134135140144152153162170171180200204206207,
%电话:210216224225230234240242252260261264270306312314
%N数N,使得N可以被(10^k-数字和(N))整除,其中k等于数字和的位数(N)。
%这个序列是无限的,因为所有数字和等于9的数字都是这个序列的一部分。
%H Pieter Post,n表,n=1..12089的a(n)</a>
%e a(1)=5,因为5除以(10-5)等于1。
%e a(7)=26,因为数字(26)=8,26除以(10-8)等于13。
%e a(20)=86,该序列的第一个成员,其中digitsum(n)>=10。数字和(86)=14,所以k=10^2-14=86,所以86是这个序列的成员。
%t fQ[n_]:=块[{d=Total@IntegerDigits@n,k},k=IntegerLength@d;
%t可分[n,10^k-d]];选择[Range@314,fQ](*或*)
%t选择[Range@314,Divisible[#,(10^(Floor[Log[10,Total@IntegerDigits@#]]+1)-Total@IntegerDigits@#)]&](*_Michael De Vlieger_,2015年8月5日*)
%o(Python)
%o定义草皮(n,m):
%o。。。。kk=0
%o。。。。当n>0时:
%o。。。。。。。。kk=kk+(n%m)
%o。。。。。。。。n=整数(n//m)
%o。。。。返回kk
%o对于范围(1,10**6)中的c:
%o。。。。k=长度(str(sod(c,10))
%o。。。。kl=10**k-sod(c,10)
%o。。。。如果c%kl==0:
%o。。。。。。。。打印(c)
%o(PARI)isok(n)=我的(sd=和数字(n),nsd=数字(sd));n%(10^nsd-sd)==0;\\_米歇尔·马库斯(Michel Marcus),2015年8月5日
%Y参考A005349、A007953、A113315。
%K nonn,基数,减去
%O 1,1号机组
%A _邮政信箱_,2015年7月23日
|