原题链接：
题意：
(已知(x+y)) % (p = b) , ((x*y)) % (p = c)(space)求x , y.若不存在则输出(-1space-1)
思路：
(通过平方差公式(x-y)^2) = ((x+y)^2) (-) (4*(x+y)) = ((b^2-4c)) % (p),
(我们就会想到二次剩余存在X^2equiv d(modspace p ) ,(d是模p的二次剩余) 来解得x-y，然后用b与(x-y)求出x,y)
(这里使用欧拉准则计算勒让德符号:)
(勒让德符号可以计算当p为素数的情况的二次剩余判别问题,而雅可比符号判定条件为a为任意整数，p为任意奇数)
(题目种给出的p为1000000007，是一个质数所以我们就用欧拉准则来做。)
(根据欧拉准则：当 a 是 modspace p 的平方剩余时有 a^frac{p-1}{2} equiv space1) ((modspace p)space，当a不是平方剩余时a^frac{p-1}{2} equiv space-1(modspace p ))
(所以我们就用快速幂来求出(b^2-4c)^frac{p-1}{2} 判断是否等于space mod-1（无解情况），以及求出最后x,y的解。)
(第一次使用markdown编辑不容易啊...)
```ruby
#include
#define ll long long
using namespace std;
const ll mod = 1000000007;
ll q_pow(ll a,ll b){
ll res=1;
while(b){
if(b&1)res=res*a % mod;
a=a*a%mod;
b>>=1;
}
return res;
}
int main(){
int T;
cin >> T;
while(T--){
ll b,c;
cin >> b >> c;
ll delta = (b * b - 4 * c + mod) % mod;
if(q_pow (delta,(mod - 1) / 2) == mod - 1) { puts("-1 -1"); continue; }
ll a = q_pow (delta,((mod + 1) / 2 + 1) / 2);
ll x = ((b - a + mod) * q_pow(2, mod - 2)) % mod;
ll y = (b + mod - x) % mod;
printf("%lld %lld
", min(x,y),max(x,y));
}
return 0;
}
```