我被卡住了。我有一个优化问题,这个问题是在Lingo上解决的,但是它变得太大了。我试着用其他语言做这个程序,我真的被卡住了。你知道我怎么用R来构建这个程序吗?
"Objective Function;
Min = x111 + x121 + x131 + x141 + x151 + x211 + x221 + x231 + x241 + x251 + x311 + x321 + x331 + x341 + x351 + x411 + x421 + x431 + x441 + x451 + x512 + x522 + x532 + x542 + x552 + x612 + x622 + x632 + x642 + x652 + x113 + x123 + x133 + x143 + x153 + x213 + x223 + x233 + x243 + x253 + x313 + x323 + x333 + x343 + x353 + x413 + x423 + x433 + x443 + x453 + x314 + x324 + x334 + x344 + x354 + x414 + x424 + x434 + x444 + x454 + x514 + x524 + x534 + x544 + x554 + x614 + x624 + x634 + x644 + x654 + x115 + x125 + x135 + x145 + x155 + x215 + x225 + x235 + x245 + x255 + x315 + x325 + x335 + x345 + x355 + x415 + x425 + x435 + x445 + x455 + x515 + x525 + x535 + x545 + x555 + x615 + x625 + x635 + x645 + x655 + x116 + x126 + x136 + x146 + x156 + x216 + x226 + x236 + x246 + x256 + x316 + x326 + x336 + x346 + x356 + x416 + x426 + x436 + x446 + x456 + x117 + x127 + x137 + x147 + x157 + x217 + x227 + x237 + x247 + x257 + x317 + x327 + x337 + x347 + x357 + x417 + x427 + x437 + x447 + x457 + x118 + x128 + x138 + x148 + x158 + x218 + x228 + x238 + x248 + x258 + x318 + x328 + x338 + x348 + x358 + x418 + x428 + x438 + x448 + x458;
Subject To
x11 + x1 + x31 + x71 + x101 + x121 + x141 <= 7 ;
x12 + x2 + x32 + x72 + x102 + x122 + x142 <= 7 ;
x13 + x3 + x33 + x73 + x103 + x123 + x143 <= 7 ;
x14 + x4 + x34 + x74 + x104 + x124 + x144 <= 7 ;
x15 + x5 + x35 + x75 + x105 + x125 + x145 <= 7 ;
x21 + x6 + x36 + x76 + x106 + x126 + x146 <= 7 ;
x22 + x7 + x37 + x77 + x107 + x127 + x147 <= 7 ;
x23 + x8 + x38 + x78 + x108 + x128 + x148 <= 7 ;
x24 + x9 + x39 + x79 + x109 + x129 + x149 <= 7 ;
x25 + x10 + x40 + x80 + x110 + x130 + x150 <= 7 ;
x31 + x11 + x41 + x51 + x81 + x111 + x131 + x151 <= 7 ;
x32 + x12 + x42 + x52 + x82 + x112 + x132 + x152 <= 7 ;
x33 + x13 + x43 + x53 + x83 + x113 + x133 + x153 <= 7 ;
x34 + x14 + x44 + x54 + x84 + x114 + x134 + x154 <= 7 ;
x35 + x15 + x45 + x55 + x85 + x115 + x135 + x155 <= 7 ;
x41 + x16 + x46 + x56 + x86 + x116 + x136 + x156 <= 7 ;
x42 + x17 + x47 + x57 + x87 + x117 + x137 + x157 <= 7 ;
x43 + x18 + x48 + x58 + x88 + x118 + x138 + x158 <= 7 ;
x44 + x19 + x49 + x59 + x89 + x119 + x139 + x159 <= 7 ;
x45 + x20 + x50 + x60 + x90 + x120 + x140 + x160 <= 7 ;
x51 + x21 + x61 + x91 <= 7 ;
x52 + x22 + x62 + x92 <= 7 ;
x53 + x23 + x63 + x93 <= 7 ;
x54 + x24 + x64 + x94 <= 7 ;
x55 + x25 + x65 + x95 <= 7 ;
x61 + x26 + x66 + x96 <= 7 ;
x62 + x27 + x67 + x97 <= 7 ;
x63 + x28 + x68 + x98 <= 7 ;
x64 + x29 + x69 + x99 <= 7 ;
x65 + x30 + x70 + x100 <= 7 ;
x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8 + x9 + x10 + x11 + x12 + x13 + x14 + x15 + x16 + x17 + x18 + x19 + x20 >= 11.8134 ;
x21 + x22 + x23 + x24 + x25 + x26 + x27 + x28 + x29 + x30 >= 15.5714 ;
x31 + x32 + x33 + x34 + x35 + x36 + x37 + x38 + x39 + x40 + x41 + x42 + x43 + x44 + x45 + x46 + x47 + x48 + x49 + x50 >= 34.2000 ;
x51 + x52 + x53 + x54 + x55 + x56 + x57 + x58 + x59 + x60 + x61 + x62 + x63 + x64 + x65 + x66 + x67 + x68 + x69 + x70 >= 0.0000 ;
x71 + x72 + x73 + x74 + x75 + x76 + x77 + x78 + x79 + x80 + x81 + x82 + x83 + x84 + x85 + x86 + x87 + x88 + x89 + x90 + x91 + x92 + x93 + x94 + x95 + x96 + x97 + x98 + x99 + x100 >= 0.0000 ;
x101 + x102 + x103 + x104 + x105 + x106 + x107 + x108 + x109 + x110 + x111 + x112 + x113 + x114 + x115 + x116 + x117 + x118 + x119 + x120 >= 0.0000 ;
x121 + x122 + x123 + x124 + x125 + x126 + x127 + x128 + x129 + x130 + x131 + x132 + x133 + x134 + x135 + x136 + x137 + x138 + x139 + x140 >= 0.0000 ;
x141 + x142 + x143 + x144 + x145 + x146 + x147 + x148 + x149 + x150 + x151 + x152 + x153 + x154 + x155 + x156 + x157 + x158 + x159 + x160 >= 4.1652 ;
x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8 + x9 + x10 + x11 + x12 + x13 + x14 + x15 + x16 + x17 + x18 + x19 + x20 + x21 + x22 + x23 + x24 + x25 + x26 + x27 + x28 + x29 + x30 + x31 + x32 + x33 + x34 + x35 + x36 + x37 + x38 + x39 + x40 + x41 + x42 + x43 + x44 + x45 + x46 + x47 + x48 + x49 + x50 + x51 + x52 + x53 + x54 + x55 + x56 + x57 + x58 + x59 + x60 + x61 + x62 + x63 + x64 + x65 + x66 + x67 + x68 + x69 + x70 + x71 + x72 + x73 + x74 + x75 + x76 + x77 + x78 + x79 + x80 + x81 + x82 + x83 + x84 + x85 + x86 + x87 + x88 + x89 + x90 + x91 + x92 + x93 + x94 + x95 + x96 + x97 + x98 + x99 + x100 + x101 + x102 + x103 + x104 + x105 + x106 + x107 + x108 + x109 + x110 + x111 + x112 + x113 + x114 + x115 + x116 + x117 + x118 + x119 + x120 + x121 + x122 + x123 + x124 + x125 + x126 + x127 + x128 + x129 + x130 + x131 + x132 + x133 + x134 + x135 + x136 + x137 + x138 + x139 + x140 + x141 + x142 + x143 + x144 + x145 + x146 + x147 + x148 + x149 + x150 + x151 + x152 + x153 + x154 + x155 + x156 + x157 + x158 + x159 + x160 >= 65.7501 ;
!x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12, x13, x14, x15, x16, x17, x18, x19, x20, x21, x22, x23, x24, x25, x26, x27, x28, x29, x30, x31, x32, x33, x34, x35, x36, x37, x38, x39, x40, x41, x42, x43, x44, x45, x46, x47, x48, x49, x50, x51, x52, x53, x54, x55, x56, x57, x58, x59, x60, x61, x62, x63, x64, x65, x66, x67, x68, x69, x70, x71, x72, x73, x74, x75, x76, x77, x78, x79, x80, x81, x82, x83, x84, x85, x86, x87, x88, x89, x90, x91, x92, x93, x94, x95, x96, x97, x98, x99, x100, x101, x102, x103, x104, x105, x106, x107, x108, x109, x110, x111, x112, x113, x114, x115, x116, x117, x118, x119, x120, x121, x122, x123, x124, x125, x126, x127, x128, x129, x130, x131, x132, x133, x134, x135, x136, x137, x138, x139, x140, x141, x142, x143, x144, x145, x146, x147, x148, x149, x150, x151, x152, x153, x154, x155, x156, x157, x158, x159, x160>=0, interger (interger);
字符串
我没有写代码的经验。有人能帮我建一个或者给予我一个模板吗?我有一个很大的矩阵,翻译不好。谢谢
1条答案
按热度按时间yrefmtwq1#
进行一系列大胆的猜测,特别是关于将目标索引转换为约束索引的猜测,并在Python/PuLP中进行演示(您提到您在其他语言中尝试了这一点,但没有指定任何关于您尝试的内容)-这可以简单到