hdu 3870 Catch the Theves 平面图转对偶图

题意：求一个平面图的最小割。

如果直接求最小割的话会超时，原图很特殊，是一个平面图，平面图的任意一个割对应对偶图的一条路径，所以即求对偶图的最短路径。

建图：

将S-T连线，将最外面的一个大面（无限大）一分为二了，一个为S*，一个为T*，原图的每个面转化为新图的一个点，对于原图中的每条边，将这条边的两面对应的点连接，再去掉S*与T*的连线。然后求S*到T*的最短路即可。

#include <cstdio>
#include <cstring>
#include <queue>
using namespace std;
#define maxn 450
#define maxe 1000000
#define INF 100000000
struct node
{
    int u,v,w,next;
}e[maxe];
int n,head[450*450],cnt,map[maxn][maxn],s,t,id[maxn][maxn];
int dist[450*450],vis[450*450];
void add(int u,int v,int w)
{
    e[cnt].u=u;
    e[cnt].v=v;
    e[cnt].w=w;
    e[cnt].next=head[u];
    head[u]=cnt++;
    e[cnt].u=v;
    e[cnt].v=u;
    e[cnt].w=w;
    e[cnt].next=head[v];
    head[v]=cnt++;
}
void spfa()
{
    queue <int> Q;
    int i,u,v;
    memset(vis,0,sizeof(vis));
    for(i=0;i<=t;i++) dist[i]=INF;
    dist[s]=0;
    Q.push(s);
    while(!Q.empty())
    {
        u=Q.front();
        Q.pop();
        vis[u]=0;
        for(i=head[u];i!=-1;i=e[i].next)
        {
            v=e[i].v;
            if(dist[v]>dist[u]+e[i].w)
            {
                dist[v]=dist[u]+e[i].w;
                if(!vis[v])
                {
                    vis[v]=1;
                    Q.push(v);
                }
            }
        }
    }
}
int main()
{
    int cas,i,j;
    scanf("%d",&cas);
    while(cas--)
    {
        scanf("%d",&n);
        for(i=1;i<=n;i++)
            for(j=1;j<=n;j++)
                scanf("%d",&map[i][j]);
        memset(head,-1,sizeof(head));
        cnt=0;
        for(i=1;i<n;i++)
            for(j=1;j<n;j++)
                id[i][j]=++cnt;
        cnt=0;
        s=0;
        t=(n-1)*(n-1)+1;
        for(i=1;i<n-1;i++)
        {
            for(j=1;j<n-1;j++)
            {
                add(id[i][j],id[i+1][j],map[i+1][j]);
                add(id[i][j],id[i][j+1],map[i][j+1]);
            }
        }
        for(i=1;i<n-1;i++)
        {
            add(id[n-1][i],id[n-1][i+1],map[n-1][i+1]);
            add(id[i][n-1],id[i+1][n-1],map[i+1][n-1]);
        }
        for(i=1;i<n;i++)
        {
            add(s,id[1][i],map[1][i]);
            add(id[n-1][i],t,map[n][i]);
        }
        for(i=1;i<n;i++)
        {
            add(s,id[i][n-1],map[i][n]);
            add(id[i][1],t,map[i][1]);
        }
        spfa();
        printf("%d\n",dist[t]);
    }
    return 0;
}