從一個集合中提取出滿足條件的元素

各位好，給定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) 等等都可以
先不說別的，上邊所有的代碼都不產生“組合”