A^x = D (mod P)

0 <= x <= M, here M is a given integer.

1 <= A, P < 2^31, 0 <= D < P, 1 <= M < 2^63

----------------------------------------------------------

裸拓展baby step giant step

先转成非拓展版,然后如果在转的过程中就出解了,则return 1,否则就找出=D的起点和循环节长度,然后直接求解。

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
const int Hash_MOD = 1000003;
typedef long long LL;
typedef unsigned long long ULL;
LL times;
struct link
{
LL link;
int val,next;
}es[100000];
struct triple
{
LL x,y,g;
triple(const LL _x = 0,const LL _y = 0,const LL _g = 0):
x(_x),y(_y),g(_g){}
};
int H[Hash_MOD * 2],ec;
triple exgcd(const LL a,const LL b)
{
if (!b) return triple(1,0,a);
const triple last(exgcd(b,a%b));
return triple(last.y,last.x - a / b * last.y,last.g);
}
LL rem_equ(const LL a,const LL b,const LL c)
{
//ax == c (mod b)
//ax mod b == c
const triple tmp(exgcd(a,b));
const LL MOD = (b / tmp.g);
return ((tmp.x * (c / tmp.g) % MOD) + MOD) % MOD;
}
LL find(const LL x)
{
for (int i = H[x%Hash_MOD]; i; i = es[i].next)
if (es[i].link == x) return es[i].val;
return -1;
}
void insert(LL x,const LL v)
{
if (find(x) != -1) return;
const LL key = x % Hash_MOD;
es[++ec].next = H[key];
H[key] = ec;
es[ec].link = x;
es[ec].val = v;
}
LL BSGS(LL A,LL P,LL D)
{
LL AA = 1 % P,x = 1 % P,MOD = P,DD = D,m = static_cast<LL>(std::ceil(std::sqrt((double)P)));
times = 0;
for (triple tmp; (tmp = exgcd(A,P)).g != 1; ) {
if (x == DD) return times++;
if (D % tmp.g) return -1;
P /= tmp.g;
D /= tmp.g;
(AA *= A / tmp.g) %= P;
++times;
(x *= A) %= MOD;
}
A %= P;
ec = 0;
memset(H,0,sizeof H);
LL tmp = 1 % P;
for (int i = 0; i < m; ++i,(tmp *= A) %= P)
insert(tmp,i);
x = 1 % P;
for (LL i = 0; i < m; ++i,(x *= tmp) %= P) {
const int j = find(rem_equ(AA*x%P,P,D));
if (j != -1) return i * m + j + times;
}
return -1;
}
LL getphi(LL x)
{
LL res = 1;
for (LL i = 2; i * i <= x; ++i)
if (x % i == 0) {
x /= i;
res *= i - 1;
while (x % i == 0) x /= i,res *= i;
}
if (x > 1) res *= x - 1;
return res;
}
LL power(LL a,LL k,LL MOD)
{
//a ^ k % MOD
LL res = 1 % MOD;
for (; k; k/=2,(a *= a) %= MOD)
if (k & 1) (res *= a) %= MOD;
return res;
}
LL calc_len(const LL start,const LL A,const LL P,const LL D)
{
//A ^ (start + *this) == A ^ start == D (mod P)
LL phi = getphi(P),res = phi;
for (LL i = 2; i * i <= phi; ++i) {
for (; phi % i == 0; phi /= i);
for (; res % i == 0 && power(A,start + res/i,P) == D; res /= i);
}
if (phi > 1)
for (; res % phi == 0 && power(A,start + res/phi,P) == D; res /= phi);
return res;
}
ULL work(const LL A,const LL P,const LL D,const ULL M)
{
LL start = BSGS(A,P,D);
//printf("%lld\n",start);
if (start == -1 || start > M) return 0;
else if (start < times) return 1;
// LL phi=getphi(P);
// if (power(A,start + phi,P) != D) return 1;
ULL len = calc_len(start,A,P,D);
return (M - start) / len + 1;
}
int main()
{
#ifndef ONLINE_JUDGE
freopen("3254.in","r",stdin);
freopen("3254.out","w",stdout);
#endif
LL A,P,D;
ULL M;
while (~scanf("%lld%lld%lld%llu",&A,&P,&D,&M))
printf("%llu\n",work(A%P,P,D,M));
}

