Script Programming(스크립트 프로그래밍)

Lanczos 알고리즘은 Stirlng 공식에 의한 알고리즘 보다 정밀하며, 십진수로 유효숫자 약 15자리 까지는 정확하게 계산해 준다.  단지 exp 함수를 이용하는 부분에서는 exp 함수의 구현에 따라 오차가 더 있을 수 있다.

* 리눅스의 Python 2.6 으로 실행한 결과

Python 2.6.5 (r265:79063, Oct  1 2012, 22:04:36)
[GCC 4.4.3] on linux2
Type "help", "copyright", "credits" or "license" for more information.
>>> x = complex(1, 2./3)
>>> x
(1+0.66666666666666663j)
>>> "{0.real:.5}+{0.imag:.5}j".format(x)
'1.0+0.66667j'
>>> "{0.real:.5f}{0.imag:+.5f}j".format(x)
'1.00000+0.66667j'
>>> print "{0.real:.5f}{0.imag:+.5f}j".format(x)
1.00000+0.66667j
>>> "{0.real:.5f}{0.imag:+.5f}j".format(complex(1, -2./3))
'1.00000-0.66667j'
>>> print "{0.real:.5f}{0.imag:+.5f}j".format(complex(1, -2./3))
1.00000-0.66667j
>>> format(1+1j, '')
'(1+1j)'
>>> print format(1+1j, '')
(1+1j)

* 윈도우의 Python 2.7 로 실행한 결과   (Python 2.6의 경우와 동일함)

Python 2.7.3 (default, Apr 10 2012, 23:31:26) [MSC v.1500 32 bit (Intel)] on win32
Type "help", "copyright", "credits" or "license" for more information.
>>> x = complex(1, 2./3)
>>> x
(1+0.6666666666666666j)
>>> "{0.real:.5}+{0.imag:.5}j".format(x)
'1.0+0.66667j'
>>> "{0.real:.5f}{0.imag:+.5f}j".format(x)
'1.00000+0.66667j'
>>> print "{0.real:.5f}{0.imag:+.5f}j".format(x)
1.00000+0.66667j
>>> "{0.real:.5f}{0.imag:+.5f}j".format(complex(1, -2./3))
'1.00000-0.66667j'
>>> print "{0.real:.5f}{0.imag:+.5f}j".format(complex(1, -2./3))
1.00000-0.66667j
>>> format(1+1j, '')
'(1+1j)'
>>> print format(1+1j, '')
(1+1j)

* 윈도우의 Python 3.2 로 실행한 결과   (결과는 print() 함수의 사용법 차이 외에는 Python 2.7의 경우와 동일함)

Python 3.2 (r32:88445, Feb 20 2011, 21:29:02) [MSC v.1500 32 bit (Intel)] on win32
Type "help", "copyright", "credits" or "license" for more information.
>>> x = complex(1, 2./3)
>>> x
(1+0.6666666666666666j)
>>> "{0.real:.5}+{0.imag:.5}j".format(x)
'1.0+0.66667j'
>>> "{0.real:.5f}{0.imag:+.5f}j".format(x)
'1.00000+0.66667j'
>>> print( "{0.real:.5f}{0.imag:+.5f}j".format(x) )
1.00000+0.66667j
>>> print( "{0.real:.5f}{0.imag:+.5f}j".format(complex(1, -2./3)) )
1.00000-0.66667j
>>> format(1+1j, '')
'(1+1j)'
>>> print( format(1+1j, '') )
(1+1j)

* pyglet 을 이용하는 간단한 Python 소스 (소스파일명: testpyglet.py)

실행 결과:

* pyglet 을 이용하여 한글을 출력하는 Python 소스 (소스파일명: testpyglet2.py)
(소스파일을 저장시 반드시 UTF-8 인코딩으로 저장해야 한다)

