poj1149

题意: 有 M 个猪圈,每个猪圈里初始时有若干头猪。一开始所有猪圈都是关闭的。依次来了 N 个顾客,每个顾客分别会打开指定的几个猪圈,从中买若干头猪。每个顾客分别都有他能够买的数量的上限。每个顾客走后,他打开的那些猪圈中的猪,都可以被任意地调换到其它开着的猪圈里,然后所有猪圈重新关上。

分析: 此题是网络流,网络流里经典的构图题。将顾客看作除源和汇以外的节点,源和每个猪圈的第一个顾客连边,边的权是开始时猪圈中猪的数目,若源和某个节点之间有重边,则将权合并,每个猪圈的前后相邻两顾客之间连边,由前一个顾客指向后一个顾客。每个顾客和汇之间连边,边的权是顾客所希望购买的猪的数目。

网络流,需要对图进行化简

化简规则:

规律 1. 如果几个节点的流量的来源完全相同,且流量为+∞,则可以把它们合并成一个。
规律 2. 如果几个节点的流量的去向完全相同,且流量为+∞,则可以把它们合并成一个。
规律 3. 如果从点 u 到点 v 有一条容量为 +∞ 的边,并且 u 是 v 的唯一流量来源,或者 v 是 u 的唯一流量去向,则可以把 u 和 v 合并成一个节点。

View Code
#include <iostream>
#include
<cstdio>
#include
<cstdlib>
#include
<cstring>
usingnamespace std;

#define maxm 1005
#define maxn 105
#define N maxn
#define E 2 * maxn * maxm
#define inf 0x3f3f3f3f

struct edge
{
int x, y, nxt, c;
} bf[E];

int ne, head[N], cur[N], ps[N], dep[N];
int n, m;
int pig[maxm];
int custom[maxm][maxn];
int num[maxm];
int buy[maxn];
int s, t;

void addedge(int x, int y, int c)
{
bf[ne].x
= x;
bf[ne].y
= y;
bf[ne].c
= c;
bf[ne].nxt
= head[x];
head[x]
= ne++;
bf[ne].x
= y;
bf[ne].y
= x;
bf[ne].c
=0;
bf[ne].nxt
= head[y];
head[y]
= ne++;
}

int flow(int n, int s, int t)
{
int tr, res =0;
int i, j, k, f, r, top;

while (1)
{
memset(dep,
-1, n *sizeof(int));
for (f = dep[ps[0] = s] =0, r =1; f != r;)
for (i = ps[f++], j = head[i]; j; j = bf[j].nxt)
{
if (bf[j].c &&-1== dep[k = bf[j].y])
{
dep[k]
= dep[i] +1;
ps[r
++] = k;
if (k == t)
{
f
= r;
break;
}
}
}
if (-1== dep[t])
break;

memcpy(cur, head, n
*sizeof(int));
for (i = s, top =0;;)
{
if (i == t)
{
for (k =0, tr = inf; k < top; ++k)
if (bf[ps[k]].c < tr)
tr
= bf[ps[f = k]].c;
for (k =0; k < top; ++k)
bf[ps[k]].c
-= tr, bf[ps[k] ^1].c += tr;
res
+= tr;
i
= bf[ps[top = f]].x;
}
for (j = cur[i]; cur[i]; j = cur[i] = bf[cur[i]].nxt)
if (bf[j].c && dep[i] +1== dep[bf[j].y])
break;
if (cur[i])
{
ps[top
++] = cur[i];
i
= bf[cur[i]].y;
}
else
{
if (0== top)
break;
dep[i]
=-1;
i
= bf[ps[--top]].x;
}
}
}
return res;
}

void input()
{
memset(num,
0, sizeof(num));
scanf(
"%d%d", &m, &n);
for (int i =0; i < m; i++)
scanf(
"%d", &pig[i]);
for (int i =0; i < n; i++)
{
int a, b;
scanf(
"%d", &a);
for (int j =0; j < a; j++)
{
scanf(
"%d", &b);
b
--;
custom[b][num[b]
++] = i;
}
scanf(
"%d", &buy[i]);
}
}

void work()
{
ne
=2;
memset(head,
0, sizeof(head));
s
= n;
t
= n +1;
for (int i =0; i < n; i++)
addedge(i, t, buy[i]);
for (int i =0; i < m; i++)
{
if (num[i] >0)
addedge(s, custom[i][
0], pig[i]);
for (int j =1; j < num[i]; j++)
addedge(custom[i][j
-1], custom[i][j], inf);
}
}

int main()
{
//freopen("t.txt", "r", stdin);
input();
work();
printf(
"%d\n", flow(n +2, s, t));
return0;
}
原文地址:https://www.cnblogs.com/rainydays/p/2155710.html