[ZOJ3254] MON 9.2009Secret Code的更多相关文章

  1. 【zoj3254】Secret Code

    题意: 给出a.p.d.m 求a^x=d(mod p) 0<=x<=m 的解的个数 题解: 今天一整天的时间大部分都在调这题Orz BSGS什么的还是太不熟了 我们可以用BSGS拓展版求出 ...

  2. Visual Studio Code 代理设置

    Visual Studio Code (简称 VS Code)是由微软研发的一款免费.开源的跨平台文本(代码)编辑器,在十多年的编程经历中,我使用过非常多的的代码编辑器(包括 IDE),例如 Fron ...

  3. 我们是怎么做Code Review的

    前几天看了<Code Review 程序员的寄望与哀伤>,想到我们团队开展Code Review也有2年了,结果还算比较满意,有些经验应该可以和大家一起分享.探讨.我们为什么要推行Code ...

  4. Code Review 程序员的寄望与哀伤

    一个程序员,他写完了代码,在测试环境通过了测试,然后他把它发布到了线上生产环境,但很快就发现在生产环境上出了问题,有潜在的 bug. 事后分析,是生产环境的一些微妙差异,使得这种 bug 场景在线下测 ...

  5. 从Script到Code Blocks、Code Behind到MVC、MVP、MVVM

    刚过去的周五(3-14)例行地主持了技术会议,主题正好是<UI层的设计模式——从Script.Code Behind到MVC.MVP.MVVM>,是前一天晚上才定的,中午花了半小时准备了下 ...

  6. 在Visual Studio Code中配置GO开发环境

    一.GO语言安装 详情查看:GO语言下载.安装.配置 二.GoLang插件介绍 对于Visual Studio Code开发工具,有一款优秀的GoLang插件,它的主页为:https://github ...

  7. 代码的坏味道(14)——重复代码(Duplicate Code)

    坏味道--重复代码(Duplicate Code) 重复代码堪称为代码坏味道之首.消除重复代码总是有利无害的. 特征 两个代码片段看上去几乎一样. 问题原因 重复代码通常发生在多个程序员同时在同一程序 ...

  8. http status code

    属于转载 http status code:200:成功,服务器已成功处理了请求,通常这表示服务器提供了请求的网页 404:未找到,服务器未找到 201-206都表示服务器成功处理了请求的状态代码,说 ...

  9. Visual Studio Code——Angular2 Hello World 之 2.0

    最近看到一篇用Visual Studio Code开发Angular2的文章,也是一篇入门教程,地址为:使用Visual Studio Code開發Angular 2專案.这里按部就班的做了一遍,感觉 ...

随机推荐

  1. IFeatureCursorProxy.flush AutomationException: 0x80041538

    添加面的时候碰到的一个问题,有些数据没问题,有些报错,后来请教一位同事说有可能是经纬度字段的数据精度问题,因为投影坐标系统不同,支持的经纬度经度不同,后来转换投影坐标系统后果然解决问题了,我一开始也怀 ...

  2. session_write_close()的作用

    简单地说,当开启session_start以后,这个session会一直开启,并且被一个用户使用.其他用户开启session的话要等待第一个session用户关闭以后才可以开启sessio,这样就造成 ...

  3. php常用函数——数学函数

    php常用函数——数学函数

  4. 每天一条linux命令(1):ls命令

    ls命令是linux下最常用的命令.ls命令就是list的缩写缺省下ls用来打印出当前目录的清单如果ls指定其他目录那么就会显示指定目录里的文件及文件夹清单. 通过ls 命令不仅可以查看linu ...

  5. Java八种基本类型

    boolean 二进制位: true ,false   byte 二进制位:8 -128 - 127   -2的7次方到2的7次方-1 char 二进制位:16 0 - 65535   short 二 ...

  6. Linux端口占用

    1.netstat netstat -anp | grep 23232 Sample: [root@BICServer 0825]# netstat -anp | grep 23232 tcp 0 0 ...

  7. python基础===100盏灯的问题

    闪存里有人这样提问这样: 第一轮操作所有电灯,第二轮操作第2盏,第4盏开关,以此类推,第三轮改变编号为3的倍数的电灯,第3盏,第6盏,如果原来那盏灯是亮的,就熄灭它,如果原来是灭的,就点亮它,以此类推 ...

  8. 健身VS不健身,完全是两种不同的人生!

    这两天一组同龄人合照 刷爆了国内健身圈, 图左是一位67岁的老人, 图右是67岁的健美运动员杨新民老师 相同年龄, 但从外观上有着强烈的距离感! 让多人不禁感叹,健身和不健身, 简直就是两种状态,两种 ...

  9. Hashtable之Properties

    properties的使用:1.Hashtable的实现类,线程安全.与HashMap不同,Hashtable不允许使用null作为key和value2.和HashMap一样,Hashtable也不能 ...

  10. Nginx常见错误及处理方法

    转载:https://www.cnblogs.com/liyongsan/p/6795851.html 404 bad request 一般原因:请求的Header过大 解决方法:配置nginx.co ...