#(1,2) (1,3,4) (5,6) (1,4,5)=> |1 1 0 0 0 0|
#                              |1 0 1 1 0 0|
#                              |0 0 0 0 1 1|
#                              |1 0 0 1 1 0|
x <- read.csv(text="id, productid
              1,1
              1,2
              2,1
              2,3
              2,4
              3,5
              3,6
              4,1
              4,4
              4,5")
a <- matrix(0, max(x$id), max(x$productid))
a[as.matrix(x[,c(1,2)])] <- 1