Prime Path-POJ3126-1尤拉篩選法+BFS-好題
阿新 • • 發佈:2018-12-13
題意: 給你兩個四位數的質數m,n(m<=n),每次變化只能改變當前數的一個數字,並且過程中的數字全是四位數的質數,求最小的步數。 思路: 先尤拉篩選法求出素數,然後從m開始BFS,每次只改變一個數位,如何只改變一個數位(這個很關鍵),這裡用一個dir陣列,dir[4] = {1,10,100,1000},具體看程式碼註釋理解。
AC code:
#include<iostream>
#include<cstdio>
#include<vector>
#include<bitset>
#include<stack>
#include<set>
#include<queue>
#include<map>
#include<cmath>
#include<string>
#include<cstring>
#include<ctime>
#include<fstream>
#include<cstdlib>
#include<algorithm>
using namespace std;
#define pii pair<int, int>
#define pb push_back
#define mem(a,b) memset(a,b,sizeof(a))
#define per(i,a,b) for(int i=a;i<=b;i++)
#define rep(i,a,b) for(int i=a;i>=b;i--)
#define all(x) x.begin(),x.end()
#define PER(i,x) for(auto i=x.begin();i!=x.end();i++)
#define PI acos(-1.0)
#define INF 0x3f3f3f3f
typedef long long LL;
const double eps=1.0e-5;
const int maxn=1e4 + 10;
int dx[4] = {0,0,-1,1};
int dy[4] = {1,-1,0,0};
int dir[4] = {1,10,100,1000};
int n = 0,m = 0;
int prime[maxn];
bool is_prime[maxn];
bool vis[maxn];
struct node{
int num,step;
node(int numm,int ss):num(numm),step(ss){
}
node(){
}
};
int ouler(int n){
fill(is_prime,is_prime+maxn,true);
is_prime[0] = is_prime[1] = false;
int cnt = 0;
per(i,2,n){
if(is_prime[i] == true){
prime[cnt++] = i;
}
for(int j = 0;j < cnt && i * prime[j] <= n;++j){
is_prime[i*prime[j]] = false;
if(i % prime[j] == 0){
break;
}
}
}
return cnt;
}
void bfs(int x,int y){
fill(vis,vis+maxn,false);
queue<node> que;
que.push(node(x,0));
vis[x] = true;
while(!que.empty()){
node u = que.front();
que.pop();
if(u.num == y){
printf("%d\n",u.step);
return ;
}
per(k,0,3){
int ty = u.num - u.num / dir[k] % 10 * dir[k];//精華所在
per(i,0,9){
if(k == 3 && i == 0){//首位不能是0
continue;
}
if(k == 0 && i % 2 == 0){//末尾不能是0
continue;
}
if(u.num / dir[k] % 10 == i){//相等
continue;
}
int tx = ty + i * dir[k];
if(vis[tx] == false && is_prime[tx]){//因為有前面的剪枝,所以這裡不需要
que.push(node(tx,u.step + 1));//再判斷 tx >= 1000 && tx <= 9999
vis[tx] = true;
}
}
}
}
printf("Impossible\n");
}
int main(){
#ifndef ONLINE_JUDGE
//freopen("a.txt","r",stdin);
#endif
ouler(10000);
int T = 0;
scanf("%d",&T);
while(T--){
scanf("%d %d",&m,&n);
bfs(m,n);
}
return 0;
}