从一个集合中提取出满足条件的元素

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

各位好，给定x_1, x_2, x_3, 集合A={(y_1, y_2, y_3): y_1+y_2<=x_1, y_1+y_3<=x_2, y_2+y_3<=x_3, y_i均为自然数}。

请问R 里面有没有一个函数或者命令直接把A中的元素拿出来？

多谢啦。 请不要用for 循环。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：各位好

for 循环 我是知道的

1. library(matlib)
   
2. A <- matrix(c(1, 1, 0,1, 0, 1,0, 1, 1), 3, 3, byrow=TRUE)  #约束条件其实是三元方程组 系数转化为矩阵
   
3. b <- c(7, 10, 15)  #假设给定的三个x是7,10,15
   
4. 
   
5. E <- echelon(A,b) #求阶梯矩阵，结果的第四列是y的最大解
   
6. 
   
7. #列出各y小于最大解的自然数
   
8. y1<-c(1:E[1,4])
   
9. y2<-c(1:E[2,4])
   
10. y3<-c(1:E[3,4])
    
11. 
    
12. #这段代码错误 不能求出全部解 且没有实现列示集合所有y组合元素 歉

复制代码

就是3元一次方程嘛。

1. x <- c(7, 8, 9)
   
2. y1 <- c(1:floor((x[1]+x[2]-x[3])/2))
   
3. y2 <- c(1:floor((x[1]+x[3]-x[2])/2))
   
4. y3 <- c(1:floor((x[2]+x[3]-x[1])/2))
   
5. A <- c(y1, y2, y3) %>%
   
6.   unique()

复制代码

drunkfish69 发表于 2020-8-27 14:55
就是3元一次方程嘛。

确实是三元方程组

drunkfish69 发表于 2020-8-27 14:55
就是3元一次方程嘛。

运行之后有错误啊，能再修改一下吗

1. library(tidyverse)
   
2. x <- c(7, 8, 9)
   
3. y1 <- c(1:floor((x[1]+x[2]-x[3])/2))
   
4. y2 <- c(1:floor((x[1]+x[3]-x[2])/2))
   
5. y3 <- c(1:floor((x[2]+x[3]-x[1])/2))
   
6. A <- c(y1, y2, y3) %>%
   
7.   unique()

复制代码

运行之后有错误啊，能再修改一下吗

1. x <- c(10, 8, 9)
   
2. x <- sort(x)
   
3. A <- c(1:floor((x[3]+x[2]-x[1])/2))

复制代码

这样就行，也不用麻烦y_1～3了

llb_321 发表于 2020-8-27 14:38

这算法肯定是有问题的

比如你算出的y1=1
其实y1=2也可以

上边的所有代码都不行
楼主要求的是所有组合
比如b <- c(7, 10, 15)  #假设给定的三个x是7,10,15
那么y = (2,1,1) （6,1,1) 等等都可以
先不说别的，上边所有的代码都不产生“组合”