POJ 4370 0 or 1 图论模型的转化

这题的优点就在于模型的转化

转化完就是个非常简单的最短路了

显然，题目给的是一个0/1规划模型。

解题的关键在于如何看出这个模型的本质。

3个条件明显在刻画未知数之间的关系，从图论的角度思考问题，容易得到下面3个结论：

1.X12+X13+...X1n=1 于是1号节点的出度为1

2..X1n+X2n+...Xn-1n=1 于是n号节点的入度为1

3.∑Xki =∑Xij 于是2~n-1号节点的入度必须等于出度

于是3个条件等价于一条从1号节点到n号节点的路径，故Xij=1表示需要经过边(i,j)，代价为Cij。Xij=0表示不经过边(i,j)。注意到Cij非负且题目要求总代价最小，因此最优答案的路径一定可以对应一条简单路径。

最终，我们直接读入边权的邻接矩阵，跑一次1到n的最短路即可.

#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
#include <cstdio>
#include <cmath>
#include <queue>
#include <map>
#include <set>
#define eps 1e-5
#define MAXN 333
#define MAXM 211111
#define INF 1000000007
using namespace std;
int n;
struct EDGE
{
    int v, next, w;
}edge[MAXM];
int head[MAXN], e;
void init()
{
    memset(head, -1, sizeof(head));
    e = 0;
}
void add(int x, int y, int w)
{
    edge[e].v = y;
    edge[e].w = w;
    edge[e].next = head[x];
    head[x] = e++;
}
int d[MAXN], vis[MAXN];
int q[MAXN * MAXN];
void spfa(int src)
{
    int h = 0, t = 0;
    for(int i = 0; i <= n; i++) d[i] = INF, vis[i] = 0;
    vis[src] = 1;
    d[src] = 0;
    q[t++] = src;
    while(h < t)
    {
        int u = q[h++];
        vis[u] = 0;
        for(int i = head[u]; i != -1; i = edge[i].next)
        {
            int v = edge[i].v;
            int w = edge[i].w;
            if(d[u] + w < d[v])
            {
                d[v] = d[u] + w;
                if(!vis[v])
                {
                    q[t++] = v;
                    vis[v] = 1;
                }
            }
        }
    }
}
int main()
{
    while(scanf("%d", &n) != EOF)
    {
        init();
        int x;
        for(int i = 1; i <= n; i++)
            for(int j = 1; j <= n; j++)
            {
                scanf("%d", &x);
                add(i, j, x);
            }
        spfa(1);
        printf("%d\n", d[n]);
    }
    return 0;
}