POJ-3126.PrimePath(欧拉筛素数打表 + BFS)

给出一篇有关素数线性筛和区间筛的博客，有兴趣的读者可以自取。

　　

　　本题大意：

　　　　给定两个四位的素数，没有前导零，每次变换其中的一位，最终使得两个素数相等，输出最小变换次数。要求变换过程中的数也都是素数。

　　本题思路：

　　　　求最小步数，用BFS，判断素数，用欧拉筛打表。

　　参考代码：

 #include <cstdio>
 #include <queue>
 #include <map>
 #include <cstring>
 using namespace std;

 int star, finish, ans, total = ;
 const int maxn =  + ;
 bool isprime[maxn];
 bool isvisited[maxn];
 int prime[];

 struct node{
     int cur, step;
 } now, Next;

 void get_prime() {//欧拉筛打表
     memset(isprime, false, sizeof(isprime));
     for(int i  = ; i < maxn; i ++) {
         if(!isprime[i])
             prime[total ++] = i;
         for(int j = ; j < total; j ++) {
             if(prime[j] * i > maxn) break;
             isprime[prime[j] * i] = true;
             if(i % prime[j] == )   break;
         }
     }
 }

 int bfs() {
     queue <node> Q;
     isvisited[star] = true;
     now.cur = star;
     now.step = ;
     Q.push(now);
     while(!Q.empty()) {
         char num[];
         now = Q.front();
         Q.pop();
         if(now.cur == finish)    return now.step;
         for(int i = ; i < ; i ++) {
             sprintf(num, "%d", now.cur);
             for(int j = ; j < ; j ++) {
                 if(j ==  && i == )    continue;
                 if(i == )
                     Next.cur = j *  + (num[] - '') *  + (num[] - '') *  + (num[] - '');
                 else if(i == )
                     Next.cur = j *  + (num[] - '') *  +(num[] - '') *  + (num[] - '');
                 else if(i == )
                     Next.cur = j *  + (num[] - '') *  + (num[] - '') *  + (num[] - '');
                 else if(i == )
                     Next.cur = j + (num[] - '') *  + (num[] - '') *  + (num[] - '') * ;
                 if(!isprime[Next.cur] && !isvisited[Next.cur]) {
                     Next.step = now.step + ;
                     isvisited[Next.cur] = ;
                     Q.push(Next);
                 }
             }
         }
     }
     return -;
 }

 int main () {
     memset(isvisited, false, sizeof(isvisited));
     get_prime();
     int t;
     scanf("%d", &t);
     while(t --) {
         memset(isvisited, false, sizeof(isvisited));
         scanf("%d %d", &star, &finish);
         ans = bfs();
         if(ans == -)   printf("Impossible\n");
         else    printf("%d\n", ans);
     }
     return ;
 }