Znn3bq.jpeg
²é¿´: 1370  |  »Ø¸´: 8
±¾Ìû²úÉú 1 ¸ö ³ÌÐòÇ¿Ìû £¬µã»÷ÕâÀï½øÐв鿴
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

wangww2011

ľ³æ (ÖøÃûдÊÖ)

[½»Á÷] Project Euler 46 Å·À­¹¤³Ì 46 Ìâ ÒÑÓÐ4È˲ÎÓë

µÂ¹úÊýѧ¼ÒChristian GoldbachÔø¾­Ìá³öÒ»¸ö²ÂÏ룺
ÈκÎÒ»¸öÆæºÏÊý¶¼ÄÜд³ÉÒ»¸öËØÊýÓëÒ»¸öƽ·½ÊýµÄ¶þ±¶µÄºÍ£¬Èç
9 = 7 + 2*1^2
15 = 7 + 2*2^2
21 = 3 + 2*3^2
25 = 7 + 2*3^2
27 = 19 + 2*2^2
33 = 31 + 2*1^2
µ«ÊǺóÀ´Ö¤ÊµÕâ¸ö²ÂÏëÊÇ´íÎóµÄ¡£
ÄÇÃ´ÍÆ·­Õâ¸ö²ÂÏëµÄ×îСµÄÆæºÏÊýÊǶàÉÙ£¿
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wangww2011

ľ³æ (ÖøÃûдÊÖ)

¡ï ¡ï
ben_ladeng: 2011-09-05 20:57:16
xzhdty(½ð±Ò+2): »¶Ó­³£À´ÌÖÂÛ 2011-09-06 12:57:05
ÁíÍâÒ»ÖÖ˼·£¬
CODE:
from math import sqrt
n=10000
primes=[2,3]
oddcomp=[]
for i in range(5,n,2):
    isprime=True
    for j in primes:
        if j*j>=i: break
        if i%j==0:
            isprime=False
            break
    if isprime:
        primes.append(i)
    else:
        oddcomp.append(i)
   
res=[p+2*j*j for p in primes for j in range(1,int(sqrt(n))+1)]
print set(oddcomp)-(set(res)|set(primes))

½á¹ûΪ
CODE:
set([5777, 5993])

PS ²»»áÖ»ÓÐÕâÁ½¸öÊý°É£¿

[ Last edited by wangww2011 on 2011-9-5 at 12:36 ]
7Â¥2011-09-05 12:28:26
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 9 ¸ö»Ø´ð

tieer

ľ³æ (ÕýʽдÊÖ)

¡ï ¡ï ¡ï
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû
xzhdty: »¶Ó­ÌÖÂÛ 2011-09-05 07:10:30
΢³¾¡¢ÃÎÏë(½ð±Ò+2): 2011-09-07 11:57:14
³õѧpython£¬²ËÄñÓ×Öɰ棬ºÃÏñÔËÐÐÁË3,4·ÖÖÓ£¬ÆÚ´ý¸ßÊָĽø
´ð°¸£¬5777£¬²»ÖªµÀÊÇ·ñÕýÈ·£¬ºÇºÇ
CODE:
# -*- coding: cp936 -*-
#Å·À­¹¤³Ì 46 Ìâ
#ÈκÎÒ»¸öÆæºÏÊý¶¼ÄÜд³ÉÒ»¸öËØÊýÓëÒ»¸öƽ·½ÊýµÄ¶þ±¶µÄºÍ,ѰÕÒÆÆÀý
#n=prime+2*m**2
from math import sqrt
def isprime(p):    #ÑéÖ¤ÊÇ·ñËØÊý£¬ËØÊý·µ»Ø±¾Éí£¬ºÏÊý·µ»ØFalse
    k=1
    for i in xrange(2,int(sqrt(p))+1):
        if p%i==0:
            k=0
            return False
            break
    if k:
        return p
n=35
while True:
    killer=1                   #ÉèÖòÂÏëµÄÅжϲÎÊý
    if not isprime(n):         #ÑéÖ¤ÊÇ·ñËØÊý£¬ËØÊýÔònµÝÔö¼ÌÐøÏÂһѭ»·£¬ºÏÊý½øÐÐÏÂÒ»²½ÑéÖ¤
        for m in xrange(1,int(sqrt((n-2)/2)+1)):
            for i in xrange(1,n-1):
                if isprime(i) and n==i+2*m**2:
                    killer=0   #·ûºÏ²ÂÏë
                    break
            if not killer:     #·ûºÏ²ÂÏë
                break
    else:
        n+=2
        continue
    if killer:                 #ÑéÖ¤Íê±Ï£¬²»·ûºÏ²ÂÏ룬¼´Îª½á¹û
        print('the number is:%d')%n
        break
    else:
        n+=2

[ Last edited by tieer on 2011-9-5 at 00:01 ]
˼¿¼£¬ÈÃÕâ¸öÊÀ½ç¸üÓÐȤ¡£
2Â¥2011-09-04 23:56:16
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

asaka

Òø³æ (³õÈëÎÄ̳)

¡ï ¡ï ¡ï
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû
΢³¾¡¢ÃÎÏë(½ð±Ò+2): 2011-09-07 11:57:24
ÎÒµÄpython°æ±¾
CODE:
from math import sqrt
primes = []
n = 3
while True :
  isp = True
  for j in primes :
    if j**2 > n : break
    if n%j == 0 :
      isp = False
      break
  if isp :
    primes.append(n)
  else :
    isf = True
    for j in reversed(primes) :
      it = (n-j)/2
      if it == int(sqrt(it))**2 :
        isf = False
        break
    if isf : break
  n = n + 2
print "The number is",n

