【TOJ1132】Square Root，二次同余方程

传送门
写在前面：ACdreamer大神太强辣！希望能达到他一半的高度啊！
思路：ACdreamer’s blog
裸题调了好久ORZ
注意：关于(a+√w）^((p+1)/2)的求法，我至今都没太想明白……
代码：

#include"bits/stdc++.h"
using namespace std;
int t,a,x,p,w;
struct os
{int u,v;
};
os mul(os a,os b)
{os c;c.u=(a.u*b.u%p+a.v*b.v%p*w%p)%p;c.v=(a.u*b.v%p+a.v*b.u%p)%p;return c;
}
int qr(int x,int y,int z)
{int ans=1;while (y){if (y&1) ans=ans*x%p;y>>=1;x=x*x%z;}return ans;
}
main()
{scanf("%d",&t);while (t--){scanf("%d%d",&x,&p);if (p==2) {printf("1\n");continue;}if (qr(x,(p-1)>>1,p)-p==-1){printf("No root\n");continue;}while (1){a=rand()%p;w=a*a-x;w%=p;while (w<0) w+=p;if (qr(w,(p-1)>>1,p)-p==-1) break;}int yy=(p+1)>>1;os ans,xx;ans.u=1;ans.v=0;xx.u=a;xx.v=1;while(yy){if (yy&1) ans=mul(ans,xx);xx=mul(xx,xx);yy>>=1;}int ans1=ans.u,ans2=p-ans.u;if (ans1==ans2) printf("%d\n",ans1);else printf("%d %d\n",min(ans1,ans2),max(ans1,ans2));}
}

【TOJ1132】Square Root，二次同余方程相关推荐

快速开平方根倒数算法(Fast inverse square root)的一点探究
文章目录一.写在前面 1. 提示 2. 背景与前情二.正文 1. 需求分析 2. 必备工具之IEEE-754浮点数表示方法 3. 同一储存单元32bits的两种不同意义 4. 公式推导 4. 本文 ...
二次剩余，二次同余方程
定义: 对于方程x2≡n(mod  p)  ,  nx^2\equiv n (mod \;p)\;,\;nx2≡n(modp),n为p的二次剩余 , x为该二次同余方程的解就如字面意思一样 , n就 ...
UVA10023 Square root【大数】
You are to determinate X by given Y , from expression X =Y\sqrt{Y}Y Input The first line is the num ...
平方根倒数速算法（Fast Inverse Square Root）
平方根倒数速算法(Fast Inverse Square Root) 一.背景平方根倒数速算法是适用于快速计算积的算术平方根的倒数(在此需取符合IEEE 754标准格式的32位浮点数)的一种算法 ...
The Square Root of Three 孤独的根号三 David B. Feinberg 大卫·范伯格
The Square Root of Three 孤独的根号三 David B. Feinberg 大卫·范伯格 I'm sure that I will always be A lonely num ...
[Swift]快速反向平方根 | Fast inverse square root
★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★ ➤微信公众号:山青咏芝(shanqingyongzhi) ➤博客园地址:山青咏芝(https://www.cnblog ...
牛顿法求解1-100的平方根python_使用牛顿-拉弗森法定义平方根函数（Newton-Raphson method Square Root Python）...
牛顿法(Newton's method)又称为牛顿-拉弗森法(Newton-Raphson method),是一种近似求解实数方程式的方法.(注:Joseph Raphson在1690年出版的< ...
欧拉工程第57题：Square root convergents
题目链接 Java程序 package projecteuler51to60;import java.math.BigInteger; import java.util.Iterator; impor ...
CodeForces - 715A Plus and Square Root(思维+构造)
题目链接:点击查看题目大意:我们在玩一个游戏,屏幕上有一个数字,我们设这个数字为x,初始值为2,我们一开始的等级是k,我们每一次可以有两种操作: 加法:可以让x加上k 开根号:可以让x开根号,并且等 ...

【TOJ1132】Square Root，二次同余方程

【TOJ1132】Square Root，二次同余方程相关推荐

最新文章

热门文章