抱歉,您的浏览器无法访问本站
本页面需要浏览器支持(启用)JavaScript
了解详情 >

传送门:洛谷


Solution

网络流 先引用一句真理:网络流最重要的就是建模 今天这道题让我深有体会 首先,观察数据范围,n=100,一般这种100-1000的图论题,很有可能是网络流. 那就直接从网络流的角度入手 考虑这样建模 kwklcV.png 建模要点如下: 1.建权值为1的边,保证每个食物和水仅用一次 2.没了 对以上的图求一个最大流,那不就是我们想要的最大的匹配数吗? 看起来是不是很OjbK? . 其实不然,这样子一头牛有可能脚踏N条食物和水,但是题目要求一头牛只能吃喝一次 反例如下: kwk0c6.png 所以说,我们要对一头牛吃的东西做一个限制,保证其只流过1 怎么限制呢? 直接把一头牛拆成两头牛,中间连一条边就OK了嘛 如下图: kwkrnO.png 接下来就可以考虑如何给点编号了 我的编号方法很直接,很暴力 源点:1 食物: 2 ~ 2+f-1 牛 2+f ~ 2+f+2n -1 水: 2+f+2n ~ 2+f+2n+d 汇点:1000 如下图所示 按照以上方法,每种物品对应的点为: 食物i : 1+i 牛i : 1+f+i 牛i的右边的分点: 1+f+n+i 水i: 1+f+2*n+i 然后dinic直接求一波最大流就可以带走啦 DINIC教程传送门: http://www.cnblogs.com/SYCstudio/p/7260613.html (个人感觉写得很清楚)

Code

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
//Luogu P2891 [USACO07OPEN]吃饭Dining
//Apr,2ed,2018
//Dinic求最大流
#include<iostream>
#include<cstdio>
#include<vector>
#include<cstring>
using namespace std;
long long read()
{
long long x=0,f=1; char c=getchar();
while(!isdigit(c)){if(c=='-') f=-1;c=getchar();}
while(isdigit(c)){x=x*10+c-'0';c=getchar();}
return x*f;
}
const int N=100+10;
const int M=N*10;
const int inf=0x3f3f3f3f;
struct line
{
int t,rev,w;
};
vector <line> e[M];
inline void AddLine(int s,int t)
{
line temp;
temp.t=t,temp.rev=e[t].size(),temp.w=1;
e[s].push_back(temp);
temp.t=s,temp.rev=e[s].size()-1,temp.w=0;
e[t].push_back(temp);
}
int n,f,d;
int dl[M],head,tail,depth[M];
bool visited[M];
bool bfs()
{
memset(visited,0,sizeof visited);
dl[1]=1,depth[1]=1;
visited[1]=true;
head=1,tail=2;
while(head<tail)
{
int now=dl[head],size=e[now].size();
for(int i=0;i<size;i++)
if(visited[e[now][i].t]==false and e[now][i].w>0)
{
visited[e[now][i].t]=true;
dl[tail++]=e[now][i].t;
depth[e[now][i].t]=depth[now]+1;
}
head++;
}
return visited[1000];
}
int dfs(int now,int f)
{
if(now==1000)
return f;
int size=e[now].size(),ans=0;
for(int i=0;i<size;i++)
if(depth[e[now][i].t]==depth[now]+1 and e[now][i].w>0)
{
int t_ans=dfs(e[now][i].t,min(f,e[now][i].w));
e[now][i].w-=t_ans;
e[e[now][i].t][e[now][i].rev].w+=t_ans;
f-=t_ans,ans+=t_ans;
if(f==0) break;
}
return ans;
}
int Dinic()
{
int ans=0;
while(bfs())
ans+=dfs(1,inf);
return ans;
}
int main()
{
n=read(),f=read(),d=read();
for(int i=1;i<=1000;i++)
e[i].reserve(16);
for(int i=1;i<=f;i++)
AddLine(1,1+i);
for(int i=1;i<=d;i++)
AddLine(1+f+2*n+i,1000);
for(int i=1;i<=n;i++)
{
int F=read(),D=read(),temp;
for(int j=1;j<=F;j++)
{
temp=read();
AddLine(1+temp,1+f+i);
}
for(int j=1;j<=D;j++)
{
temp=read();
AddLine(1+f+n+i,1+f+2*n+temp);
}
}
for(int i=1;i<=n;i++)
AddLine(1+f+i,1+f+n+i);

printf("%d",Dinic());
return 0;
}

C++

后记

写的时候发现自己真的有点生疏了,这种比较靠记忆的算法还是久不久搞一下来恢复记忆为妙 网络流的时间复杂度真玄学

评论