4356: Ceoi2014 Wall 最短路

Description
给出一个N*M的网格图，有一些方格里面存在城市，其中首都位于网格图的左上角。你可以沿着网络的边界走，要求你走的路线是一个环并且所有城市都要被你走出来的环圈起来，即想从方格图的外面走到任意一个城市一定要和你走的路线相交。你沿着方格的边界走是需要费用的，不同的边界费用可能不同，求最小代价。
1<=N,M<=400,走过边界的代价为正整数且不超过10^9

题解：

先求出左上角到每个城市左上角的最短路，然后有个结论：最优方案一定会包着这些最短路，也就是说不会穿过这些最短路。有了这个之后就好办了，可以把每个格点拆成四个点，拆出来的四个点分别跟他们相邻的点连边，若两点穿过最短路，那么权值为inf，否则为0，对于原图中的边，拆出来的对应的点连就好了，注意还要删掉那些在城市中的点，这样跑一次最短路就可以了。
我犯的错误：
1、make_pair的第一维写成int
2、mark只搞了一半
3、输入$n\times m$的矩阵，搞成$m\times m$
以后要注意啊，比赛时尽量用结构体，不要用make_pair写dijkstra

代码：

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<LL,int>
const int Maxn=410;
const int Maxnm=650100;
const LL inf=(1LL<<50);
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
    return x*f;
}
int n,m,a[Maxn][Maxn],d1[Maxn][Maxn],d2[Maxn][  Maxn];
bool mark[Maxnm][4],del[Maxnm];//0123上下左右 
struct Edge{int y,next;LL d;}e[Maxnm<<4];
int last[Maxnm],len=0;
void ins(int x,int y,LL d){int t=++len;e[t].y=y;e[t].d=d;e[t].next=last[x];last[x]=t;}
void addedge(int x,int y,LL d){if(del[x]||del[y])return;ins(x,y,d),ins(y,x,d);}
LL f[Maxnm];bool vis[Maxnm];
int pre[Maxnm];
void dijkstra(int st)
{
    priority_queue<pa,vector<pa>,greater<pa> >q;
    memset(f,63,sizeof(f));f[st]=0;
    memset(vis,false,sizeof(vis));
    q.push(make_pair(0,st));
    while(!q.empty())
    {
        pa t=q.top();q.pop();
        int x=t.second;
        if(vis[x])continue;
        vis[x]=true;
        for(int i=last[x];i;i=e[i].next)
        {
            int y=e[i].y;
            if(f[x]+e[i].d<f[y])
            {
                f[y]=f[x]+e[i].d;
                pre[y]=x;
                q.push(make_pair(f[y],y));
            }
        }
    }
}
int o(int x,int y){return (x-1)*(m+1)+y;}
bool Vis[Maxnm];
void work(int x)
{
    if(!x||Vis[x])return;
    Vis[x]=true;
    int t=pre[x];
         if(t+1==x)mark[t][3]=true,mark[x][2]=true;
    else if(t-1==x)mark[t][2]=true,mark[x][3]=true;
    else if(t+m+1==x)mark[t][1]=true,mark[x][0]=true;
    else if(t-m-1==x)mark[t][0]=true,mark[x][1]=true;
    work(t);
}
int main()
{
    n=read(),m=read();
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
    a[i][j]=read();
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m+1;j++)
    addedge(o(i,j),o(i+1,j),d1[i][j]=read());
    for(int i=1;i<=n+1;i++)
    for(int j=1;j<=m;j++)
    addedge(o(i,j),o(i,j+1),d2[i][j]=read());
    dijkstra(1);
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
    if(a[i][j])
    {
        work(o(i,j));
        del[((o(i,j)-1)<<2)+4]=del[((o(i,j+1)-1)<<2)+3]=
        del[((o(i+1,j)-1)<<2)+2]=del[((o(i+1,j+1)-1)<<2)+1]=true;
    }
    len=0;
    for(int i=1;i<=o(n+1,m+1);i++)last[i]=0;
    mark[1][0]=mark[1][2]=true;
    for(int i=1;i<=n+1;i++)
    for(int j=1;j<=m+1;j++)
    {
        int x=o(i,j);
        addedge(((x-1)<<2)+1,((x-1)<<2)+2,mark[x][0]?inf:0);
        addedge(((x-1)<<2)+3,((x-1)<<2)+4,mark[x][1]?inf:0);
        addedge(((x-1)<<2)+1,((x-1)<<2)+3,mark[x][2]?inf:0);
        addedge(((x-1)<<2)+2,((x-1)<<2)+4,mark[x][3]?inf:0);
    }
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m+1;j++)
    {
        int x=o(i,j),y=o(i+1,j);
        addedge(((x-1)<<2)+3,((y-1)<<2)+1,d1[i][j]);
        addedge(((x-1)<<2)+4,((y-1)<<2)+2,d1[i][j]);
    }
    for(int i=1;i<=n+1;i++)
    for(int j=1;j<=m;j++)
    {
        int x=o(i,j),y=o(i,j+1);
        addedge(((x-1)<<2)+2,((y-1)<<2)+1,d2[i][j]);
        addedge(((x-1)<<2)+4,((y-1)<<2)+3,d2[i][j]);
    }
    dijkstra(2);
    printf("%lld",f[3]);
}