Shifted bilinear sums of Salié sums and the distribution of modular square roots of shifted primes

We establish various upper bounds on Type-I and Type-II shifted bilinear sums with Salié sums modulo a large prime $q$. We use these bounds to study, for fixed integers $a,b\equiv ̸0q$, the distribution ofsolutions to the congruence $x^2\equiv ap+bq$, over primes $p\le P$. This is similar to the recently studied case of $b=0$, however the case $b\equiv ̸0q$ exhibits some new difficulties.