51Nod1123 X^A Mod B 数论 中国剩余定理 原根 BSGS
原文链接https://www.cnblogs.com/zhouzhendong/p/51Nod1123.html
题目传送门 - 51Nod1123
题意
$T$ 组数据。
给定 $A,B,C$,求出使得 $x^A \equiv C \pmod B$ 的所有 $x$,保证解的个数不超过 $\sqrt B$ 。
$T\leq 100,1\leq A,B,C \leq 10^9$
题解
先记一下写这一题的感受:
1. 写的过程中代码长度峰值达到过 300 行,好久没写码农题了,感到自己码力大减。
2. 我貌似只会推 exCRT 了,互质的 CRT 我居然都记错了,而且还推错了,浪费了好多时间。
3. BZOJ2219 的数据太水,我在做这题的时候找到了 hack 我那一题的数据……
4. 数论真有(du)趣(liu)。
UPD(2018-09-10):
详见数论总结。
传送门 - https://www.cnblogs.com/zhouzhendong/p/Number-theory-Residue-System.html
代码
#include <bits/stdc++.h>
using namespace std;
const int N=100005;
int read(){
int x=0;
char ch=getchar();
while (!isdigit(ch))
ch=getchar();
while (isdigit(ch))
x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
return x;
}
int pcnt,f[N],Prime[N];
void Get_Prime(int n){
memset(f,0,sizeof f);
pcnt=0;
for (int i=2;i<=n;i++){
if (f[i])
continue;
Prime[++pcnt]=i;
for (int j=i+i;j<=n;j+=i)
f[j]=1;
}
}
void Divide(int x,int *p,int *q,int &cnt){
cnt=0;
for (int i=1;i<=pcnt&&Prime[i]*Prime[i]<=x;i++){
if (x%Prime[i])
continue;
p[++cnt]=Prime[i],q[cnt]=0;
while (x%p[cnt]==0)
x/=p[cnt],q[cnt]++;
}
if (x>1)
p[++cnt]=x,q[cnt]=1;
}
int Pow(int x,int y,int mod){
int ans=1;
for (;y;y>>=1,x=1LL*x*x%mod)
if (y&1)
ans=1LL*ans*x%mod;
return ans;
}
int Pow(int x,int y){
return Pow(x,y,2e9);
}
int gcd(int a,int b){
return b?gcd(b,a%b):a;
}
int Fac[50],Fac_cnt=0;
bool Get_g_Check(int P,int C,int x){
int phi=Pow(P,C-1)*(P-1),pw=Pow(P,C);
if (C>1&&Pow(x,phi/P,pw)==1)
return 0;
for (int i=1;i<=Fac_cnt;i++)
if (Pow(x,phi/Fac[i],pw)==1)
return 0;
return 1;
}
int Get_g(int P,int C){
int v=P-1;
Fac_cnt=0;
for (int i=1;i<=pcnt&&Prime[i]*Prime[i]<=v;i++)
if (v%Prime[i]==0){
Fac[++Fac_cnt]=Prime[i];
while (v%Prime[i]==0)
v/=Prime[i];
}
if (v>1)
Fac[++Fac_cnt]=v;
for (int i=2;;i++)
if (Get_g_Check(P,C,i))
return i;
return -1;
}
struct hash_map{
static const int Ti=233,mod=1<<16;
int cnt,k[mod+1],v[mod+1],nxt[mod+1],fst[mod+1];
int Hash(int x){
int v=x&(mod-1);
return v==0?mod:v;
}
void clear(){
cnt=0;
memset(fst,0,sizeof fst);
}
void update(int x,int a){
int y=Hash(x);
for (int p=fst[y];p;p=nxt[p])
if (k[p]==x){
v[p]=a;
return;
}
k[++cnt]=x,nxt[cnt]=fst[y],fst[y]=cnt,v[cnt]=a;
return;
}
int find(int x){
int y=Hash(x);
for (int p=fst[y];p;p=nxt[p])
if (k[p]==x)
return v[p];
return 0;
}
int &operator [] (int x){
int y=Hash(x);
for (int p=fst[y];p;p=nxt[p])
if (k[p]==x)
return v[p];
k[++cnt]=x,nxt[cnt]=fst[y],fst[y]=cnt;
return v[cnt]=0;
}
}Map;
int BSGS(int A,int B,int P){
int M=max((int)(0.8*sqrt(1.0*P)),1),AM=Pow(A,M,P);
Map.clear();
for (int b=0,pw=B;b<M;b++,pw=1LL*pw*A%P)
Map.update(pw,b+1);
for (int a=M,pw=AM;a-M<P;a+=M,pw=1LL*pw*AM%P){
int v=Map.find(pw);
if (v)
return a-(v-1);
}
return -1;
}
int exgcd(int a,int b,int &x,int &y){
if (!b){
x=1,y=0;
return a;
}
int res=exgcd(b,a%b,y,x);
y-=(a/b)*x;
return res;
}
int RHD(vector <int> &v,int A,int B,int P,int C,int flag){
// printf("RHD :: %d %d %d %d\n",A,B,P,C);
int g=Get_g(P,C);
int t=BSGS(g,B,Pow(P,C));
int mod=(P-1)*Pow(P,C-1),pw=mod/(P-1)*P;
int GCD=gcd(mod,gcd(A,t));
if (!flag)
return gcd(A,mod)>GCD?0:GCD;
v.clear();
if (gcd(A,mod)>GCD)
return 0;
int _A=A/GCD,_t=t/GCD,_mod=mod/GCD;
int x,y;
exgcd(_A,_mod,y,x);
y=(y%_mod+_mod)%_mod;
y=1LL*y*_t%_mod;
for (int i=y;i<mod;i+=_mod)
v.push_back(Pow(g,i,pw));
return GCD;
}
vector <int> tmp;
void dfs(vector <int> &v,int A,int B,int C,int t,int d){
if (Pow(t,A,1<<d)!=B%(1<<d))
return;
if (C==d)
return v.push_back(t);
dfs(v,A,B,C,t,d+1);
dfs(v,A,B,C,t|(1<<d),d+1);
}
int solve2(vector <int> &v,int A,int B,int C){
v.clear();
B%=1<<C;
dfs(v,A,B,C,0,0);
return (int)v.size();
}
int solve(vector <int> &v,int A,int B,int P,int C,int flag){
if (P==2)
return solve2(v,A,B,C);
if (flag)
v.clear();
int pw=Pow(P,C),Phi=(P-1)*Pow(P,C-1);
B%=pw;
if (B==0){
if (!flag)
return Pow(P,C-((C+A-1)/A));
int base=Pow(P,((C+A-1)/A));
for (int i=0;i<pw;i+=base)
v.push_back(i);
return (int)v.size();
}
int g=gcd(B,pw),Q=0;
B/=g;
while (g>1)
g/=P,Q++;
if (!flag)
return Pow(P,Q-Q/A)*((Q%A)?0:RHD(v,A,B,P,C-Q,0));
if (Q%A)
return 0;
RHD(tmp,A,B,P,C-Q,1);
int base=Pow(P,C-Q+Q/A),tg=Pow(P,Q/A);
for (vector <int> :: iterator i=tmp.begin();i!=tmp.end();i++)
for (int j=*i=1LL*(*i)*tg%base;j<pw;j+=base)
v.push_back(j);
return (int)v.size();
}
vector <int> res[50];
set <int> Ans;
int crtk[50];
void dfsCRT(int x,int d,int dcnt,int *p,int *q,int Fac){
if (d>dcnt){
Ans.insert(x);
return;
}
for (int j=0;j<res[d].size();j++)
dfsCRT((1LL*res[d][j]*crtk[d]+x)%Fac,d+1,dcnt,p,q,Fac);
}
int main(){
Get_Prime(1e5);
int T=read();
while (T--){
int A=read(),P=read(),B=read();
int cnt,p[50],q[50];
Divide(P,p,q,cnt);
swap(p[1],p[cnt]);
swap(q[1],q[cnt]);
int ans=1;
for (int i=1;i<=cnt&&ans>0;i++)
ans*=solve(res[i],A,B,p[i],q[i],0);
if (!ans){
puts("No Solution");
continue;
}
for (int i=1;i<=cnt;i++)
solve(res[i],A,B,p[i],q[i],1);
for (int i=1;i<=cnt;i++){
int v1=Pow(p[i],q[i]),v2=P/Pow(p[i],q[i]),x,y;
exgcd(v2,v1,x,y);
x=(x%v1+v1)%v1;
crtk[i]=1LL*v2*x%P;
}
Ans.clear();
dfsCRT(0,1,cnt,p,q,P);
for (set <int> :: iterator i=Ans.begin();i!=Ans.end();i++)
printf("%d ",*i);
puts("");
}
return 0;
}
51Nod1123 X^A Mod B 数论 中国剩余定理 原根 BSGS的更多相关文章
- BZOJ2219 数论之神 数论 中国剩余定理 原根 BSGS
原文链接https://www.cnblogs.com/zhouzhendong/p/BZOJ2219.html 题目传送门 - BZOJ2219 题意 求同余方程 $x^A\equiv B \pmo ...
- BZOJ2219数论之神——BSGS+中国剩余定理+原根与指标+欧拉定理+exgcd
题目描述 在ACM_DIY群中,有一位叫做“傻崽”的同学由于在数论方面造诣很高,被称为数轮之神!对于任何数论问题,他都能瞬间秒杀!一天他在群里面问了一个神题: 对于给定的3个非负整数 A,B,K 求出 ...
- 数学--数论--中国剩余定理 拓展 HDU 1788
再次进行中国余数定理 问题描述 我知道部分同学最近在看中国剩余定理,就这个定理本身,还是比较简单的: 假设m1,m2,-,mk两两互素,则下面同余方程组: x≡a1(mod m1) x≡ a2(mod ...
- 【bzoj1951】: [Sdoi2010]古代猪文 数论-中国剩余定理-Lucas定理
[bzoj1951]: [Sdoi2010]古代猪文 因为999911659是个素数 欧拉定理得 然后指数上中国剩余定理 然后分别lucas定理就好了 注意G==P的时候的特判 /* http://w ...
- POJ 1006 Biorhythms (数论-中国剩余定理)
Biorhythms Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 111285 Accepted: 34638 Des ...
- 【中国剩余定理】【容斥原理】【快速乘法】【数论】HDU 5768 Lucky7
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5768 题目大意: T组数据,求L~R中满足:1.是7的倍数,2.对n个素数有 %pi!=ai 的数 ...
- 数论F - Strange Way to Express Integers(不互素的的中国剩余定理)
F - Strange Way to Express Integers Time Limit:1000MS Memory Limit:131072KB 64bit IO Format: ...
- 中国剩余定理(CRT)与欧拉函数[数论]
中国剩余定理 ——!x^n+y^n=z^n 想必大家都听过同余方程这种玩意,但是可能对于中国剩余定理有诸多不解,作为一个MOer&OIer,在此具体说明. 对于同余方程: x≡c1(mod m ...
- BZOJ1951 [Sdoi2010]古代猪文 中国剩余定理 快速幂 数论
原文链接http://www.cnblogs.com/zhouzhendong/p/8109156.html 题目传送门 - BZOJ1951 题意概括 求 GM mod 999911659 M=∑i ...
随机推荐
- 时间日期date/cal
命令: date 作用: 查看下系统时间 使用: date 命令: cal 对应英文: calendar 作用: 查看日历 选项: -y:可查看一年的日历 使用: cal cal -y
- Struts2配置拦截器
<package name="loginaction" namespace="/" extends="struts-default"& ...
- Java希尔排序算法
希尔排序就是对直接插入排序的一个优化.比如有这么一种情况:对一个无序数组进行从小到大的排序,但是数组的最后一个位置的数是最小的,我们要把它挪到第一个位置,其他位置的都要往后移动,要是这个数组非常大,那 ...
- ACM-ICPC 2018 焦作赛区网络预赛 G Give Candies
There are NNN children in kindergarten. Miss Li bought them NNN candies. To make the process more in ...
- Webform中<%%>
其实<%%>很早之前就见过了,只是会用一点功能,其它的不甚了解.今天偶尔见到了它的庐山真面目,现在共享给大家. 语法 代码块呈现(<%%>)定义了当呈现页时执行的内联代码或内联 ...
- Confluence 6 属性的一个示例
下面是有关 Confluence 页面被调用的前几行的访问概述. [344ms] - /display/ds/Confluence+Overview [313ms] - SiteMesh: parse ...
- 访问 Confluence 6 的计划任务配置
希望访问 Confluence 计划任务配置界面: 进入 > 基本配置(General Configuration) > 计划任务(Scheduled Jobs) 所有的计划任务将会按照 ...
- Confluence 6 "Duplicate Key" 相关问题解决
如果你遇到了下面的错误信息,例如: could not insert: [bucket.user.propertyset.BucketPropertySetItem#bucket.user.prope ...
- 基于多进程的Tcp套接字服务器
服务端 import socketfrom multiprocessing import Process def task(c): print('顾客吃点啥') while True: data = ...
- Python计算器实操
要求: 开发一个简单的python计算器 实现加减乘除及拓号优先级解析 用户输入 1 - 2 * ( (60-30 +(-40/5) * (9-2*5/3 + 7 /3*99/4*2998 +10 * ...