poj 3126 prime path 简单广搜

  1 #include <algorithm>
  2 #include <iostream>
  3 #include <sstream>
  4 #include <cstdlib>
  5 #include <cstring>
  6 #include <cstdio>
  7 #include <string>
  8 #include <bitset>
  9 #include <vector>
 10 #include <queue>
 11 #include <stack>
 12 #include <cmath>
 13 #include <list>
 14 #include <map>
 15 #include <set>
 16 using namespace std;
 17 /*10^8-----1s*/
 18 /***************************************/
 19 typedef vector<int> VI;
 20 typedef vector<char> VC;
 21 typedef vector<string> VS;
 22 typedef set<int> SI;
 23 typedef set<string> SS;
 24 typedef map<int ,int> MII;
 25 typedef map<string,int> MSI;
 26 typedef pair<int,int> PII;
 27 typedef vector<PII> VII;
 28 typedef vector<VI > VVI;
 29 /***************************************/
 30 #define min(a,b) (a>b?b:a)
 31 #define max(a,b) (a>b?a:b)
 32 
 33 #define clr(a,b) memset(a,b,sizeof(a))
 34 #define all(x)    (x).begin(), (x).end()
 35 #define sz(x) ((int)(x).size())
 36 #define ll long long
 37 #define int64 __int64
 38 #define pb push_back
 39 #define mp make_pair
 40 #define LL(x) ((x)<<1)
 41 #define RR(x) ((x)<<1|1)
 42 #define ri(x) scanf("%d",&x)
 43 #define rii(x,y) scanf("%d%d",&x,&y)
 44 #define rd(x) scanf("%lf",&x)
 45 #define rdd(x,y) scanf("%lf%lf",&x,&y)
 46 #define rs(x) scanf("%s",x)
 47 #define pi(x) printf("%d",x)
 48 #define pin(x) printf("%d\n",x)
 49 #define ps(x) printf("%s",x)
 50 #define pn()  printf("\n")
 51 #define sqr(x) ((x)*(x))
 52 #define rep(i,a,b)  for(int i=(a);i<(b);i++)
 53 #define repu(i,a,b) for(int i=(a);i<=(b);i++)
 54 #define repd(i,a,b) for(int i=(a);i>=(b);i--)
 55 #define repc(i,a,c) for(int i=(a);(c);i++)
 56 /***************************************/
 57 const int INF = 0x7f7f7f7f;
 58 const double eps = 1e-8;
 59 const double PIE=acos(-1.0);
 60 const int dx[]= {0,-1,0,1};
 61 const int dy[]= {1,0,-1,0};
 62 const int fx[]= {-1,-1,-1,0,0,1,1,1};
 63 const int fy[]= {-1,0,1,-1,1,-1,0,1};
 64 /***************************************/
 65 void openfile()
 66 {
 67     freopen("data.in","rb",stdin);
 68     freopen("data.out","wb",stdout);
 69 }
 70 /**********************The End OF The Template*****************/
 71 
 72 int primes[10004];
 73 int visit[10004];//标记 同时 记录转换的次数
 74 int a,b;
 75 int w[5];
 76 int flag;
 77 
 78 
 79 /*打表  筛选法求素数（2是最小的素数）*/
 80 void isprime()
 81 {
 82     int i,j;
 83     primes[0]=0;//等于0表示非素数
 84     primes[1]=0;
 85     for(i=2; i<=10002; i++)
 86         primes[i]=1;
 87     for(i=2; i<=10002; i++)
 88     {
 89         if(primes[i])
 90         {
 91             for(j=i+i; j<=10002; j+=i)
 92                 primes[j]=0;
 93         }
 94     }
 95 }
 96 
 97 /*广搜找最小需要的转换次数*/
 98 int bfs(int x)
 99 {
100     int xx,k,t,i,j,l;
101     queue<int >q;
102     q.push(x);
103     flag=0;
104     while(!q.empty())
105     {
106         xx=q.front();
107         q.pop();
108         if(xx==b)
109         {
110             flag=1;
111             return visit[xx];
112         }
113         w[1]=xx%10;
114         w[2]=(xx/10)%10;
115         w[3]=(xx/100)%10;
116         w[4]=xx/1000;
117 
118         for(i=1; i<=4; i++) // 1 个 2 十
119             for(j=0; j<10; j++)
120             {
121                 if(i==4&&j==0)//最高位不能为0
122                     continue;
123                 if(i==1&&j%2==0)//个位为偶数跳过
124                     continue;
125                 k=0;
126                 t=w[i];
127                 w[i]=j;
128                 for(l=4; l>0; l--)
129                 {
130                     k=k*10+w[l];
131                 }
132                 w[i]=t;
133 
134                 if(primes[k]&&!visit[k])
135                 {
136                     visit[k]=visit[xx]+1;
137                     q.push(k);
138                 }
139 
140             }
141     }
142 }
143 
144 
145 
146 int main()
147 {
148     int t;
149     isprime();
150     scanf("%d",&t);
151     while(t--)
152     {
153         memset(visit,0,sizeof(visit));
154         scanf("%d %d",&a,&b);
155         visit[a]=1;
156         int m=bfs(a);
157         if(flag)
158             printf("%d\n",m-1);
159         else
160             printf("Impossible\n");
161     }
162     return 0;
163 }