百度之星初赛B 度度熊的交易计划(最小费用可行流)

原创 2017年08月14日 16:38:42

这道题说实话挺简单的,但是建图上太蠢了导致当时一直TLE,现在才补出来。

当时用的是把一个点拆成三个点的做法,导致复杂度不知道扩大了多少倍。。。现在来总结一下,当时建图的方法按照样例来看是:


现在看来这个图是完全可以优化的,因为源点连出去的边全都是流量为INF,费用为0的边,那么此时源点连出去的点是无用的,直接合并到源点即可。汇点也是。

所以就可以得到下面这种建图:


这样瞬间点数少了四个,边数也少了四条,这样跑费用流时候的复杂度整整下降了一半以上。到了点数边数十分多的时候效果更加明显。

这道题由于是要求最大费用,所以我们直接将原本是对答案贡献为正的边的费用建为负值,对答案贡献为负的费用建为正值。最后得到答案再取一次相反数即可。

由于我们是要最大收益,所以这道题最后得到的网络不一定是最大流,我们只需要取出费用为负的增广路即可。所以这题是可行流而不是最大流。

代码如下:

#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
#include<utility>
#include<stack>
#include<algorithm>
#include<cstring>
#include<string>
#include<cmath>
using namespace std;
typedef pair<int, int> pii;
const int INF = 0x3f3f3f3f;
const int maxn = 505 + 5;
const int maxm = 40000 + 5;
int n, m;
int head[maxn], to[maxm], front[maxm], flow[maxm], cost[maxm], ppp;
int dis[maxn], minflow[maxn];
bool flag[maxn];
pair<int, int> par[maxn];
int a[maxn], b[maxn], c[maxn], d[maxn];
int G[maxn][maxn];

struct MIN_COST_MAX_FLOW{
	
	void init() {
		memset(head, -1, sizeof(head));
		ppp = 0;
	}
	
	bool spfa(int s, int t)
	{
		int u, v;
		for(int i = 0; i <= t; i++)
			dis[i] = INF;
		memset(flag, 0, sizeof(flag));
		dis[s] = 0;
		minflow[s] = INF;
		queue <int> q;
		q.push(s);
		while(!q.empty())
		{
			u = q.front();
			q.pop();
			flag[u] = 0;
			for(int i = head[u]; ~i; i = front[i])
			{
				v = to[i];
				if(flow[i] && dis[v] > dis[u] + cost[i])
				{
					dis[v] = dis[u] + cost[i];
					par[v] = (make_pair(u, i));
					minflow[v] = min(minflow[u], flow[i]);
					if(!flag[v])
					{
						flag[v] = 1;
						q.push(v);
					}
				}
			}
		}
		return dis[t] != INF;
	}
	
	int slove(int s, int t)
	{
		int ans = 0, p;
		while(spfa(s, t))
		{
			if(dis[t] > 0)
				break;
			p = t;
			while(p != s)
			{
				flow[par[p].second] -= minflow[t];
				flow[par[p].second^1] += minflow[t];
				p = par[p].first;
			}
			ans += dis[t] * minflow[t];  
		}
		return -ans;
	}
	
	void add_edge(int u, int v, int f, int c)
	{
		to[ppp] = v, front[ppp] = head[u], flow[ppp] = f, cost[ppp] = c, head[u] = ppp++;
		to[ppp] = u, front[ppp] = head[v], flow[ppp] = 0, cost[ppp] = -c, head[v] = ppp++;
	}
}mcmf;



int main() {
#ifndef ONLINE_JUDGE
	freopen("in.txt", "r", stdin);
#endif
	while(cin >> n >> m) {
		memset(head, -1, sizeof(head));
		ppp = 0;
		memset(G, 0x3f, sizeof(G));
		int s = 0, t = n + 1;
		for(int i = 1; i <= n; i++) {
			scanf("%d%d%d%d", &a[i], &b[i], &c[i], &d[i]);
			mcmf.add_edge(s, i, b[i], a[i]);
			mcmf.add_edge(i, t, d[i], -c[i]);
		}
		while(m--) {
			int u, v, val;
			scanf("%d%d%d", &u, &v, &val);
			if(u == v)
				continue;
			G[u][v] = G[v][u] = min(val, G[u][v]);
		}
		for(int i = 1; i <= n; i++) {
			for(int j = i + 1; j <= n; j++) {
				if(G[i][j] != INF) {
					mcmf.add_edge(i, j, INF, G[i][j]);
					mcmf.add_edge(j, i, INF, G[i][j]);
				}
			}
		}
		int ans = mcmf.slove(s, t);
		printf("%d\n", ans);
	}
	return 0;
}