1. # coding: MS949
3. #  Filename: testSyntheticDivision2.py
4. #
5. #  Purpose:  Find the quotient and remainder when some polynomial is
6. #            divided by a monic polynomial of the first degree.
7. #
8. #  Execute:  python testSyntheticDivision2.py 5 -4 7 8 6 8
9. #        Or  jython testSyntheticDivision2.py 5 -4 7 8 6 8
11. import sys
13. # 사용법 표시
14. def printUsage():
15.     print("사용법: python testSyntheticDivision2.py [제식의 계수] [피제식의 계수들]")
16.     print("조립제법(synthetic method)에 의한 다항식 나눗셈 결과를 보여준다.")
18. # 부동소수점수의 표현이 .0 으로 끝나는 경우 이를 잘라낸다.
19. # 전체 문자열 표시 너비는 매개변수 width 로 전달받아 처리한다.
20. def simplify(v, width=None):
21.     t = "%g" % v
22.     tlen = len(t)
23.     if t[tlen-2:tlen-1] == ".0":
24.         t = t[0, -2]
25.     if width != None:
26.         if tlen < width:
27.             t = "              "[0: width - tlen] + t
28.     return t
30. # "/수" 표시하기
31. # 부동소수점수의 표현이 .0 으로 끝나는 경우 이를 잘라낸다.
32. # 전체 문자열 표시 너비는 매개변수 width 로 전달받아 처리한다.
33. def simplifyReciprocal(v, width=None):
34.     t = "/%g" % v
35.     tlen = len(t)
36.     if t[tlen-2:tlen-1] == ".0":
37.         t = t[0, -2]
38.     if width != None:
39.         if tlen < width:
40.             t = "              "[0: width - tlen] + t
41.     return t
43. # 다항식을 내림차순의 스트링 표현으로 반환
44. def toPolyString(c):
45.     t = ""
46.     sc0 = simplify(c[0])
47.     if len(c) > 2:
48.         if sc0 == "1":
49.             t += "x^%d" % (len(c)-1)
50.         elif sc0 == "-1":
51.             t += "-x^%d" % (len(c)-1)
52.         else:
53.             t += sc0 + " x^%d" % (len(c)-1)
54.     elif len(c) == 2:
55.         if sc0 == "1":
56.             t += "x"
57.         elif sc0 == "-1":
58.             t += "-x"
59.         else:
60.             t += sc0 + " x"
61.     elif len(c) == 1:
62.         t += sc0
64.     for i in range(1, len(c)):
65.         k = len(c) - 1 - i
66.         sc = simplify(c[i])
67.         if k > 1:
68.             if c[i] > 0.0:
69.                 if sc == "1":
70.                     t += " + " + "x^%d" % k
71.                 else:
72.                     t += " + " + sc + " x^%d" % k
73.             elif c[i] < 0.0:
74.                 if sc == "-1":
75.                     t += " - " + "x^%d" % k
76.                 else:
77.                     t += " - " + simplify(abs(c[i])) + " x^%d" % k
78.         elif k == 1:
79.             if c[i] > 0.0:
80.                 if sc == "1":
81.                     t += " + " + "x"
82.                 else:
83.                     t += " + " + sc + " x"
84.             elif c[i] < 0.0:
85.                 if sc == "-1":
86.                     t += " - " + "x"
87.                 else:
88.                     t += " - " + simplify(abs(c[i])) + " x"
89.         elif k == 0:
90.             if c[i] > 0.0:
91.                 t += " + " + sc
92.             elif c[i] < 0.0:
93.                 t += " - " + simplify(abs(c[i]))
94.     return t
96. # 다항식 나눗셈 결과를
97. #     (피제식) = (제식)(몫) + (나머지)
98. # 형태로 출력
99. def printDivisionResult(a, c, d, b):
100.     strLine = "  " + toPolyString(c)
101.     print( strLine )
103.     strLine = "    = ( " + toPolyString( [a[0], -a[1]] ) + " )"
104.     tmp = [0.0] * (len(b) - 1)
105.     for i in range(0, len(tmp)):
106.         tmp[i] = b[i]
107.     strLine += "( " + toPolyString(tmp) + " )"
108.     r = d[len(d) - 1]
109.     if r > 0.0:
110.         strLine += " + " + simplify(r)
111.     elif r < 0.0:
112.         strLine += " - " + simplify(abs(r))
114.     print( strLine )
116. # 조립제법 계산표 출력 함수
117. def printSyntheticTable(a, c, s, d, q):
118.     strLine = "         | "
119.     strLine += simplify(c[0], 8)
120.     for i in range(1, len(c)):
121.         strLine += "  " + simplify(c[i], 8)
122.     print( strLine )
124.     strLine = simplify(a[1], 8) + " |"
125.     strLine += "           "
126.     strLine += simplify(s[1], 8)
127.     for i in range(2, len(s)):
128.         strLine += "  " + simplify(s[i], 8)
129.     print( strLine )
131.     strLine = "         |"
132.     for i in range(0, len(q)):
133.         strLine += "----------"
134.     print( strLine )
136.     strLine = "         | "
137.     strLine += simplify(d[0], 8)
138.     for i in range(1, len(d)):
139.         strLine += "  " + simplify(d[i], 8)
140.     print( strLine )
142.     strLine = simplifyReciprocal(a[0], 8) + " | "
143.     strLine += simplify(q[0], 8)
144.     for i in range(1, len(q) - 1):
145.         strLine += "  " + simplify(q[i], 8)
146.     print( strLine )
148. # 실행 시작 지점
149. if len(sys.argv) < 3:
150.     printUsage()
151.     sys.exit(1)
153. ######################################################
154. # 피제식은 c_0 x^n +  c_1 x^(n -1) + ... + c_n
155. # 제식은 a x -  b
156. a = [float(sys.argv[1]), - float(sys.argv[2])]
157. c = [0.0]*(len(sys.argv) - 3)
158. s = [0.0]*(len(sys.argv) - 3)
159. d = [0.0]*(len(sys.argv) - 3)
160. b = [0.0]*(len(sys.argv) - 3)
162. for i in range(0, len(c)):
163.     c[i] = float(sys.argv[i + 3])
165. ######################################################
166. # 조립제법의 주요 부분
167. s[0] = 0.0
168. d[0] = c[0]
169. b[0] = c[0] / a[0]
170. for i in range(1, len(c)):
171.     s[i] = b[i-1]*a[1]
172.     d[i] = c[i] + s[i]
173.     b[i] = d[i] / a[0]
175. ######################################################
176. # 몫의 계수와 나머지를 출력한다.
177. print("몫의 계수는"),
178. for i in range(0, len(b) - 2):
179.     print(simplify(b[i], None) + ", " ),
180. print(simplify(b[len(b) - 2], None)),
181. print("이고, 나머지는 " + simplify(d[len(d) - 1], None) + " 이다.")
182. print
184. ######################################################
185. # 조립제법 표를 출력한다.
186. printSyntheticTable(a, c, s, d, b)
187. print
189. ######################################################
190. # (피제식) = (제식) x (몫) + (나머지)
191. printDivisionResult(a, c, d, b)

