2520 is the smallest number that can be divided by each of the numbers from 1 to 10 without any remainder. What is the smallest positive number that is evenly divisible by all of the numbers from 1 to 20?
计算最小的能被1-20所有数整除的数,换个说法就是求{2-20}的最小公倍数。
根据题目给出的内容,2520是{2-10}的最小公倍数,那么{2-20}的最小公倍数应当可以表示为2520X。
其次,可以由因数分解算法得出,2520=(2^3)*(3^2)*5*7。整除2520的数同时整除了2*2*3=12, 2*2*5=20, 2*7=14, 3*5=15, 3*3*2=18.
所以,在求取{2-20}的最小公倍数2520X,可以不再考虑{1-10, 12, 14, 15, 18, 20}的倍数,仅仅测试{11, 13, 16, 17, 19}
def isOK(n): for i in [11, 13, 16, 17, 19]: if n % i: return False return True def smallestDivisibleNumber(): i = 2520 while not isOK(i): i += 2520 return i
这个解法需要先用另一个算法求出将2520分解,不好不好。
设N为2到K最小公倍数,根据“任何数可以分解为有限个素数的乘积”,N的因子应当覆盖所有2-K范围内的素数,N为这些数的乘积。
例如:
K=2, N=2
K=3, N=2*3=6 # 2,3的最小公倍数是6
K=4, N=2*3*2=12 # 2,3,4的最小公倍数是12, 4=2*2.
K=5, N=2*3*2*5=60 # 2,3,4,5的最小公倍数是60
k=6, N=2*3*2*5=60 # 2,3,4,5,6的最小公倍数是60, 6=2*3
……
k=20, N=2*3*2*5*7*2*3*11*13*2*17*19
用手就算出来了~~~
总结规律:
1、N必须被小于N的所有素数整除。
2、N必须被小于N的所有合数整除。
3、这些合数可以以小于N的素数的指数形式表示。
如上:
K=4, N=(2^2)*3
K=5, N=(2^2)*3*5
……
K=60, N=(2^4)*(3^2)*5*7*11*13*17*19
4、计算这些素数的指数部分时,对于20来说:
(1)、2^4=15, 2^5=25
(2)、3^2=9, 3^3=27
(3)、5^2=25, 7^2=25 ……
将素数表示为p, 其指数部分为a,则p^a <= K。可由公式log(p^a)=a*log(p) = K, a = log(k) / log(p), 对a进行地板取整便是合适的值。而对于那些大于sqrt(K)的素数,它们的指数部分必定为1,因为(sqrt(p)^2)>K了。这样求出的公倍数就不只是2-K的了,还包含了超出K的部分。
def SmallestDivisibleNum(): k = 20 n = i = 1 flag = True limit = math.sqrt(k) for p in [2, 3, 5, 7, 11, 13, 17, 19]: a = 1 if flag: if pi <= limit: a = math.floor(math.log(k)/math.log(p)) else: flag = False n *= p**a return n