POJ3126 Prime Path【数论】【BFS】
2015-10-31 18:40
375 查看
题目链接:
http://poj.org/problem?id=3126
题目大意:
给你两个有四位数字的素数 A、B,问:每次只改变一个数字,且改变前后的数都是
素数,那么从 A 变到 B,最少需要多少次。
解题思路:
用 BFS 来做。判断素数用筛法求素数打表预处理一下,不过注意 1000 以下的数要
当非素数看待。
每次改变一位数字,并且如果改变后的数仍为素数的话入队,并用 cnt[] 来记录步数。
直到得到目标数的时候,返回结果。
AC代码:
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
const int MAXN = 10100;
bool Prime[MAXN];
void IsPrime()
{
Prime[0] = Prime[1] = false;
for(int i = 2; i < MAXN; ++i)
Prime[i] = true;
for(int i = 2; i < MAXN; ++i)
{
if(Prime[i])
{
for(int j = i+i; j < MAXN; j+=i)
Prime[j] = false;
}
}
for(int i = 2; i < 1000; ++i)
Prime[i] = false;
}
bool vis[MAXN];
int t[5],cnt[MAXN];
int Bfs(int S,int T)
{
queue<int> q;
memset(vis,false,sizeof(vis));
memset(cnt,0,sizeof(cnt));
q.push(S);
vis[S] = true;
while(!q.empty())
{
int v = q.front();
q.pop();
//存取各个位上数字
t[0] = v/1000;
t[1] = v/100%10;
t[2] = v/10%10;
t[3] = v%10; t[0]=v/1000;
// printf("%d %d %d %d\n",t[0],t[1],t[2],t[3]);
for(int j = 0; j < 4; ++j)
{
int temp = t[j]; //记录第 j 个数字
for(int i = 0; i <= 9; ++i) //将 t[j] 改变为 i
{
if(i != temp)
{
t[j] = i;
int vtemp = t[0]*1000 + t[1]*100 + t[2]*10 + t[3];
if(!vis[vtemp] && Prime[vtemp])
{
cnt[vtemp] = cnt[v] + 1;
vis[vtemp] = true;
q.push(vtemp);
}
if(vtemp == T)
return cnt[vtemp];
}
}
t[j] = temp; //变回原本数字,继续改变另外数字
}
if(v == T)
return cnt[v];
}
return -1;
}
int main()
{
int T,A,B;
IsPrime();
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&A,&B);
int Ans = Bfs(A,B);
if(Ans == -1)
printf("Impossible\n");
else
printf("%d\n",Ans);
}
return 0;
}
http://poj.org/problem?id=3126
题目大意:
给你两个有四位数字的素数 A、B,问:每次只改变一个数字,且改变前后的数都是
素数,那么从 A 变到 B,最少需要多少次。
解题思路:
用 BFS 来做。判断素数用筛法求素数打表预处理一下,不过注意 1000 以下的数要
当非素数看待。
每次改变一位数字,并且如果改变后的数仍为素数的话入队,并用 cnt[] 来记录步数。
直到得到目标数的时候,返回结果。
AC代码:
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
const int MAXN = 10100;
bool Prime[MAXN];
void IsPrime()
{
Prime[0] = Prime[1] = false;
for(int i = 2; i < MAXN; ++i)
Prime[i] = true;
for(int i = 2; i < MAXN; ++i)
{
if(Prime[i])
{
for(int j = i+i; j < MAXN; j+=i)
Prime[j] = false;
}
}
for(int i = 2; i < 1000; ++i)
Prime[i] = false;
}
bool vis[MAXN];
int t[5],cnt[MAXN];
int Bfs(int S,int T)
{
queue<int> q;
memset(vis,false,sizeof(vis));
memset(cnt,0,sizeof(cnt));
q.push(S);
vis[S] = true;
while(!q.empty())
{
int v = q.front();
q.pop();
//存取各个位上数字
t[0] = v/1000;
t[1] = v/100%10;
t[2] = v/10%10;
t[3] = v%10; t[0]=v/1000;
// printf("%d %d %d %d\n",t[0],t[1],t[2],t[3]);
for(int j = 0; j < 4; ++j)
{
int temp = t[j]; //记录第 j 个数字
for(int i = 0; i <= 9; ++i) //将 t[j] 改变为 i
{
if(i != temp)
{
t[j] = i;
int vtemp = t[0]*1000 + t[1]*100 + t[2]*10 + t[3];
if(!vis[vtemp] && Prime[vtemp])
{
cnt[vtemp] = cnt[v] + 1;
vis[vtemp] = true;
q.push(vtemp);
}
if(vtemp == T)
return cnt[vtemp];
}
}
t[j] = temp; //变回原本数字,继续改变另外数字
}
if(v == T)
return cnt[v];
}
return -1;
}
int main()
{
int T,A,B;
IsPrime();
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&A,&B);
int Ans = Bfs(A,B);
if(Ans == -1)
printf("Impossible\n");
else
printf("%d\n",Ans);
}
return 0;
}
相关文章推荐
- 面向对象的4要素:抽象,封装,继承,多态
- CodeForces 337C 找规律 ,等比数列求和
- Python 创建XML
- weka –Apriori算法 关联规则挖掘实验
- jquery技巧总结
- javascript常用数组算法总结
- 【HDU5521 2015沈阳赛区M】【拆点最短路 dijkstra+heap】Meeting 集合内距离相同
- c++ list的用法(整理)
- 蓝桥杯——01字串
- Android Studio自定义组合控件
- 安装两次ORACLE客户端,PL-SQL登录报错 ORA-12154
- 使用Weka进行数据挖掘
- 设计模式之简单工厂模式
- linux awk命令详解
- Android手机适配
- Android OpenGL ES 开发教程 从入门到精通
- 四元数转欧拉角
- java学习笔记.04——线程池
- AndroidWheel 省市区三级联动(高仿IOS版本)
- iOS-CALayer && CAAnimation