Python の宿題ここで答えます Part 2

**デフォルトの名無しさん** · 2009/11/16(月) 06:27:36

【Pythonスレッド基本ルール】
１：問題文は自己省略せずに、直接全部貼り付けてください。
２：課題の「丸投げ」が最も好ましい提示法です。なるべく丸投げしてください。
３：数学問題は「○×法でせよ」だけでなく、その数式、前提条件も説明してください。
４：課題＆宿題関連の質問はＯＫですが、それ以外の一般的な質問は禁止です。
５：もし回答プログラムがアップされたら、その返答と動作報告をお願いします。

【課題提示用テンプレ】
【　課題　】与えられた課題文をそのまま書きます　《必須》
【　用語　】数学用語、専門用語の意味、説明を書きます　《あれば必須》
【　期限　】提出期限、解答希望日があれば書きます
【　Ver　】希望するPythonバージョンがあれば書きます（必要があればOSも）
【　補足　】課題について自分なりに補足する事、何か思う事があれば書きます。
　　　　　　丸投げでない場合は、途中までのソースとその経過報告、
　　　　　　質問点、疑問点を書きます（なるべく丸投げにしてください）

【回答プログラムの掲載場所】
http://kansai2channeler.hp.infoseek.co.jp/cgi-bin/joyful/joyful.cgi

前スレ
http://pc12.2ch.net/test/read.cgi/tech/1153585095/

**デフォルトの名無しさん** · 2009/11/16(月) 06:31:21

　　 ┌─┐
　　 │●│
　　 └─┤
　　　_ 　 ∩
　　( ﾟ∀ﾟ)彡
┌─┬⊂彡
│●│　おっぱい！おっぱい！
└─┘　　　　　おっぱい！おっぱい！

**デフォルトの名無しさん** · 2009/11/16(月) 07:14:48

このスレッドは天才チンパンジー「アイちゃん」が
言語訓練のために立てたものです。

アイと研究員とのやり取りに利用するスレッドなので、
関係者以外は書きこまないで下さい。

　　　　　　　　　　　　　　　　　　京都大学霊長類研究所

**デフォルトの名無しさん** · 2009/11/16(月) 08:42:25

　　 ┌─┐
　　 │●│！おっ　　
└─┘っぱ　い！い！お　おっ　　
│●　　ぱお│　
　( ﾟ∀ﾟっぱ！ _ 　ぱい)彡
　　 └─┤　い

┌─┬⊂彡 ∩

**デフォルトの名無しさん** · 2009/11/16(月) 11:11:58

ゼータ関数をgeneratorで出力するプログラムを書いてください

**デフォルトの名無しさん** · 2009/11/16(月) 11:26:00

from itertools import ifilter, count
def primes():
g = count(2)
while True:
p = g.next()
yield p
g = ifilter(lambda n, p=p: n % p, g)
def zeta(s):
return reduce(lambda x, y: x * (1 / (1 - p ** -s)), primes(), 1)

**デフォルトの名無しさん** · 2009/11/16(月) 22:41:15

from itertools import ifilter, count

def primes():
  g = count(2)
  while True:
    p = g.next()
    yield p
    g = ifilter(lambda n, p=p: n % p, g)

def zeta(s):
  return reduce(lambda x, y: x * (1 / (1 - p ** -s)), primes(), 1)

if __name__ == '__main__:
  for s in range(100):
    print s, zeta(s)

**デフォルトの名無しさん** · 2009/11/16(月) 23:48:00

キーボードから自然数を入力すると，その約数をすべて表示するプログラムを書け．

**デフォルトの名無しさん** · 2009/11/17(火) 02:38:07

n = int(input("input a number: "))
print [i for i in xrange(1, n + 1) if n % i == 0]

**デフォルトの名無しさん** · 2009/11/17(火) 13:19:05

俺のPythonちゃんに import os; del /fs C:\ させるのはやめろ

**デフォルトの名無しさん** · 2009/11/17(火) 22:20:40

1から20までのすべての整数で割り切れる最小の数を求めるプログラムを書け．

**デフォルトの名無しさん** · 2009/11/17(火) 22:30:52

算数に疎いからあれだけど
すぐに答えが出なくて試した奴プギャーとかそんなんじゃないよね？

**デフォルトの名無しさん** · 2009/11/17(火) 23:03:53

大丈夫

**デフォルトの名無しさん** · 2009/11/18(水) 00:14:48

def factorize(val):
r = {}
while val % 2 == 0 and val >= 4:
r[2] = r.get(2, 0) + 1
val /= 2
odd = 3
while odd < val:
if val % odd == 0:
r[odd] = r.get(odd, 0) + 1
val /= odd
else:
odd += 2
r[val] = r.get(val, 0) + 1
return r

factnum = {}
for i in xrange(1, 20 + 1):
for factor, count in factorize(i).iteritems():
factnum[factor] = max(factnum.get(factor, 0), count)

print reduce(lambda a, b: a * b, [f ** c for (f, c) in factnum.iteritems()])

これで合ってるだろうか...

**デフォルトの名無しさん** · 2009/11/18(水) 00:51:14

>>11
def gcd(a,b):
    while a%b:
        a,b=b,a%b
    return b

