SPOJ Problem Set (classical) 11582. A Famous Grid

//螺旋数，很多时候不是你不会，是对未知的东西充满了恐惧//其实静下心来，这个螺旋数是完全可以写出来的，剩下的就是简单的bfs了//PS：图是无限大的，虽然x,y给出的范围是10000，也就是100*100的图，但是可以绕着走，所以要稍微开大一点儿#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>#include<queue>using namespace std;const int maxn=104;int a,b;int map[maxn+10][maxn+10],tt;struct in{int x,y;int cnt;}p[maxn*maxn+100],que[maxn*maxn+100],start,p0,p1;int prim[maxn*maxn+100]={0};int vis[maxn+10][maxn+10]={0};int path[maxn+10][maxn+10];void deal()//螺旋数的构造，从外向内，思维很直观，所以思考很重要，不能因为陌生而放弃 {  int l,r,u,d,now=maxn*maxn;  l=1,r=maxn,u=1,d=maxn;  while(now>0)  {      for(int i=l;i<=r;i++)      {          map[u][i]=now;          p[now].x=u;          p[now].y=i;          now--;      }      u++;      for(int i=u;i<=d;i++)      {          map[i][r]=now;          p[now].x=i;          p[now].y=r;          now--;      }      r--;      for(int i=r;i>=l;i--)      {          map[d][i]=now;          p[now].x=d;          p[now].y=i;          now--;      }      d--;      for(int i=d;i>=u;i--)      {          map[i][l]=now;          p[now].x=i;          p[now].y=l;          now--;      }      l++;  }}void deal2(){int i,j,n=maxn*maxn;prim[1]=1;for(i=2;i<=n;i++)if(prim[i]==0){for(j=i+i;j<=n;j+=i)prim[j]=1;}}int findx[][2]={{0,-1},{0,1},{-1,0},{1,0}};int f;int bfs(int t,int s){    if(prim[t]==0||prim[s]==0)    return -1;    else if(s==t)    return 0;    else    {        queue <in> que;        que.push(start);        while(!que.empty())        {            p0=que.front();            que.pop();            for(int i=0;i<4;i++)            {            //cout<<tt<<endl;                int nextx=p0.x+findx[i][0];                int nexty=p0.y+findx[i][1];                int xx=map[nextx][nexty];                if(nextx<1 || nextx>maxn || nexty<1 || nexty>maxn) continue;  if(vis[nextx][nexty]==1) continue;   if(prim[xx]==0) continue;     vis[nextx][nexty]=1;                    p1.x=nextx;                    p1.y=nexty;                    p1.cnt=p0.cnt+1;                    if(p1.x==a&&p1.y==b)                    {                    f=1;                        return p1.cnt;                    }                    else                    {                                       que.push(p1);                    }                               }      }//cout<<f<<endl;        if(f==0)        return -1;     }                }int main(){int s,t,sum;tt=1;deal();deal2();while(~scanf("%d%d",&s,&t)){memset(vis,0,sizeof(vis));//每次初始化vis，以免之前的数据影响以后的，如果想省下初始化的时间，可以将vis每次的赋值为tt printf("Case %d: ",tt++);start.x=p[s].x;        start.y=p[s].y;        start.cnt=0;        vis[start.x][start.y]=1;         a=p[t].x,b=p[t].y;         f=0;sum=bfs(t,s);if(sum==-1) puts("impossible");else printf("%d\n",sum);}return 0;}