二次剩余模板

。。历时两周终于把它看懂了，太不容易了qaq（当然，是断断续续看的）

这是自己认为讲的非常详细的两篇博客了
here1 and here2

自用模板：

ll q_pow(ll a,ll b,ll p)
{
	ll ans=1;
	while(b)
	{
		if(b&1) ans=ans*a%p;
		a=a*a%p;
		b>>=1;
	}
	return ans;
}
ll Quadratic_residue(ll n,ll mod)
{
	if(n==0)	
		return 0;
	if(mod==2) return 1;
	if(q_pow(n,(mod-1)>>1,mod)==mod-1) return -1;
	ll b=(rand()<<14^rand())%mod;
	while(q_pow(b,(mod-1)>>1,mod)!=mod-1) 
		b=(rand()<<14^rand())%mod;
	ll s=mod-1,t=0,x,inv=q_pow(n,mod-2,mod),f=1;
	while(s%2==0) s>>=1,t++,f<<=1;
	t--,x=q_pow(n,(s+1)>>1,mod),f>>=1;
	while(t)
	{
		f>>=1;
		if(q_pow(inv*x%mod*x%mod,f,mod)!=1) x=x*q_pow(b,s,mod)%mod;
		t--,s<<=1;
	}
	return min(x,mod-x);
}