如何在两个向量的交叉组合上应用双变量函数以得到R中的矩阵？

Question

说我有两个向量：

a <- 1:4
b <- 1:2

和双变量函数：

f <- function(x,y) x**y

我想得到一个简单有效的方法（一个单行？）（对于这个具体的例子）：

[,1] [,2]
[1,]    1    1
[2,]    2    4
[3,]    3    9
[4,]    4   16

我能做到：

res <- matrix(nrow=length(a), ncol=length(b))
for (i in 1:length(b)){
  res[,i] <- mapply(f, a , b[i])
}

但我想避免循环。