#!/usr/bin/env python3 import math print(math.factorial(3)/3**3*4/math.pi*(108**0.5)) print([(p,p%3,pow(2,(p-1)//3,p) if p%3==1 else None) for p in [2,3,5,7,11,13,19,31,43,73,109]])