数学:二次剩余与n次剩余
二次剩余求的是这个东西
如果给定x,再给定若干个大的质数p,如果结果a相同,那么x是完全平方数?
/*poj 1808
题意:
判断平方剩余,即判断(x^2)%p=a是否有解。
限制:
|a| <= 1e9 && a % p !=0; 2 < p < 1e9 && p为奇素数。
思路:
用欧拉准则计算勒让德符号(用来判断平方剩余)
*/
#include<iostream>
#include<cstdio>
using namespace std;
#define LL __int64
LL a_b_MOD_c(LL a,LL b,LL mod){
LL ret = ;
a %= mod;
while(b){
if(b & ) ret = ret * a % mod;
a = a * a % mod;
b >>= ;
}
return ret;
}
//(x^2)%n=a 求平方剩余,n必须是奇素数
//注意:如果a为负,则看题意,是否要化为a=(a%n+n)%n
int modsqr(int a,int n){
int b,k,i,x;
if(n==) return a%n;
if(a_b_MOD_c(a,(n-)/,n)==){
if(n%==)
x=a_b_MOD_c(a,(n+)/,n);
else{
for(b=;a_b_MOD_c(b,(n-)/,n)==;b++){
i=(n-)/;
k=;
}
do{
i/=;
k/=;
if((a_b_MOD_c(a,i,n)*a_b_MOD_c(b,k,n)+)%n==)
k+=(n-)/;
}
while(i%==);
x=(a_b_MOD_c(a,(i+)/,n)*a_b_MOD_c(b,k/,n))%n;
}
if(x*>n)
x=n-x;
return x;
}
return -;
}
//用欧拉准则计算勒让德符号(用来判断平方剩余)
//表示为(a|p) a为整数,p为奇素数(所以m=2不适用勒让德符号),有三种情况。
//1. (a|p)=0, if(a%p==0)
//2. (a|p)=1, if(a%p!=0 && (x^2)%p=a 有整数解)
//3. (a|p)=-1,if((x^2)%p=a 无整数解)
//注意:如果a为负,则看题意,是否要化为a=(a%p+p)%p
int lrd(LL a,LL p){
LL ret=a_b_MOD_c(a,(p-)>>,p);
if(ret==)
return ;
return -;
}
int main(){
int T,cas=;
int a,n;
scanf("%d",&T);
while(T--){
scanf("%d%d",&a,&n);
a=(a+n)%n; //以后注意给出余数的时候,要注意它是不是负的 //求平方剩余
//cout<<modsqr(a,n)<<endl;
printf("Scenario #%d:\n%d\n\n",++cas,lrd(a,n));
}
return ;
}
然后是n次剩余
/*hdu 3930
题意:
给定newx, k, m, 方程 (x^k)%m=newx, 求在模m意义下的所有解x。
限制:
0 <= newx, m, k <= 1.5*10^15; m是素数。
思路:
N次剩余
*/
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <vector>
#include <algorithm>
using namespace std;
#define LL __int64
#define PB push_back
LL mul(LL a,LL b,LL m){
LL ret = ;
a %= m;
while(b){
if(b & ) ret = (ret + a) % m;
a = (a + a) % m;
b >>= ;
}
return ret;
}
LL a_b_MOD_c(LL a,LL b,LL m){
LL ret = ;
a %= m;
while(b){
if(b&) ret = mul(ret,a,m);
a = mul(a,a,m);
b >>= ;
}
return ret;
} LL ext_gcd(LL a,LL b,LL &x,LL &y){
if(b==) { x=, y=; return a; }
LL ret= ext_gcd(b,a%b,y,x);
y-= a/b*x;
return ret;
}
vector<LL> a;
bool g_test(LL g,LL p){
for(LL i=;i<a.size();++i)
if(a_b_MOD_c(g,(p-)/a[i],p)==)
return ;
return ;
}
LL pri_root(LL p){
a.clear();
LL tmp=p-;
for(LL i=;i<=tmp/i;++i)
if(tmp%i==){
a.push_back(i);
while(tmp%i==)
tmp/=i;
}
if(tmp!=)
a.push_back(tmp);
LL g=;
while(true){
if(g_test(g,p))
return g;
++g;
}
}
const int HASH_MOD=;
LL key[HASH_MOD], val[HASH_MOD];
int head[HASH_MOD], next[HASH_MOD];
struct Hash{
int tot;
void init(){
memset(head, -, sizeof(head));
tot = ;
}
LL insert(LL x, LL y){
int k = x % HASH_MOD;
key[tot] = x;
val[tot] = y;
next[tot] = head[k];
head[k] = tot++;
}
LL find(LL x){
int k = x % HASH_MOD;
for(int i = head[k]; i != -; i = next[i])
if(key[i] == x)
return val[i];
return -;
}
}hs;
//求解模方程a^x=b(mod m),n为素数,无解返回-1
//注意:要求0 < a < m; 0 <= b < m; 否则按题意自己转化。
//复杂度O(sqrt(m))
LL log_mod(LL a, LL b, LL m){
hs.init();
LL s = ceil(sqrt(m + 0.5));
LL cur = ;
for (int i = ; i < s; ++i){
if(hs.find(cur)==-) hs.insert(cur,i); //记得先判重,在插入
cur = cur * a % m;
} LL v = a_b_MOD_c(a, (m - s - + m) % m, m);
for(int i = ; i < s; ++i){
LL tmp = hs.find(b);
if(tmp!=-)
return s * i + tmp;
b=b*v%m;
}
return -;
}
/*n次剩余
任务:
给定N, a, p, 求出(x^N)%p=a 在模p意义下的所有解x。
说明:
令g为p的原根,因为p为素数,所以phi(p)=p-1。
由原根的性质得:
如果g为p的原根,则:g^i mod p != g^j mod p (p为素数), 其中i != j且i, j介於1至(p-1)之间
所以,可以设g^y=x, g^t=a,则有:
g^(y*N)%p=g^t
又由原根的性质:
g^(y*N)%p=g^t -> (y*N)%(p-1)=t (此方程可以由拓展欧几里得解)
另外g^t=a可以由离散对数求出
*/
vector<LL> residue(LL p, LL N, LL a){
LL g = pri_root(p);
g %= p;
LL m = log_mod(g, a, p);
vector<LL> ret;
if(a == ){
ret.PB();
return ret;
}
if(m == -)
return ret;
LL A = N, B = p - , C = m, x, y;
LL d = ext_gcd(A, B, x, y);
if(C % d != ) return ret;
x = x * (C / d) % B;
LL delta = B / d;
for(int i = ; i < d; ++i){
x = ((x + delta) % B + B) % B;
ret.PB(a_b_MOD_c(g, x, p));
}
sort(ret.begin(), ret.end());
ret.erase(unique(ret.begin(), ret.end()), ret.end());
return ret;
}
int main(){
int cas = ;
LL k,m,newx;
while(scanf("%I64d%I64d%I64d",&k, &m, &newx)!=EOF){
vector<LL> ans;
ans = residue(m,k,newx);
printf("case%d:\n",++cas);
if(ans.size()==) puts("-1");
for(int i = ; i < ans.size(); ++i)
printf("%I64d\n",ans[i]);
}
return ;
}
数学:二次剩余与n次剩余的更多相关文章
- 二次剩余、三次剩余、k次剩余
今天研究了一下这块内容...首先是板子 #include <iostream> #include <stdio.h> #include <math.h> #incl ...
- N次剩余和二次剩余
N次剩余 给定 \(N,a,P\),且 \(P\) 最好为质数 可以算出 \(x^N\equiv a(mod~p)\) 的解 首先可以算出 \(P\) 的原根 \(g\) 解方程 \(g^y\equi ...
- 数学杂烩总结(多项式/形式幂级数+FWT+特征多项式+生成函数+斯特林数+二次剩余+单位根反演+置换群)
数学杂烩总结(多项式/形式幂级数+FWT+特征多项式+生成函数+斯特林数+二次剩余+单位根反演+置换群) 因为不会做目录所以请善用ctrl+F 本来想的是笔记之类的,写着写着就变成了资源整理 一些有的 ...
- ACM数学
1.burnside定理,polya计数法 这个专题我单独写了个小结,大家可以简单参考一下:polya 计数法,burnside定理小结 2.置换,置换的运算 置换的概念还是比较好理解的,< ...
- OI数学汇总
最前面:\(\LaTeX\)可能需要加载一会,请耐心等待o~ 前言 数学在\(\text{OI}\)中十分重要.其中大多都是数论. 什么是数论? \[ 研究整数的理论 --zzq \] 本文包含所有侧 ...
- RSA加密数学原理
RSA加密数学原理 */--> *///--> *///--> UP | HOME RSA加密数学原理 Table of Contents 1 引言 2 RSA加密解密过程 2.1 ...
- codeforces 687B - Remainders Game 数学相关(互质中国剩余定理)
题意:给你x%ci=bi(x未知),是否能确定x%k的值(k已知) ——数学相关知识: 首先:我们知道一些事情,对于k,假设有ci%k==0,那么一定能确定x%k的值,比如k=5和ci=20,知道x% ...
- [HNOI 2013] 旅行 (数学)
感觉此题难啊,数学还是太渣了,看了半天的题解才算明白了点儿. 题目大意 给一个长度为n且仅由1和-1组成的序列ai, i = 1, 2, ..., n,每个位置都有另一个值vi,要求用某种方案将序列划 ...
- POJ 2359 Questions(约瑟夫环——数学解法)
题目链接: http://poj.org/problem?id=2359 题意描述: 输入一个字符串 按照下面的规则,如果剩下的最后一个字符是'?',输出"Yes",如果剩下的最后 ...
随机推荐
- Linux下安装maven(mvn命令)
Maven(mvn)是基于项目对象模型(POM project object model),可以通过一小段描述信息(配置)来管理项目的构建,报告和文档的软件项目管理工具(百度百科) 简单理解为一个打包 ...
- Unity接SDK通用方法总结
第一篇博客,回顾接UnitySDK的坑 接SDK主要有两种方式,第三方SDK和手动接(我用的android studio) 首先接触到的SDK就是Facebook-Unity的SDK,主要就是face ...
- UE4中Timeline的使用
UE4中经常需要一些和时间相联系的功能,例如在一段时间内完成一个动作,播放一段动画,或者只是单纯的延迟函数的执行时间,即调整事件的执行顺序.在UE4的蓝图自带函数中有一个很好用的函数可以完美地解决这些 ...
- spring boot 2.0 源码分析(二)
在上一章学习了spring boot 2.0启动的大概流程以后,今天我们来深挖一下SpringApplication实例变量的run函数. 先把这段run函数的代码贴出来: /** * Run the ...
- ejs 模板使用方法
http://embeddedjs.com/ Embedded JS Templates Embedded JS(EJS) 来源于ERB模板,且与ERB有很多相似之处.它有着与ERB相同的Tag,且包 ...
- (Alpha)Let's-展示博客
Let's Alpha 项目答辩 ·选题由来 手机端——用户相对较多,使用环境限制相对宽松 手机游戏?校园p2p应用?线下交流!(滴滴打水?) 模式的选择:发起——加入活动 ...
- 老李的blog使用日记(2)
寥寥数语结束一个不曾期待的遇见,可还是剧情不会这样结束,他也会在我的时间里注册自己的专属账号,无论什么时候,他会时而需要被注视着,为了达到目的,即使不择手段,只为一次擦肩而过的邂逅,极短的一段时间,相 ...
- ElasticSearch搜索实例含高亮显示及搜索的特殊字符过滤
应用说明见代码注解. 1.简单搜索实例展示: public void search() throws IOException { // 自定义集群结点名称 String clusterName = & ...
- 第九周PSP&进度条
PSP 一.表格: D日期 C类型 C内容 S开始时间 E结束时间 I时间间隔 T净时间(mins) 预计花费时间(mins) 11月11号 讨论 讨论beta发布 09:00 09:54 1 ...
- PHP的魔术方法
PHP中把以两个下划线__开头的方法称为魔术方法(Magic methods) 魔术方法包括: __construct(),类的构造函数 __destruct(),类的析构函数 __call(),在对 ...