Aerodynamic shape optimization of two-dimensional airfoils in inviscid compressible flows is performed using a continuous adjoint formulation on unstructured meshes. Accurate evaluation of the gradient is achieved by using a reconstruction scheme based on the Laplacian averaging. A least-square method with extended stencil is used for flow gradient calculations. Proper convergence criterion is studied on Euler and adjoint equations for efficient design. The present method has been applied to RAE2822 and NACA0012 airfoils such that wave drag can be minimized by removing the shock wave. An inverse design is also performed to recover the shock wave on the designed RAE2822 airfoil.