bzoj3122 [SDOI2013]随机数生成器
给定一个递推式, \(X_i=(aX_{i-1}+b)\mod P\)
求满足 \(X_k=t\) 的最小整数解,无解输出 \(-1\)
\(0\leq a,\ b,\ t,\ P\leq10^9,\ P\) 为质数
BSGS
首先化式子,推得
\]
因此
\]
所以上 \(BSGS\)
然而这题特判很恶心,不加特判 \(0\text{pts}\)
特判如下:
- \(x=t:ans=1\)
- \(a=1\)
- \(b=0:ans=-1\)
- \(b\neq0:ans=\frac{t-x}{b}+1\)
- \(a=0\)
- \(b=t:ans=2\)
- \(b\neq t:ans=-1\)
时间复杂度 \(O(T\sqrt P)\)
代码
#include <bits/stdc++.h>
using namespace std;
int P;
int qp(int a, int k) {
int res = bool(a);
for (; k; k >>= 1, a = 1ll * a * a % P) {
if (k & 1) res = 1ll * res * a % P;
}
return res % P;
}
int bsgs(int a, int b) {
if (!a && b) return -1;
map <int, int> s;
int sz = sqrt(P), inv_a = qp(a, P - 2), pw = qp(a, sz), cur = 1;
for (int i = 0; i <= sz; i++) {
s.insert(make_pair(1ll * b * cur % P, i)), cur = 1ll * cur * inv_a % P;
}
cur = 1;
map <int, int> :: iterator it;
for (int i = 0; i <= sz; i++, cur = 1ll * cur * pw % P) {
if ((it = s.find(cur)) != s.end()) {
return i * sz + (it -> second);
}
}
return -1;
}
int main() {
int Tests, a, b, x, t, A, B;
scanf("%d", &Tests);
while (Tests--) {
scanf("%d %d %d %d %d", &P, &a, &b, &x, &t);
a %= P, b %= P, x %= P, t %= P;
if (x == t) {
puts("1"); continue;
} else if (a == 1) {
if (!b) {
puts("-1"); continue;
}
printf("%d\n", 1ll * (t - x + P) * qp(b, P - 2) % P + 1);
continue;
} else if (!a) {
puts(b == t ? "2" : "-1");
continue;
}
A = a, B = 1ll * (1ll * a * t - t + b + P) % P * qp((b - x + 1ll * a * x + P) % P, P - 2) % P;
int ans = bsgs(A, B);
printf("%d\n", ~ans ? ans + 1 : ans);
}
return 0;
}
bzoj3122 [SDOI2013]随机数生成器的更多相关文章
- BZOJ3122: [Sdoi2013]随机数生成器(BSGS)
题意 题目链接 Sol 这题也比较休闲. 直接把\(X_{i+1} = (aX_i + b) \pmod P\)展开,推到最后会得到这么个玩意儿 \[ a^{i-1} (x_1 + \frac{b}{ ...
- bzoj千题计划259:bzoj3122: [Sdoi2013]随机数生成器
http://www.lydsy.com/JudgeOnline/problem.php?id=3122 等比数列求和公式+BSGS #include<map> #include<c ...
- [bzoj3122][SDOI2013]随机数生成器 ——BSGS,数列
题目大意 给定递推序列: F[i] = a*F[i-1] + b (mod c) 求一个最小的i使得F[i] == t 题解 我们首先要化简这个数列,作为一个学渣,我查阅了一些资料: http://d ...
- BZOJ3122 [Sdoi2013]随机数生成器 【BSGS】
题目 输入格式 输入含有多组数据,第一行一个正整数T,表示这个测试点内的数据组数. 接下来T行,每行有五个整数p,a,b,X1,t,表示一组数据.保证X1和t都是合法的页码. 注意:P一定为质数 输出 ...
- 【BZOJ3122】[Sdoi2013]随机数生成器 BSGS+exgcd+特判
[BZOJ3122][Sdoi2013]随机数生成器 Description Input 输入含有多组数据,第一行一个正整数T,表示这个测试点内的数据组数. 接下来T行,每行有五个整数p,a,b, ...
- 【bzoj3122】: [Sdoi2013]随机数生成器 数论-BSGS
[bzoj3122]: [Sdoi2013]随机数生成器 当a>=2 化简得 然后 BSGS 求解 其他的特判 : 当 x=t n=1 当 a=1 当 a=0 判断b==t /* http: ...
- 【BZOJ-3122】随机数生成器 BSGS
3122: [Sdoi2013]随机数生成器 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 1362 Solved: 531[Submit][Sta ...
- 【BZOJ 3122】 [Sdoi2013]随机数生成器 (BSGS)
3122: [Sdoi2013]随机数生成器 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 1442 Solved: 552 Description ...
- 洛咕 P3306 [SDOI2013]随机数生成器
洛咕 P3306 [SDOI2013]随机数生成器 大力推式子??? \(X_{i}=\underbrace{a(a(\cdots(a(a}_{i-1个a}X_1+b)))\cdots)\) \(=b ...
随机推荐
- centos中安装虚拟机
我这里选用的虚拟机软件为virtualbox 下载wget http://download.virtualbox.org/virtualbox/5.2.2/VirtualBox-5.2-5.2.2_1 ...
- leetcode — remove-duplicates-from-sorted-array
import java.util.Arrays; /** * Source : https://oj.leetcode.com/problems/remove-duplicates-from-sort ...
- NET快速信息化系统开发框架 V3.2 -> “用户管理”主界面使用多表头展示、增加打印功能
RDIFrameowrk.NET 用户管理是使用非常频繁的模块,由于需要展示的字段比较多,以前的展示方式显得不是太规范,现3.2版本用户管理主界面进行了全新的设计,数据列表展示使用了Dev家族全新的G ...
- python实战学习之matplotlib绘图续
学习完matplotlib绘图可以设置的属性,还需要学习一下除了折线图以外其他类型的图如直方图,条形图,散点图等,matplotlib还支持更多的图,具体细节可以参考官方文档:https://matp ...
- OpenCV在C#中应用—OpenCVSharp
1.什么是OpenCVSharp 之前一直是基于OpenCV开发视觉算法,但C++语言对于GUI的开发相对于C#来说确实很不方便,之前就了解到C#下使用OpenCV可以使用EmguCV,这段时间 ...
- spark问题
使用IDEA运行spark程序,除了需要导入spark的一些依赖包之外,还需要注意的是 当启动spark报找不到可执行的hadoop winutils.exe 可已下载相应版本的winutils.ex ...
- Java开发笔记(十二)布尔变量论道与或非
在编程语言的设计之初,它们除了可以进行数学计算,还常常用于逻辑推理和条件判断.为了实现逻辑判断的功能,Java引入了一种布尔类型boolean,用来表示“真”和“假”.该类型的变量只允许两个取值,即t ...
- Android开发过程中的坑及解决方法收录(六)
1. file.listFiles 空指针异常 最近在弄个小项目,类似一个文件管理器,需要获得手机存储里的目录之后显示,但是运行过程中出现错误,搜索了资料,得出了以下的解决办法 问题产生的原因: an ...
- Serialization
title: Serialization date: 2018-03-26 15:18:26 tags: [JAVA,Read] categories: other --- 概述 程序运行时,变量的内 ...
- Java Socket网络编程学习笔记(一)
0.前言 其实大概半年前就已经看过网络编程Socket的知识了(传统IO),但是因为长时间的不使用导致忘的一干二净,最近正好准备校招,又重新看了网络编程这一章, 是传统IO(BIO)相关的内容,故在此 ...