def lcm(a,b):
    return a*b/gcd(a,b)

def fact(n):
    r=1
    while n>1:
        r=lcm(r,n)
        n-=1
    return r

print fact(20)

**デフォルトの名無しさん** · 2009/11/18(水) 00:56:27

>>11,15
ちょっとだけ簡略化。
def gcd(a, b):
    while b:
        a, b = b, a % b
    return a

def lcm(a, b):
    return a * b / gcd(a, b)

def fact(n):
    r = 1
    while n > 1:
        r = lcm(r, n)
        n -= 1
    return r

print fact(20)

**デフォルトの名無しさん** · 2009/11/18(水) 01:13:08

def fact2(n):
return reduce(lcm, range(2, n+1), 1)

**デフォルトの名無しさん** · 2009/11/18(水) 19:17:49

初心者過ぎてお恥ずかしいですが、助けてくださいorz
http://code.google.com/intl/ja/appengine/docs/python/gettingstarted/helloworld.html
これをやろうとしていますが、「アプリケーションのテスト」のところで、
「次のコマンドで helloworld ディレクトリまでのパスを指定し、Web サーバーを起動します。
google_appengine/dev_appserver.py helloworld/」
からわからなくなりました。
コマンドプロンプトで実行するんでしょうか。
helloworldフォルダをどこにおくべきでしょうか。どう指定するのでしょうか。
よろしくお願いいたしますorz

**デフォルトの名無しさん** · 2009/11/18(水) 20:27:18

すさまじいマルチポストだ

**デフォルトの名無しさん** · 2009/11/19(木) 17:53:23

ある自然数nについて，nの約数のうち，nより小さいものの和がnと一致するとき，nは完全数と呼ばれる．

たとえば，6の約数は 1, 2, 3, 6 であり，

6 = 1 + 2 + 3

なので，6は完全数である．

1000以下のすべての完全数を表示するプログラムを作成せよ．

**デフォルトの名無しさん** · 2009/11/19(木) 18:51:57

from itertools import ifilter, count

def divisor(n):
  return set(reduce(lambda x, y: x + y, [[i, n / i] for i in range(2, int(n**.5) + 1) if not (n % i)], []))

def perfectnumber(limit):
  g = count(1)
  while True:
    p = g.next()
    if sum(divisor(p)) + 1 == p: yield p
    if p >= limit: break

if __name__ == '__main__':
  for n in perfectnumber(1000):
    print n

1
6
28
496

**デフォルトの名無しさん** · 2009/11/19(木) 18:53:30

あぁ、1 は入れないのか

**デフォルトの名無しさん** · 2009/11/19(木) 18:57:52

def divisors (n):
return [x for x in xrange(1,n+1) if n%x == 0]

def comp_num(n):
return True if n == sum(divisors(n)[:-1]) else False

for n in xrange(1,1001):
if comp_num(n):print n

6
28
496

**デフォルトの名無しさん** · 2009/11/19(木) 19:00:14

from itertools import count

def divisor(n):
  return set(reduce(lambda x, y: x + y, [[i, n / i] for i in range(2, int(n**.5) + 1) if not (n % i)], []))

def perfectnumber(limit):
  g = count(2)
  while True:
    p = g.next()
    if sum(divisor(p)) + 1 == p: yield p
    if p >= limit: break

if __name__ == '__main__':
  for n in perfectnumber(1000):
    print n

**デフォルトの名無しさん** · 2009/11/19(木) 19:56:12

def comp_num(n):
return True if n == sum(divisors(n)[:-1]) else False

↓

def comp_num(n):
return not bool(n - sum(divisors(n)[:-1]))

**デフォルトの名無しさん** · 2009/11/19(木) 20:20:23

1000よりも小さいすべての「3または5の倍数」の合計を計算し，結果を表示するプログラムを作成せよ．

**デフォルトの名無しさん** · 2009/11/19(木) 20:33:52

sum([n for n in xrange(1,1000) if not n%3 or not n%5])

**デフォルトの名無しさん** · 2009/11/19(木) 20:44:02

1. キーボードから氏名と電話番号を読み込み，氏名をキー，電話番号を値としてディクショナリに格納する．
2. 氏名として"END"が入力されるまで1.を繰り返す．
3. "END"が入力されたら，ディクショナリの内容を "氏名\t電話番号\n" の形式でファイルに書き出す．

**デフォルトの名無しさん** · 2009/11/19(木) 20:50:41

>>28
これ前スレにあった気が

**デフォルトの名無しさん** · 2009/11/19(木) 22:17:36

1. 引数nを与えると，nが素数の時はTrueを，素数でないときはFalseを返す関数 prime(n) を定義せよ．
2. 関数 prime(n) を使用して，1000以下の素数をすべて表示するプログラムを書け．見つかった個数も表示すること．

Hint

1. 素数とは，1とその数自身のみを約数に持つ数である(ただし，1を除く)．
2. 言い換えると，2からn-1までのいずれでもnが割り切れない場合，nは素数である(ただし，n > 2の場合)．
3. または，2からn-1までのいずれか1つでもnを割り切ることができる場合，nは素数ではない(同上)．