In the present paper, we focus on the bright-dark solitons and interaction behavior associated with a dimensionally reduced p -gBKP equation. New test functions are constructed by using the Bilinear Neural Network Method (BNNM). These test functions are composed of arbitrary functions in ‘3-2-2-2’ neural network model or specific activation functions in single-layer model. The bright-dark solitons and interaction solutions between lump and double exponential function or hyperbolic tangent function are obtained by choosing appropriate values and activation functions. Via various three-dimensional plots, 2D plots, contour plots and density plots, the evolution characteristics of these waves are exhibited well.