[SDOI2009]晨跑

看懂题意就很简单了233以前语文差的看不懂题意

将每个点拆成出点入点,流量为1费用为0,控制每个点经过一次,然后最小费用最大流

Code:

#include <cstdio>
#include <queue>
#include <cstring>
#include <bitset>
using namespace std;
#define M 100001
#define N 10001
#define INF (0x7FFFFFFF)
#define S(x) (x<<1)
#define T(x) ((x<<1)-1)
struct E
{
	int to, f, c, next, op;
}G[M];
int point[N], cnt, S, T, f, n, m, cost, D[N], path[N], prev[N];
bitset<N> vis;
queue<int> Q;
inline void add(int u, int v, int f, int c)
{
	cnt++, G[cnt] = (E){v, f, c, point[u], cnt+1}, point[u] = cnt;
	cnt++, G[cnt] = (E){u, 0, -c, point[v], cnt-1}, point[v] = cnt;
}
inline bool spfa()
{
	for (int i = S;i <= T; ++i)D[i] = INF;
	D[S] = 0;
	vis[S] = 1;
	prev[S] = -1;
	Q.push(S);
	while (!Q.empty())
	{
		int u = Q.front();
		Q.pop();
		for (int i = point[u];i;i = G[i].next)
		{
			int v = G[i].to;
			if (G[i].f > 0 && D[v] > D[u] + G[i].c)
			{
				D[v] = D[u] + G[i].c;
				path[v] = i;
				prev[v] = u;
				if (!vis[v])
					vis[v] = 1, Q.push(v);
			}
		}
		vis[u] = 0;
	}
	return D[T] != INF;
}
inline int argument()
{
	int delta(INF), i, j, res(0);
	for (i = T;prev[i] != -1;i = prev[i])
	{
		j = path[i];
		delta = min(delta, G[j].f);
	}
	f += delta;
	for (i = T;prev[i] != -1;i = prev[i])
	{
		j = path[i];
		G[j].f -= delta;
		G[G[j].op].f += delta;
		res += delta * G[j].c;
	}
	return res;
}
int main()
{
	scanf("%d %d", &n, &m);
	S = T(1), T = S(n);
	for (int i = 1;i <= m; ++i)
	{
		int a, b, c, x(1), y(n);
		scanf("%d %d %d", &a, &b, &c);
		x = T(a), y = S(b);
		add(x, y, 1, c);
	}
	for (int i = 1;i <= n; ++i)
		if (i != 1 && i != n)
			add(S(i), T(i), 1, 0);
	while (spfa())
		cost += argument();
	printf("%d %d\n", f, cost);
}	


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值