moved Pk(x, a) to Pk.cpp