This application is based upon and claims the benefit of priority of the prior Japanese Patent Application No. 2023-147780, filed on Sep. 12, 2023, the entire contents of which are incorporated herein by reference.
The embodiment discussed herein is related to a non-transitory computer-readable recording medium storing an information processing program and the like for solving a vehicle routing problem.
A vehicle routing problem (hereinafter referred to as VRP) is a problem of obtaining routes for delivering packages from a specific point called a depot by a plurality of mobile bodies, for example, transport vehicles, and returning to the depot. The VRP is one of representative problems of a combinatorial optimization problem, and is actively studied from a classical VRP to a VRP in the real world. It is not uncommon for the VRP in the real world to include problem setting specific to each delivery company, and there are various variations.
Since before, in order to cope with the various variations of the VRP, work of developing a solution system based on specific problem setting preset in the VRP is performed for each VRP.
Japanese Laid-open Patent Publication No. 2021-165196 is disclosed as related art.
According to an aspect of the embodiments, there is provided a non-transitory computer-readable recording medium storing an information processing program for causing a computer to execute processing including: generating routes in a case where a plurality of mobile bodies that have started from a specific node passes through a plurality of nodes and returns to the specific node, which satisfy a plurality of newly set conditions related to the routes, from a set of the plurality of routes that indicates order in which the plurality of mobile bodies visit the nodes, by applying a column generation method to the set of the plurality of routes; generating a route selection problem by using the generated routes; solving the generated route selection problem; and outputting a solution result of the route selection problem.
The object and advantages of the invention will be realized and attained by means of the elements and combinations particularly pointed out in the claims.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory and are not restrictive of the invention.
Although studies on solution systems for various vehicle routing problems have been conducted around the world, the solution systems that have been studied are each specialized in a particular object VRP and are difficult to apply to other variations of the VRP.
For example, in a case where a business that provides a VRP solution system to a delivery company that desires to solve a VRP is assumed, even when the solution system may be developed, a high level of human resources and a sufficient period are needed each time, which is not realistic from a viewpoint of business scale.
In one aspect, an object of an embodiment is to provide an information processing program, an information processing method, and an information processing device capable of solving even a VRP for which unknown problem setting is made.
Hereinafter, an embodiment of an information processing program, an information processing method, and an information processing device disclosed in the present application will be described in detail with reference to the drawings. Note that the present disclosure is not limited by the embodiment.
Before describing the present embodiment, representative problem setting of a vehicle routing problem (VRP) in the real world will be described. For example, a simplest definition of the VRP is indicated as a definition A1.
Definition A1: A plurality of transport vehicles that has started from a specific point called a depot returns to the depot again via customers. Order of passage of the transport vehicle is called a route, and a requested package is delivered to each customer. An optimization problem of obtaining the routes that satisfy all requests and minimize a total travel distance of the respective routes.
Even in a case where it is simply referred to as the VRP, the VRP may include a capacity restriction of a transport vehicle and an upper limit constraint on the number of transport vehicles, and an objective function (a function desired to be minimized) may be a function related to a travel time or the number of transport vehicles. In the present embodiment, problem setting such as a capacity constraint is appropriately referred to as an “attribute”. The VRP in the real world includes a wide variety of attributes. Hereinafter, examples of representative attributes will be described.
Capacity (C): Each transport vehicle needs to have a total size of loads equal to or less than a maximum load capacity.
Heterogeneous Fleet (HF): A plurality of types of transport vehicles having different maximum load capacities and the like is handled.
Backhauls (B): Packages are delivered from a depot to customers, packages are picked up from customers on a way back to the depot, and a return to the depot is made.
Time-Windows (TW): A time frame for delivery or pickup is given, and delivery or pickup needs to be performed within the time frame.
Multiple Depots (MD): There are multiple depots, and a transport vehicle may start from any one of the multiple depots.
Multi-Trips (MT): A temporary return to the depot and a re-start with packages may be made.
Split Deliveries (SD): Split delivery may be performed for one customer.
Pickup and Deliveries (PD): A package is picked up at a visit to a customer and is delivered to another customer.
Time Dependent (TD): In order to consider road congestion, a travel time between points depends on a time zone of the day.
In the following description, for example, the attribute “Time-Windows” may be abbreviated as the attribute “TW”. In addition to a travel distance and a travel time of a transport vehicle, a fixed cost that varies depending on a type of the transport vehicle may be considered as the frequent objective function.
Subsequently, a VRP general form will be described. In an information processing device according to the present embodiment, an attribute of a VRP to be solved and basic information such as a travel time between the respective nodes are given as input. The node corresponds to a point such as a depot or a customer. The VRP solved by the information processing device is internally called a VRP general form. In order to define the VRP general form, each of a “task”, a “route constraint”, a “transport vehicle number constraint”, a “route cost”, and a “travel time calculation method” is defined. In the following description, the information processing device according to the present embodiment is referred to as an “information processing device 100”.
A definition of the “task” will be described. The task is represented as a set of information on work of picking up and delivering a package. For example, the task includes the following information. *Pickup node: a node where a package is picked up, *Delivery node: a node where a package is delivered, *Priority: a priority related to the task, *Load amount vector: a vector related to a load amount, indicating a type, an amount, or the like of a package, *Pickup time frame: a time frame related to pickup of a package, *Delivery time frame: a time frame related to delivery of a package, and *Possibility of split: whether split delivery is possible.
In the simplest VRP, all packages are picked up at a depot and delivered to customers. For example, in a case where the attribute PD is considered, a package is delivered from a customer to a different customer. In the case of the attribute B, there are two types of packages: a package delivered from a depot to a customer, and a package picked up from a customer on the way before a return to the depot. In a VRP in practical work, there are cases where a plurality of attributes is mixed, and one of objects of introducing the concept of the task is to express them flexibly as optimization problems in the same framework.
Next, a definition of the “route constraint” will be described. The route constraint is a constraint condition for each route, which does not depend on another route. Note that the route commonly refers to visit order of a transport vehicle, and here, it is assumed that the route also includes constraint conditions related to information such as start and arrival times and a load amount accompanying the route. The route constraint in practical work includes, for example, a constraint condition related to a time at which a return to a depot needs to be made, a constraint condition related to a travel time different for each type of a transport vehicle, a constraint condition related to a node that needs to be visited before a return to a depot, and the like.
Next, a definition of the “transport vehicle number constraint” will be described. For example, when a positive integer k and a condition related to a transport vehicle are given, the transport vehicle number constraint is a constraint condition that the number of transport vehicles satisfying the condition is equal to or less than k or equal to k. For example, the following constraint conditions correspond to the transport vehicle number constraint. *Constraint condition related to the number of usable specific transport vehicles, *Constraint condition related to an upper limit of the total number of transport vehicles, and *Constraint condition related to the number of transport vehicles capable of long-distance travel.
Next, a definition of the “route cost” will be described. A cost that may be defined for one route is called the route cost. For example, the route cost includes a fixed cost depending on a travel time and a travel distance of a route, a type of a transport vehicle, and the like.
Next, a definition of the “travel time calculation method” will be described. The travel time calculation method is a calculation method related to a travel time between optional two nodes (edge). In the simplest VRP, a travel time of each edge is uniquely determined. In a case where road congestion is considered, the travel time calculation method according to a start time of a start point node of an edge is determined.
Next, a definition of the VRP general form will be described. For example, the definition of the VRP general form is indicated as a definition B1.
Definition B1: The following optimization problem is called the VRP general form. *A plurality of transport vehicles that has started from a specific node called a depot passes through another node and returns to the depot again (depot start-arrival constraint), *One or more types of transport vehicles having different sizes, start nodes (depots), and the like are available, *Packages of all tasks are picked up and delivered, and *Each task is assigned to a transport vehicle, and the following constraint conditions (task constraints) b-1, b-2, b-3, and b-4 related to tasks are satisfied.
Note that the representative attributes described with the simplest VRP may be expressed as follows in the VRP general form.
The information processing device 100 solves the VRP general form. Since the VRP general form may increase attributes corresponding to the route constraint, the transport vehicle number constraint, and the route cost, there is a possibility that the VRP general form may be applied to various variations of the VRP.
Subsequently, a route generation approach of the information processing device 100 will be described. The information processing device 100 adopts a column generation method to generate routes satisfying the depot start-arrival constraint, the task constraint, and all the route constraints based on the VRP general form. Furthermore, the information processing device 100 selects, from the generated routes, routes that minimize an objective function so as to satisfy other constraint conditions.
Before describing the column generation method, a route selection problem used in the present embodiment is formulated. First, the following are defined.
In the present embodiment, an integer variable xr indicating the number of adoptions of a route (the number of vehicles using the corresponding route) is introduced for each route r∈R. At this time, the route selection problem is formulated by Expression (5). Note that it is assumed that the xr included in Expression (5) satisfies Expressions (5a), (5b), (5c), (5d), and (5e).
When solving the route selection problem, the information processing device 100 searches for a value of the xr that minimizes a value of Expression (5) under the condition that the xr satisfies Expressions (5a) to (5e). Note that, in a case where the value of the xr is 0, the corresponding route is excluded from routes to be adopted. For example, in a case where a value of xr1 is 0, it indicates that a route r1 is not adopted.
Subsequently, the column generation method for the VRP general form will be described. In a case where it is assumed that all the routes satisfying the constraint conditions may be listed, an optimal solution for the VRP general form may be obtained by solving the route selection problem. However, listing all the routes is not realistic from a viewpoint that a memory of the information processing device 100 and the route selection problem become large in scale.
In practice, since an inefficient route is also included in all the routes, an object of the column generation method is to generate only useful routes. A goal of the column generation method is to obtain an optimal value of a linear relaxation problem of the route selection problem in a case where all the routes are listed. It is empirically known that a set of routes generated for obtaining an optimal solution of a linear programming problem is excellent in configuring an optimal solution of an original route selection problem.
Here, an example of a processing procedure of the column generation method executed by the information processing device 100 will be described.
The information processing device 100 solves the reduced cost minimization problem (step S12). In a case where an objective function value of an obtained solution is negative (Step S13, Yes), the information processing device 100 adds a route corresponding to the obtained solution to the route set R (step S14), and proceeds to step S11.
On the other hand, in a case where the objective function value of the obtained solution is not negative (Step S13, No), the information processing device 100 outputs the route set R (step S15).
Since the linear relaxation problem in step S11 of
Therefore, the information processing device 100 applies a local search to the reduced cost minimization problem to obtain the solution. For example, the information processing device 100 executes processing of
The information processing device 100 sets a route set R′ in which one task is assigned to the route r (step S21). The information processing device 100 selects, from among routes included in the route set R′, a first route being a route which satisfies the constraint condition C and minimizes a value of an objective function g( ) of the reduced cost minimization problem, and sets the selected first route as a route r′ (i.e., r′∈R′) (step S22).
The information processing device 100 selects, from among the nodes, a node v other than a depot (step S23). The information processing device 100 sets a route rv obtained by adding the node v to a visit node of the route r (step S24). The information processing device 100 sets a route set R′v in which one task is assigned to the route rv (step S25).
The information processing device 100 selects, from among routes included in the route set R′v, a second route being a route which satisfies the constraint condition C and minimizes a value of the objective function g( ) of the reduced cost minimization problem, and sets the selected second route as a route r′v (i.e., r′v∈R′) (step S26).
In a case where g(r′)>g(r′v) is satisfied, the information processing device 100 sets r′:=r′v (step S27). For example, in the step S27, the route r′v is set into a value of a variable indicating the route r′. In a case where all the nodes other than the depot have not been selected (step S28, No), the information processing device 100 proceeds to step S23. On the other hand, in a case where all the nodes other than the depot have been selected (step S28, Yes), the information processing device 100 proceeds to step S29.
In a case where g(r)<g(r′) is not satisfied (step S29, No), the information processing device 100 sets r:=r′ (step S30), and proceeds to step S21. On the other hand, in a case where g(r)<g(r′) is satisfied (step S29, Yes), the information processing device 100 outputs the route r (step S31).
Here, the information processing device 100 includes a route constraint determination unit, a transport vehicle number condition determination unit, a route cost calculation unit, and a travel time calculation unit, and copes with a VRP including an unexpected attribute using the route constraint determination unit, the transport vehicle number condition determination unit, the route cost calculation unit, and the travel time calculation unit.
The route constraint determination unit checks a constraint condition violation of the route constraint. For example, the route constraint determination unit is used in the processing described in step S22 of
The transport vehicle number condition determination unit performs condition determination of the transport vehicle number constraint. For example, the transport vehicle number condition determination unit is used in the processing described in steps S22, S26, S27, and S29 of
The route cost calculation unit calculates a cost of a route. For example, the route cost calculation unit is used in the processing described in steps S22, S26, and S27 of
The travel time calculation unit calculates a travel time of an optional edge (optional two nodes) from information regarding order of visit nodes, a task in charge, and a type of a transport vehicle. For example, in a case where a travel time between two nodes is calculated by the travel time calculation unit, a travel time when travel is performed among all the nodes (in the following description, referred to as a travel time for a route) is determined for each route. The travel time for a route may be used for a check of a constraint condition violation of the route constraint, condition determination of the transport vehicle number constraint, and calculation of the route cost.
A relationship between the route constraint and the travel time for a route will be described. For example, in an attribute “Route Duration Constraints”, an upper limit is imposed on the travel time for a route. Therefore, the travel time for a route is used in a case where the route constraint determination unit determines whether or not the route constraint is violated. For example, in a case where a vehicle may not travel within the travel time, it may be determined that the route constraint is violated.
A relationship between the transport vehicle number constraint and the travel time for a route will be described. For example, the travel time for a route is used in a case where the transport vehicle number condition determination unit performs condition determination of the transport vehicle number constraint such as “the number of transport vehicles whose travel time for a route is allowed to be a predetermined time (8 hours) or more is equal to or less than three”.
A relationship between the route cost and the travel time for a route will be described. The route cost calculation unit uses the travel time for a route in a case where the route cost depends on the travel time for a route. For example, when the route cost is the travel time for a route, the route cost depends on the travel time for a route. When a driver of a transport vehicle is paid by the hour, the route cost depends on the travel time for a route. For example, when the travel time for a route is 1.5 hours, a labor cost of the driver for two hours is the route cost.
Subsequently, an example of a processing procedure of the information processing device 100 according to the present embodiment will be described.
The information processing device 100 generates a VRP general form based on the input VRP data (step S102). The information processing device 100 generates routes by applying the column generation method to the VRP general form (step S103). For example, the processing procedure in step S103 of
The information processing device 100 generates a route selection problem using the generated routes (corresponding to the route set R in step S15 of
The information processing device 100 solves the generated route selection problem (step S105). The information processing device 100 converts a solution result of the route selection problem into output data (step S106). The information processing device 100 outputs the output data (step S107).
Next, a configuration example of the information processing device 100 that executes the processing described with reference to
The communication unit 110 executes data communication with an external device or the like owned by a delivery company via a network. The communication unit 110 is a network interface card (NIC) or the like. For example, the communication unit 110 may receive VRP data 141 from the external device or the like.
The input unit 120 is an input device that inputs various types of information to the control unit 150 of the information processing device 100. For example, the input unit 120 corresponds to a keyboard, a mouse, a touch panel, or the like. A user may input a solution parameter 143 or the like by operating the input unit 120.
The display unit 130 is a display device that displays information output from the control unit 150. For example, the display unit 130 displays output data of a solution result of a route selection problem.
The storage unit 140 includes the VRP data 141, a VRP general form database (DB) 142, and the solution parameter 143. The storage unit 140 is a memory or the like.
The VRP data 141 includes various types of data used in a VRP. For example, the VRP data 141 has data of each task, a travel time between nodes, the number of transport vehicles, and a route cost coefficient. Furthermore, the VRP data 141 includes a travel distance between nodes, a maximum load capacity, a problem setting flag, and the like. The VRP data 141 may include an add-in program for executing a route constraint determination unit 153, a transport vehicle number condition determination unit 154, a route cost calculation unit 155, and a travel time calculation unit 156 described later. For example, the VRP data 141 includes conditions related to a specific route of a delivery company. The conditions related to the specific route of the delivery company include a condition for determining whether or not the corresponding route satisfies the route constraint, a condition for determining whether or not a transport vehicle number condition is satisfied, a definition of the route cost, and data for calculating a travel time.
The VRP general form DB 142 holds data of a VRP general form (optimization problem) generated based on the VRP data 141. The VRP general form is defined in the definition B1 as described above.
The solution parameter 143 includes a parameter such as a time limit. The information processing device 100 ends solution in a case where a time needed for the solution exceeds the time limit set in the solution parameter 143.
Subsequently, description proceeds to description of the control unit 150. The control unit 150 includes an acquisition unit 151, a generation unit 152, the route constraint determination unit 153, the transport vehicle number condition determination unit 154, the route cost calculation unit 155, the travel time calculation unit 156, and a solution unit 157. The control unit 150 is a central processing unit (CPU), a graphics processing unit (GPU), or the like.
The acquisition unit 151 acquires the VRP data 141 from an external device or the like via the communication unit 110, and stores the acquired VRP data 141 in the storage unit 140. For example, the acquisition unit 151 receives a digital signal that indicates the VRP data 141 from the external device via the communication unit 110, and acquires the VRP data from the received digital signal by decoding the received digital signal. The acquisition unit 151 acquires the solution parameter 143 from the input unit 120 or the like, and stores the solution parameter in the storage unit 140.
The generation unit 152 generates a VRP general form based on the VRP data 141. For example, the generation unit 152 changes initial values of the VRP general form DB 142 based on the VRP data 141.
In the example illustrated in
For example, the generation unit 152 updates the initial values of the VRP general form DB 142 as illustrated in
For example, it is assumed that the generation unit 152 converts the initial values of the VRP general form DB 142 based on a preset conversion table (a rule for converting a value for a name of the VRP data 141 into a value corresponding to a name of the VRP general form DB 142). Here, the case has been described where the generation unit 152 generates the VRP general form DB 142 based on the VRP data 141. However, the generation unit 152 may acquire data of the VRP general form having set values from the input unit 120 or the like, and update the initial values of the VRP general form DB 142 based on the acquired data of the VRP general form.
The route constraint determination unit 153 checks a constraint condition violation of the route constraint. When the solution unit 157 solves a reduced cost minimization problem by a local search, the route constraint determination unit 153 acquires the object router from the solution unit 157 and determines whether or not the route r satisfies the route constraint set in the VRP general form DB 142. For example, the route constraint determination unit 153 is used in the processing described in step S22 of
The transport vehicle number condition determination unit 154 performs condition determination of the transport vehicle number constraint. When the solution unit 157 calculates a value of an objective function g( ) of a reduced cost minimization problem, the transport vehicle number condition determination unit 154 determines whether or not the value set to the objective function g( ) satisfies the transport vehicle number constraint set in the VRP general form DB 142. For example, the transport vehicle number condition determination unit 154 is used in the processing described in steps S22, S26, S27, and S29 of
The route cost calculation unit 155 calculates a cost of a route. When the solution unit 157 calculates a value of an objective function g( ) of a reduced cost minimization problem, the route cost calculation unit 155 calculates the route cost based on the route cost list, the route cost coefficient list, and the like set in the VRP general form DB 142. For example, the route cost calculation unit 155 is used in the processing described in steps S22, S26, and S27 of
The travel time calculation unit 156 calculates a travel time of an optional edge (optional two nodes) from information regarding order of visit nodes, a task in charge, and a type of a transport vehicle, which is set in the VRP general form DB 142. For example, in a case where a travel time between two nodes is calculated by the travel time calculation unit 156, a travel time when travel is performed among all nodes (a travel time for a route) is determined for each route. The travel time for the route calculated by the travel time calculation unit 156 is used for a check of a constraint condition violation of the route constraint by the route constraint determination unit 153, condition determination of the transport vehicle number constraint by the transport vehicle number condition determination unit 154, and calculation of a route cost by the route cost calculation unit 155. The travel time calculation unit 156 outputs data of the travel time corresponding to the route to the route constraint determination unit 153, the transport vehicle number condition determination unit 154, and the route cost calculation unit 155.
The solution unit 157 specifies routes by executing, at a processor circuit such as a central processing unit, the following processing to solve a VRP general form. First, the solution unit 157 generates routes by applying the column generation method to data of the VRP general form DB 142. For example, the solution unit 157 generates routes (route set R) by executing the processing described with reference to
The solution unit 157 executes the local search described with reference to
The solution unit 157 generates routes by the column generation method, and then generates a route selection problem based on the generated routes. For example, the solution unit 157 formulates the route selection problem by Expression (5).
The solution unit 157 solves the generated route selection problem. For example, in the case of solving the route selection problem, the solution unit 157 searches for a value of the integer variable xr that minimizes a value of Expression (5) under the condition that the xr satisfies Expressions (5a) to (5e). The solution unit 157 generates output data by setting, in a route data list, a route in which the value of the xr does not become 0 among the plurality of routes generated by the column generation method. For example, in a case where the value of the xr1 is 0, the route r1 is not added to the route data list. Note that data of the routes included in the route data list includes a type of an assigned transport vehicle, node visit order, a node visit time, task load amount data, and the like.
The solution unit 157 outputs the route data list (output data) to the display unit 130. The solution unit 157 may notify an external device or the like of the route data list via the communication unit 110. Note that, in the case of generating the route selection problem and in the case of solving the route selection problem, the solution unit 157 ends solution in a case where the time limit set in the solution parameter 143 is exceeded.
Next, effects of the information processing device 100 according to the present embodiment will be described. The information processing device 100 applies the column generation method to the VRP general form DB 142, and generates, from a set of a plurality of routes, routes that satisfy newly set conditions of the routes and minimize a value of an objective function. The information processing device 100 generates a route selection problem using the routes generated by applying the column generation method, solves the route selection problem, and outputs a solution result. As a result, it is possible to solve even a VRP for which unknown problem setting is made.
In a case where the routes are generated by applying the column generation method, the information processing device 100 generates, by using the route constraint determination unit 153, the routes that satisfy a route constraint and minimize the value of the objective function. As a result, even when a new route constraint is set in the VRP data 141, the routes satisfying the new route constraint may be generated.
In a case where the routes are generated by applying the column generation method, the information processing device 100 generates, by using the transport vehicle number condition determination unit 154, the routes that minimize the value of the objective function within a range satisfying a transport vehicle number constraint. As a result, even when a new transport vehicle number condition is set in the VRP data 141, the routes satisfying the new transport vehicle number condition may be generated.
In a case where the routes are generated by applying the column generation method, the information processing device 100 generates, by using the route cost calculation unit 155 and a route cost calculated by the route cost calculation unit 155, the routes that minimize the value of the objective function. As a result, even when a new definition of the route cost is set in the VRP data 141, optimal routes may be generated.
In a case where the routes are generated by applying the column generation method, the information processing device 100 performs, by using the travel time calculation unit 156 and a travel time calculated by the travel time calculation unit 156, determination as to whether or not the route constraint is satisfied, determination as to whether or not the transport vehicle number constraint is satisfied, and calculation of the route cost. As a result, even when a new definition of the travel time is set in the VRP data 141, optimal routes may be generated.
Next, an example of a hardware configuration of a computer that implements functions similar to those of the information processing device 100 described above will be described.
As illustrated in
The hard disk device 307 includes an acquisition program 307a, a generation program 307b, a determination and calculation program 307c, and a solution program 307d. The CPU 301 reads each of the programs 307a to 307d and develops the programs 307a to 307d in the RAM 306.
The acquisition program 307a functions as an acquisition process 306a. The generation program 307b functions as a generation process 306b. The determination and calculation program 307c functions as a determination and calculation process 306c. The solution program 307d functions as a solution process 306d.
Processing of the acquisition process 306a corresponds to the processing of the acquisition unit 151. Processing of the generation process 306b corresponds to the processing of the generation unit 152. Processing of the determination and calculation process 306c corresponds to the processing of the route constraint determination unit 153, the transport vehicle number condition determination unit 154, and the route cost calculation unit 155. Processing of the solution process 306d corresponds to the processing of the solution unit 157.
Note that each of the programs 307a to 307d does not necessarily have to be stored in the hard disk device 307 beforehand. For example, each of the programs is stored in a “portable physical medium” to be inserted in the computer 300, such as a flexible disk (FD), a compact disc read only memory (CD-ROM), a digital versatile disc (DVD), a magneto-optical disk, or an integrated circuit (IC) card. Then, the computer 300 may read and execute each of the programs 307a to 307d.
All examples and conditional language provided herein are intended for the pedagogical purposes of aiding the reader in understanding the invention and the concepts contributed by the inventor to further the art, and are not to be construed as limitations to such specifically recited examples and conditions, nor does the organization of such examples in the specification relate to a showing of the superiority and inferiority of the invention. Although one or more embodiments of the present invention have been described in detail, it should be understood that the various changes, substitutions, and alterations could be made hereto without departing from the spirit and scope of the invention.
| Number | Date | Country | Kind |
|---|---|---|---|
| 2023-147780 | Sep 2023 | JP | national |