解析数论导轮中的数学实验(python)

最新推荐文章于 2024-03-24 11:29:53 发布

buck

最新推荐文章于 2024-03-24 11:29:53 发布

阅读量510

点赞数

分类专栏：解析数论导论数学文章标签： python 数学

本文链接：https://blog.csdn.net/buck84/article/details/65938009

版权

数学同时被 2 个专栏收录

17 篇文章 2 订阅

订阅专栏

解析数论导论

4 篇文章 0 订阅

订阅专栏

前置代码

import math

def isprime(n):
    sqrtn = math.ceil(math.sqrt(n))
    for i in range(2, int(sqrtn)):
        if n % i == 0:
            return False
    return True

def primes(n):
    ret = [2]
    for i in range(3, n+1):
        if isprime(i):
            ret.append(i)
    return ret
print("calc primes")
ps = primes(100000)
print("calc primes end")

第2章 Arithmetical Functions and Dirichlet Multiplication

这里写代码片

第3章 Averages of Arithmetical Functions

3.11


for i in x:
    lx = 0
    for j in range(1,i+1):
        lx = lx + math.log(j)

    mx = 0
    for j in range(1,i+1):
        mx = mx + lambdaFunc(j)*math.floor(i/j)
    print "info", i, lx, mx