data/math/sierpinski


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26

# Sierpinski triangles
import sys
size = 3
modulus = 2

lines = modulus**size

vector = [1]
for i in range(1,lines+1):
  vector.insert(0,0)
  vector.append(0)

for i in range(0,lines):
  newvector = vector[:]
  for j in range(0,len(vector)-1):
    if (newvector[j] == 0):
      sys.stdout.write(" ")
    else:
      remainder = newvector[j] % modulus
      if (remainder == 0):
        sys.stdout.write("O")
      else:
        sys.stdout.write(".")
    newvector[j] = vector[j-1] + vector[j+1]
  print
  vector = newvector[:]