题目链接(洛谷)

题目大意

给定两个数 \(u\) , \(v\) 。有三种操作:

  1. \(u=u+1(mod\) \(p)\) 。
  2. \(u=u+p−1(mod\) \(p)\) 。
  3. \(u=u^{p−2}(mod\) \(p)\) 。

思路

BFS

状态太多导致队列装不下。

迭代加深

\(TLE\) ,浪费了太多时间在每一层上,没有记忆化且状态很多。

IDA*

不行,无法得出乐观股价函数。

双向BFS

这样会将步数很为两半,状态相较于普通的 \(BFS\) 会少很多。

先来看操作一和操作二,他们的关系是可以互逆的。一个对于原数 \(+1\) ,另一个对于原数 \(-1\) 。

操作三和操作三是互逆的,由费马小定理可知:若 \(p\) 为质数,则 \(a^{p-1}≡1(mod\) \(p)\)。

可得出:\((u^{p-2})^{p-2}≡u^{(p-2)(p-2)}≡u^{(p-1)(p-3)+1}≡(u^{p-1})^{p-3}u≡u(mod\) \(p)\)

那么就分别由开始状态与结束状态来向中间推进。

Code

#include <map>
#include <queue>
#include <cstdio>
#include <iostream>
using namespace std;
struct Status {//队列中保存的状态
int step, number, flag;//分别是:步数,当前状态的数,正向或者反向
Status() {}
Status(int S, int N, int F) {
step = S;
number = N;
flag = F;
}
};
const int MAXN = 1e6 + 5;
queue<Status> q;
map<int, int> real;
bool vis[2][MAXN];//是否访问过
int dis[2][MAXN];//步数
pair<int, int> pre[2][MAXN];//first记录前一个数的哈希值,second记录操作的序号
int u, v, p;
int tot;
int Quick_Pow(int fpx, int fpy) {//快速幂
long long res = 1;
long long x = fpx;
long long y = fpy;
while(y) {
if(y & 1)
res = (res * x) % p;
x = (x * x) % p;
y >>= 1;
}
int ans = res;
return ans;
}
int Get_Hash(int x) {//map映射假哈希
map<int, int>::iterator it = real.find(x);
if(it != real.end())
return (*it).second;
real[x] = ++tot;
return tot;
}
void Print(int x, int y) {//输出路径:记录每个前缀
if(y == -1)
return;
if(!x) {//前半部分倒着输出
if(pre[x][y].first != -1) {
Print(x, pre[x][y].first);
printf("%d ", pre[x][y].second);
}
}
else {//后半部分正着输出
if(pre[x][y].first != -1) {
printf("%d ", pre[x][y].second);
Print(x, pre[x][y].first);
}
}
}
void DoubleBfs() {
int tmp;
q.push(Status(0, u, 0));//初始化两个状态
q.push(Status(0, v, 1));
tmp = Get_Hash(u);
vis[0][tmp] = 1;
pre[0][tmp].first = -1;
tmp = Get_Hash(v);
vis[1][tmp] = 1;
pre[1][tmp].first = -1;
while(!q.empty()) {
Status now = q.front();
q.pop();
int skt = Get_Hash(now.number);
if(vis[!now.flag][skt]) {//碰头了输出并跳出
printf("%d\n", dis[!now.flag][skt] + dis[now.flag][skt]);
if(pre[0][skt].first != -1) {
Print(0, pre[0][skt].first);
printf("%d ", pre[0][skt].second);
}
if(pre[1][skt].first != -1) {
printf("%d ", pre[1][skt].second);
Print(1, pre[1][skt].first);
}
return;
}
Status next = now;
next.step++;
next.number = (next.number + 1) % p;
tmp = Get_Hash(next.number);
if(!vis[now.flag][tmp]) {//没有被访问则访问
vis[now.flag][tmp] = 1;
dis[now.flag][tmp] = next.step;
pre[now.flag][tmp].first = skt;
if(now.flag)
pre[now.flag][tmp].second = 2;//若是倒着的,则该操作为1
else
pre[now.flag][tmp].second = 1;//若是正着的,则该操作为2
q.push(next);
}
next = now;
next.step++;
next.number = (next.number + p - 1) % p;
tmp = Get_Hash(next.number);
if(!vis[now.flag][tmp]) {//同上
vis[now.flag][tmp] = 1;
dis[now.flag][tmp] = next.step;
pre[now.flag][tmp].first = skt;
if(now.flag)
pre[now.flag][tmp].second = 1;
else
pre[now.flag][tmp].second = 2;
q.push(next);
}
next = now;
next.step++;
next.number = Quick_Pow(next.number, p - 2) % p;
tmp = Get_Hash(next.number);
if(!vis[now.flag][tmp]) {//同上
vis[now.flag][tmp] = 1;
dis[now.flag][tmp] = next.step;
pre[now.flag][tmp].first = skt;
pre[now.flag][tmp].second = 3;//自己的逆操作就是自己
q.push(next);
}
}
}
int main() {
scanf("%d %d %d", &u, &v, &p);
DoubleBfs();
return 0;
}

