http://poj.org/problem?id=3126
题意：给你两个四位数a,b，从a开始 每次只能改变上一次数的其中一位，问至少需要几步才能得到b
分析：求最小路 典型的广搜 表面上是 40入口的bfs 但是除去有的数不是素数 入口数远小于40
可以写一个 判断一个数是否为素数的函数 ， 每次去 调用 判断一个数是否要进队列
也可以 事先打一个素数表 这样会快点
注意：output ：either with a number stating the minimal cost or containing the word Impossible
用c++才交对的 g++ not ok i am confused
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 }
时间： 2024-12-26 23:15:04