%I#29 2023年2月24日18:37:08
%S 1,3,9,27,3710130390924391011098912987158732564127027,
%电话:303033703747619769239090910101136986330303039090909,
%电话:1234567927027370370101010101243902439303030390909091010101109890109891298701298715873015873
%N正整数N(具有k位数字),如果具有k+1位数字的正整数m可以被N整除,那么m的所有旋转都可以被N除尽。
%C·约翰·D·库克的博客(见下面的链接)证明了“如果一个三位数可以被37整除,那么如果你旋转它的数字,它仍然可以被37除尽。”
%H John D.Cook,<a href=“https://www.johndcook.com/blog/2023/02/12/rotating-multiples-of-37/“>37的旋转倍数。
%e对于a(4)=27,405是27的3位数倍数,405的两次旋转(即54和540)也是27的倍数。
%e对于a(5)=37,185是37的3位数倍数,185的两次旋转(即851和518)也是37的倍数。
%e对于a(9)=2439,12195是2439的5位数倍数,12195的四次旋转(即21951、19512、95121和51219)也是2439的倍数。
%o(Python)
%o定义旋转(str):
%o first_char=字符串[0:1]
%o remaining_chars=str[1:]
%o返回(remaining_chars+first_char)
%o定义get_rrotations(n):
%o n_as_str=字符串(n)
%o旋转=[]
%o旋转_ as_str=旋转_ as_str
%o对于范围内的i(len(n_as_str)-1):
%o rotation_as_str=旋转(rotationaas_str)
%o rotations.append(int(rotation_as_str))
%o返回旋转
%o序列=[]
%o最大值=9999999
%o表示范围(1,max_n+1)内的n:
%o n_len=长度(str(n))
%o系数=2
%o为True时:
%o prod=n*系数
%o prod_len=长度(str(prod))
%o如果prod_len<n_len+1:
%o系数=系数+1
%o elif prod_len>n_len+1:
%o顺序追加(n)
%o中断
%o其他:
%o#prod_len==n_len+1
%o旋转=获取旋转(prod)
%o如果全部(旋转%n==0表示旋转):
%o系数=系数+1
%o其他:
%o中断
%o打印(seq)
%Y参考A034089、A066484。
%K nonn,基础
%O 1,2号机组
%A _罗贝尔特·里昂,2023年2月14日
%E a(25)-a(35),来自恰瓦乌,2023年2月24日
|