规模稍大 GLPK 就跪，请问是多大规模会跪啊...

等于没说。。。这也是教程？骗谁的？

Quote: 引用 5 楼 lenxx 的回复:

规模稍大 GLPK 就跪，请问是多大规模会跪啊...

我都是用最大规模来调试的，gplk最大规模的一个都过不了，lpsolve有的会超时，cbc相对来说比较稳定，基本能在10s出解甚至最优

Quote: 引用 6 楼 open_source_library 的回复:

Quote: 引用 5 楼 lenxx 的回复:

规模稍大 GLPK 就跪，请问是多大规模会跪啊...

我都是用最大规模来调试的，gplk最大规模的一个都过不了，lpsolve有的会超时，cbc相对来说比较稳定，基本能在10s出解甚至最优

请查收私信，交流一下。

针对这个题目，我们内部有一套数学模型可以保证找出最优解，但是速度是分钟级别的。如果楼主能够找到一个最优解情况下速度又很快的方案，我们也很期待。

针对这个题目，我们内部有一套数学模型可以保证找出最优解，但是速度是分钟级别的。如果楼主能够找到一个确定找到最优解、速度又很快的方案，我们也很期待。

首先申明本人不用，也不去用开源库的。只是为了帮大家谋福利，激激楼主。。。。还有就是严重怀疑楼主，说求解器非常耗时，但是能解出最优解，而且在很耗时的基础上，能进赛区前三。这不太可能吧(本人江苏赛区的，所以对一些前几名的解的状况还是有所了解的)。。。

楼主可否推荐几篇建模论文参考？

请叫自己抄模型而不是写模型

请叫自己抄模型而不是写模型

Quote: 引用 11 楼 frostbolt 的回复:

针对这个题目，我们内部有一套数学模型可以保证找出最优解，但是速度是分钟级别的。如果楼主能够找到一个最优解情况下速度又很快的方案，我们也很期待。

我们所有用例都在10S内算出最优解了，可是复赛不让用开源库了，您的肯定是对我们最好的鼓励！

我不知道求解器的速度如何，但是我们组没有用，我们用了一篇论文里的算法解决该整数规划问题，复杂度为为O(n^3)，对于600个点的，只需要0.04ms，求出来的路径可能是成环的。求解器是针对一般情况的，对于这个问题，存在有更高效率的算法。PS：我们的时间主要花在去除环路上

我不知道求解器的速度如何，但是我们组没有用，我们用了一篇论文里的算法解决该整数规划问题，复杂度为为O(n^3)，对于600个点的，只需要0.04ms，求出来的路径可能是成环的。求解器是针对一般情况的，对于这个问题，存在有更高效率的算法。PS：我们的时间主要花在去除环路上

Quote: 引用 11 楼 frostbolt 的回复:

针对这个题目，我们内部有一套数学模型可以保证找出最优解，但是速度是分钟级别的。如果楼主能够找到一个确定找到最优解、速度又很快的方案，我们也很期待。

谢谢版主支持！当听到复赛不能用求解器的时候，我就没修改代码了。我有一个自认为很好的方案，是结合搜索和模型的，可能的话，我希望比赛结束后好好实现和优化，希望有机会和赛题组的老师们交流一下。

好啊，我们也希望多和大家聊聊天，你的队伍名称方便发私信给我下吧。
另外，我们本次比赛结束以后会有一个悬赏题，奖金更高，题目是真实存在的的、我们内部还没有完全解决的问题，对解题方式不做任何限制，有兴趣请继续关注，谢谢。。

Quote: 引用 13 楼 open_source_library 的回复:

Quote: 引用 11 楼 frostbolt 的回复:

针对这个题目，我们内部有一套数学模型可以保证找出最优解，但是速度是分钟级别的。如果楼主能够找到一个确定找到最优解、速度又很快的方案，我们也很期待。

谢谢版主支持！当听到复赛不能用求解器的时候，我就没修改代码了。我有一个自认为很好的方案，是结合搜索和模型的，可能的话，我希望比赛结束后好好实现和优化，希望有机会和赛题组的老师们交流一下。

好啊，我们也希望多和大家聊聊天，你的队伍名称方便发私信给我下吧。
另外，我们本次比赛结束以后会有一个悬赏题，奖金更高，题目是真实存在的的、我们内部还没有完全解决的问题，对解题方式不做任何限制，有兴趣请继续关注，谢谢。。

