POJ 3126 Prime Path

BFS。

一个质数至少要变多少次才能到另外一个质数。

中间变得数也必须是质数。每次变换一个位置。

Sample Input

3
1033 8179
1373 8017
1033 1033

Sample Output

6
7
0

记得判断找出来的数是否是质数，最好打印1000——10000的质数表（很多办法）。

然后记得标记数是否被搜索到。

#include<stdio.h>
#include<string.h>
#include<math.h>
int a,b,ans,top,end;
short q[10001];
short flag[10001];
short prime[10001];
int isprime()
{
    int i,j,k;
    for(i=1001;i<10000;i+=2)
    {
        k=sqrt(i*1.0);
        for(j=2;j<=k;j++)
        if(i%j==0)
            break;
        if(j>=k+1)
            flag[i]=1;
    }
}
int bfs()
{
    int i,j;
    q[end]=a;
    while(1)
    {
        short ok=0;
        for(i=end;i<=top;i++)
        {
            if(q[i]==b)return 1;
            prime[q[i]]=1;
        }
        for(i--;end<=i;end++)
        {
            int t;
            for(j=1;j<=9;j++)
            {
                t=q[end]%1000+j*1000;
                if(flag[t]&&!prime[t])
                q[++top]=t,ok=1;
            }
            for(j=0;j<=9;j++)
            {
                t=q[end]/1000*1000+q[end]%100+j*100;
                if(flag[t]&&!prime[t])
                q[++top]=t,ok=1;
            }
            for(j=0;j<=9;j++)
            {
                t=q[end]/100*100+q[end]%10+j*10;
                if(flag[t]&&!prime[t])
                q[++top]=t,ok=1;
            }
            for(j=0;j<=9;j++)
            {
                t=q[end]/10*10+j;
                if(flag[t]&&!prime[t])
                q[++top]=t,ok=1;
            }
        }
        if(!ok)return 0;
        ans++;
    }
}
int main()
{
    int n;
    memset(flag,0,sizeof(flag));
    isprime();
    while(~scanf("%d",&n))
    {
        while(n-->0)
        {
            memset(prime,0,sizeof(prime));
            scanf("%d%d",&a,&b);
            ans=end=0,top=1;
            if(bfs())
            printf("%d\n",ans);
            else
            printf("Impossible\n");
        }
    }
}