POJ_1808_Quadratic Residues

最新推荐文章于 2019-10-26 13:08:30 发布

edward_108

最新推荐文章于 2019-10-26 13:08:30 发布

阅读量170

点赞数

分类专栏：随笔杂记

本文链接：https://blog.csdn.net/edward_108/article/details/75452745

版权

随笔杂记专栏收录该内容

37 篇文章 0 订阅

订阅专栏

Quadratic Residues

Time Limit: 1000MS		Memory Limit: 30000K
Total Submissions: 1614		Accepted: 542

Description

Background
In 1801, Carl Friedrich Gauss (1777-1855) published his "Disquisitiones Arithmeticae", which basically created modern number theory and is still being sold today. One of the many topics treated in his book was the problem of quadratic residues.
Consider a prime number p and an integer a !≡ 0 (mod p). Then a is called a quadratic residue mod p if there is an integer x such that
x² ≡ a (mod p),
and a quadratic non residue otherwise. Lagrange (1752-1833) introduced the following notation, called the "Legendre symbol":

For the calculation of these symbol there are the following rules, valid only for distinct odd prime numbers p, q and integers a, b not divisible by p:

The statements 1. to 3. are obvious from the definition, 4. is called the Completion Theorem, and 5. is the famous Law of Quadratic Reciprocity for which Gauss himself gave no less than six different proofs in the "Disquisitiones Arithmeticae". Knowing these facts, one can calculate all possible Legendre symbols as in the following example:

Input

The first line contains the number of scenarios.
For each scenario, there is one line containing the integers a and p separated by a single blank, where 2 < p < 10⁹ is an odd prime, and a satisfies both a !≡ 0 (mod p) and |a| <= 10⁹.

Output

Start the output for every scenario with a line containing "Scenario #i:", where i is the number of the scenario starting at 1. Then print a single line containing (a/p), followed by a blank line.

Sample Input

Sample Output

Scenario #1:
-1

Scenario #2:
-1

Scenario #3:
1

Source

TUD Programming Contest 2003, Darmstadt, Germany

判定平方剩余是否有解

#include <iostream>
#include <string>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <climits>
#include <cmath>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>
using namespace std;
typedef long long           LL ;
typedef unsigned long long ULL ;
const int    maxn = 1000 + 10  ;
const int    inf  = 0x3f3f3f3f ;
const int    npos = -1         ;
const int    mod  = 1e9 + 7    ;
const int    mxx  = 100 + 5    ;
const double eps  = 1e-6       ;

LL pow_mod(LL a, LL k, LL m){
	LL r=1%m;
	while(k){
		if(1&k)r=(r*a)%m;
		a=(a*a)%m;
		k>>=1;
	}
	return r%m;
}
LL mod_sqr(LL a, LL n){
	if(1==pow_mod(a,(n-1)/2,n))
		return 1;
	return -1;
}
int T, kase, a, p;
int main(){
	// freopen("in.txt","r",stdin);
	// freopen("out.txt","w",stdout);
	while(~scanf("%d",&T)){
		for(kase=1;kase<=T;kase++){
			scanf("%lld %lld",&a,&p);
			printf("Scenario #%d:\n%d\n\n",
				kase,mod_sqr((a%p+p)%p,p));
		}
	}
	return 0;
}