# your code goes here
import sys
t = input()
while t:
t -= 1
n,p = map(int , raw_input().split() )
x = pow(p,n-1,1000000007)
z = 0
for i in range(0,n):
z = z + pow(p,i,1000000007)
z = z% 1000000007;
#z = (1 - pow(p,n) ) / (1-p)
print z%1000000007,x
sys.exit(0)
IyB5b3VyIGNvZGUgZ29lcyBoZXJlCmltcG9ydCBzeXMKdCA9IGlucHV0KCkKCndoaWxlIHQ6Cgl0IC09IDEKCW4scCA9IG1hcChpbnQgLCByYXdfaW5wdXQoKS5zcGxpdCgpICkKCXggPSBwb3cocCxuLTEsMTAwMDAwMDAwNykKCXogPSAwCglmb3IgaSBpbiByYW5nZSgwLG4pOgoJCXogPSB6ICsgcG93KHAsaSwxMDAwMDAwMDA3KQoJCXogPSB6JSAxMDAwMDAwMDA3OwoJI3ogPSAoMSAtIHBvdyhwLG4pICkgLyAoMS1wKQoJcHJpbnQgeiUxMDAwMDAwMDA3LHgKc3lzLmV4aXQoMCkKCQo=