실행> python testSyntheticDivision2.py 5 -4 7 8 6 8
몫의 계수는 1.4,  2.72,  3.376 이고, 나머지는 21.504 이다.

         |        7         8         6         8
       4 |                5.6     10.88    13.504
         |----------------------------------------
         |        7      13.6     16.88    21.504
      /5 |      1.4      2.72     3.376

  7 x^3 + 8 x^2 + 6 x + 8
    = ( 5 x - 4 )( 1.4 x^2 + 2.72 x + 3.376 ) + 21.504




이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.

Python 언어로 숫자 맞추기 게임을 작성해 보았다.
input()  함수가 내부적으로 Python 2.x와 Python3.x에서 다르게 동작한다.

    In Python2.*, input([prompt]) was equivalent to eval(raw_input([prompt])).
    In Python3.*, input([prompt]) of Python2.* was removed and raw_input([prompt])
                        was renamed as input([prompt]).

아래 소스의 16째 줄 guess = int(sbuf) 는 스트링을 정수로 타입변환하는 과정이다.



소스 파일명: guessNumber01.py

1. #!/usr/bin/env python
2. #
3. #   Filename: guessNumber01.py
4. #   Purpose:  Interatice game guessing a given number.
5. #                 if CONDITION:
6. #                     ......
7. #                 else:
8. #                     ......
9. #   Execute: python guessNumber01.py
11. def doGuessing(num):
12.     print("Enter your guess:")
13.     #sbuf = raw_input()     # for Python 2.x
14.     sbuf = input()          # for Python 3.x
15.     guess = int(sbuf)
16.     if guess == num:
17.         print("You win!")
18.         return
20.     # we won't get here if guess == num
21.     if guess < num:
22.         print("Too low!")
23.         doGuessing(num)
24.     else:
25.         print("Too high!")
26.         doGuessing(num)
28. doGuessing(123)

