#!/usr/bin/env python
# coding: utf-8

# [Oregon Curriculum Network](http://www.4dsolutions.net/ocn) <br />
# [Discovering Math with Python](Introduction.ipynb)
# 
# 
# Crystal Ball Sequence
# ============
# The face-centered cubic (FCC) lattice is not always presented in this simplest form, ditto the cubic close packing (CCP), which amounts to the same thing.  A nuclear ball is surrounded by a layer of twelve, all touching it, and adjacent neighbors.  The shape so formed is not a cube, but a cuboctahedron, with eight triangular faces and six square.  This is where I can type stuff.
# 
# As the cuboctahedral packing continues to expand outward, layer by layer, the cumulative number of balls or points forms the Crystal Ball Sequence.
# 
# cubocta(), a generator, yields the number of balls in each successive layer of the cuboctahedron, according to a simple formula derived by R. Buckminster Fuller, a prolific inventor and philosopher [1].  cummulative( ) delegates to cubocta( ) while accumulating the number in each layer to provide a running total.
# 
# ![Image of Cubocta](http://www.4dsolutions.net/ocn/graphics/cubanim.gif)

# In[1]:


from itertools import accumulate, islice

def cubocta():
    """
    Classic Generator:  Cuboctahedral / Icosahedral #s
    https://oeis.org/A005901
    """
    yield 1  # nuclear ball
    f = 1
    while True:
        elem = 10 * f * f + 2  # f for frequency
        yield elem    # <--- pause / resume here
        f += 1
        
def cummulative(n):
    """
    https://oeis.org/A005902 (crystal ball sequence)
    """
    yield from islice(accumulate(cubocta()),0,n)

print("{:=^30}".format(" Crystal Ball Sequence "))
print("{:^10} {:^10}".format("Layers", "Points"))
for f, out in enumerate(cummulative(30),start=1):
    print("{:>10} {:>10}".format(f, out))
   


# Octet Truss
# ======
# 
# When adjacent CCP ball centers interconnect, what do you get?  Why the octet truss of course, a well known space frame, used a lot in architecture.  Alexander Graham Bell was fascinated by this construction.[2]
# 
# <a data-flickr-embed="true"  href="https://www.flickr.com/photos/kirbyurner/23636692173/in/album-72157664250599655/" title="Business Accelerator Building"><img src="https://farm2.staticflickr.com/1584/23636692173_103b411737.jpg" width="500" height="375" alt="Business Accelerator Building"></a><script async src="//embedr.flickr.com/assets/client-code.js" charset="utf-8"></script>

# [1] Siobahn Roberts. *King of Infinite Space*. New York: Walker & Company (2006). pp 179-180.
# 
# "Coxeter sent back a letter saying that one equation would be 'a remarkable discovery, justifying Bucky's evident pride,' if only it weren't too good to be true.  The next day, Coxeter called: 'On further reflection, I see that it *is* true'. Coxeter told Fuller how impressed he was with his formula -- on the cubic close-packing of balls."
# 
# [2]  http://worldgame.blogspot.com/2006/02/octet-truss.html (additional info on the octet truss)

# Pascal's Triangle
# =========
# 
# Pascal's Triangle connects to the Binomial Theorem (originally proved by Sir Isaac Newton) and to numerous topics in probability theory.  The triangular and tetrahedral number sequences may be discovered lurking in its columns.
# 
# pascal(), a generator, yields successive rows of Pascal's Triangle.  By prepending and appending a zero element and adding vertically, a next row is obtained.  For example, from [1] we get [0, 1] + [1, 0] = [1, 1].  From [1, 1] we get [0, 1, 1] + [1, 1, 0] = [1, 2, 1] and so on.
# 
# Notice the triangular numbers (1, 3, 6, 10...) and tetrahedral number sequences (1, 4, 10, 20...) appear in the slanted columns. [3]
# 
# ![Image of PT](http://www.4dsolutions.net/ocn/graphics/randtrianim.gif)

# In[2]:


from itertools import islice

def pascal():
    row = [1]
    while True:
        yield row
        row = [i+j for i,j in zip([0]+row, row+[0])]

print("{0:=^60}".format(" Pascal's Triangle "))
print()
for r in islice(pascal(),0,11):
    print("{:^60}".format("".join(map(lambda n: "{:>5}".format(n), r))))


# Each number in Pascal's Triangle may be understood as the number of unique pathways to that position, were falling balls introduced through the top and allowed to fall left or right to the next row down.  This apparatus is sometimes called a Galton Board.
# 
# For example, a ball could reach the 6 in the middle of the 5th row going 1,1,2,3,6 in four ways (counting left and right mirrors), or 1,1,1,3,6 in two ways.  The likely pattern when many balls fall through this maze will be a bell curve, as shown in the simulation below.

# In[1]:


from IPython.display import YouTubeVideo
YouTubeVideo("9xUBhhM4vbM")


# [3] for more on both Pascal's Triangle and Cuboctahedral Numbers see:
# http://www.4dsolutions.net/ocn/numeracy0.html
#