2017-10-24离线赛总结

失分小结：
估分：260
实际分数：160
第一题爆零。。。原因是题目没有看清楚就开始写了
又痛失100分 与day1加起来只有460 幸好day1AK，勉强省一

第一题O(n)被我强行O(nlogn)分治强搞
第二题是近几年最简单的一道dp题，数据结构优化就好了，也有贪心的O(n)做法
第三题华容道是到贼恶心的题目
正解是很早就想到了，但实现上却困难重重整场考试只写这道也未必能调出来
可以发现目标方块的移动一定要靠空格实现，
那么可以预处理空格的移动
得到目标方块移动的最小花费
然后就跑一边最短路，就可以得到答案了
复杂度为$O(n^4+q*X)$ X为最短路复杂度

#include<bits/stdc++.h>
using namespace std;
#define FOR(i,x,y) for(int i=(x);i<=(y);i++)
#define DOR(i,x,y) for(int i=(x);i>=(y);i--)
#define pb push_back
#define M 35

int cx[]={-1,0,0,1};
int cy[]={0,1,-1,0};
struct name{int x,y,step;};
struct Node{int i,x,y;};
queue<name>Q;
queue<Node>ac;
bool mark[M][M],mark1[4005];
int d[4][M][M],dis[M][M][4][4],Dis[M][M],A[M][M];
int a3,b3,n,m,k;

bool check(int x,int y){return (x<0||y<0||x>=n||y>=m||!A[x][y]);}

void BFS(int x,int y){
    memset(mark,0,sizeof(mark));
    memset(Dis,-1,sizeof(Dis));
    Q.push((name){x,y,0});
    mark[x][y]=1;
    while(!Q.empty()){
        name now=Q.front();Q.pop();
        Dis[now.x][now.y]=now.step;
        FOR(i,0,3){
            name nxt=now;
            nxt.step++,nxt.x+=cx[i],nxt.y+=cy[i];
            if(check(nxt.x,nxt.y)||mark[nxt.x][nxt.y])continue;
            mark[nxt.x][nxt.y]=1;
            Q.push(nxt);
        }
    }
}

int Hash(int k,int a2,int b2){return k*900+a2*30+b2;}

void spfa(int a1,int b1,int a2,int b2){
    memset(d,-1,sizeof(d));
    memset(mark1,0,sizeof(mark1));
    A[a2][b2]=0;
    BFS(a1,b1);
    A[a2][b2]=1;
    FOR(i,0,3){
        int x=a2+cx[i],y=b2+cy[i];
        if(check(x,y))continue;
        int way=Dis[x][y];
        if(way==-1)continue;
        ac.push((Node){3-i,x,y});
        d[3-i][x][y]=way+1;
        mark1[Hash(3-i,x,y)]=1;
    }

    while(!ac.empty()){
        Node now=ac.front();ac.pop();
        mark1[Hash(now.i,now.x,now.y)]=0;
        int gg=d[now.i][now.x][now.y];

        for(int i=0;i<4;i++){
            Node nxt=now;
            nxt.x+=cx[i],nxt.y+=cy[i];
            if(check(nxt.x,nxt.y))continue;
            int cost=dis[now.x][now.y][now.i][i]+1;
            if(cost==0)continue;
            nxt.i=3-i;
            if(cost+gg<d[nxt.i][nxt.x][nxt.y]||d[nxt.i][nxt.x][nxt.y]==-1){
                d[nxt.i][nxt.x][nxt.y]=cost+gg;
                int v=Hash(nxt.i,nxt.x,nxt.y);
                if(!mark1[v]){
                    mark1[v]=1;
                    ac.push(nxt);
                }
            }
        }
    }
}
int main(){
    scanf("%d%d%d",&n,&m,&k);
    FOR(i,0,n-1)FOR(j,0,m-1)scanf("%d",&A[i][j]);

    memset(dis,-1,sizeof(dis));
    FOR(i,0,n-1)FOR(j,0,m-1){
        if(A[i][j])FOR(k,0,3){
            int ii=i+cx[k],jj=j+cy[k];
            if(check(ii,jj))continue;
            A[i][j]=0;
            BFS(ii,jj);
            A[i][j]=1;
            FOR(h,0,3)dis[i][j][k][h]=Dis[i+cx[h]][j+cy[h]];
        }
    }

    FOR(i,1,k){
        int a1,a2,b1,b2;
        scanf("%d%d%d%d%d%d",&a1,&b1,&a2,&b2,&a3,&b3);
        a1--,a2--,a3--,b1--,b2--,b3--;
        if(a2==a3&&b2==b3){puts("0");continue;}

        spfa(a1,b1,a2,b2);
        int ans=-1;
        for(int i=0;i<4;i++)
            if(d[i][a3][b3]!=-1)
                if(d[i][a3][b3]<ans||ans==-1)
                    ans=d[i][a3][b3];
        printf("%d\n",ans);
    }
    return 0;
}