题目连接:http://acm.hdu.edu.cn/showproblem.php?pid=3572
题目大意:有m台机器,给你n个任务,对于每个任务都有三个值, 做完任务所需天数p, 第s天之后开始做, 必须在第e天之前结束。每个任务可以在任何时间打断,每个任务可以能被所有机器做, 但同一天每个任务只能被一台机器做, 问你所给的任务能否在限定的时间内被全部做完。
这是网络流的判满流。 网络流最关键的是建图。在我看完这道题后,我根本就没想到给网络流有什么关系, 看别人的博客才知道这是网络流中的一道经典题,时间区间网络流,就是把时间区间里的每个点都连上一条边的建图方式。
这道题我们把0看成超级源点, 把源点和每个任务都连一条边, 容量为做完任务所需的天数, 再把每个任务和其对应的时间区间的每个点都连上一条边, 容量为1;最后再找一个超级汇点, 把所有的时间点都与超级汇点连上一条边, 容量为机器的个数, 应为同一天内最多做m个任务。最后求出网络流中的最大流。判断最大流和做完所有任务所需的总天数是否相等即可。
这道题的边很多,为稠密图, EK算法的O(n*m^2)的复杂度是过不了的, Dinic算法的O(m*n^2)的复杂度是可以过的,其实所有网络流的题EK算法一般都是过不了的,建议学习Dinic 或SAP算法。
#include<stdio.h>
#include<string.h>
#include<queue>
#include<vector>
#include<algorithm>
using namespace std;
const int maxn = 1000 + 10;
const int INF = 0x3f3f3f3f;
struct Edge
{
int from, to, cap, flow;
Edge(int u, int v, int c, int f) : from(u), to(v), cap(c), flow(f){}
};
struct Dinic
{
int n, m, s, t;
vector<int> G[maxn];
vector<Edge> edges;
bool vis[maxn];
int d[maxn], cur[maxn];
void inin(int n)
{
this->n = n;
for(int i = 0; i <= n; i++)
G[i].clear();
edges.clear();
}
void AddEdge(int from, int to, int cap)
{
edges.push_back(Edge(from, to, cap, 0));
edges.push_back(Edge(to, from, 0, 0));
m = edges.size();
G[from].push_back(m-2);
G[to].push_back(m-1);
}
bool bfs()
{
queue<int> Q;
memset(vis, false, sizeof(vis));
d[s] = 0;
Q.push(s);
vis[s] = true;
while(!Q.empty())
{
int x = Q.front();
Q.pop();
for(int i = 0; i < G[x].size(); i++)
{
Edge& e = edges[G[x][i]];
if(!vis[e.to] && e.cap>e.flow)
{
d[e.to] = d[x] + 1;
Q.push(e.to);
vis[e.to] = true;
}
}
}
return vis[t];
}
int dfs(int x, int a)
{
if(x==t || a==0)
return a;
int f, flow = 0;
for(int& i = cur[x]; i < G[x].size(); i++)
{
Edge& e = edges[G[x][i]];
if(d[e.to]==d[x]+1 && (f=dfs(e.to, min(a, e.cap-e.flow)))>0)
{
e.flow += f;
edges[G[x][i]^1].flow -= f;
flow += f;
a -= f;
if(a == 0)
break;
}
}
return flow;
}
int max_flow(int s, int t)
{
this->s = s;
this->t = t;
int flow = 0;
while(bfs())
{
memset(cur, 0, sizeof(cur));
flow += dfs(s, INF);
}
return flow;
}
};
int main()
{
int T, n, m;
Dinic solve;
scanf("%d", &T);
for(int cnt = 1; cnt <= T; cnt++)
{
scanf("%d%d", &n, &m);
int max_day = 0, sum = 0;
for(int i = 1; i <= n; i++)
{
int a, b, c;
scanf("%d%d%d", &a, &b, &c);
sum += a;
max_day = max(max_day, c+n);
solve.AddEdge(0, i, a);
for(int j = b; j <= c; j++)
solve.AddEdge(i, j+n, 1);
}
for(int i = n+1; i <= max_day; i++)
solve.AddEdge(i, max_day+1, m);
int k = solve.max_flow(0, max_day+1);
printf("Case %d: ", cnt);
if(k == sum)
puts("Yes");
else
puts("No");
puts("");
solve.inin(max_day+1);
}
return 0;
}