【NEU1440】The minimum square sum (平方剩余和欧拉准则)

最新推荐文章于 2019-09-20 18:26:03 发布

「已注销」

最新推荐文章于 2019-09-20 18:26:03 发布

阅读量363

点赞数

分类专栏：数学定理

本文链接：https://blog.csdn.net/loveyan0213/article/details/45648789

版权

数学定理专栏收录该内容

2 篇文章 0 订阅

订阅专栏

1440: The minimum square sum

时间限制: 1 Sec 内存限制: 256 MB
提交: 101 解决: 28
[ 提交][ 状态][ 讨论版]

题目描述

　　Given a prime p (p<10^8),you are to find min{x^2+y^2},where x and y belongs to positive integer, so that x^2+y^2=0 (mod p).

输入

　　Every line is a p. No more than 10001 test cases.

输出

　　The minimum square sum as described above.

样例输入

2
3
5

样例输出

2
18
5

若p=2或p=4*k+1 则p可以表成两平方数的和的形式 (欧拉和费马已证明，并且有求的方法) 所以答案是p
若p=4*k+3 设a^2=n(mod p) (n!=0)  可以证明不存在b,b^2=p-n(mod p) 即若n是p的平方剩余 则p-n不是p的平方剩余
证明:因为a^2=n(mod p) 所以由欧拉准则 得n^((p-1)/2)=1(mod p)
若b^2=-n(mod p) 那么(-n)^((p-1)/2)=1(mod p)
左边把符号提出来 得(-1)^((p-1)/2)*n^((p-1)/2)
因为p是4*k+3型的 所以(p-1)/2是奇数 所以左边不可能等于1 假设与事实矛盾所以
因此a^2=0(mod p) b^2=0(mod p) 所以答案是2*p*p
代码如下：
#include<cstdio>
using namespace std;

int main()
{
    long long p;
    while(scanf("%lld",&p))
    {
        if(p==2 || p%4==1)
            printf("%lld\n",p);
        else
            printf("%lld\n",p*p*2);
    }
    return 0;
}