没问题！在R里完全可以不用for循环实现两个列表所有元素对的穷尽交集，我给你分享两种简洁高效的方案，兼容base R和tidyverse生态：

首先先还原你的示例数据（保证结果可复现）：

set.seed(7)
x <- 1:25
x1 <- sample(x, 11)
x2 <- sample(x, 7)
x3 <- sample(x, 15)
list1 <- list(x1 = x1, x2 = x2, x3 = x3)
y <- 1:20
y1 <- sample(y, 13)
y2 <- sample(y, 10)
list2 <- list(y1 = y1, y2 = y2)

方法1：使用tidyverse的purrr包（更直观易读） #

利用expand.grid生成所有元素对的组合，再用map2批量计算交集：

library(purrr)

# 生成list1和list2所有元素的笛卡尔积组合
element_pairs <- expand.grid(
  list1_item = names(list1),
  list2_item = names(list2),
  stringsAsFactors = FALSE
)

# 对每一对元素计算交集
intersection_results <- map2(
  list1[element_pairs$list1_item],
  list2[element_pairs$list2_item],
  intersect
)

# 给结果命名，方便识别对应关系
names(intersection_results) <- paste(
  element_pairs$list1_item,
  element_pairs$list2_item,
  sep = "_vs_"
)

# 查看最终结果
intersection_results

这个方法的优势是逻辑清晰：先明确所有需要计算的元素对，再批量执行交集操作，结果的命名也能让你一眼看出是哪两个向量的交集。

方法2：Base R原生实现（无需额外安装包） #

用outer函数直接对两个列表的所有元素对应用交集计算：

# 自定义一个返回列表的交集函数，适配outer的矩阵输出
calc_intersect <- function(vec_a, vec_b) {
  list(intersect(vec_a, vec_b))
}

# 生成所有元素对的交集结果（矩阵格式的列表）
intersection_matrix <- outer(list1, list2, FUN = calc_intersect)

# 转成扁平化的命名列表（可选，方便查看）
flattened_results <- unlist(intersection_matrix, recursive = FALSE)
names(flattened_results) <- paste(
  rep(names(list1), length(list2)),
  rep(names(list2), each = length(list1)),
  sep = "_vs_"
)

# 查看结果
flattened_results

这个方法完全依赖base R，不需要安装任何扩展包，适合追求轻量解决方案的场景。

验证结果 #

你可以手动验证其中一组的结果，比如计算x1和y1的交集：

intersect(list1$x1, list2$y1)

会发现和两种方法输出的x1_vs_y1结果完全一致，说明方案是正确的。

内容的提问来源于stack exchange，提问作者SteveM