@misc{TateiwaShinanoYasudaetal.2021, author = {Tateiwa, Nariaki and Shinano, Yuji and Yasuda, Masaya and Kaji, Shizuo and Yamamura, Keiichiro and Fujisawa, Katsuki}, title = {Massively parallel sharing lattice basis reduction}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-85209}, year = {2021}, abstract = {For cryptanalysis in lattice-based schemes, the performance evaluation of lattice basis reduction using high-performance computers is becoming increasingly important for the determination of the security level. We propose a distributed and asynchronous parallel reduction algorithm based on randomization and DeepBKZ, which is an improved variant of the block Korkine-Zolotarev (BKZ) reduction algorithm. Randomized copies of a lattice basis are distributed to up to 103,680 cores and independently reduced in parallel, while some basis vectors are shared asynchronously among all processes via MPI. There is a trade-off between randomization and information sharing; if a substantial amount of information is shared, all processes will work on the same problem, thereby diminishing the benefit of parallelization. To monitor this balance between randomness and sharing, we propose a metric to quantify the variety of lattice bases. We empirically find an optimal parameter of sharing for high-dimensional lattices. We demonstrate the efficacy of our proposed parallel algorithm and implementation with respect to both performance and scalability through our experiments.}, language = {en} } @misc{FujiiKimKojimaetal.2023, author = {Fujii, Koichi and Kim, Sunyoung and Kojima, Masakazu and Mittelmann, Hans D. and Shinano, Yuji}, title = {An Exceptionally Difficult Binary Quadratic Optimization Problem with Symmetry: a Challenge for The Largest Unsolved QAP Instance Tai256c}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-93072}, year = {2023}, abstract = {Tai256c is the largest unsolved quadratic assignment problem (QAP) instance in QAPLIB. It is known that QAP tai256c 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. As the BQOP is much simpler than the original QAP, the conversion increases the possibility to solve the QAP. Solving exactly the BQOP, however, is still very difficult. Indeed, a 1.48\% gap remains between the best known upper bound (UB) and lower bound (LB) of the unknown optimal value. This paper shows that the BQOP admits a nontrivial symmetry, a property that makes the BQOP very hard to solve. The symmetry induces equivalent subproblems in branch and bound (BB) methods. To effectively improve the LB, we propose an efficient BB method that incorporates a doubly nonnegative relaxation, the standard orbit branching and a technique to prune equivalent subproblems. With this BB method, a new LB with 1.25\% gap is successfully obtained, and computing an LB with 1.0\% gap is shown to be still quite difficult.}, language = {en} } @article{YokoyamaKamadaShinanoetal.2021, author = {Yokoyama, Ryohei and Kamada, Hiroki and Shinano, Yuji and Wakui, Tetsuya}, title = {A hierarchical optimization approach to robust design of energy supply systems based on a mixed-integer linear model}, volume = {229}, journal = {Energy}, doi = {https://doi.org/10.1016/j.energy.2021.120343}, year = {2021}, abstract = {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 hierarchical mixed-integer linear programming (MILP) method is proposed to maximize the robustness of energy supply 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 criterion, and considering relationships among integer design variables, uncertain energy demands, and integer 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. Different types of optimization problems are solved by applying a hierarchical MILP method developed for ordinary optimal design problems without and 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.}, language = {en} } @article{FujiiKimKojimaetal.2024, author = {Fujii, Koichi and Kim, Sunyoung and Kojima, Masakazu and Mittelmann, Hans D. and Shinano, Yuji}, title = {An exceptionally difficult binary quadratic optimization problem with symmetry: a challenge for the largest unsolved QAP instance Tai256c}, journal = {Optimization Letters}, publisher = {Springer Science and Business Media LLC}, issn = {1862-4472}, doi = {10.1007/s11590-024-02157-2}, year = {2024}, abstract = {Tai256c is the largest unsolved quadratic assignment problem (QAP) instance in QAPLIB. It is known that QAP tai256c 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. As the BQOP is much simpler than the original QAP, the conversion increases the possibility to solve the QAP. Solving exactly the BQOP, however, is still very difficult. Indeed, a 1.48\% gap remains between the best known upper bound (UB) and lower bound (LB) of the unknown optimal value. This paper shows that the BQOP admits a nontrivial symmetry, a property that makes the BQOP very hard to solve. Despite this difficulty, it is imperative to decrease the gap in order to ultimately solve the BQOP exactly. To effectively improve the LB, we propose an efficient BB method that incorporates a doubly nonnegative relaxation, the orbit branching and the isomorphism pruning. With this BB method, a new LB with 1.25\% gap is successfully obtained, and computing an LB with gap is shown to be still quite difficult.}, language = {en} } @article{HosodaMaherShinanoetal.2024, author = {Hosoda, Junko and Maher, Stephen J. and Shinano, Yuji and Villumsen, Jonas Christoffer}, title = {A parallel branch-and-bound heuristic for the integrated long-haul and local vehicle routing problem on an adaptive transportation network}, volume = {165}, journal = {Computers \& Operations Research}, publisher = {Elsevier BV}, issn = {0305-0548}, doi = {10.1016/j.cor.2024.106570}, year = {2024}, abstract = {Consolidation of commodities and coordination of vehicle routes are fundamental features of supply chain management problems. While locations for consolidation and coordination are typically known a priori, in adaptive transportation networks this is not the case. The identification of such consolidation locations forms part of the decision making process. Supply chain management problems integrating the designation of consolidation locations with the coordination of long haul and local vehicle routing is not only challenging to solve, but also very difficult to formulate mathematically. In this paper, the first mathematical model integrating location clustering with long haul and local vehicle routing is proposed. This mathematical formulation is used to develop algorithms to find high quality solutions. A novel parallel framework is developed that combines exact and heuristic methods to improve the search for high quality solutions and provide valid bounds. The results demonstrate that using exact methods to guide heuristic search is an effective approach to find high quality solutions for difficult supply chain management problems.}, language = {en} } @article{YokoyamaTakeuchiShinanoetal.2021, author = {Yokoyama, Ryohei and Takeuchi, Kotaro and Shinano, Yuji and Wakui, Tetsuya}, title = {Effect of model reduction by time aggregation in multiobjective optimal design of energy supply systems by a hierarchical MILP method}, volume = {228}, journal = {Energy}, doi = {https://doi.org/10.1016/j.energy.2021.120505}, year = {2021}, abstract = {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 efficiently. In addition, a method of reducing model by time aggregation has been proposed to search design candidates accurately and efficiently at the upper level. In this paper, the hierarchical MILP method and model reduction by time aggregation are applied to the multiobjective optimal design. The methods of clustering periods by the order of time series, by the k-medoids method, and based on an operational strategy are applied for the model reduction. As a case study, the multiobjective optimal design of a gas turbine cogeneration system is investigated by adopting the annual total cost and primary energy consumption as the objective functions, 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, but that the model reduction only by the k-medoids method is effective very limitedly when importance is given to minimizing the second objective function.}, language = {en} } @article{BertholdKochShinano2021, author = {Berthold, Timo and Koch, Thorsten and Shinano, Yuji}, title = {MILP. Try. Repeat.}, volume = {2}, journal = {Proceedings of the 13th International Conference on the Practice and Theory of Automated Timetabling - PATAT 2021}, year = {2021}, language = {en} } @misc{TateiwaShinanoYamamuraetal.2021, author = {Tateiwa, Nariaki and Shinano, Yuji and Yamamura, Keiichiro and Yoshida, Akihiro and Kaji, Shizuo and Yasuda, Masaya and Fujisawa, Katsuki}, title = {CMAP-LAP: Configurable Massively Parallel Solver for Lattice Problems}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-82802}, year = {2021}, abstract = {Lattice problems are a class of optimization problems that are notably hard. There are no classical or quantum algorithms known to solve these problems efficiently. Their hardness has made lattices a major cryptographic primitive for post-quantum cryptography. Several different approaches have been used for lattice problems with different computational profiles; some suffer from super-exponential time, and others require exponential space. This motivated us to develop a novel lattice problem solver, CMAP-LAP, based on the clever coordination of different algorithms that run massively in parallel. With our flexible framework, heterogeneous modules run asynchronously in parallel on a large-scale distributed system while exchanging information, which drastically boosts the overall performance. We also implement full checkpoint-and-restart functionality, which is vital to high-dimensional lattice problems. Through numerical experiments with up to 103,680 cores, we evaluated the performance and stability of our system and demonstrated its high capability for future massive-scale experiments.}, language = {en} } @inproceedings{TateiwaShinanoYamamuraetal.2021, author = {Tateiwa, Nariaki and Shinano, Yuji and Yamamura, Keiichiro and Yoshida, Akihiro and Kaji, Shizuo and Yasuda, Masaya and Fujisawa, Katsuki}, title = {CMAP-LAP: Configurable Massively Parallel Solver for Lattice Problems}, booktitle = {HiPC 2021 proceedings}, year = {2021}, abstract = {Lattice problems are a class of optimization problems that are notably hard. There are no classical or quantum algorithms known to solve these problems efficiently. Their hardness has made lattices a major cryptographic primitive for post-quantum cryptography. Several different approaches have been used for lattice problems with different computational profiles; some suffer from super-exponential time, and others require exponential space. This motivated us to develop a novel lattice problem solver, CMAP-LAP, based on the clever coordination of different algorithms that run massively in parallel. With our flexible framework, heterogeneous modules run asynchronously in parallel on a large-scale distributed system while exchanging information, which drastically boosts the overall performance. We also implement full checkpoint-and-restart functionality, which is vital to high-dimensional lattice problems. Through numerical experiments with up to 103,680 cores, we evaluated the performance and stability of our system and demonstrated its high capability for future massive-scale experiments.}, language = {en} } @article{WakayamaYokoyamaShinanoetal.2018, author = {Wakayama, Yuki and Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {エネルギー供給システムの設計および運用の階層的関係を考慮した最適化(時間集約によるモデル縮約とその効果)}, journal = {エネルギー・資源学会第34回エネルギーシステム・経済・環境コンファレンス講演論文集}, pages = {401 -- 404}, year = {2018}, abstract = {A hierarchical mixed-integer linear programming method to solve optimal design problems of energy supply systems efficiently has been proposed. The original problem is solved by dividing it into a relaxed optimal design problem at the upper level and optimal operation problems which are independent of each other at the lower level. In this paper, a method of reducing model by time aggregation is proposed to search design solution candidates efficiently in the relaxed optimal design problem at the upper level. This method is realized only by clustering periods and averaging energy demands for clustered periods, while it guarantees to derive the optimal solution. On one hand, the method may decrease the number of design variables and costraints at the upper level, and thus the computation time at the upper level. On the other hand, it may increase the numbers of generated design solution candidates and solved optimal operation problems, and thus the computation time at both the levels. Though a case study on the optimal design of a cogeneration system, it is clarified how the model reduction is effective to enhance the computation efficiency.}, language = {ja} } @article{TaniguchiYokoyamaShinanoetal.2015, author = {Taniguchi, Shusuke and Yokoyama, Ryohei and Shinano, Yuji and Ohkura, Masashi and Wakui, Tetsuya}, title = {エネルギー供給システムの設計および運用の階層的関係を考慮した最適化(多目的最適化におけるK-ベスト解の探索)}, journal = {エネルギー・資源学会第31回エネルギーシステム・経済・環境コンファレンス講演論文集}, pages = {267 -- 272}, year = {2015}, abstract = {A mixed-integer linear programming method utilizing the hierarchical relationship between design and operation variables proposed to solve the optimal design problem of energy supply systems efficiently is extended to search K-best solutions in multiobjective optimization: At the upper level, the optimal values of design variables are searched by the branch and bound method with operation variables relaxed to continuous ones; At the lower level, the values of operation variables are optimized independently at the respective periods set for variations in energy demands by the branch and bound method with the values of design variables given tentatively during the search at the upper level. A weighting method is employed for multiobjective optimization, and a weighted sum of the annual total cost and primary energy consumption is adopted as the objective function to be minimized. A practical case study on the optimal design of a cogeneration system is conducted. The validity and effectiveness of the proposed method are clarified, and the trade-off relationship between the annual total cost and primary energy consumption is confirmed.}, language = {ja} } @article{TakeuchiYokoyamaShinanoetal.2019, author = {Takeuchi, Kotaro and Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {エネルギー供給システムの 設計と運用の階層的関係を考慮した最適化 (運用を考慮した期間クラスタリングによるモデル縮約)}, journal = {エネルギー・資源学会第35回エネルギーシステム・経済・環境コンファレンス講演論文集}, pages = {169 -- 174}, year = {2019}, abstract = {To attain the highest performance of energy supply systems, it is necessary to determine design specifications optimally in consideration of operational strategies corresponding to seasonal and hourly variations in energy demands. Mixed-integer linear programming (MILP) methods have been applied widely to such multi-period optimal design problems. A hierarchical MILP method has been proposed to solve the problems very efficiently. In addition, by utilizing features of the hierarchical MILP method, a method of reducing model by clustering periods has also been proposed to search design solution candidates efficiently in the relaxed optimal design problem at the upper level. In this paper, by utilizing features of the hierarchical MILP method, a method of clustering periods is proposed based on the optimal operational strategies of energy supply systems obtained by solving the relaxed optimal design problem. As a case study, the method is applied to the optimal design of a gas turbine cogeneration system, and it is clarified that the method is effective to enhance the computation efficiency in comparison with a conventional method of clustering periods regularly.}, language = {ja} } @article{KamadaYokoyamaShinanoetal.2019, author = {Kamada, Hiroki and Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {混合整数線形モデルによる エネルギー供給システムのロバスト最適設計 (階層的最適化手法の適用)}, journal = {エネルギー・資源学会第35回エネルギーシステム・経済・環境コンファレンス講演論文集}, pages = {163 -- 168}, year = {2019}, abstract = {A robust optimal design method of energy supply systems under uncertain energy demands has been proposed using a mixed- integer linear model for constituent equipment. A robust optimal design problem has been formulated as a three-level min-max- min optimization one by expressing uncertain energy demands by intervals, evaluating the robustness in a performance criterion based on the minimax regret criterion, and considering hierarchical relationships among design variables, uncertain energy demands, and operation variables. However, this method takes a long computation time, and thus it can be applied only to small-scale problems. In this paper, mixed-integer linear programming method in consideration of the hierarchical relationship between design and operation variables is applied to parts of the robust optimal design method which take long computation times to solve problems efficiently. In a case study, this revised method is applied to the robust optimal design of a cogeneration system with a simple configuration, and the validity and effectiveness of the method are ascertained.}, language = {ja} } @article{YokoyamaShinanoTaniguchietal.2014, author = {Yokoyama, Ryohei and Shinano, Yuji and Taniguchi, Syusuke and Ohkura, Masashi and Wakui, Tetsuya}, title = {エネルギー供給システムの設計および運用の階層的関係を考慮した最適化(手法の提案)}, journal = {第33回エネルギー・資源学会研究発表会講演論文要旨集・講演論文集}, pages = {183 -- 188}, year = {2014}, abstract = {A mixed-integer linear programming (MILP) method utilizing the hierarchical relationship between design and op- eration variables is proposed to solve the optimal design problem of energy supply systems efficiently: At the upper level, the optimal values of design variables are searched by the branch and bound method; At the lower level, the values of operation variables are optimized by the branch and bound method under the values of design variables given tentatively during the search at the upper level; Lower bounds for the optimal value of the objective function are evaluated, and are utilized for the bounding operations at both the levels. This method is implemented into open and commercial MILP solvers. Illustrative and practical case studies on the optimal design of cogeneration systems are conducted, and the validity and effectiveness of the proposed method are clarified.}, language = {ja} } @article{YokoyamaShinanoTaniguchietal.2014, author = {Yokoyama, Ryohei and Shinano, Yuji and Taniguchi, Syusuke and Ohkura, Masashi and Wakui, Tetsuya}, title = {エネルギー供給システムの設計および運用階層的関係を考慮した最適化(K-ベスト解の探索)}, journal = {日本機械学会第11回最適化シンポジウム2014講演論文集}, pages = {1 -- 4}, year = {2014}, abstract = {A mixed-integer linear programming methodutilizing the hierarchical relationship between design and operation variables proposed ot solve the optimal design problem of energy supply systems efficiently is extended ot search K-best solutions: At the upper level, the optimal values ofdesign variables are searched with operation variables relaxed to continuous ones; At the lower level, the values of operation variables are optimized with the values of design variables given tentatively; The obtained solution is used to renew K-best incumbent solutions, and the upper bound for the value of the objective function for K-best solutions is replaced correspondingly. A practical case study is conducted, and the validity and effectiveness of the proposed method are clarified.}, language = {ja} } @article{WakayamaYokoyamaShinanoetal.2018, author = {Wakayama, Yuki and Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {階層的最適化によるエネルギ ー供給システムの最適設計・運用 (期間のクラスタリングに よるモデル縮約とその効果)}, journal = {日本機械学会関西支部第93期定時総会講演会講演論文集}, pages = {259 -- 260}, year = {2018}, abstract = {Ahierarchical mixed-integer linear programmingmethod has been proposed to solve optimal design problems of energy supply systems efficiently. In this paper, a method of reducing model by clustering periods is proposed to search design solution candidates efficiently in the relaxed optimal design problem at the upper level. This method is realized only by clustering periods and averaging energy demands for clustered periods, while it guarantees to derive the optimal solution. Through acase study on the optimaldesign of a cogeneration system, ti is clarified how the model reduction si effective ot enhance the computation efficiency.}, language = {ja} } @article{YokoyamaShinanoTakeuchietal.2019, author = {Yokoyama, Ryohei and Shinano, Yuji and Takeuchi, Kotaro and Wakui, Tetsuya}, title = {エネルギー供給システムの設計および運用の 階層的関係を考慮した最適化 (モデル縮約のための期間クラスタリング手法の比較)}, journal = {第38回エネルギー・資源学会研究発表会講演論文集}, pages = {109 -- 114}, year = {2019}, abstract = {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, three clustering methods are applied to time aggregation and compared with one another in terms of the computation efficiency. Especially, the k-medoids method is applied newly in addition to the time-series and operation-based methods applied previously. A case study is conducted on the optimal design of a gas turbine cogeneration system for district energy supply. Through the study, it turns out the k-medoids method is effective to shorten the computation time as compared with the time-series method, although it is necessary to set the number of clusters artifically in both the methods. It also turns out that the operation-based method is more effective than the k-medoids method in terms of the computation efficiency even with the number of clusters set automatically.}, language = {ja} } @article{TakeuchiYokoyamaShinanoetal.2020, author = {Takeuchi, Kotaro and Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {エネルギー供給システムの 設計および運用の階層的関係を考慮した最適化 (多目的最適設計へのモデル縮約の適用)}, journal = {エネルギー・資源学会第36回エネルギーシステム・経済・環境コンファレンス講演論文集}, pages = {724 -- 729}, year = {2020}, abstract = {To attain the highest performance of energy supply systems, it is necessary to determine design specifications optimally in consideration of operational strategies corresponding to seasonal and hourly variations in energy demands. Mixed-integer linear programming (MILP) methods have been applied widely to such optimal design problems. A hierarchical MILP method has been proposed to solve the problems very efficiently. In addition, by utilizing features of the hierarchical MILP method, a method of reducing model by clustering periods based on the optimal operational strategies of equipment has been proposed to search design solution candidates efficiently in the relaxed optimal design problem at the upper level. In this paper, these methods are applied to the multiobjective optimal design of a cogeneration system by considering the annual total cost and primary energy consumption as the objective functions to be minimized. Through a case study, it turns out that the model reduction by the operation-based time-period clustering is effective in terms of the computation efficiency when importance is given to the first objective function, while it is not when importance is given to the second objective function.}, language = {ja} } @article{YokoyamaShinanoWakui2021, author = {Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {エネルギー供給システムの設計および運用の 階層的関係を考慮した最適化 (蓄エネルギー機器を有するシステムへの適用)}, journal = {第40回エネルギー・資源学会研究発表会講演論文集}, pages = {398 -- 403}, year = {2021}, abstract = {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 has been extendedly applied to the optimal design of energy supply systems with storage units. Especially, the method of re- ducing model is extended by aggregating representative days and sampling times differently in consideration of the characteristics of storage units. A case study is conducted on the optimal design of a gas turbine cogeneration system with a thermal storage unit for district energy supply. Through the study, it turns out the hierarchical MILP method is effective to derive the optimal solution as compared with a conventional method. It also turns out that the model reduction with the special time aggregation is effective to shorten the computation time as compared with that without time aggregation in case that the number of candidates for equipment capacities is relatively small.}, language = {ja} } @article{YokoyamaShinanoWakui2022, author = {Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {時間集約によるエネルギー供給システムの 近似最適設計解の導出および評価 (階層的最適化を援用した設計解の評価)}, journal = {エネルギー・資源学会第38回エネルギーシステム・経済・環境コンファレンス講演論文集}, pages = {468 -- 473}, year = {2022}, abstract = {For the purpose of attaining the highest performance of energy supply systems, it is important to design the systems optimally in consideration of their operational strategies for seasonal and hourly variations in energy demands. An ap- proach to solve such an optimal design problem with a large number of periods efficiently is to derive an approximate optimal design solution by aggregating periods with a clustering method. However, such an approach does not provide any information on the accuracy for the optimal value of the objective function. The purpose of this paper is to provide a time aggregation method for deriving aprroximate optimal design solutions and evaluting their values of the objective function. Especially, a method of evaluating design solutions is presented here using both methods of evaluating the robustness under uncertain energy demands and solving optimal design problems by a hierarchical approach. A case study is conducted for a cogeneration system with a practical configuration, and it turns out that the proposed approach enables one to evaluate effective lower bounds for the optimal value of the objective function as compared with those obtained by a conventional approach.}, language = {ja} } @article{YokoyamaShinanoWakui2022, author = {Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {時間集約によるエネルギー供給システムの 近似最適設計解の導出および評価 (階層的最適化を援用した近似最適設計)}, journal = {第41回エネルギー・資源学会研究発表会講演論文集}, pages = {144 -- 148}, year = {2022}, abstract = {For the purpose of attaining the highest performance of energy supply systems, it is important to design the systems optimally in consideration of their operational strategies for seasonal and hourly variations in energy demands. An ap- proach to efficiently solve such an optimal design problem with a large number of periods for variations in energy de- mands is to derive an approximate optimal design solution by aggregating periods with a clustering method. However, such an approach does not provide any information on the accuracy for the optimal value of the objective function. The purpose of this paper is to propose a time aggregation approach for deriving suitable aprroximate optimal design solutions and evaluting their values of the objective function accurately. This time aggregation approach is realized by combining a robust optimal design method under uncertain energy demands and a hierarchical approach for solving large scale optimal design problems. A case study is conducted for a cogeneration system with a practical configuration, and it turns out that the proposed approach enables one to evaluate effective upper and lower bounds for the optimal value of the objective function as compared with those obtained by a conventional approach.}, language = {ja} } @article{YokoyamaShinanoWakayamaetal.2019, author = {Yokoyama, Ryohei and Shinano, Yuji and Wakayama, Yuki and Wakui, Tetsuya}, title = {Model reduction by time aggregation for optimal design of energy supply systems by an MILP hierarchical branch and bound method}, volume = {181}, journal = {Energy}, pages = {782 -- 792}, year = {2019}, abstract = {Mixed-integer linear programming (MILP) methods have been applied widely to optimal design of en- ergy supply systems in consideration of multi-period operation. A hierarchical MILP method has been proposed to solve such optimal design problems efficiently. An original problem has been solved by dividing it into a relaxed optimal design problem at the upper level and optimal operation problems which are independent of one another at the lower level. In addition, some strategies have been pro- posed to enhance the computation efficiency furthermore. In this paper, a method of reducing model by time aggregation is proposed as a novel strategy to search design candidates efficiently in the relaxed optimal design problem at the upper level. In addition, the previous strategies are modified in accor- dance with the novel strategy. This method is realized only by clustering periods and averaging energy demands for clustered periods, while it guarantees to derive the optimal solution. Thus, it may 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 model reduction is effective to enhance the computation efficiency in comparison and combination with the modified previous strategies.}, language = {en} } @article{KamadaYokoyamaShinanoetal.2020, author = {Kamada, Hiroki and Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {階層的最適化手法の援用による エネルギー供給システムのロバスト最適設計}, journal = {エネルギー・資源学会第36回エネルギーシステム・経済・環境コンファレンス講演論文集}, pages = {730 -- 735}, year = {2020}, abstract = {A robust optimal design method of energy supply systems under uncertain energy demands has been proposed using a mixed- integer linear model for constituent equipment. A robust optimal design problem has been formulated as a three-level min-max- min optimization one by expressing uncertain energy demands by intervals, evaluating the robustness in a performance criterion based on the minimax regret criterion, and considering hierarchical relationships among design variables, uncertain energy demands, and operation variables. Since this problem must be solved by a special algorithm and is too difficult to solve even using a commercial solver, a hierarchical optimization approach has been applied to solve the problem but its application is limited only to small scale toy problems. In this paper, some strategies are introduced into the hierarchical optimization approach to enhance the computation efficiency for the purpose of applying the approach to large scale practical problems. In a case study, the proposed approach is applied to the robust optimal design of a cogeneration system with a complex configuration, and the validity and effectiveness of the method are ascertained.}, language = {ja} } @article{TakeuchiYokoyamaShinanoetal.2020, author = {Takeuchi, Kotaro and Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {エネルギー供給システムの設計および運用の階層的関係を考慮した最適化 (k-メドイド法に基づく期間クラスタリングの適用)}, journal = {日本機械学会関西支部第95期定時総会講演会講演論文集}, pages = {1 -- 4}, year = {2020}, abstract = {To attain the highest performance of energy supply systems, it is necessary to determine design specifications optimally in consideration of operational strategies corresponding to seasonal and hourly variations in energy demands. A hierarchical mixed-integer linear programming method has been proposed to solve such an optimal design problem efficiently. In this paper, a method of reducing model by clustering periods with the k-medoids method is applied to the relaxed optimal design problem at the upper level. Through a case study, it is clarified how the proposed method is effective to enhance the computation efficiency in a large scale optimal design problem.}, language = {ja} } @article{KamadaYokoyamaShinanoetal.2020, author = {Kamada, Hiroki and Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {階層的最適化手法を用いた エネルギー供給システムのロバスト性評価}, journal = {日本機械学会関西支部第95期定時総会講演会講演論文集}, pages = {1 -- 4}, year = {2020}, abstract = {A robust optimal design method of energy supply systems under uncertain energy demands has been proposed using a mixed-integer linear model for constituent equipment. However, this method takes a long computation time, and thus it can be applied only to small-scale problems. In this paper, a hierarchical optimization method is applied to two types of optimization problems for evaluating robustness to solve them efficiently. In a case study, the proposed method is applied to a cogeneration system with a complex configuration, and the validity and effectiveness of the method are ascertained.}, language = {ja} } @article{YokoyamaShinanoWakui2019, author = {Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {Evaluation of Design Alternatives for a Gas Turbine Cogeneration Plant Based on Multiobjective K-Best Solutions}, journal = {International Conference on Power Engineering-2019}, pages = {1112 -- 1117}, year = {2019}, abstract = {To realize the best performance in energy supply, it is an important task to determine types, capacities, and numbers of energy conversion equipment appropriately. It is also necessary to take account of seasonal and hourly variations in energy demands and corresponding operational strategies of equipment. The mixed-integer linear programming (MILP) method has been utilized widely to derive the optimal solution for such a design problem. From the design viewpoint, it is important to generate not only the optimal solution but also suboptimal ones which follow the optimal one without omission, what are called K-best solutions. In this paper, an MILP method utilizing the hierarchical relationship between design and operation variables is applied along with some strategies to efficiently drive K-best solutions of a multiobjective optimal design problem of a gas turbine cogeneration system for district energy supply. The annual total cost and primary energy consumption are adopted as the objective functions to be minimized. A case study is conducted, and K-best solutions with an allowable increase in the value of the combined objective function are derived for each value of the weight for the annual total cost. Through the study, it is clarified how the weight affects the computation time, the number of the K-best solutions, and the design specifications and the values of the objective functions for the K-best solutions.}, language = {en} } @article{YokoyamaHiramatsuShinanoetal.2022, author = {Yokoyama, Ryohei and Hiramatsu, Yuto and Shinano, Yuji and Wakui, Tetsuya}, title = {Evaluation of robustness in multiple performance criteria for designing energy supply systems based on a mixed-integer linear model}, journal = {The 35th International Conference on Efficiency, Cost, Optimization, Simulation and Environmental Impact of Energy Systems (ECOS 2022)}, pages = {1937 -- 1948}, year = {2022}, abstract = {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. As a result, it turns out how the values of these performance criteria are close to the optimal ones in relation to the uncertainty in energy demands. This work will be an important step toward developing a multiobjective robust optimal design method.}, language = {en} } @article{YokoyamaShinanoTaniguchietal.2019, author = {Yokoyama, Ryohei and Shinano, Yuji and Taniguchi, Syusuke and Wakui, Tetsuya}, title = {Search for K-best solutions in optimal design of energy supply systems by an extended MILP hierarchical branch and bound method}, volume = {184}, journal = {Energy}, doi = {https://doi.org/10.1016/j.energy.2018.02.077}, pages = {45 -- 57}, year = {2019}, abstract = {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.}, language = {en} } @article{YokoyamaHiramatsuShinanoetal.2022, author = {Yokoyama, Ryohei and Hiramatsu, Yuto and Shinano, Yuji and Wakui, Tetsuya}, title = {Evaluation of robustness in multiple performance criteria for designing energy supply systems based on a mixed-integer linear model}, journal = {Proceedings of the 35th International Conference on Efficiency, Cost, Optimization, Simulation and Environmental Impact of Energy Systems (ECOS 2022)}, pages = {1937 -- 1948}, year = {2022}, abstract = {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.}, language = {en} } @article{YokoyamaShinanoTakeuchietal.2019, author = {Yokoyama, Ryohei and Shinano, Yuji and Takeuchi, Kotaro and Wakui, Tetsuya}, title = {Operation-based time-period clustering for optimal design of energy supply systems by a hierarchical MILP method}, journal = {THE 32ND INTERNATIONAL CONFERENCE ON EFFICIENCY, COST, OPTIMIZATION, SIMULATION AND ENVIRONMENTAL IMPACT OF ENERGY SYSTEMS(ECOS 2019)}, pages = {527 -- 539}, year = {2019}, abstract = {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.}, language = {en} } @article{RehfeldtKochShinano2023, author = {Rehfeldt, Daniel and Koch, Thorsten and Shinano, Yuji}, title = {Faster exact solution of sparse MaxCut and QUBO problems}, volume = {15}, journal = {Mathematical Programming Computation}, doi = {10.1007/s12532-023-00236-6}, pages = {445 -- 470}, year = {2023}, abstract = {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.}, language = {en} } @article{YokoyamaShinanoWakui2022, author = {Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {Day and time aggregations for optimal design of energy supply systems with storage units by a hierarchical MILP method}, journal = {The 35th International Conference on Efficiency, Cost, Optimization, Simulation and Environmental Impact of Energy Systems (ECOS 2022)}, pages = {1891 -- 1902}, year = {2022}, abstract = {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.}, language = {en} } @article{YokoyamaShinanoWakui2022, author = {Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {Day and time aggregations for optimal design of energy supply systems with storage units by a hierarchical MILP method}, journal = {PROCEEDINGS OF ECOS 2022 - THE 35 TH INTERNATIONAL CONFERENCE ON EFFICIENCY, COST, OPTIMIZATION, SIMULATION AND ENVIRONMENTAL IMPACT OF ENERGY SYSTEMS}, pages = {1891 -- 1902}, year = {2022}, abstract = {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.}, language = {en} } @misc{RehfeldtFujisawaKochetal.2022, author = {Rehfeldt, Daniel and Fujisawa, Katsuki and Koch, Thorsten and Nakao, Masahiro and Shinano, Yuji}, title = {Computing single-source shortest paths on graphs with over 8 trillion edges}, issn = {1438-0064}, doi = {10.12752/8818}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-88180}, year = {2022}, abstract = {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.}, language = {en} } @misc{ShinanoVigerske2025, author = {Shinano, Yuji and Vigerske, Stefan}, title = {Smoothie: Mixing the strongest MIP solvers to solve hard MIP instances on supercomputers - Phase I development}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-101906}, year = {2025}, abstract = {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.}, language = {en} } @article{ShinanoFujii2025, author = {Shinano, Yuji and Fujii, Koichi}, title = {Introduction to Quadratic Programming Solvers}, volume = {64}, journal = {Journal of The Society of Instrument and Control Engineers}, number = {8}, doi = {10.11499/sicejl.64.532}, pages = {532 -- 540}, year = {2025}, language = {en} } @article{YokoyamaShinanoWakui2023, author = {Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {An effective approach for deriving and evaluating approximate optimal design solutions of energy supply systems by time series aggregation}, volume = {11}, journal = {Frontiers in Energy Research}, doi = {10.3389/fenrg.2023.1128681}, pages = {1 -- 14}, year = {2023}, abstract = {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.}, language = {en} } @article{TateiwaShinanoYasudaetal.2023, author = {Tateiwa, Nariaki and Shinano, Yuji and Yasuda, Masaya and Kaji, Shizuo and Yamamura, Keiichiro and Fujisawa, Katsuki}, title = {Development and analysis of massive parallelization of a lattice basis reduction algorithm}, volume = {40}, journal = {Japan Journal of Industrial and Applied Mathematics}, number = {1}, doi = {10.1007/s13160-023-00580-z}, pages = {13 -- 56}, year = {2023}, abstract = {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.}, language = {en} } @misc{HosodaMaherShinano2024, author = {Hosoda, Junko and Maher, Stephen J. and Shinano, Yuji}, title = {Combining Benders' decomposition and parallelisation to solve large-scale adaptive network supply chain management problems}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-98744}, year = {2024}, abstract = {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.}, language = {en} }