hdu3870——平面图最小割

2008OI论文原题～周冬《两极相通——浅析最大—最小定理在信息学竞赛中的应用》

可以转化为最短路径问题～

方法讲的很详细，就不再写了～注意这种方法只能针对s, t在平面图最外层边上的情况～

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
usingnamespace std;

constint N =405;

struct node
{
    int u, v, w;
}e[(2*N*N+4*N)*2];

int n, s, t, tot, h[N*N], nxt[(2*N*N+4*N)*2], a[N][N], dis[N*N];
bool visit[N*N];

queue <int> q;

void add(int a, int b, int c)
{
    e[tot].u = a;
    e[tot].v = b;
    e[tot].w = c;
    nxt[tot] = h[a];
    h[a] = tot++;
    e[tot].u = b;
    e[tot].v = a;
    e[tot].w = c;
    nxt[tot] = h[b];
    h[b] = tot++;
}

void spfa()
{
    int u, i, v;
    while(!q.empty())  q.pop();
    q.push(s);
    memset(dis, -1, sizeof(dis));
    memset(visit, false, sizeof(visit));
    dis[s] =0;
    while(!q.empty())
    {
        u = q.front();
        q.pop();
        visit[u] =false;
        for(i = h[u]; i !=-1; i = nxt[i])
        {
            v = e[i].v;
            if(dis[v]==-1|| (dis[v]>dis[u]+e[i].w))
            {
                dis[v] = dis[u] + e[i].w;
                if(!visit[v])
                {
                    visit[v] =true;
                    q.push(v);
                }
            }
        }
    }
}

int main()
{
    int T, i, j;
    scanf("%d", &T);
    while(T--)
    {
        scanf("%d", &n);
        s =0;
        t = (n-1)*(n-1)+1;
        for(i =1; i <= n; i++)
            for(j =1; j <= n; j++)
                scanf("%d", &a[i][j]);
        memset(h, -1, sizeof(h));
        tot =0;
        for(i =1; i < n; i++)
            for(j =1; j < n; j++)
            {
                if(i==1|| j==n-1)
                {
                    if(i ==1)  add((i-1)*(n-1)+j, t, a[i][j]);
                    if(j == n-1)  add((i-1)*(n-1)+j, t, a[i][j+1]);
                    if(i==1&& j<n-1)  add((i-1)*(n-1)+j, (i-1)*(n-1)+j+1, a[i][j+1]);
                    if(i>1&& j==n-1)  add((i-1)*(n-1)+j, (i-2)*(n-1)+j, a[i][j]);
                }
                else
                {
                    add((i-1)*(n-1)+j, (i-2)*(n-1)+j, a[i][j]);
                    add((i-1)*(n-1)+j, (i-1)*(n-1)+j+1, a[i][j+1]);
                }
                if(j==1|| i==n-1)
                {
                    if(j ==1)  add(s, (i-1)*(n-1)+j, a[i][j]);
                    if(i == n-1)  add(s, (i-1)*(n-1)+j, a[i+1][j]);
                    if(j==1&& i>1)  add((i-1)*(n-1)+j, (i-2)*(n-1)+j, a[i][j]);
                    if(i==n-1&& j<n-1)  add((i-1)*(n-1)+j, (i-1)*(n-1)+j+1, a[i][j+1]);
                }
            }
        spfa();
        printf("%d\n", dis[t]);
    }
    return0;
}