poj3126prime path

最开始就跳过的题居然是最简单的==而且是1A!!!太激动了有木有!而且这个题是这次周赛里通过的最少的!!  光辉的印记

简单说一下思路,咳咳咳……

正常搜索就是一个二维或是三维的空间里用二重循环或是三重循环是吧,这个题问的是给已知一个1000-9999的素数让你每次改变一个数字变成一个素数,再变一个,再变一个……中间过程要求都是1000-9999的素数,最后导到最终的素数==问你需要几步,不就是把四个方向换成0-9嘛

#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
bool prime[100006];
int p[10000],k=0;
bool vis[100000];
int t,m,n;
struct node
{
    int x,step;
}q[100000];
void isprime()
{
    long long i,j;
    k=0;
    memset(prime,true,sizeof(prime));
    for(i=2;i<=100000;i++)
    {
        if(prime[i])
        {
            p[k++]=i;
            for(j=i*i;j<100000;j+=i)
            prime[j]=false;
        }
    }
    //for(int i=999;i<=9999;i++)if(prime[i]) cout<<i<<" ";
}
int fun(int n)
{
    int tmp=1;
    while(n--) tmp*=10;
    return tmp;
}
int bfs(int sx)
{
    int l=0,r=0;
    q[r].x=sx,q[r++].step=0;
    vis[sx]=1;
    while(l<r)
    {
        node cur=q[l++],next;
        for(int i=0;i<4;i++) for(int j=0;j<=9;j++)
        {
            if((cur.x/fun(i))%10==j) continue;
            if(i==3&&j==0) continue;
            next.x=cur.x-(cur.x/fun(i))%10*fun(i)+j*fun(i);
            next.step=cur.step+1;
            if(prime[next.x]&&!vis[next.x])
            {
                if(next.x==m) return next.step;
                vis[next.x]=1;
                q[r++]=next;
            }
        }
    }
}
int main()
{
    //freopen("cin.txt","r",stdin);
    isprime();
    while(~scanf("%d",&t))
    {
        while(t--)
        {
            scanf("%d%d",&n,&m);
            if(n==m) {cout<<0<<endl;continue;}
            memset(vis,0,sizeof(vis));
            cout<<bfs(n)<<endl;
        }
    }
    return 0;
}



你可能感兴趣的:(算法,数论,ACM,poj,广搜)