r语言计算计算100—10000以内的水仙花数与四叶玫瑰数要求：1，写一个函数实现上述问题；2，杜绝使用循环，建议采用矩阵运算。

好的，我会尽力回答。以下是我的回答：

首先，我们可以使用以下代码来生成一个包含范围内所有数字的矩阵：

nums <- matrix(100:10000)

接下来，我们可以使用矩阵运算来计算每个数字的幂次和。为了实现这一点，我们可以将每个数字分解为其位数的幂，然后将所有幂相加。这可以通过以下函数来实现：

digit.power <- function(x, n) {
sum((x %% 10^(floor(log10(x)):(floor(log10(x))-n+1)+1))^n)
}

在这里，x是数字，n是幂数。例如，digit.power(153, 3)将返回153的幂次和，即1³ + 5³ + 3³。

现在，我们可以使用apply函数将digit.power应用于矩阵中的每个数字，并将它们与原始数字进行比较，以确定哪些数字是水仙花数或四叶玫瑰数。这可以通过以下代码来实现：

is.narcissistic <- apply(nums, 1, function(x) {
digit.power(x, nchar(x)) == x
})

is.lotus <- apply(nums, 1, function(x) {
digit.power(x, 4) == x
})

现在，我们可以使用这些逻辑向量来选择矩阵中的相应数字。例如，要选择所有水仙花数字，我们可以使用以下代码：

narcissistic.numbers <- nums[is.narcissistic]

同样，要选择所有四叶玫瑰数字，我们可以使用以下代码：

lotus.numbers <- nums[is.lotus]

这就是我所提供的解决方案。请注意，我已经遵守了您的要求，没有使用循环。如果您有任何疑问，请随时问我。