BZOJ 2007: [Noi2010]海拔

同1001一样,对偶图最小割转最短路

又被卡spfa= =

话说为什么总是有些人总喜欢卡spfa,又好写一般情况又快,为何老是要逼人写迪杰斯特拉= =

要就题目直接说会卡spfa嘛

好吧其实是前几天考试被spfa卡掉一半点发牢骚的= =不要介意= =

话说最近真的有些背,各种卡时卡空间(差2MB就过了的惨痛经历啊QAQ)

好吧不说那么多了= =
CODE:

#include<cstdio>

#include<iostream>

#include<cstring>

#include<algorithm>

#include<queue>

using namespace std;

#define maxk 510

#define maxn maxk*maxk

#define maxm maxn*4

#define inf 0x7fffffff

struct edges{

    int to,next,dist;

}edge[maxm];

struct node{

    int u,v;

};

bool operator < (node x,node y) {

return x.v>y.v;

}


int l,next[maxn];

void addedge(int x,int y,int z){

    edge[++l]=(edges){y,next[x],z};next[x]=l;

}

int dist[maxn],s,t;

bool b[maxn];

priority_queue<node > q;

int spfa(){

    for (int i=1;i<=t;i++) dist[i]=inf;

    dist[s]=0;

    q.push((node){s,0});

    while (!q.empty()) {

node x=q.top();q.pop();

int u=x.u;

if (b[u]) continue;

b[u]=1;

for (int i=next[u];i;i=edge[i].next) 

   if (dist[edge[i].to]>dist[u]+edge[i].dist) {

dist[edge[i].to]=dist[u]+edge[i].dist;

q.push((node){edge[i].to,dist[edge[i].to]});

   }

    }

    return dist[t];

}

int n,x,id[maxk][maxk],cnt;

int main(){

    scanf("%d",&n);

    n++;

    for (int i=1;i<n;i++) 

for (int j=1;j<n;j++) id[i][j]=++cnt;

    s=++cnt,t=++cnt;

    for (int i=1;i<=n;i++) 

for (int j=1;j<n;j++) {

   scanf("%d",&x);

   addedge(i==1?s:id[i-1][j],i==n?t:id[i][j],x);

}

    for (int i=1;i<n;i++) 

for (int j=1;j<=n;j++) {

   scanf("%d", &x);

   addedge(j==n?s:id[i][j],j==1?t:id[i][j-1],x);

}

    for (int i=1;i<=n;i++) 

for (int j=1;j<n;j++) {

   scanf("%d",&x);

   addedge(i==n?t:id[i][j],i==1?s:id[i-1][j],x);

}

    for (int i=1;i<n;i++) 

for (int j=1;j<=n;j++) {

   scanf("%d",&x);

   addedge(j==1?t:id[i][j-1],j==n?s:id[i][j],x);

}

    printf("%d ",spfa());

    return 0;

}


原文地址:https://www.cnblogs.com/New-Godess/p/4348901.html