洛谷·[网络流24题]餐巾计划问题

初见安~咕咕咕了这么久终于开始写网络流了呢。这里是传送门:洛谷P1251 餐巾计划问题

题解

这个题跟洛谷P2223几乎一模一样,只是对于【过几天】这个边界有略微的逻辑上的区别。

比较明显是个最小费用流。首先大概的边界我们可以构建出来:从S向每天的点连边,又从每天的点向T连边,并且容量都是该天要用的纸巾的数量。然后就要连代表新买,从快洗过来的和从慢洗过来的纸巾的边,容量都是INF,费用是各自的费用。

现在我们细化。因为从S和T都要分别连出去,而各个点内部又有支配流量的关系,所以我们把每天拆成两个点,第一个点用来送餐巾去洗,与S相连,第二个点用来接收来自快洗慢洗和新买的纸巾,与T相连。这样一来我们对图的建立就很清晰了——

将S连向第i天的第一个点,并提且往第i+m和i+n天的第二个点供餐巾,将第二个点连向T。以及将S连向第i天的第二个点表示新购买餐巾。这样看起来就完了,但是你会发现你WA了。因为有的时候你明明不用再买一条,但你买了。比如m=2的时候,第一个点洗了餐巾给第三天,但是第三天用不到这么多餐巾,所以就屯在那里了,也就是说第四天第五天都是可以用这个餐巾的。所以我们对于每个i还要向i+1匀一条边。是i的第一个点连向i+1的第一个点还是i的第二个点连向i+1的第二个点都是可以的,一个表示延期送去洗,一个表示洗了延期用。

上代码啦——

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<queue>
#define maxn 5000
#define maxm 50000
using namespace std;
typedef long long ll;
const int INF = 0x3f3f3f3f;
int read() {
	int x = 0, f = 1, ch = getchar();
	while(!isdigit(ch)) {if(ch == '-') f = -1; ch = getchar();}
	while(isdigit(ch)) x = (x << 1) + (x << 3) + ch - '0', ch = getchar();
	return x * f;
}

struct edge {int to, w, f, nxt;} e[maxm];
int head[maxn], k = 0;
void add(int u, int v, int w, int f) {
	e[k] = {v, w, f, head[u]}; head[u] = k++;
	e[k] = {u, 0, -f, head[v]}; head[v] = k++;
} 

int n, A, B, fa, fb, F, r[maxn];
int S, T;
int dis[maxn], pre[maxn], incf[maxn];
bool vis[maxn];
bool spfa() {
	memset(dis, 0x3f, sizeof dis); dis[S] = 0;
	memset(vis, 0, sizeof vis); vis[S] = true;
	queue<int> q; incf[S] = INF; q.push(S);
	while(q.size()) {
		register int u = q.front(), v; q.pop(); vis[u] = false;
		for(int i = head[u]; ~i; i = e[i].nxt) {
			v = e[i].to; if(e[i].w && dis[u] + e[i].f < dis[v]) {
				dis[v] = dis[u] + e[i].f;
				incf[v] = min(incf[u], e[i].w); pre[v] = i;
				if(!vis[v]) q.push(v), vis[v] = true;
			}
		}
	}
	if(dis[T] == INF) return false; return true;
}

ll ans = 0;
void update() {
	register int u = T, i;
	while(u != S) {
		i = pre[u];
		e[i].w -= incf[T], e[i ^ 1].w += incf[T];
		u = e[i ^ 1].to;
	}
	ans += 1ll * incf[T] * dis[T];
}//上面都是最小费用流基本操作

signed main() {
	n = read(); memset(head, -1, sizeof head);
	for(int i = 1; i <= n; i++) r[i] = read();
	F = read(), A = read(), fa = read(), B = read(), fb = read();
	
	S = 0, T = n * 2 + 1;
	for(int i = 1; i <= n; i++) {
		add(S, i, r[i], 0), add(i + n, T, r[i], 0);//如题解所言建图
		add(S, i + n, INF, F);
		if(i + A <= n) add(i, i + A + n, INF, fa);
		if(i + B <= n) add(i, i + B + n, INF, fb);
		if(i < n) add(i , i + 1, INF, 0);
	}
	
	while(spfa()) update();
	printf("%lld\n", ans);
	return 0;
}

迎评:)
——End——

发布了158 篇原创文章 · 获赞 23 · 访问量 2万+

猜你喜欢

转载自blog.csdn.net/qq_43326267/article/details/104188773