import time
import math

def P(number):
    return (number*(3*number - 1))/2

def isP(number):
    a = 3
    b = -1
    c = -2*number
    z = math.sqrt(b**2 - 4*a*c)
    if z>b:
        x = (-b+z)/(2*a)
        return x==int(x)
    else:
        return false

def solve():    
    pentagonals = []
    n = 1
    while True:
        newP = P(n)
        for oldP in pentagonals:
            if isP(newP+oldP) and (newP-oldP) in pentagonals: 
                print(newP,oldP)
                print(newP-oldP)
                return
            
        pentagonals.append(newP)
        n+=1

t = time.clock()    
solve()
print(time.clock()-t)
