とりあえず(pr^qr)で割って、わかるのは
a'=(Π[k=1,r](1+pk+…+pk^qk))/(pr^qr)
ここまで。

a'=(Π[k=1,r-1](1+pk+…+pk^qk))
にどうやってたどり着くのかな。