Randomized incomplete $U$-statistics in high dimensions