CF995E Number Clicker (双向BFS)的更多相关文章

  1. CF995E Number Clicker 解题报告

    CF995E Number Clicker 题目描述 Allen is playing Number Clicker on his phone. He starts with an integer u ...

  2. CF995E Number Clicker

    题目分析 首先,我们必须明白,操作都是互逆的,\(1,2\)之间是可以互相转化的,这是不需证明的,对于操作\(3\),实际上,是求当前数的逆元,我们知道,逆元就是求当前数在模另一个数下的倒数,那么,逆 ...

  3. Number Clicker CodeForces - 995E(双向bfs)

    双向bfs  注意数很大  用map来存 然后各种难受....

  4. CodeForces - 995E Number Clicker (双向BFS)

    题意:给出u,v,p,对u可以进行三种变化: 1.u=(u+1)%p ; 2.u = (u+p-1)%p;  3.u = 模p下的逆元.问通过几步可以使u变成v,并且给出每一步的操作. 分析:朴素的b ...

  5. HDU 3085 Nightmare Ⅱ (双向BFS)

    Nightmare Ⅱ Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Tota ...

  6. Hdu1401-Solitaire(双向bfs)

    Solitaire is a game played on a chessboard 8x8. The rows and columns of the chessboard are numbered ...

  7. UVA1601-The Morning after Halloween(双向BFS)

    Problem UVA1601-The Morning after Halloween Accept: 289 Submit: 3136 Time Limit: 12000 mSec  Problem ...

  8. Eight (HDU - 1043|POJ - 1077)(A* | 双向bfs+康拓展开)

    The 15-puzzle has been around for over 100 years; even if you don't know it by that name, you've see ...

  9. HDU3085(双向BFS+曼哈顿距离)题解

    Nightmare Ⅱ Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Tota ...

随机推荐

  1. Selenium刚玩一会儿,就感受了私人秘书的体验

    学习python的过程中,少不了接触第三方库,毕竟作为胶水语言python的强大之处也就是第三方库体量庞大,无疑大大增强了python的战斗力. 有时候想完成网页自动化操作,这时候Selenium进入 ...

  2. spring再学习之设计模式

    今天我们来聊一聊,spring中常用到的设计模式,在spring中常用的设计模式达到九种. 第一种:简单工厂 三种工厂模式:https://blog.csdn.net/xiaoddt/article/ ...

  3. Linux系统启动过程内核文件丢失解决方法

    一.问题描述 公司近期因机房断电,导致服务器重启后,引导进入不了操作系统.经过检查发现启动文件缺失,导致系统启动失败,网上搜了好多资料,解决都比较零散,现结合实际处理经验和网友的建议整理接方案. 二. ...

  4. python3 anaconda 安装pyhook3 pythoncom(pywin32)

    为什么不安装pyhook 1.pyhook不支持python3 2.网络上有一些方法下载pyhook的whl然后pip安装到python3,可以运行,但是会因为编码问题导致移动到窗口标题含有非ASCI ...

  5. ysoserial-URLDNS学习

    简述 ysoserial很强大,花时间好好研究研究其中的利用链对于了解java语言的一些特性很有帮助,也方便打好学习java安全的基础,刚学反序列化时就分析过commoncollections,但是是 ...

  6. Deep Learning Specialization 笔记

    1. numpy中的几种矩阵相乘: # x1: axn, x2:nxb np.dot(x1, x2): axn * nxb np.outer(x1, x2): nx1*1xn # 实质为: np.ra ...

  7. windows10 浏览器跑分对比!

    2015-12-12 windows10 浏览器跑分对比! YOUR BROWSER SCORES!                MaxScore=555http://html5test.com/i ...

  8. cnblogs & 502 Bad Gateway

    cnblogs & 502 Bad Gateway 博客园 502 Bad Gateway 服务器发生了一些错误,请联系 contact@cnblogs.com 可以查看,不可以编辑 HTTP ...

  9. position: absolute; not work

    position: absolute; not work https://stackoverflow.com/questions/11928294/css-position-absolute-with ...

  10. Flutter: moor_flutter库,简化sqlite操作

    入门 video moor_flutter 示例代码仓库 install dependencies: ... moor_flutter: dev_dependencies: ... moor_gene ...