TY - JOUR A1 - Yokoyama, Ryohei A1 - Shinano, Yuji A1 - Taniguchi, Syusuke A1 - Wakui, Tetsuya T1 - Search for K-best solutions in optimal design of energy supply systems by an extended MILP hierarchical branch and bound method JF - Energy N2 - For the purpose of making a decision in the optimal design of an energy supply system, it is important to investigate not only the optimal solution but also suboptimal ones which follow the optimal one without any omissions, what are called K-best solutions. In this paper, a mixed-integer linear programming method utilizing the hierarchical relationship between design and operation variables proposed previ- ously is extended to search the K-best solutions very efficiently. In addition, methods for updating the incumbents are incorporated into the extended method for three options for the criterion set newly in deriving the K-best solutions. This extended method is implemented into open and commercial MILP solvers, and is applied to illustrative and practical case studies, respectively, on the optimal design of cogeneration systems. Through the studies, it turns out that the proposed method is much superior in terms of solution optimality and computation efficiency to a conventional method, and that the computation efficiency to derive one of the K-best solutions by the proposed method increases with the number of K-best solutions. In addition, features of the K-best solutions in the value of objective function are clarified. Y1 - 2019 U6 - https://doi.org/https://doi.org/10.1016/j.energy.2018.02.077 VL - 184 SP - 45 EP - 57 ER - TY - JOUR A1 - Yokoyama, Ryohei A1 - Hiramatsu, Yuto A1 - Shinano, Yuji A1 - Wakui, Tetsuya T1 - Evaluation of robustness in multiple performance criteria for designing energy supply systems based on a mixed-integer linear model JF - Proceedings of the 35th International Conference on Efficiency, Cost, Optimization, Simulation and Environmental Impact of Energy Systems (ECOS 2022) N2 - In designing energy supply systems, it is important to consider the uncertainty in energy demands, evaluate the robustness in some performance criteria, and heighten the robustness. A robust optimal design method has been previously proposed to maximize the robustness in a single performance criterion against the uncertainty in energy demands based on a mixed-integer linear model. In this paper, as a preliminary step toward multiobjective robust optimal design, a method of evaluating robustness in multiple performance criteria against the uncertainty in energy demands is proposed based on a mixed-integer linear model. The problems of evaluating the robustness in the performance criteria are formulated as bilevel mixed-integer linear programming (MILP) ones. They are solved by evaluating lower and upper bounds for the maximum regrets in the performance criteria alternately and repeatedly with the aid of a solution method based on reformulation and decomposition. In addition, for the purpose of applying the proposed method to practical problems, a hierarchical MILP method is used to efficiently solve some MILP problems in the solution process. Through a case study on a cogeneration system, the robustness in the annual total cost and primary energy consumption is evaluated and its trade-off relationship is clarified. Y1 - 2022 SP - 1937 EP - 1948 ER - TY - JOUR A1 - Yokoyama, Ryohei A1 - Shinano, Yuji A1 - Takeuchi, Kotaro A1 - Wakui, Tetsuya T1 - Operation-based time-period clustering for optimal design of energy supply systems by a hierarchical MILP method JF - THE 32ND INTERNATIONAL CONFERENCE ON EFFICIENCY, COST, OPTIMIZATION, SIMULATION AND ENVIRONMENTAL IMPACT OF ENERGY SYSTEMS(ECOS 2019) N2 - Mixed-integer linear programming (MILP) methods have been applied widely to optimal design of energy supply systems. A hierarchical MILP method has been proposed to solve such optimal design problems efficiently. In addition, some strategies have been proposed to enhance the computation efficiency furthermore. As one of the strategies, a method of reducing model by time aggregation has been proposed to search design candidates efficiently in the relaxed optimal design problem at the upper level. In this paper, a method of clustering periods has been proposed based on the optimal operational strategies of the systems to avoid a large decrease in the lower bound for the optimal value of the objective function by model reduction. This method has been realized only by solving the relaxed optimal design problem at the upper level in advance. The method can decrease the number of operation variables and constraints at the upper level, and thus can decrease the computation time at the upper level. Through a case study on the optimal design of a gas turbine cogeneration system, it is clarified how the proposed clustering method is effective to enhance the computation efficiency in comparison with the conventional one which clusters periods regularly in time series. Y1 - 2019 SP - 527 EP - 539 ER - TY - JOUR A1 - Rehfeldt, Daniel A1 - Koch, Thorsten A1 - Shinano, Yuji T1 - Faster exact solution of sparse MaxCut and QUBO problems JF - Mathematical Programming Computation N2 - The maximum-cut problem is one of the fundamental problems in combinatorial optimization. With the advent of quantum computers, both the maximum-cut and the equivalent quadratic unconstrained binary optimization problem have experienced much interest in recent years. This article aims to advance the state of the art in the exact solution of both problems—by using mathematical programming techniques. The main focus lies on sparse problem instances, although also dense ones can be solved. We enhance several algorithmic components such as reduction techniques and cutting-plane separation algorithms, and combine them in an exact branch-and-cut solver. Furthermore, we provide a parallel implementation. The new solver is shown to significantly outperform existing state-of-the-art software for sparse maximum-cut and quadratic unconstrained binary optimization instances. Furthermore, we improve the best known bounds for several instances from the 7th DIMACS Challenge and the QPLIB, and solve some of them (for the first time) to optimality. Y1 - 2023 U6 - https://doi.org/10.1007/s12532-023-00236-6 VL - 15 SP - 445 EP - 470 ER - TY - JOUR A1 - Yokoyama, Ryohei A1 - Shinano, Yuji A1 - Wakui, Tetsuya T1 - Day and time aggregations for optimal design of energy supply systems with storage units by a hierarchical MILP method JF - The 35th International Conference on Efficiency, Cost, Optimization, Simulation and Environmental Impact of Energy Systems (ECOS 2022) N2 - Mixed-integer linear programming (MILP) methods have been applied widely to optimal design of energy supply systems. A hierarchical MILP method has been proposed to solve such optimal design problems efficiently. In addition, some strategies have been proposed to enhance the computation efficiency furthermore. As one of the strategies, a method of reducing model by time aggregation has been proposed to search design candidates efficiently in the relaxed optimal design problem at the upper level. In this paper, the hierarchical MILP method with the strategies is extendedly applied to the optimal design of energy supply systems with storage units. Especially, the method of reducing model is extended by aggregating both representative days and sampling times separately in consideration of the characteristics of energy storage units. A case study is conducted on the optimal design of a cogeneration system with a thermal storage tank. Through the study, it turns out the hierarchical MILP method is effective to derive the optimal solutions in short computation times. It also turns out that the model reduction with day and time aggregations is effec- tive to shorten the computation times furthermore when the number of candidates for equipment capacities is relatively small. Y1 - 2022 SP - 1891 EP - 1902 ER - TY - JOUR A1 - Yokoyama, Ryohei A1 - Shinano, Yuji A1 - Wakui, Tetsuya T1 - Day and time aggregations for optimal design of energy supply systems with storage units by a hierarchical MILP method JF - PROCEEDINGS OF ECOS 2022 - THE 35 TH INTERNATIONAL CONFERENCE ON EFFICIENCY, COST, OPTIMIZATION, SIMULATION AND ENVIRONMENTAL IMPACT OF ENERGY SYSTEMS N2 - Mixed-integer linear programming (MILP) methods have been applied widely to optimal design of energy supply systems. A hierarchical MILP method has been proposed to solve such optimal design problems efficiently. In addition, some strategies have been proposed to enhance the computation efficiency furthermore. As one of the strategies, a method of reducing model by time aggregation has been proposed to search design candidates efficiently in the relaxed optimal design problem at the upper level. In this paper, the hierarchical MILP method with the strategies is extendedly applied to the optimal design of energy supply systems with storage units. Especially, the method of reducing model is extended by aggregating both representative days and sampling times separately in consideration of the characteristics of energy storage units. A case study is conducted on the optimal design of a cogeneration system with a thermal storage tank. Through the study, it turns out the hierarchical MILP method is effective to derive the optimal solutions in short computation times. It also turns out that the model reduction with day and time aggregations is effective to shorten the computation times furthermore when the number of candidates for equipment capacities is relatively small. Y1 - 2022 SP - 1891 EP - 1902 ER - TY - GEN A1 - Rehfeldt, Daniel A1 - Fujisawa, Katsuki A1 - Koch, Thorsten A1 - Nakao, Masahiro A1 - Shinano, Yuji T1 - Computing single-source shortest paths on graphs with over 8 trillion edges N2 - This paper introduces an implementation for solving the single-source shortest path problem on distributed-memory machines. It is tailored to power-law graphs and scales to trillions of edges. The new implementation reached 2nd and 10th place in the latest Graph500 benchmark in June 2022 and handled the largest and second-largest graphs among all participants. T3 - ZIB-Report - 22-22 Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-88180 SN - 1438-0064 ER - TY - GEN A1 - Shinano, Yuji A1 - Vigerske, Stefan T1 - Smoothie: Mixing the strongest MIP solvers to solve hard MIP instances on supercomputers - Phase I development N2 - Mixed-Integer Linear Programming (MIP) is applicable to such a wide range of real-world decision problems that the competition for the best code to solve such problems has lead to tremendous progress over the last decades. While current solvers can solve some of the problems that seemed completely out-of-reach just 10 years ago, there are always relevant MIP problems that currently cannot be solved. With the Smoothie solver we intend to solve extremely hard MIP problems by building on the many years that went into the development of several state-of-the-art MIP solvers and by utilizing some of the largest computing resources available. The high-level task parallelization framework UG (Ubiquity Generator) is used and extended by Smoothie to build a solver that uses large-scale parallelization to distribute the solution of a single MIP on a shared- or distributed-memory computing infrastructure, thereby employing several established MIP solvers simultaneously. For the first development phase, which is the topic of this report, both FICO Xpress and Gurobi are used in concurrent mode on a single machine, while information on incumbent solutions and explored branch-and-bound subtrees is exchanged. A dynamic restarting mechanism ensures that solver configurations are selected that promise most suitable for the MIP to be solved. We report on initial findings using this early version of Smoothie on unsolved problems from MIPLIB 2017. T3 - ZIB-Report - 25-18 KW - MIPLIB, MIP, UG, Parallelization Y1 - 2025 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-101906 SN - 1438-0064 ER - TY - JOUR A1 - Shinano, Yuji A1 - Fujii, Koichi T1 - Introduction to Quadratic Programming Solvers JF - Journal of The Society of Instrument and Control Engineers Y1 - 2025 U6 - https://doi.org/10.11499/sicejl.64.532 VL - 64 IS - 8 SP - 532 EP - 540 ER - TY - JOUR A1 - Yokoyama, Ryohei A1 - Shinano, Yuji A1 - Wakui, Tetsuya T1 - An effective approach for deriving and evaluating approximate optimal design solutions of energy supply systems by time series aggregation JF - Frontiers in Energy Research N2 - It is important to design multi-energy supply systems optimally in consideration of their operations for variations in energy demands. An approach for efficiently solving such an optimal design problem with a large number of periods for variations in energy demands is to derive an approximate optimal design solution by time series aggregation. However, such an approach does not provide any information on the accuracy for the optimal value of the objective function. In this paper, an effective approach for time series aggregation is proposed to derive an approximate optimal design solution and evaluate a proper gap between the upper and lower bounds for the optimal value of the objective function based on a mixed-integer linear model. In accordance with aggregation, energy demands are relaxed to uncertain parameters and the problem for deriving an approximate optimal design solution and evaluating it is transformed to a three-level optimization problem, and it is solved by applying both the robust and hierarchical optimization methods. A case study is conducted on a cogeneration system with a practical configuration, and it turns out that the proposed approach enables one to derive much smaller gaps as compared with those obtained by a conventional approach. Y1 - 2023 U6 - https://doi.org/10.3389/fenrg.2023.1128681 VL - 11 SP - 1 EP - 14 ER - TY - JOUR A1 - Tateiwa, Nariaki A1 - Shinano, Yuji A1 - Yasuda, Masaya A1 - Kaji, Shizuo A1 - Yamamura, Keiichiro A1 - Fujisawa, Katsuki T1 - Development and analysis of massive parallelization of a lattice basis reduction algorithm JF - Japan Journal of Industrial and Applied Mathematics N2 - The security of lattice-based cryptography relies on the hardness of solving lattice problems. Lattice basis reduction is a strong tool for solving lattice problems, and the block Korkine–Zolotarev (BKZ) reduction algorithm is the de facto standard in cryptanalysis. We propose a parallel algorithm of BKZ-type reduction based on randomization. Randomized copies of an input lattice basis are independently reduced in parallel, while several basis vectors are shared asynchronously among all processes. There is a trade-off between randomization and information sharing; if a substantial amount of information is shared, all processes might work on the same problem, which diminishes the benefit of parallelization. To monitor the balance between randomness and sharing, we propose a new metric to quantify the variety of lattice bases, and we empirically find an optimal parameter of sharing for high-dimensional lattices. We also demonstrate the effectiveness of our parallel algorithm and metric through experiments from multiple perspectives. Y1 - 2023 U6 - https://doi.org/10.1007/s13160-023-00580-z VL - 40 IS - 1 SP - 13 EP - 56 ER - TY - CHAP A1 - Yokoyama, Ryohei A1 - Kamada, Hiroki A1 - Shinano, Yuji A1 - Wakui, Tetsuya T1 - A Hierarchical Optimization Approach to Robust Design of Energy Supply Systems Based on a Mixed-Integer Linear Model T2 - Conference N2 - In designing energy supply systems, designers should heighten the robustness in performance criteria against the uncertainty in energy demands. In this paper, a robust optimal design method using a hierarchi- cal mixed-integer linear programming (MILP) method is proposed to maximize the robustness of energy sup- ply systems under uncertain energy demands based on a mixed-integer linear model. A robust optimal design problem is formulated as a three-level min-max-min MILP one by expressing uncertain energy demands by intervals, evaluating the robustness in a performance criterion based on the minimax regret cri- terion, and considering relationships among integer design variables, uncertain energy demands, and inte- ger and continuous operation variables. This problem is solved by evaluating upper and lower bounds for the minimum of the maximum regret of the performance criterion repeatedly outside, and evaluating lower and upper bounds for the maximum regret repeatedly inside. Since these different types of optimization problems are difficult to solve even using commercial MILP solvers, they are solved by applying a hierarchi- cal MILP method developed for ordinary optimal design problems with its modifications. In a case study, the proposed approach is applied to the robust optimal design of a cogeneration system. Through the study, its validity and effectiveness are ascertained, and some features of the obtained robust designs are clarified. Y1 - 2020 SP - 601 EP - 613 ER - TY - CHAP A1 - Yokoyama, Ryohei A1 - Takeuchi, Kotaro A1 - Shinano, Yuji A1 - Wakui, Tetsuya T1 - Effect of Model Reduction by Time Aggregation in Multiobjective Optimal Design of Energy Supply Systems by a Hierarchical MILP Method T2 - Conference N2 - The mixed-integer linear programming (MILP) method has been applied widely to optimal design of energy supply systems. A hierarchical MILP method has been proposed to solve such optimal design problems effi- ciently. As one of the strategies to enhance the computation efficiency furthermore, a method of reducing model by time aggregation has been proposed to search design candidates accurately and efficiently in the relaxed optimal design problem at the upper level. In this paper, the hierarchical MILP method and model reduction by time aggregation are applied to the multiobjective optimal design. In applying the model reduc- tion, the methods of clustering periods by the order of time series, based on an operational strategy, and by the k-medoids method are applied. As a case study, the multiobjective optimal design of a gas turbine cogeneration system with a practical configuration is investigated by adopting the annual total cost and pri- mary energy consumption as the objective functions to be minimized simultaneously, and the clustering methods are compared with one another in terms of the computation efficiency. It turns out that the model reduction by any clustering method is effective to enhance the computation efficiency when importance is given to minimizing the first objective function. It also turns out that the model reduction only by the k- medoids method is effective very limitedly when importance is given to minimizing the second objective function. Y1 - 2020 SP - 627 EP - 639 ER - TY - CHAP A1 - Koch, Thorsten A1 - Rehfeldt, Daniel A1 - Shinano, Yuji T1 - An exact high performance solver for Steiner tree problems in graphs and related problems T2 - Modeling, Simulation and Optimization of Complex Processes HPSC 2018 Y1 - 2020 PB - Springer ER - TY - GEN A1 - Shinano, Yuji A1 - Achterberg, Tobias A1 - Berthold, Timo A1 - Heinz, Stefan A1 - Koch, Thorsten A1 - Winkler, Michael T1 - Solving Previously Unsolved MIP Instances with ParaSCIP on Supercomputers by using up to 80,000 Cores N2 - Mixed-integer programming (MIP) problem is arguably among the hardest classes of optimization problems. This paper describes how we solved 21 previously unsolved MIP instances from the MIPLIB benchmark sets. To achieve these results we used an enhanced version of ParaSCIP, setting a new record for the largest scale MIP computation: up to 80,000 cores in parallel on the Titan supercomputer. In this paper, we describe the basic parallelization mechanism of ParaSCIP, improvements of the dynamic load balancing and novel techniques to exploit the power of parallelization for MIP solving. We give a detailed overview of computing times and statistics for solving open MIPLIB instances. T3 - ZIB-Report - 20-16 KW - Mixed Integer Programming, Parallel processing, Node merging, Racing, ParaSCIP, Ubiquity Generator Framework, MIPLIB Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-78393 SN - 1438-0064 ER - TY - JOUR A1 - Hosada, Junko A1 - Maher, Stephen A1 - Shinano, Yuji A1 - Villumsen, Jonas Christoffer T1 - Location, transshipment and routing: An adaptive transportation network integrating long-haul and local vehicle routing JF - EURO Journal on Transportation and Logistics Y1 - 2022 U6 - https://doi.org/10.1016/j.ejtl.2022.100091 VL - 11 PB - Elsevier ER - TY - GEN A1 - Hosoda, Junko A1 - Maher, Stephen J. A1 - Shinano, Yuji T1 - Combining Benders’ decomposition and parallelisation to solve large-scale adaptive network supply chain management problems N2 - Integrated supply chain network design—combining warehouse clustering with inter- and intra-cluster vehicle routing—is a fundamental feature of supply chain management problems. In this paper, a Benders’ decomposition-based parallel heuristic approach for the integrated long-haul and local vehicle routing problem is proposed. The proposed algorithm comprises lower- and upper-bound search algorithms running in parallel to find high-quality solutions for challenging supply chain management problems. A novel Benders’ decomposition approach is proposed to aid the improvement in the lower bound for large-scale problems. The results demonstrate that the Benders’ decomposition approach is effective in finding strong lower bounds and is beneficial for obtaining high-quality solutions. T3 - ZIB-Report - 24-12 Y1 - 2024 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-98744 SN - 1438-0064 ER - TY - GEN A1 - Fujii, Koichi A1 - Kim, Sunyoung A1 - Kojima, Masakazu A1 - Mittelmann, Hans D. A1 - Shinano, Yuji T1 - The Largest Unsolved QAP Instance Tai256c Can Be Converted into A 256-dimensional Simple BQOP with A Single Cardinality Constraint N2 - Tai256c is the largest unsolved quadratic assignment problem (QAP) instance in QAPLIB; a 1.48% gap remains between the best known feasible objective value and lower bound of the unknown optimal value. This paper shows that the instance can be converted into a 256 dimensional binary quadratic optimization problem (BQOP) with a single cardinality constraint which requires the sum of the binary variables to be 92.The converted BQOP is much simpler than the original QAP tai256c and it also inherits some of the symmetry properties. However, it is still very difficult to solve. We present an efficient branch and bound method for improving the lower bound effectively. A new lower bound with 1.36% gap is also provided. T3 - ZIB-Report - 22-18 Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-88086 SN - 1438-0064 ER - TY - GEN A1 - Rehfeldt, Daniel A1 - Koch, Thorsten A1 - Shinano, Yuji T1 - Faster exact solution of sparse MaxCut and QUBO problems T3 - ZIB-Report - 22-02 Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-85715 SN - 1438-0064 ER - TY - GEN A1 - Fujisawa, Katsuki A1 - Shinano, Yuji A1 - Waki, Hayato T1 - Optimization in the Real World T2 - Mathematics for Industry (MFI) Y1 - 2015 SN - 978-4-431-55419-6 U6 - https://doi.org/10.1007/978-4-431-55420-2 VL - 13 PB - Springer Tokyo Heidelberg New York Dordrecht London ER - TY - GEN A1 - Ralphs, Ted A1 - Shinano, Yuji A1 - Berthold, Timo A1 - Koch, Thorsten ED - Hamadi, Youssef T1 - Parallel Solvers for Mixed Integer Linear Optimization T2 - Handbook of Parallel Constraint Reasoning Y1 - 2018 U6 - https://doi.org/10.1007/978-3-319-63516-3_8 SP - 283 EP - 336 PB - Springer Nature ER -