We develop a mathematical model considering behavioral changes and underreporting to describe the first major COVID-19 wave in Metro Manila, Philippines. Key parameters are fitted to the cumulative cases in the capital from March to September 2020. A bi-objective optimization problem is formulated that allows for the easing of restrictions at an earlier time and minimizes the number of additional beds ensuring sufficient capacity in healthcare facilities. The well-posedness of the model and stability of the disease-free equilibria are established. Simulations show that if the behavior was changed one to four weeks earlier before the easing of restrictions, cumulative cases can be reduced by up to 55% and the peak delayed by up to four weeks. If reporting is increased threefold in the first three months of the estimation period, cumulative cases can be reduced by 61% by September 2020. Among the Pareto optimal solutions, the peak of cases is lowest if strict restrictions were eased on May 20, 2020 and with at least 56 additional beds per day.

Without vaccines and medicine, non-pharmaceutical interventions (NPIs) such as social distancing, have been the main strategy in controlling the spread of COVID-19. Strict social distancing policies may lead to heavy economic losses, while relaxed social distancing policies can threaten public health systems. We formulate optimization problems that minimize the stringency of NPIs during the prevaccination and vaccination phases and guarantee that cases requiring hospitalization will not exceed the number of available hospital beds. The approach utilizes an SEIQR model that separates mild from severe cases and includes a parameter $\mu$ that quantifies NPIs. Payoff constraints ensure that daily cases are decreasing at the end of the prevaccination phase and cases are minimal at the end of the vaccination phase. Using a penalty method, the constrained minimization is transformed into a non-convex, multi-modal unconstrained optimization problem. We solve this problem using the improved multi-operator differential evolution, which fared well when compared with other optimization algorithms. We apply the framework to determine optimal social distancing strategies in the Republic of Korea given different amounts and types of antiviral drugs. The model considers variants, booster shots, and waning of immunity. The optimal $\mu$ values show that fast administration of vaccines is as important as using highly effective vaccines. The initial number of infections and daily imported cases should be kept minimum especially if the bed capacity is low. In Korea, a gradual easing of NPIs without exceeding the bed capacity is possible if there are at least seven million antiviral drugs and the effectiveness of the drug in reducing severity is at least 86%. Model parameters can be adapted to a specific region or country, or other infectious diseases. The framework can be used as a decision support tool in planning economic policies, especially in countries with limited healthcare resources.