POJ 3126 Prime Path（筛法，双向搜索）

题意：一个4位的素数每次变动一个数位，中间过程也要上素数，问变成另一个的最小步数。

线性筛一遍以后bfs就好。我写的双向，其实没有必要。

#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<queue>
#include<vector>
#include<stack>
#include<vector>
#include<map>
#include<set>
#include<algorithm>
#include<cmath>
//#include<bits/stdc++.h>
using namespace std;

typedef long long ll;

const int maxn = 1e5;
bool isNot[maxn];

void seive(int n = maxn-)
{
    int m = sqrt(n+0.5);
    for(int i = ; i <= m; i++){
        if(!isNot[i])
        for(int j = i*i; j <= n; j += i){
            isNot[j] = true;
        }
    }
    isNot[] = isNot[] = true;
}

int d[maxn];
int vis[][maxn],clk;
#define PB push_back
#define PS push
const int wei[] = {,,,};

void doubleBfs(int s,int t)
{
    if(s == t) { puts(""); return; }
    queue<int> Q[];
    queue<int> *q1 = Q, *q2 = Q+, *nxt = Q+;
    int *v1 = *vis, *v2 = vis[];
    q1->PS(s); q2->PS(t);
    v1[s] = ++clk; v2[t] = clk;
    d[s] = ; d[t] = ;
    while(q1->size() && q2->size()){
        if(q1->size()>q2->size()){
            swap(q1,q2); swap(v1,v2);
        }
        while(q1->size()){
            int u = q1->front(); q1->pop();
            int dig[];
            for(int i = ; i < ; i++){
                dig[i] = u/wei[i]%;
            }
            for(int i = ; i < ; i++){
                int t = u - dig[i]*wei[i];
                for(int j = i == ; j < ; j++){
                    if(j != dig[i]){
                        int v = t+j*wei[i];
                        if(!isNot[v] && v1[v] != clk){
                            if(v2[v] == clk) {
                                printf("%d\n",d[v]+d[u]+);
                                return;
                            }
                            d[v] = d[u]+;
                            nxt->PS(v);
                            v1[v] = clk;
                        }
                    }
                }
            }
        }
        swap(q1,nxt);
    }
    puts("Impossible");
}

//#define LOCAL
int main()
{
#ifdef LOCAL
    freopen("in.txt","r",stdin);
#endif
    seive();
    int T; scanf("%d",&T);
    while(T--){
        int a,b; scanf("%d%d",&a,&b);
        doubleBfs(a,b);
    }
    return ;
}