Seite 1 von 2
Mehr Hilfe gibt es nicht ...
Verfasst: Freitag 12. Januar 2007, 00:09
von sunmountain
Code: Alles auswählen
def sieb_des_eratosthenes(min,max):
zahlen = range(min,max)
for i in xrange(2,max):
entferne_vielfache(i,zahlen)
return zahlen
def entferne_vielfache(i,zahlen):
for j in xrange(i**2,zahlen[-1],i):
if j in zahlen:
zahlen.remove(j)
return zahlen
MIN=10
MAX=20
primzahlen = sieb_des_eratosthenes(MIN,MAX)
primepaare = [ (primzahlen[i],primzahlen[i+1]) for i in xrange(len(primzahlen)-1) if primzahlen[i]+2==primzahlen[i+1] ]
print "Primepaare von %i bis %i" % (MIN,MAX)
for p in primepaare:
print p[0]
print p[1]
print
Ich denke der Code ist so einfach,
das ich mir Kommentare spare.
[/code]
Verfasst: Freitag 12. Januar 2007, 08:48
von Luzandro
Cthulhu hat geschrieben:ich möchte jetzt bitte hilfe zu meinem quelltext haben und nicht wieder auf das sieb des was weiß ich was abgelenkt werden!
Und darf ich dich dennoch mit der Frage nach dem Grund nerven? Dein Code ist höchst umständlich und schwer verständlich geschrieben, weshalb du ja auch selber nicht weißt, was er richtig oder falsch tut. Die Hinweise sollten keine "Ablenkung" sein, sondern dich auf einen sinnvolleren und einfacheren Weg bringen
@sunmountain
Du gehst hier bei der ersten Funktion ALLE Werte von 2 bis max durch, obwohl ein Großteil davon selbst schon Vielfache von Primzahlen sein werden.
Code: Alles auswählen
def primzahl_zwillinge(a,b):
bereich = range(2,b)
primzahlen = []
while len(bereich) > 0:
primzahl = bereich.pop(0)
primzahlen.append(primzahl)
bereich = [ zahl for zahl in bereich if zahl % primzahl != 0 ]
return [ (x,y) for x in primzahlen for y in primzahlen if x >= a and x+2 == y ]
Auch nicht optimal, weil die Listcomprehension am Ende im Gegensatz zu deinem Code unnötigerweise 2 Loops enthält, aber ich finde es lest sich so irgendwie schön funktional...
edit:
wobei sich folgendes auch nicht schlechter liest
Code: Alles auswählen
[ (x,x+2) for x in primzahlen if x >= a and x+2 in primzahlen ]
3 Experten ...
Verfasst: Freitag 12. Januar 2007, 09:17
von sunmountain
... 4 Meinungen
.
Ich denke mal, jeder hat seine Vorlieben beim Codieren.
Ich bevorzuge den Weg, mir Gedanken über Redudanzen zu machen
wenn Sie die Lesbarkeit beeinträchtigen oder das Laufzeitverhalten.
Solange die Aufgabe nicht lautet, die Primzahlen von 1 bis 1E10 zu suchen,
sollte jede halbwegs "normal" programmierte Version in endlicher
Zeit zu den Ergebnissen kommen.
Verfasst: Freitag 12. Januar 2007, 09:24
von Luzandro
Ja schon klar - wollts nur erwähnen, v.a. auch weil ja die Idee vom Sieb darauf beruht, dass man die Vielfachen der Primzahlen entfernt (und ich es irgendwie witzig finde, dass du dann trotzdem xrange verwendest
)
Wirklich bitter
Verfasst: Freitag 12. Januar 2007, 09:49
von sunmountain
fand ich allerdings die Tatsache, das der junge Freund den Hinweisen
der anderen auf die wikipedia nicht gefolgt ist - dort gibt es ein gut leserliches Pascal Programm, das man leicht portieren kann.
(So hab' ich es dann gemacht
...)
Wer lesen kann ist klar im Vorteil ...
Verfasst: Freitag 12. Januar 2007, 12:09
von helmut
Hallo Cthulhu,
wenn Du unbedingt den/die Fehler in deinem Code haben möchtest, dann such sie. Als Hilfe für die Suche habe ich Deinen Code so mit einigen "print"-Anweisungen ergänzt, dass Du schrittweise den Ablauf des Programmes verfolgen kannst. Setze beim Suchen kreuz und quer gerade und ungerade Zahlen für die Grenzen ein. Sieh auch mal in die Dokumentation zur Funktion "range(von,
bis, schritt)".
Code: Alles auswählen
#Primzahlzwillinge
print "Um alle Priemzahlzwillinge in einem bestimmten Bereich rauszufinden, bitte den Bereich Definieren:"
a = input("Von: ")
b = input(" bis: ")
for p in range(a,b):
indent = ""
print "In for-Schleife. Laufvariable p = %i" % p
i=1
j=1
h=p/2
print "vor 1. while-Schleife"
print "Werte: i = %i, j = %i, h = %i" % (i,j,h)
while j!=0 and i<h:
indent = " "
print indent,"in 1. while-Schleife"
print indent,"Werte vor Anweisungsblock: i = %i, j = %i, h = %i" % (i, j, h)
i=i+1
j=p%i
print indent,"Werte nach Anweisungsblock: i = %i, j = %i, h = %i" % (i, j, h)
dummy = raw_input()
indent = ""
print "\nNach Durchlaufen/Ueberspringen der 1. while-Schleife"
print "Vor Pruefung 'if j!=0 and p!=a and p!=1:'"
print "Werte vor if: j = %i, p = %i, a = %i" % (j, p, a)
if j!=0 and p!=a and p!=1:
indent = " "
print indent,"Werte vor Anweisungsblock: p = %i" % (p)
z=p
print indent, "z = %i" % z
t=1
print "t = %i" % i
k=z+2
print "k = %i" % z
print indent, "Werte nach Anweisungsblock: z = %i, p = %i, t = %i, k = %i" % (z, p, t, k)
print indent, "vor 2. while-Schleife"
print indent, "Werte: t = %i, i = %i" % (t,i),"\n"
while t!=0 and i!=h+2:
indent = " "
print indent,"in 2. while-Schleife"
print indent,"Werte vor Anweisungsblock: i = %i, t = %i, k = %i" % (i, t, k)
i=i+1
t=k%i
print indent,"Werte nach Anweisungsblock: i = %i, t = %i," % (i, t)
dummy = raw_input
indent = " "
print indent, "\nNach Durchlaufen/Ueberspringen der 2. while-Schleife"
print indent,"Werte : t = %i, k = %i, a = %i, p = %i" % (t, k, a, p)
if t>0 and k!=a and p!=1:
print "\n****************"
print p
print k
print "****************\n"
dummy = raw_input
Viel Erfolg, Helmut
Verfasst: Freitag 12. Januar 2007, 20:20
von BlackJack
@sunmointain: Das `remove()` ist unschön. Damit müssen ziemlich viele Werte unnötig in der Gegend herum kopiert werden.
Luzandro hat geschrieben:Code: Alles auswählen
def primzahl_zwillinge(a,b):
bereich = range(2,b)
primzahlen = []
while len(bereich) > 0:
primzahl = bereich.pop(0)
primzahlen.append(primzahl)
bereich = [ zahl for zahl in bereich if zahl % primzahl != 0 ]
return [ (x,y) for x in primzahlen for y in primzahlen if x >= a and x+2 == y ]
Hier gibt's ein ähnliches "Problem".
Auch nicht optimal, weil die Listcomprehension am Ende im Gegensatz zu deinem Code unnötigerweise 2 Loops enthält, aber ich finde es lest sich so irgendwie schön funktional...
Das habe ich jetzt als Herausforderung aufgefasst. Schön funktional liest's sich so:
Code: Alles auswählen
from __future__ import division
import math
from itertools import ifilter, islice, izip, tee
def primes(start, end):
sieve = [True] * end
sieve[0] = sieve[1] = False
for number in xrange(2, int(math.sqrt(end))):
if sieve[number]:
for not_prime in xrange(2 * number, end, number):
sieve[not_prime] = False
for number in xrange(start, end):
if sieve[number]:
yield number
def tail(iterable):
return islice(iterable, 1, None)
def pairing(iterable):
iterator_a, iterator_b = tee(iterable)
return izip(iterator_a, tail(iterator_b))
def twin_primes(start, end):
def is_twin((a, b)):
return b - a == 2
return ifilter(is_twin, pairing(primes(start, end)))
def main():
print list(twin_primes(10, 100))
Verfasst: Samstag 13. Januar 2007, 08:35
von Luzandro
Ich sollte nicht so leichtfertig das Wort "funktional" verwenden
Also ich korrigiere: es lest sich so schön wie ausführbarer Pseudocode
Aber wirklich eine wunderschön funktionale und effiziente Methode hier - ich bin sicher Cthulhu wird auch sofort denken, wieso er nicht gleich darauf gekommen ist
(ich hab im ersten Moment genauso darauf gestarrt, wie auf Cthulhus Code...)
Was mir noch aufgefallen ist - das Stop beim xrange ist exclusiv, d.h. es würde noch überall ein +1 dazu gehören.
Noch die ausgebesserte "nicht funktionale Lösung" ohne dem Problem
Code: Alles auswählen
def primzahl_zwillinge(a, b):
stop = b+1
ist_prim = [True] * stop
ist_prim[0] = ist_prim[1] = False
for zahl in xrange(2, int(math.sqrt(b))+1):
if ist_prim[zahl]:
for vielfaches in xrange(zahl * 2, stop, zahl):
ist_prim[vielfaches] = False
return [ (x, x+2) for x in xrange(a, stop-2) if ist_prim[x] and ist_prim[x+2] ]
Verfasst: Samstag 13. Januar 2007, 12:43
von EyDu
So richtig funktional liest sich aber noch keine Lösung. Dazu müsste man wenigstens die Primzahlberechnung noch rekursiv machen. Was natürlich den Stack innerhalb von Bruchteilen einer Sekunde vernichtet
Daher möchte ich doch mal eine lesbare, und semi-funktionale Lösung vorschlagen. Wer will, darf natürlich die letzte Zeile noch in eine eigene Funktion schieben.
Code: Alles auswählen
def twins(a, b):
r = []
for i in range(2, b+1):
for p in r:
if i % p == 0:
break
else:
r.append(i)
return filter(lambda (x,y): x+2==y and x>=a, zip(r, r[1:]))
Verfasst: Samstag 20. Januar 2007, 13:49
von Cthulhu
also wie glaube ich schonmal gesagt kommt es nicht sehr gut biem lehrer an wenn ich hier irgendwelhe sachen benutze von denen er selber noch nie was gehört hat und er sicherlich merkt dass ich dass vorgesagt bekommen habe!
ich habe jetzt wie von blackjack empfohlen mal ein bischen verändert!
Aber ein problem bleibt immenroch! Es funktioniert ja alles einwandfrei bis auf die stelle an der das progr. über prüfen soll ob die primzahl +2 auch eine primzahl ist!
und an der stelle bräuchte ich dann hilfe bitte!
hier wieder der quellcode:
Code: Alles auswählen
#Primzahlzwillinge
print"Um alle Priemzahlzwillinge in einem bestimmten Bereich rauszufinden, bitte den Bereich Definieren:"
a=input("Von ")
b=input(" bis ")
for p in range(a,b):
i=1
j=1
h=p/2
while j!=0 and i<h:
i=i+1
j=p%i
if j!=0 and p!=a and p!=1:
t=1
k=p+2
while t!=0 and i!=h+2:
i=i+1
t=k%i
if t>0 :
print p
print k
print ""
Hoffe dass ich es dann endlich geschafft kriege!
danke
Verfasst: Samstag 20. Januar 2007, 15:34
von EyDu
Du gehst das Problem glaube ich von der falschen Seite an. An Stelle zu prüfen, ob die gefundene Primzahl+2 ebenfalls eine Primzahl ist, würde ich mir die letzte Primzahl merken, und dann prüfen, ob die aktuelle Primzahl-2 dieser entspricht.
Ich hoffe mal, dass das so als Denkanstoß genügt.
Verfasst: Samstag 20. Januar 2007, 17:44
von Gnushi
Hallo zusammen!
Hier ist meine Lösung. Die Untergrenze interessierte mich wenig und die Primzahlberechnung erfolgt auch nicht mit dem Sieb:
Code: Alles auswählen
def primzahlen(Grenze = 10):
"""
Erzeugt eine Liste von Primzahlen.
Grenze: Integerwert > 2
return: Liste
"""
Liste = [2]
def ist_prim(a):
"""
a ist eine Primzahl, wenn sie durch keine vorhergehende Primzahl
teilbar ist
return: True, wenn a prim, False sonst.
"""
tmp_ist_prim = True
for i in Liste:
if a % i == 0:
tmp_ist_prim = False
break
return tmp_ist_prim
# Grenze im richtigen Bereich?
assert(isinstance(Grenze, int))
assert(Grenze > 2)
# Alle Primzahlen auflisten
for Zahl in xrange(3, Grenze + 1):
if ist_prim(Zahl):
Liste.append(Zahl)
return Liste
Primzahlen = primzahlen(100)
Zwillinge = [(x,y) for x in Primzahlen for y in Primzahlen if x + 2 == y]
print Zwillinge
Liebe Grüße
Gnushi
Sieb++
Verfasst: Samstag 20. Januar 2007, 19:45
von sunmountain
Gnushi hat geschrieben:Hallo zusammen!
Hier ist meine Lösung. Die Untergrenze interessierte mich wenig und die Primzahlberechnung erfolgt auch nicht mit dem Sieb:
Code: Alles auswählen
def primzahlen(Grenze = 10):
"""
Erzeugt eine Liste von Primzahlen.
Grenze: Integerwert > 2
return: Liste
"""
Liste = [2]
def ist_prim(a):
"""
a ist eine Primzahl, wenn sie durch keine vorhergehende Primzahl
teilbar ist
return: True, wenn a prim, False sonst.
"""
tmp_ist_prim = True
for i in Liste:
if a % i == 0:
tmp_ist_prim = False
break
return tmp_ist_prim
# Grenze im richtigen Bereich?
assert(isinstance(Grenze, int))
assert(Grenze > 2)
# Alle Primzahlen auflisten
for Zahl in xrange(3, Grenze + 1):
if ist_prim(Zahl):
Liste.append(Zahl)
return Liste
Primzahlen = primzahlen(100)
Zwillinge = [(x,y) for x in Primzahlen for y in Primzahlen if x + 2 == y]
print Zwillinge
Liebe Grüße
Gnushi
Naja. die Ähnlichkeit zum Sieb ist aber schon da ...
Verfasst: Samstag 20. Januar 2007, 23:32
von EyDu
Den Block:
Code: Alles auswählen
tmp_ist_prim = True
for i in Liste:
if a % i == 0:
tmp_ist_prim = False
break
return tmp_ist_prim
kannst du übrigens noch vereinfachen. Ist halt so ein typsicher "Anfängerfehler".
Code: Alles auswählen
for i in Liste:
if a % i == 0:
return False
return True
Verfasst: Sonntag 21. Januar 2007, 01:40
von sape
Edit: War falsch was ich geschrieben hatte.
Verfasst: Sonntag 21. Januar 2007, 02:04
von sape
Hier mal mein eben schnell zusammenhackter versuch.
Der Primzahl test ist von Dookie aus diesem Thread:
http://www.python-forum.de/post-8344.html#8344
EDIT: Funktionierender Code -> siehe unten.
Verfasst: Sonntag 21. Januar 2007, 07:08
von sape
Ups, da habe ich einen Fehler gemacht.
Dieser code sollte nun funktionieren.
Code: Alles auswählen
#!/usr/bin/env python
# -*- coding: utf-8 -*-
import math
def isprim(n, plist=[2,3]):
#liste der Primzahlen bis sqrt(n) ergänzen falls nötig
for i in xrange(plist[-1], int(math.sqrt(n))+1,2):
for j in plist:
if i % j == 0:
break
else:
plist.append(i)
# n testen
if n in plist:
return True
else:
for i in plist:
if n % i == 0:
return False
else:
return True
def make_prims(from_, to):
return (x for x in xrange(from_, to) if isprim(x))
def get_prim_twins(list_):
for x in list_:
for y in list_:
if x + 2 == y:
yield x, y
elif y > x + 2:
break
if __name__ == "__main__":
print tuple(get_prim_twins(tuple(make_prims(0, 100))))
Code: Alles auswählen
((1, 3), (3, 5), (5, 7), (11, 13), (17, 19), (29, 31), (41, 43), (59, 61), (71, 73))
Verfasst: Sonntag 21. Januar 2007, 08:32
von BlackJack
`get_prim_twins()` ist aber nicht besonders effizient, da werden Unmengen an unnötigen vergleichen gemacht.
Verfasst: Sonntag 21. Januar 2007, 09:51
von sape
Jepp. Und zwar der bereich ``for y in list_:`` ist ineffizient. Es wird immer komplet die Liste durchlaufen, obwohl die "abgearbeiteten" längst irrelevant sind.
Verfasst: Sonntag 21. Januar 2007, 09:58
von sape
Code: Alles auswählen
def get_prim_twins(list_):
list2 = list_
for x in list_:
for idx, y in enumerate(list2):
if x + 2 == y:
yield x, y
list2 = list2[idx:]
break
Nun wird in der zweiten for schleife nicht die liste immer komplett von vorne durchgangen, sondern ab dem "letzten" Punkt.
EDIT:
Oder ist die Variante besser von der Lesbarkeit?
Code: Alles auswählen
def get_prim_twins(list_):
list2 = list_
for x in list_:
for idx in xrange(len(list2)):
if x + 2 == list2[idx]:
yield x, list2[idx]
list2 = list2[idx:]
break
EDIT:
Code: Alles auswählen
def get_prim_twins(list_):
last_idx = 0
for x in list_:
for idx in xrange(len(list_[last_idx:])):
if x + 2 == list_[last_idx+idx]:
last_idx += idx
yield x, list_[last_idx]
break