一个神经病的日志

摘要： 在Mathematica里，内置函数一般拥有属性Protected，这意味着你不能随便内置函数的定义进行修改。当然，对于绝大部分函数（除去那些有Locked属性的）都可以Unprotect了之。但是这样做是极为危险的，有时你都忘记了自己对内置函数做了哪些修改，虽然这些修改在重启Kernel之后都会恢复，但还是不建议这样做。今天介绍一种暂时修改内置函数的方法，这个方法用到了没有文档的内部函数，Internal`InheritedBlock。此函数的用法和Block类似，它能将内置函数复制一份并进行任意的修改且不影响源函数。举个例子来介绍其使用方法。假如我们希望编写一个函数WithOrderedP 阅读全文

摘要： 找出200万以下的质数之和 def primeQ(n): for i in range(2, int (n**0.5) + 1): if n % i == 0: return False return Trueprint sum(filter(primeQ, range(2, 2000001)))>>>142913828922惨不忍睹！>< 阅读全文

摘要： 已知存在并且只存在一个毕达哥拉斯三元组满足条件a + b + c = 1000。 找出该三元组中abc的乘积。 for c in range(1, 500): for b in range(1, c): if (1000 - b - c)**2 + b**2 == c**2: print (1000 - b - c) * b * c ... 阅读全文

摘要： 找出以下这个1000位的整数中连续5个数字的最大乘积。（例如前五个数字的乘积是7*3*1*6*7=882) string = '''7316717653133062491922511967442657474235534919493496983520312774506326239578318016984801869478851843858615607891129494954595017379... 阅读全文

摘要： 第10001个质数是多少？ def primeQ(n): for i in range(2, int (n**0.5) + 1): if n % i == 0: return False return Truex = 1n = 1while n < 10001: x = x + 2 if primeQ(x): n =... 阅读全文

摘要： 找出前一百个自然数的平方和与和平方的差。 rang = range(1, 101)print (sum(rang))**2 - sum(map(lambda x: x**2, rang))>>>25164150map和lambda即为Mathematica中的Map和Function。 阅读全文

摘要： 2520是最小的能被1-10中每个数字整除的正整数。 最小的能被1-20中每个数整除的正整数是多少？ def gcd(a, b): product = a * b while b != 0: t = b b = a % b a = t return product / atemp = gcd(1, 2)for i in rang... 阅读全文

摘要： 找出最大的有由个三位数乘积构成的回文数。 def huiwenQ(num): string = str(num) if string == string[::-1]: return True else: return Falseresult = 0for i in range(999, 99, -1): for j in range(99... 阅读全文

摘要： 13195的质数因子有5,7,13和29. 600851475143的最大质数因子是多少？ #求num的最大质因数num = 600851475143sqrt = int (num**0.5)def primeQ(n): for i in range(2, int (n**0.5) + 1): if n % i == 0: return Fals... 阅读全文