lt

lt

文章
389
评论
1641
推荐
260
收藏
3
社区会龄
8 年
个人网站
--
个人简介
还没有输入个人简介……

发表了

评论了

  • lt
    今年 06-21 09:39
    lt 评论了文章 欧拉计划721:无理数高次幂

    快速幂取模算法 https://www.cnblogs.com/wkfvawl/p/9125224.html

  • lt
    今年 06-21 09:25
    lt 评论了文章 欧拉计划721:无理数高次幂

    测试代码 import math def f(a,n): return math.floor((math.ceil(math.sqrt(a))+math.sqrt(a))**n) >>> f(5,2) 27 >>> f(5,5) 3935 >>> def G(n): ... s=0 ... for a in range(1,n+1) :s+=f(a,a*a) ... return s ... >>> G(10) 31997447160056074891755545250568015145391367309910918290177259841659637836252530747180

  • lt
    今年 06-11 16:02
    lt 评论了文章 欧拉计划719:拆分数

    最后的递归程序 在麒麟980的m6平板上，用termux的python 3.8.3 84秒 在i3 4010上 python 3.8.3 149秒 pypy 7.3 28秒

  • lt
    今年 06-11 15:36
    lt 评论了文章 欧拉计划719:拆分数

    import time def t(n,x): #print(n,x,-1) if n==int(x):return 1 if n>int(x):return 0 for i in range(1,len(x)): r= t(n-int(x[-i:]),x[:-i]) if r==1: return r return 0 x='123' y=x[:-2] z=x[-2:] #print(z,y) s=0 n=10**6 tm=time.time() for i in range(9,n+1,9): for j in [0,1]: i+=j z=t(i,str(i**2))*i**2 #print(z) s+=z print(time.time()-tm) print(s)

  • lt
    今年 06-07 21:28
    lt 评论了文章 欧拉计划719:拆分数

    拆成7个数只有大于10w才需要 import time def test(N): t=time.time() s=0 for i in range(1,N+1): if (adds(i)%9==0 or adds(i)%9==1) and (is_sn(i)==1 or is_sn2(i)==1 or is_sn3(i)==1 or is_sn4(i)==1 or (i>10**4 and is_sn5(i)==1) or (i>10**5 and is_sn6(i)==1)): s+=i*i #print(i,i*i,is_sn(i),is_sn2(i),is_sn3(i),is_sn4(i),is_sn5(i),is_sn6(i),adds(i)) print(time.time()-t) return s