【BZOJ】2242: [SDOI2011]计算器
http://www.lydsy.com/JudgeOnline/problem.php?id=2242
题意:(前两个问略...)第三个问是,求$a^x \equiv b \pmod{p}$最小的$x$,或者输出无解,它们范围都是$10^9$哒= =
#include <bits/stdc++.h>
using namespace std; typedef long long ll;
int mpow(ll a, ll b, ll p) {
ll r=1; a%=p;
while(b) { if(b&1) r=((ll)r*a)%p; a=((ll)a*a)%p; b>>=1; }
return r;
}
void gcd(ll a, ll b, ll &d, ll &x, ll &y) {
if(!b) { d=a; x=1; y=0; return; }
gcd(b, a%b, d, y, x); y-=a/b*x;
}
void ni(ll a, ll b, ll p) {
ll d, x, y, t;
gcd(a, p, d, x, y); if(b%d) { puts("Orz, I cannot find x!"); return; }
t=p/d;
while(x<0) x+=t;
while(x>=t) x-=t;
printf("%lld\n", (x*b)%p);
}
map<int, int> s;
void bsgs(ll y, ll z, ll p) {
y%=p; z%=p;
if(z==1) { puts("0"); return; }
if(!y && !z) { puts("1"); return; }
if(!y) { puts("Orz, I cannot find x!"); return; }
s.clear();
int m=sqrt(p+0.5), t=1, w=y;
for(int i=0; i<m; ++i) s[((ll)z*t)%p]=i, t=((ll)t*w)%p;
w=mpow(y, m, p); y=1; t=(p-1)/m+1; bool flag=1;
for(int i=0; i<=t; ++i) if(s.count(y)) { printf("%lld\n", (ll)m*i-s[y]); flag=0; break; } else y=((ll)y*w)%p;
if(flag) puts("Orz, I cannot find x!");
}
int main() {
int z, y, p, c, T;
scanf("%d%d", &T, &c);
while(T--) {
scanf("%d%d%d", &y, &z, &p);
if(c==1) printf("%d\n", mpow(y, z, p));
else if(c==2) ni(y, z, p);
else bsgs(y, z, p);
}
return 0;
}
bsgs裸题....其实就是一种分块思想..(为啥有那么牛的名字呢= =其实是我不想加分类了= =)即小块暴力然后大块就解决的思想,相信你们都能秒懂= =
要求
$$a^x \equiv b \pmod{p}$$
的最小的$x$,那么
首先我们随便选一个$m$,使得$x=km-t, 0<=t<m$,(这虽然有点区别于取余,但是这是为了方便= =)
然后推得
$$a^{km} \equiv ba^t \pmod{p}$$
然后就是右边暴力预处理,左边枚举$k$...由于枚举$k$复杂度是$O(n/m)$,显然取$m=\sqrt{n}$最优= =...由于懒,开个set记录右边= =于是总复杂度是$O(\sqrt{n}log(\sqrt{n}))$
哦最后忘记一件事= =要特判啊= =比如说$b=1$显然$x=0$啊,而且取模了$a, b$后注意特判$a=0$的情况啊= =关于一些情况没特判的同学= =我要hack你们...比如数据
1 3
4 1 2
【BZOJ】2242: [SDOI2011]计算器的更多相关文章
- bzoj 2242: [SDOI2011]计算器 BSGS+快速幂+扩展欧几里德
2242: [SDOI2011]计算器 Time Limit: 10 Sec Memory Limit: 512 MB[Submit][Status][Discuss] Description 你被 ...
- BZOJ 2242: [SDOI2011]计算器( 快速幂 + 扩展欧几里德 + BSGS )
没什么好说的... --------------------------------------------------------------------- #include<cstdio&g ...
- BZOJ 2242: [SDOI2011]计算器 [快速幂 BSGS]
2242: [SDOI2011]计算器 题意:求\(a^b \mod p,\ ax \equiv b \mod p,\ a^x \equiv b \mod p\),p是质数 这种裸题我竟然WA了好多次 ...
- bzoj 2242 [SDOI2011]计算器(数论知识)
Description 你被要求设计一个计算器完成以下三项任务: 1.给定y,z,p,计算Y^Z Mod P 的值: 2.给定y,z,p,计算满足xy≡ Z ( mod P )的最小非负整数: 3.给 ...
- BZOJ.2242.[SDOI2011]计算器(扩展欧几里得 BSGS)
同余方程都不会写了..还一直爆int /* 2.关于同余方程ax ≡b(mod p),可以用Exgcd做,但注意到p为质数,y一定有逆元 首先a%p=0时 仅当b=0时有解:然后有x ≡b*a^-1( ...
- BZOJ 2242 [SDOI2011]计算器(快速幂+Exgcd+BSGS)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=2242 [题目大意] 给出T和K 对于K=1,计算 Y^Z Mod P 的值 对于K=2 ...
- bzoj 2242 [SDOI2011]计算器——BSGS模板
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2242 第一道BSGS! 咳咳,我到底改了些什么?…… 感觉和自己的第一版写的差不多……可能是 ...
- BZOJ 2242 [SDOI2011]计算器 BSGS+高速幂+EXGCD
题意:id=2242">链接 方法: BSGS+高速幂+EXGCD 解析: BSGS- 题解同上.. 代码: #include <cmath> #include <c ...
- bzoj 2242: [SDOI2011]计算器
#include<cstdio> #include<iostream> #include<map> #include<cmath> #define ll ...
- [原博客] BZOJ 2242 [SDOI2011] 计算器
题目链接 noip级数论模版题了吧.让求三个东西: 给定y,z,p,计算`Y^Z Mod P` 的值. 给定y,z,p,计算满足`xy≡ Z ( mod P )`的最小非负整数. 给定y,z,p,计算 ...
随机推荐
- Delphi多线程开发注意事项
Q1: 多线程中需避免多个线程同时向全局变量进行写入操作,导致访问冲突问题. A1: 可以通过使用加锁机制(比如:临界区.互斥.信号量)解决此问题. Q2:多线程中对于结构体和CLASS类型的全局变 ...
- MVC中Form表单的提交
概述 Web页面进行Form表单提交是数据提交的一种,在MVC中Form表单提交到服务器.服务端接受Form表单的方式有多种,如果一个Form有2个submit按钮,那后台如何判断是哪个按钮提交的数据 ...
- 【译】DotNet 5.4或者说 .NET平台标准
原文:DotNet 5.4 or .NET Platform Standards 是的,你看到的确实是"DotNet 5.4".使用Visual Studio 2015 RC1 u ...
- php获取文件夹下面的文件列表和文件夹列表
function getDir($dir) { $dirArray[] = NULL; if (false != ($handle = opendir( $dir ))) { $i=0; while ...
- PHP计算程序运行时间的类
<?php class runTime{ private $starTime; private $stopTime; private function getMicTime(){ $mictim ...
- hdu 4753 2013南京赛区网络赛 记忆化搜索 ****
看到范围基本可以想到dp了,处理起来有点麻烦 #include<iostream> #include<cstdio> #include<cstring> #incl ...
- 第二篇:SOUI源码的获取及编译
源代码的获取 SOUI的源码采用SVN管理. SVN:http://code.taobao.org/svn/soui2 这里主要包含两个目录:trunk 及 third-part. trunk目录保存 ...
- SpringHttpInvoker解析2-服务端实现
主要的配置文件 <!-- 在Spring的httpInvoker服务 --> <bean id="httpInvokerUserService" class=&q ...
- 2016.9.13 JavaScript入门之六基础函数
1.Math.random()函数: 2.math.floor(x)返回小于参数x的最大整数,即对浮点数向下取整. 例如:random本身只产生(0~1)之间的小数,random()*10 意思是产生 ...
- DSP using MATLAB 示例 Example3.19
代码: % Analog Signal Dt = 0.00005; t = -0.005:Dt:0.005; xa = exp(-1000*abs(t)); % Discrete-time Signa ...