Eneste egentlige trick her er at kende pakken gmp, og de ekstra datatyper til meget store tal den giver adgang til.

library(gmp)
t <- as.bigz(1:100)

answer <- sum(t)^2 - sum(t^2)