lt

  • 关注
  • 短消息
文章
404
评论
1672
推荐
270
收藏
3
社区会龄
9 年
个人网站
--
个人简介
还没有输入个人简介……
  • lt 推荐

    欧拉计划751:串联重合

    一个实数c通过以下过程,可以从任何正实值生成非递减整数序列: b[1]=c b[n]=floor(b[n-1])*(floor(b[n-1])-b[n] + 1),对n>=2 a[n]=floor(b[n]) 例如,c=2.956938891377988...生成斐…...

  • lt 1推荐

    欧拉计划749:近似幂和

    一个正整数n,其10进制表达的每位的正整数k次方的和,如果等于n-1或n+1,称为近似幂和。比如,35是近似幂和,因为3^2 + 5^2 =34。 设s(n)为小于或等于n位的近似幂和的总和。 已知s(2)=100,s(6)=2562701。 求s(16) 测试程序 imp…...

  • lt 1推荐

    数论概论中的费马降阶法求一个质数分解成2个数的平方和与穷举法的效率比较

    源代码改写自《HDU 3542 (费马降阶法)》原文链接:https://blog.csdn.net/Feynman1999/article/details/81903899 import random as rd import time def ran(n): if (n…...

  • lt 1推荐

    欧拉计划748:倒丢蕃多方程

    对于丢蕃多方程的倒写 1/x2 +1/y2=13/z2 设(x,y,z)是一组正整数解,其中gcd(x,y,z)=1,为原始解, 设S(n)是满足1<=x,y,z<=n,x<=y的原始解之和。 对于n=100,有2组原始解,各是(2,3,6)和(5,90,18…...

  • lt 1推荐

    欧拉计划747:三角披萨

    https://projecteuler.net/problem=747 妈妈烤了一个三角披萨。她想把披萨饼切成小块。她首先在三角形披萨的内部(而不是边界)选择一个点P,然后进行切割,所有切割都从P开始并直接延伸到披萨的边界,这样所有的片都是三角形,并且具有相同的面积。 让我…...

评论了

  • 今年 03-14 19:18

    把浮点数改成整数,通过 Congratulations, the answer you gave to problem 751 is correct. You are the 172nd person to have solved this problem. Return to Problems page.

  • 今年 03-14 07:50

    如果把232改成100怎么解?

  • 今年 03-13 18:22

    利用https://www.ituring.com.cn/article/516443中的代码,用费马降阶法求w*13对应的x和y >>>> c1=0 >>>> t=time.time();print(f6(2*10**7,13));print(time.time()-t,c1) 7312375200 0.37439990043640137 839 >>>> c1=0 >>>> t=time.time();print(f7(2*10**7,13));print(time.time()-t,c1) 7312375200 0.296400785446167 445

  • 今年 03-11 20:42

    求任意多个4n+1型质数之积的分解 def e(l,s2): l2=[] for s1 in l: u,v=s1[0],s1[1] a,b=s2[0],s2[1] l2.append([abs(u*a+v*b),abs(v*a-u*b)]) l2.append([abs(u*a-v*b),abs(v*a+u*b)]) return l2 def test(pl): sl=[solve(p) for p in pl] l1=[sl[0]] for i in range(1,len(sl)): l1=e(l1,sl[i]) return l1 >>> test([5,13,17]) [[12, 31], [4, 33], [32, 9], [24, 23]] >>> test([5,13,17,29]) [[179, 2], [131, 122], [173, 46], [157, 86], [109, 142], [19, 178], [163, 74], [67, 166]] >>> test([5,13]) [[8, 1], [4, 7]]

  • 今年 02-28 18:28

    第138个,把k的循环改为 for k in range(int(math.log((10**d) /d,9)),int(math.log(10**d,2))+1): #次方