BZOJ1834 [ZJOI2010] network 网络扩容

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1834

Description

给定一张有向图，每条边都有一个容量C和一个扩容费用W。这里扩容费用是指将容量扩大1所需的费用。求： 1、 在不扩容的情况下，1到N的最大流； 2、 将1到N的最大流增加K所需的最小扩容费用。

Input

输入文件的第一行包含三个整数N,M,K，表示有向图的点数、边数以及所需要增加的流量。 接下来的M行每行包含四个整数u,v,C,W，表示一条从u到v，容量为C，扩容费用为W的边。

Output

输出文件一行包含两个整数，分别表示问题1和问题2的答案。
就是想学一个费用流模版，交到BZOJ上无限TLE

下载到了数据，发现数据是“错误”的，然后才发现有的是Windows/DOS格式有的是UNIX格式……

然后调，发现程序在读入一半的时候就会崩掉，开调试把我带进了头文件，以为自己遇到了OI生涯最大的玄学

然后在来学校的车上猛然想起自己记录边的数组开小了，maxm打成maxn了……

真是“玄学”，原来是低级错误在作怪……

似乎我也荒废挺久了

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <queue>
#define rep(i,l,r) for(int i=l; i<=r; i++)
#define clr(x,y) memset(x,y,sizeof(x))
#define travel(x) for(Edge *p=last[x]; p; p=p->pre)
using namespace std;
const int INF = 0x3f3f3f3f;
const int maxn = 1010;
const int maxm = 5010;
inline int read(){
int ans = 0, f = 1;
char c = getchar();
for(; !isdigit(c); c = getchar())
if (c == '-') f = -1;
for(; isdigit(c); c = getchar())
ans = ans * 10 + c - '0';
return ans * f;
}
struct Edge{
Edge *pre,*rev; int to,cap,cost;
}edge[maxm<<2],*last[maxn],*cur[maxn],*pre[maxn],*pt;
int n,m,k,u[maxm],v[maxm],c[maxm],x,d[maxn],S,T;
bool isin[maxn];
queue <int> q;
inline void init(){
clr(last,0); pt = edge;
}
inline void add(int x,int y,int z,int w){
pt->pre = last[x]; pt->to = y; pt->cap = z; pt->cost = w; last[x] = pt++;
pt->pre = last[y]; pt->to = x; pt->cap = 0; pt->cost = -w; last[y] = pt++;
last[x]->rev = last[y]; last[y]->rev = last[x];
}
bool bfs(){
while (!q.empty()) q.pop();
clr(d,-1); d[S] = 0; q.push(S);
while (!q.empty()){
int now = q.front(); q.pop();
travel(now){
if (p->cap > 0 && d[p->to] == -1){
d[p->to] = d[now] + 1;
q.push(p->to);
if (p->to == T) return 1;
}
}
}
return 0;
}
int dfs(int x,int flow){
if (x == T || (!flow)) return flow; int w = 0;
for(Edge *p = cur[x]; p && w < flow; p = p->pre){
if (d[p->to] == d[x] + 1 && p->cap > 0){
int delta = dfs(p->to,min(p->cap,flow-w));
p->cap -= delta; p->rev->cap += delta; w += delta;
if (p->cap) cur[x] = p;
}
}
if (w < flow) d[x] = -1;
return w;
}
int maxflow(){
int ans = 0;
while (bfs()){
rep(i,1,n) cur[i] = last[i];
ans += dfs(S,0x7fffffff);
}
return ans;
}
bool spfa(){
while (!q.empty()) q.pop();
clr(d,INF); clr(isin,0);
d[S] = 0; isin[S] = 1; q.push(S);
while (!q.empty()){
int now = q.front(); q.pop(); isin[now] = 0;
travel(now){
if (p->cap && d[p->to] > d[now] + p->cost){
d[p->to] = d[now] + p->cost;
pre[p->to] = p;
if (!isin[p->to]){
isin[p->to] = 1; q.push(p->to);
}
}
}
}
return d[T] != INF;
}
int mincost(){
int cost = 0; int x = INF;
while (spfa()){
for(Edge *p = pre[T]; p; p = pre[p->rev->to]) x = min(x,p->cap);
for(Edge *p = pre[T]; p; p = pre[p->rev->to]){
cost += x * p->cost; p->cap -= x; p->rev->cap += x;
}
}
return cost;
}
int main(){
n = read(); m = read(); k = read(); init();
rep(i,1,m){
u[i] = read(); v[i] = read(); x = read(); c[i] = read();
add(u[i],v[i],x,0);
}
S = 1; T = n; printf("%d ",maxflow());
S = 0; T = n + 1;
rep(i,1,m) add(u[i],v[i],INF,c[i]);
add(0,1,k,0); add(n,n+1,k,0);
printf("%d\n",mincost());
return 0;
}

View Code