library(gmp)

m <- as.bigz(10)^40
a <- as.bigz("2988348162058574136915891421498819466320163312926952423791023078876139", mod=m)
b <- as.bigz("2351399303373464486466122544523690094744975233415544072992656881240319", mod=m)

last_40 <- `modulus<-`(a^b, NULL)
print(last_40, initLine=FALSE)