Quote: 引用 29 楼 frostbolt 的回复:

好啊，我们也希望多和大家聊聊天，你的队伍名称方便发私信给我下吧。
另外，我们本次比赛结束以后会有一个悬赏题，奖金更高，题目是真实存在的的、我们内部还没有完全解决的问题，对解题方式不做任何限制，有兴趣请继续关注，谢谢。。

版主，能有机会接触下这个悬赏题目吗？

针对这个题目，我们内部有一套数学模型可以保证找出最优解，但是速度是分钟级别的。如果楼主能够找到一个确定找到最优解、速度又很快的方案，我们也很期待。

Quote: 引用 11 楼 frostbolt 的回复:

针对这个题目，我们内部有一套数学模型可以保证找出最优解，但是速度是分钟级别的。如果楼主能够找到一个确定找到最优解、速度又很快的方案，我们也很期待。

我只想问，如果都用楼主这个方案 是不是该当作弊处理？ 如果不当作弊处理，那么我们这些自己写算法的还有一丝机会吗？ 如果楼主是为了分享知识，大可以在复活赛后分享，这个时候分享算什么？ 让我们不用开源库的全部被淘汰吗？ 你把模型都给出来了，请问 使用开源库的人还需要什么技能？ 也请大赛组委想想，不使用开源库的人还有一丝机会吗？这也算公平竞技吗？ 

Quote: 引用 7 楼 uniquew1 的回复:

等于没说。。。这也是教程？骗谁的？

我已经告诉了你最好的求解器是cbc，也告诉怎么建立一个整数规划模型，你还需要什么？直接上代码？这个不是很简单的事情？

所有的开源库，都会有example的，照壶画瓢就行啊。。

比如谷歌搜“coiner cbc example”，

#include "CbcModel.hpp"

// Using CLP as the solver
#include "OsiClpSolverInterface.hpp"

int main (int argc, const char *argv[])
{
  OsiClpSolverInterface solver1;

  // Read in example model in MPS file format
  // and assert that it is a clean model
  int numMpsReadErrors = solver1.readMps("../../Mps/Sample/p0033.mps","");
  assert(numMpsReadErrors==0);

  // Pass the solver with the problem to be solved to CbcModel 
  CbcModel model(solver1);

  // Do complete search
  model.branchAndBound();

  /* Print the solution.  CbcModel clones the solver so we
     need to get current copy from the CbcModel */
  int numberColumns = model.solver()->getNumCols();
    
  const double * solution = model.bestSolution();
    
  for (int iColumn=0;iColumn<numberColumns;iColumn++) {
    double value=solution[iColumn];
    if (fabs(value)>1.0e-7&&model.solver()->isInteger(iColumn)) 
      printf("%d has value %g\n",iColumn,value);
   }
  return 0;
}    


这些就是。。如果这些都不是教程的话，我觉得可能只能等比赛结束直接奉上我的代码才算教程了

Quote: 引用 7 楼 uniquew1 的回复:

等于没说。。。这也是教程？骗谁的？

我已经告诉了你最好的求解器是cbc，也告诉怎么建立一个整数规划模型，你还需要什么？直接上代码？这个不是很简单的事情？

所有的开源库，都会有example的，照壶画瓢就行啊。。

比如谷歌搜“coiner cbc example”，

#include "CbcModel.hpp"

// Using CLP as the solver
#include "OsiClpSolverInterface.hpp"

int main (int argc, const char *argv[])
{
  OsiClpSolverInterface solver1;

  // Read in example model in MPS file format
  // and assert that it is a clean model
  int numMpsReadErrors = solver1.readMps("../../Mps/Sample/p0033.mps","");
  assert(numMpsReadErrors==0);

  // Pass the solver with the problem to be solved to CbcModel 
  CbcModel model(solver1);

  // Do complete search
  model.branchAndBound();

  /* Print the solution.  CbcModel clones the solver so we
     need to get current copy from the CbcModel */
  int numberColumns = model.solver()->getNumCols();
    
  const double * solution = model.bestSolution();
    
  for (int iColumn=0;iColumn<numberColumns;iColumn++) {
    double value=solution[iColumn];
    if (fabs(value)>1.0e-7&&model.solver()->isInteger(iColumn)) 
      printf("%d has value %g\n",iColumn,value);
   }
  return 0;
}    


这些就是。。如果这些都不是教程的话，我觉得可能只能等比赛结束直接奉上我的代码才算教程了