版权声明:本文为博主原创文章,未经博主允许不得转载。

2017百度之星资格赛题解

1003度度熊与邪恶大魔王 Accepts: 1503 Submissions: 9026 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: ...
  • winter2121
  • winter2121
  • 2017年08月05日 17:49
  • 2233

Deep learning:四十(龙星计划2013深度学习课程小总结)

头脑一热,坐几十个小时的硬座北上去天津大学去听了门4天的深度学习课程,课程预先的计划内容见:http://cs.tju.edu.cn/web/courseIntro.html。上课老师为微软研究院的大...
  • mydear_11000
  • mydear_11000
  • 2016年03月13日 15:12
  • 1422

百度地图绘制实时路线以及最短线路规划

如何使用百度地图绘制实时路线以及最短线路规划最近在做百度地图的实时路线绘制,发现一些问题,比如由于定位漂移带来的路线绘制偏差,还有由于定位漂移,导致人未走动时,也会绘制路线等。百度鹰眼的线路纠偏个人感...
  • zero172
  • zero172
  • 2016年05月03日 13:36
  • 6398

百度之星初赛1006(计算几何:能包含凸包的最小矩形面积)

矩形面积    Accepts: 717    Submissions: 1619  Time Limit: 2000/1000 MS (Java/Others)    Memory...
  • u013445530
  • u013445530
  • 2015年05月30日 20:25
  • 1097

百度之星2012初赛Day1 A - 度度熊就是要第一个出场

开始还很担心的~~今天出来结果~~过初赛了~~Day1晋级的~~第一场我是做了两道~~骗了一道(差不多过不了几个点..因为我C题原意是用旋转卡壳骗点..结果赛后发现有个地方写错了..囧.)..结果还好...
  • kk303
  • kk303
  • 2012年06月08日 20:25
  • 898

2015年百度之星初赛(1) --- B 找连续数

找连续数 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submi...
  • u013371163
  • u013371163
  • 2017年03月05日 17:17
  • 138

2016"百度之星" 初赛(Astar Round2B) 1001 区间的价值

题意:我们定义“区间的价值”为一段区间的最大值*最小值。 一个区间左端点在L,右端点在R,那么该区间的长度为(R−L+1)。 现在聪明的杰西想要知道,对于长度为k的区间,最大价值的区间价值是多少。...
  • u014664226
  • u014664226
  • 2016年05月22日 20:43
  • 471

2016"百度之星" - 初赛(Astar Round2B)1003 瞬间移动(HDU5698)(逆元+快速幂)

2016”百度之星” - 初赛(Astar Round2B)1003 瞬间移动(HDU5698)(逆元+快速幂)Time Limit: 4000/2000 MS (Java/Others) ...
  • ping1632743560
  • ping1632743560
  • 2016年05月28日 22:47
  • 394

2016"百度之星" - 初赛(Astar Round2B)题解

1001我们定义“区间的价值”为一段区间的最大值*最小值。一个区间左端点在LL,右端点在RR,那么该区间的长度为(R-L+1)(R−L+1)。我们想要知道的是,对于长度为1-n的区间,最大价值的区间价...
  • Miracle_ma
  • Miracle_ma
  • 2016年05月23日 10:52
  • 915

2016"百度之星" - 初赛(Astar Round2B)

5696 区间的价值 对于一个区间,考虑最小的数做出的贡献。对于最小的数,一定是和最大的数组成的对数才能对答案产生影响,否则对答案是没有作用的。所以就可以从最小值的位置,把整个区间分裂成两段,然后再...
  • u010697167
  • u010697167
  • 2016年05月22日 21:34
  • 622
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:百度之星初赛B 度度熊的交易计划(最小费用可行流)
举报原因:
原因补充:

(最多只允许输入30个字)