BZOJ 2007: [Noi2010]海拔

这道题挺不错的。。不难发现都是0和1 原题还有小数迷惑人2333（想想 如果这条边很优 为什么不干脆弄到1？）
然后很快就知道 其实就等于把图分成了割成了两块
那么这就是平面图的最小割了。。
我太菜了 后面就不会了
然后去研究了一下 其实就是平面图转对偶图求最短路 太强了 好劲啊
资料: http://blog.sina.com.cn/s/blog_60707c0f01011fnn.html
其实倒也没有卡spfa 只不过对spfa很不友善
于是我打了个dijkstra 好快啊。。

#include<bits/stdc++.h>#define me(a,x) memset(a,x,sizeof a)using namespace std;const int N=251000;inline int read(){    char ch=getchar(); int x=0,f=1;    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0'; ch=getchar();}    return x*f;}struct node{    int d,x;    friend bool operator <(const node &x,const node &y){return x.d>y.d;}};struct edge{int y,c,next;}a[N<<2]; int len,first[N];void ins(int x,int y,int c){a[++len]=(edge){y,c,first[x]},first[x]=len;}priority_queue<node>q; int d[N],s,t;void dij(){    me(d,63); d[s]=0;    q.push((node){0,s});    while(!q.empty())    {        node x=q.top(); q.pop();        if(d[x.x]<x.d)continue;        for(int k=first[x.x];k;k=a[k].next)        {            int y=a[k].y;            if(d[y]>d[x.x]+a[k].c)            {                d[y]=d[x.x]+a[k].c;                q.push((node){d[y],y});            }        }    }}int main(){    int n=read(),i,j,c ;s=n*n+1,t=s+1;    for(i=1;i<=n;i++)c=read(),ins(i,t,c);    for(i=1;i<n;i++)for(j=1;j<=n;j++)c=read(),ins(i*n+j,(i-1)*n+j,c);    for(i=1;i<=n;i++)c=read(),ins(s,n*(n-1)+i,c);    for(i=1;i<=n;i++)    {        c=read(),ins(s,(i-1)*n+1,c);        for(j=1;j<n;j++)c=read(),ins((i-1)*n+j,(i-1)*n+j+1,c);        c=read(),ins(i*n,t,c);    }    for(i=1;i<=n;i++)c=read(),ins(t,i,c);    for(i=1;i<n;i++)for(j=1;j<=n;j++)c=read(),ins((i-1)*n+j,i*n+j,c);    for(i=1;i<=n;i++)c=read(),ins(n*(n-1)+i,s,c);    for(i=1;i<=n;i++)    {        c=read(),ins((i-1)*n+1,s,c);        for(j=1;j<n;j++)c=read(),ins((i-1)*n+j+1,(i-1)*n+j,c);        c=read(),ins(t,i*n,c);    }    dij();    printf("%d\n",d[t]);    return 0;}