The Big Bang

import os
from random import randrange

universe_temp = 1000000000
matter = True
matter_count = 0
anti_matter = False
anti_matter_count = 0
photons = 0
protons = 0
neutrons = 0
electrons = 0
hydrogen_atom = 1
helium_atom = 1
lithium_atom = 1
at_nu_complete = {"Hydrogen": hydrogen_atom, "Helium": helium_atom, "Lithium": lithium_atom}


class Atomic_nuclei():
    def __init__(self, protons, neutrons, electrons):
        self.protons = protons
        self.neutrons = neutrons
        self.electrons = electrons


        self.electrons += 1
        if protons == 1 and neutrons == 0:
            at_nu_complete["Hydrogen"] += 1
        elif protons == 2 and neutrons == 2:
            at_nu_complete["Helium"] += 1
        elif protons == 3 and neutrons == 3:
            at_nu_complete["Lithium"] += 1

def display():
    print"Universe Temp =  %d" % universe_temp
    print"Matter        =  %d" % matter_count
    print"Anti-Matter   =  %d" % anti_matter_count
    print"Photons       =  %s" % photons
    print(at_nu_complete)
    os.system("cls")

def at_nu_gen():
    protons = randrange(0, 4)
    neutrons = randrange(0, 4)
    at_nu = Atomic_nuclei(protons, neutrons, electrons)
    return at_nu

while True:
    if universe_temp >= 1000000000:
        matter_count += 1
        anti_matter_count += 1
        photons += 1000000000
        if anti_matter_count > 1000000000:
             matter_count += 1

        at_nu = at_nu_gen()
        universe_temp += -1
        display()

    if universe_temp <= 999999999:
        electo_magentic_force = True
        weak_force = True
        if matter > 1:
            matter_count += -1000000
            anti_matter_count += -1000000

        if photons > 1000000000:
            matter_count = 0
            a = 1000000000
            for a in photons:
                matter_count += 1

        at_nu = at_nu_gen()
        universe_temp += -1
        display()