ÔËÐнá¹û
The number is 5777
real    0m0.038s

[ Last edited by asaka on 2011-9-5 at 04:37 ]
3Â¥2011-09-05 03:26:21
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

asaka

Òø³æ (³õÈëÎÄ̳)

¡ï ¡ï ¡ï
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû
xzhdty(½ð±Ò+2): »¶Ó­³£À´ 2011-09-05 07:10:04
Áí¸½£ºÑÛ»¨çÔÂÒÖ®---FORTRAN goto °æ
CODE:
      integer p(800)
      j=0;n=1
   1  n=n+2
      do 2 i=1,j
        if(p(i)**2.gt.n) goto 3
   2    if(mod(n,p(i)).eq.0) goto 4
   3  j=j+1
      p(j)=n;goto 1
   4  do 5 i=j,1,-1
        k=(n-p(i))/2
   5    if(k.eq.int(sqrt(float(k)))**2) goto 1
      write(6,*) "The number is ",n
      End

ÔËÐнá¹û£º
The number is         5777
real    0m0.003s
4Â¥2011-09-05 04:10:13
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[Óлú½»Á÷] ¶þ±½¼×ͪËáÀàÑÜÉúÎï 50+3 С°×°®Ö÷ÈË 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[¿¼ÑÐ] Ò»Ö¾Ô¸Öпƴó²ÄÁÏÓ뻯¹¤£¬353·Ö»¹Óе÷¼ÁѧУÂð +10 ·ñ¼«Ì©À´2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[½Ìʦ֮¼Ò] ɽ¶«Ë«·ÇԺУ¿¼ºË³¬¼¶ÎÞµ×Ïߣ¬Áìµ¼ÐÒÔÖÀÖ»ö£¬½ÌʦÔâÑê¿Ö +4 qut2026 2026-04-11 8/400 2026-04-17 16:10 by »á·ÉµÄÖí157
[ÂÛÎÄͶ¸å] ÓÐûÓнÓÊձȽϿìµÄsciÆÚ¿¯Ñ½£¬×îºÃÔÚÒ»¸öÔÂÖ®Äڵģ¬ÑÐÈýº¢×ÓÇó±ÏÒµ 20+4 Ö®»¤×Å 2026-04-16 5/250 2026-04-17 10:02 by bobvan
[¿¼ÑÐ] 297£¬¹¤¿Æµ÷¼Á? +4 ºÓÄÏũҵ´óѧ-ÄÜ 2026-04-14 4/200 2026-04-16 22:52 by wulijun2012
[¿¼ÑÐ] 279ѧ˶ʳƷרҵÇóµ÷¼ÁԺУ 20+7 ¹Â¶ÀµÄÀǰ®³ÔÑò 2026-04-12 29/1450 2026-04-16 09:00 by screening
[¿¼ÑÐ] Çóµ÷¼ÁÍÆ¼ö +8 СÄô°®Ñ§Ï° 2026-04-14 8/400 2026-04-16 07:22 by ѧԱJpLReM
[¿¼ÑÐ] ÇóÖúµ÷¼Á£¬¿çµ÷ +19 XÊ®¸¦´çY 2026-04-11 20/1000 2026-04-15 21:18 by cuisz
[¿¼ÑÐ] 085404 22408 309·ÖÇóµ÷¼Á +9 lzmk 2026-04-14 10/500 2026-04-15 20:02 by ѧԱJpLReM
[¿¼ÑÐ] ҩѧÇóµ÷¼Á +11 RussHu 2026-04-12 13/650 2026-04-15 19:07 by zhuwenxu
[¿¼ÑÐ] 297¹¤¿Æµ÷¼Á? +14 ºÓÄÏũҵ´óѧ-ÄÜ 2026-04-13 15/750 2026-04-15 13:25 by ºÚ¿Æ¼¼¿óÒµ
[¿¼ÑÐ] ҩѧ305Çóµ÷¼Á +7 Â꿨°Í¿¨boom 2026-04-11 7/350 2026-04-15 13:21 by Î÷±±Íû¡ª·çɳ
[¿¼ÑÐ] µ÷¼Á +12 ÔÂ@163.com 2026-04-11 12/600 2026-04-14 15:37 by zs92450
[¿¼ÑÐ] 105500ҩѧÇóµ÷¼Á +4 x_skys 2026-04-12 4/200 2026-04-14 13:37 by rndfc
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤329·ÖÇóµ÷¼Á +24 Ò¶zilin 2026-04-13 25/1250 2026-04-14 09:20 by ÊÔ¹ÜÆÆÁÑ
[¿¼ÑÐ] 293Çóµ÷¼Á +16 ÎÒ°®¸ßÊý¸ßÊý°®Î 2026-04-12 18/900 2026-04-13 21:47 by ѧԱJpLReM
[¿¼ÑÐ] Ò»Ö¾Ô¸085802 323·ÖÇóµ÷¼Á +13 drizzle_9 2026-04-12 14/700 2026-04-13 10:26 by Faiz5552
[¿¼ÑÐ] 346·Ö£¬¹¤¿Æ0854Çóµ÷¼Á£¬×¨Ë¶ +6 moser233 2026-04-12 7/350 2026-04-12 22:11 by fqwang
[¿¼ÑÐ] 331Çóµ÷¼Á +5 Íõ¹ú˧ 2026-04-11 5/250 2026-04-11 22:56 by Ϫ½§Á÷Ë®
[¿¼ÑÐ] Ò»Ö¾Ô¸985»úеѧ˶380Çóµ÷¼Á +5 ¹Ø¹ØöÂð¯10 2026-04-11 5/250 2026-04-11 10:10 by ÖªÄî¡£A
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û