In actual flight, aircraft rarely operate under a single design condition; multiple flight states must be considered to meet performance requirements. With the push for green and low-carbon aviation, there is growing demand for high-performance, fuel-efficient aircraft. This study focuses on the Blended Wing Body (BWB) configuration. To address large-scale design variables and multiple constraints, a discrete adjoint-based aerodynamic optimization method is developed, improving computational efficiency and reducing cost.The optimization results show reduced drag coefficients across various flight conditions and enhanced drag divergence performance. The robustness of the multi-point optimization approach is validated, confirming its ability to improve aircraft performance across different states. The proposed method is practical and provides an effective reference for aerodynamic design of BWB aircraft.