실행> python guessNumber01.py
Enter your guess:
111
Too low!
Enter your guess:
222
Too high!
Enter your guess:
123
You win!



위의 소스는 Jython 2.2.1이나 Jython 2.5.1로 실행해도 똑같이 실행된다.

실행> jython guessNumber01.py
Enter your guess:
111
Too low!
Enter your guess:
222
Too high!
Enter your guess:
123
You win!





이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.

[파일명:  testStringFindInList.py]------------------------------------------------
# coding=ms949

import sys

def find(arr, s):
    for i in range(0, len(arr)):
     if arr[i].find(s) >= 0:
      return i
    return -1;

def printArray(arr):
    sys.stdout.write("[")
    for i in range(0, len(arr) - 1):
         sys.stdout.write(arr[i] + ", ")
    if len(arr) > 0:
        sys.stdout.write(arr[len(arr) - 1])
    print("]")

words = ["하나", "둘", "셋", "넷", "다섯", "여섯"]

sys.stdout.write("list: ")
printArray(words)
where = find(words, "셋")
if where > 0:
    sys.stdout.write("발견!  ")
    print("Next word of 셋 in list: " + words[where+1])

print("Sorting...")
words.sort()

sys.stdout.write("list: ")
printArray(words)
where = find(words, "셋")
if where > 0:
    sys.stdout.write("발견!  ")
    print("Next word of 셋 in list: " + words[where+1])
------------------------------------------------


실행> python testStringFindInList.py
list: [하나, 둘, 셋, 넷, 다섯, 여섯]
발견!  Next word of 셋 in list: 넷
Sorting...
list: [넷, 다섯, 둘, 셋, 여섯, 하나]
발견!  Next word of 셋 in list: 여섯

크리에이티브 커먼즈 라이선스

이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-동일조건변경허락 2.0 대한민국 라이선스에 따라 이용하실 수 있습니다.

[파일명:  testSort.py]------------------------------------------------
import sys

def printArray(a):
    sys.stdout.write("[")
    for i in range(0, len(a) - 1):
        sys.stdout.write("%s, " % a[i])
    if len(a) > 0:
        sys.stdout.write("%s" % a[-1])
    sys.stdout.write("]\n")

list = sys.argv[1::]
list.sort()
printArray(list)
------------------------------------------------


실행> python testSort.py one two three four five
[five, four, one, three, two]

실행> python testSort.py 하나 둘 셋 넷 다섯
[넷, 다섯, 둘, 셋, 하나]


IronPython으로 실행하기:
실행> ipy testSort.py one two three four five
[five, four, one, three, two]

실행> ipy testSort.py 하나 둘 셋 넷 다섯
[넷, 다섯, 둘, 셋, 하나]


Jython으로 실행하기:
실행> jython testSort.py one two three four five
[five, four, one, three, two]

( Jython으로 실행하는 경우, )
( 버전 2.2.1이든 2.5b3이든 명령행 옵션으로 넣어준 한글을 잘 처리하지 못한다. )
실행> jython testSort.py 하나 둘 셋 넷 다섯
[?, ??, ?, ?, ??]

크리에이티브 커먼즈 라이선스

이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-동일조건변경허락 2.0 대한민국 라이선스에 따라 이용하실 수 있습니다.