BZOJ3438 小M的作物(和拓展)
题目链接:戳我
我们如果要选择一种种植情况的话,一定是其他的选择都不可行了。这种决策问题用最小割来处理最好不过。
建图方式——A为源点,B为汇点。然后将每个点分别向A,B连边,边权为种植它的价值。组合的话,我们考虑新建两个节点,一个连A,对应集合中的每个数,连INF;一个连B,从集合中的每个数向它连INF。
最后的答案是总点权-最小割。割边表示不选。
如果选择组合,我们肯定是先割掉了集合中节点一侧的所有边,然后发现组合对应节点还没有割掉,我们把它割掉就行了,而割掉的这个边也肯定是同一侧的。
如果不选择组合,我们肯定是割掉了集合中一部分节点的一侧边,一部分节点的另一侧边。这时候我们发现组合对应的这两个节点,都可以通过组合中节点没有割掉的边形成增广路。所以我们不得不把组合对应的两个点的边都割掉,以表示不选这个集合。
代码如下:
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<queue>
#define INF 0x3f3f3f3f
#define MAXN 100010
#define S 0
#define T n+m*2+1
using namespace std;
int n,m,t=1,ans;
int head[MAXN],dep[MAXN],a[MAXN],b[MAXN],cur[MAXN];
struct Edge{int nxt,to,dis;}edge[MAXN<<4];
inline void add(int from,int to,int dis)
{
edge[++t].nxt=head[from],edge[t].to=to,edge[t].dis=dis,head[from]=t;
edge[++t].nxt=head[to],edge[t].to=from,edge[t].dis=0,head[to]=t;
}
inline bool bfs()
{
queue<int>q;
memset(dep,0x3f,sizeof(dep));
memcpy(cur,head,sizeof(cur));
q.push(S);dep[S]=0;
while(!q.empty())
{
int u=q.front();q.pop();
for(int i=head[u];i;i=edge[i].nxt)
{
int v=edge[i].to;
if(dep[v]==0x3f3f3f3f&&edge[i].dis)
dep[v]=dep[u]+1,q.push(v);
}
}
if(dep[T]==0x3f3f3f3f) return false;
return true;
}
inline int dfs(int x,int f)
{
if(!f||x==T) return f;
int w,used=0;
for(int i=cur[x];i;i=edge[i].nxt)
{
int v=edge[i].to;
cur[x]=i;
if(dep[v]==dep[x]+1&&(w=dfs(v,min(edge[i].dis,f))))
{
edge[i].dis-=w,edge[i^1].dis+=w;
f-=w,used+=w;
if(!f) break;
}
}
return used;
}
inline int dinic()
{
int cur_ans=0;
while(bfs()) cur_ans+=dfs(S,INF);
return cur_ans;
}
int main()
{
#ifndef ONLINE_JUDGE
freopen("ce.in","r",stdin);
#endif
scanf("%d",&n);
for(int i=1;i<=n;i++) scanf("%d",&a[i]);
for(int i=1;i<=n;i++) scanf("%d",&b[i]);
scanf("%d",&m);
for(int i=1;i<=m;i++)
{
int k,suma,sumb;
scanf("%d%d%d",&k,&suma,&sumb);
ans+=suma+sumb;
add(S,i+n,suma),add(i+m+n,T,sumb);
//printf("[%d %d] %d\n",S,i+n,suma);
//printf("[%d %d] %d\n",i+m+n,T,sumb);
for(int j=1;j<=k;j++)
{
int cur;
scanf("%d",&cur);
add(i+n,cur,INF),add(cur,i+n+m,INF);
//printf("[%d %d] %d\n",i+n,cur,INF);
//printf("[%d %d] %d\n",cur,i+n+m,INF);
}
}
for(int i=1;i<=n;i++) add(S,i,a[i]);//rintf("[%d %d] %d\n",S,i,a[i]);
for(int i=1;i<=n;i++) add(i,T,b[i]);//printf("[%d %d] %d\n",i,T,b[i]);
for(int i=1;i<=n;i++) ans+=a[i]+b[i];
printf("%d\n",ans-dinic());
return 0;
}
其实类似的模型还有这道题
这道题不仅是把组合分给A,B的代价不同,而且不要这个组合的话,还需要倒扣钱。
但是思路还是一样的,我们先把A,B的得分加上,然后组合的a,b,c的得分也都加上。
然后我们这样子连边——
add(S,id_a(i),A[i]);
add(id_b(i),T,B[i]);
add(S,id_a(i),EA+EC);
add(id_a(i),u,INF);
add(id_a(i),v,INF);
add(id_b(i),T,EB+EC);
add(u,id_b(i),INF);
add(v,id_b(i),INF);
BZOJ3438 小M的作物(和拓展)的更多相关文章
- bzoj3438: 小M的作物(那年花开最小割)
3438: 小M的作物 题目:传送门 题解: 最小割标准水题(做了几天的最小割之后表示是真的水) 为什么水:博主已经做过两道基本一样的题目了... 详情参考:bzoj3894 代码: #include ...
- BZOJ3438 小M的作物(最小割)
题目 Source http://www.lydsy.com/JudgeOnline/problem.php?id=3438 Description 小M在MC里开辟了两块巨大的耕地A和B(你可以认为 ...
- BZOJ3438小M的作物——最小割
题目描述 小M在MC里开辟了两块巨大的耕地A和B(你可以认为容量是无穷),现在,小P有n中作物的种子,每种作物的种子 有1个(就是可以种一棵作物)(用1...n编号),现在,第i种作物种植在A中种植可 ...
- BZOJ3438 小M的作物
AC通道:http://www.lydsy.com/JudgeOnline/problem.php?id=3438 这题觉得和上题有点类似吧. 如果没有联合在一起的收成,可以比较好做[我们将属于A的表 ...
- 【BZOJ3438】小M的作物 最小割
[BZOJ3438]小M的作物 Description 小M在MC里开辟了两块巨大的耕地A和B(你可以认为容量是无穷),现在,小P有n中作物的种子,每种作物的种子 有1个(就是可以种一棵作物)(用1. ...
- 【BZOJ-3438】小M的作物 最小割 + 最大权闭合图
3438: 小M的作物 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 825 Solved: 368[Submit][Status][Discuss ...
- 「BZOJ3438」小M的作物(最小割
3438: 小M的作物 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 1891 Solved: 801[Submit][Status][Discus ...
- BZOJ 3438: 小M的作物( 最小割 )
orz出题人云神... 放上官方题解... 转成最小割然后建图跑最大流就行了... ---------------------------------------------------------- ...
- P1361 小M的作物
P1361 小M的作物 题目描述 小M在MC里开辟了两块巨大的耕地A和B(你可以认为容量是无穷),现在,小P有n中作物的种子,每种作物的种子有1个(就是可以种一棵作物)(用1...n编号). 现在,第 ...
随机推荐
- springboot - 应用实践(1)认识springboot
1.为什么要推出springboot springboot设计的目的是用来简化新spring应用的初始搭建以及开发过程.springboot遵循“约定优于配置”原则. 2.springboot默认的配 ...
- 重装java后hadoop配置文件的修改
1.删除hdfs-site.xml中dfs.namenode.name.dir目录和dfs.datanode.data.dir目录 然后 hdfs namenode -format 不然将无法启动na ...
- C中的异常处理
1,C 语言崇尚简洁高效,因此语言本身并没有异常处理的相关语法规则,但是异常处理在 C 语言中 是存在的,我们有必要从 C 语言开始先看一看 C 语言中的异常处理是怎样, 然后对比 C++ 里面的异常 ...
- python-day18(正式学习)
目录 numpy模块 numpy简介 为什么要用numpy 创建numpy数组 numpy数组的基本属性 获取numpy数组的行列数 切割numpy数组 numpy数组元素替换 numpy数组的合并 ...
- springboot导包spring-boot-starter-parent出现错误
<parent> <groupId>org.springframework.boot</groupId> <artifactId>spring-boot ...
- openstack docker build error
1. _proto_tcp = socket.getprotobyname('tcp') OSError: protocol not found you should have a /etc/prot ...
- Flask开发系列之模板
Flask开发系列之模板 本文对<FlaskWeb开发:基于python的Web应用开发实战>模板一节做的总结. Jinja2模板引擎 模板 模板是一个包含响应文本的文件,其中包含用占位变 ...
- Shiro多Realm验证
在单Realm的基础上,进行如下内容的修改 原有的ShiroRealm.java: package com.atguigu.shiro.realms; import org.apache.shiro. ...
- 从零开始学MySQL(二)
鉴于上节篇幅以安装为主,因此对于调用mysql所需要使用的“命令”只是略微提及.随之而来就会带给读者诸多不解了,因为你会思考,这串长长的字符到底有什么特殊的含义呢?聪明的你可能早就抱着好奇心去“摆渡” ...
- VMware虚拟机中CentOS/redhat设置固定IP
你的笔记本中的VMware中redhat或centOS系统,如果想在上面建站,而又如果你需要在家里和公司都能访问该站(至少希望你自己的笔记本能访问),那么就需要将虚拟机IP设置为固定IP了.以下介绍两 ...