BZOJ 5421: 收藏家
直接搞很复杂,考虑转化问题
题目只要求第1个人最多能获得的物品数量
所以如果一种物品拥有多个和一个是没区别的
那么考虑每种物品对第1个人怎样贡献
显然要经过一些交换最终到达第一个人那里
发现很像一个流,那么考虑建立网络流模型
建一个源点向每个点连一条最大流量为1的边,相当于初始每个点有1个物品
点1向汇点连一条 $a_1$ 的边,因为点1最多能放 $a_1$ 个物品
因为一个点 i 同时最多只能有 $a_i$ 种物品,所以也把其他个点拆成两个,之间连一条流量为 $a_i$ 的边
考虑点与点之间的连接
因为交换是按时间顺序的,所以不可能直接连
同样拆点
把一个点 i 分成 m 个点,表示 i 在 1~m 的时间点上的状态,显然这 m 个点之间的流量为 $a_i$(一个点 i 最多同时有 $a_i$ 种不同的物品)
这样就可以把不同的点连起来了
具体说来就是:如果在时间 i , a 和 b 发生了交换,那么在 a 拆出来的第 i 个点和 b 拆出来的第 i 个点之间连双向边
显然边权均为 1 (一次只能换一个物品)
可以参照下面的丑图:

但是这样有$nm$个点,显然不行
但是可以发现只有向其他点连接的点是有用的,即一串 m 个点只有向外面其他点有连接的点是有用的,所以我们动态地加点连边,只有对于需要的点我们才要连边
这样因为边数是 m 所以要加的点数就是 2m
要注意一些细节
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<queue>
using namespace std;
typedef long long ll;
inline int read()
{
int x=,f=; char ch=getchar();
while(ch<''||ch>'') {if(ch=='-') f=-; ch=getchar(); }
while(ch>=''&&ch<='') { x=(x<<)+(x<<)+(ch^); ch=getchar(); }
return x*f;
}
const int N=3e5+,INF=1e9+;
int fir[N],from[N<<],to[N<<],val[N<<],cntt=,Fir[N];
inline void add(int a,int b,int c)
{
from[++cntt]=fir[a]; fir[a]=cntt;
to[cntt]=b; val[cntt]=c;
from[++cntt]=fir[b]; fir[b]=cntt;
to[cntt]=a; val[cntt]=;
}
queue <int> q;
int n,m,tot,ans;
int S,T,p[N],pre[N],dep[N];
//p是每个点的容量,pre[i]是为了动态加点,存点i的上个时间的点的编号
//以下为Dinic
bool BFS()
{
memset(dep,,sizeof(dep));
q.push(S); dep[S]=; int x;
while(!q.empty())
{
x=q.front(); q.pop();
for(int i=fir[x];i;i=from[i])
{
int &v=to[i]; if(dep[v]||!val[i]) continue;
dep[v]=dep[x]+; q.push(v);
}
}
for(int i=;i<=tot;i++) Fir[i]=fir[i];
return dep[T] ? : ;
}
int DFS(int x,int mif)
{
if(!mif||x==T) return mif;
int fl=,res=;
for(int i=Fir[x];i;i=from[i])
{
Fir[x]=i; int &v=to[i]; if(dep[v]!=dep[x]+) continue;
if( res=DFS(v, min(mif,val[i]) ) )
{
fl+=res; mif-=res;
val[i]-=res; val[i^]+=res;
if(!mif) break;
}
}
return fl;
}
//以上为Dinic
void slove()//处理读入
{
int a,b; ans=tot=; cntt=;//多组数据记得初始化
memset(fir,,sizeof(fir));
n=read(); m=read();
S=++tot; T=++tot;//建源点和汇点
for(int i=;i<=n;i++)
{
add(S,pre[i]=++tot,);//源点向每个初始点连边
p[i]=read();
}
while(m--)
{
a=read(); b=read();
int na=++tot,nb=++tot;//新的时间的点
add(pre[a],na,p[a]); add(pre[b],nb,p[b]);//和上一时间的点相连
add(na,nb,); add(nb,na,);//之间连双向边
pre[a]=na; pre[b]=nb;//更新pre
}
add(pre[],T,p[]);//最后点1连向汇点,容量为p[1]
while(BFS()) ans+=DFS(S,INF);
printf("%d\n",ans);
}
int main()
{
//freopen("collection.in","r",stdin);
//freopen("collection.out","w",stdout);
int T=read();
while(T--) slove();
return ;
}
BZOJ 5421: 收藏家的更多相关文章
- BZOJ 2127: happiness [最小割]
2127: happiness Time Limit: 51 Sec Memory Limit: 259 MBSubmit: 1815 Solved: 878[Submit][Status][Di ...
- BZOJ 3275: Number
3275: Number Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 874 Solved: 371[Submit][Status][Discus ...
- BZOJ 2879: [Noi2012]美食节
2879: [Noi2012]美食节 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 1834 Solved: 969[Submit][Status] ...
- bzoj 4610 Ceiling Functi
bzoj 4610 Ceiling Functi Description bzoj上的描述有问题 给出\(n\)个长度为\(k\)的数列,将每个数列构成一个二叉搜索树,问有多少颗形态不同的树. Inp ...
- BZOJ 题目整理
bzoj 500题纪念 总结一发题目吧,挑几道题整理一下,(方便拖板子) 1039:每条线段与前一条线段之间的长度的比例和夹角不会因平移.旋转.放缩而改变,所以将每条轨迹改为比例和夹角的序列,复制一份 ...
- 【sdoi2013】森林 BZOJ 3123
Input 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20. 第二行包含三个整数N,M,T,分别表示节点数.初始边数.操作数.第三行包含N个非负整数 ...
- 【清华集训】楼房重建 BZOJ 2957
Description 小A的楼房外有一大片施工工地,工地上有N栋待建的楼房.每天,这片工地上的房子拆了又建.建了又拆.他经常无聊地看着窗外发呆,数自己能够看到多少栋房子. 为了简化问题,我们考虑这些 ...
- 【splay】文艺平衡树 BZOJ 3223
Description 您需要写一种数据结构(可参考题目标题),来维护一个有序数列,其中需要提供以下操作:翻转一个区间,例如原有序序列是5 4 3 2 1,翻转区间是[2,4]的话,结果是5 2 3 ...
- bzoj 刷水
bzoj 3856: Monster 虽然是sb题,,但是要注意h可能<=a,,,开始忘记判了WA得很开心. #include <iostream> #include <cst ...
随机推荐
- mask rcnn训练自己的数据集参考文章(推荐)
最近用Mask_RCNN训练模型,下面几篇文章提供了不少帮助,汇总出来,方便以后查找,并向几位博主老师表示感谢 https://blog.csdn.net/qq_29462849/article/de ...
- 面试题:Java开发中的23种设计模式详解(转)
设计模式(Design Patterns) ——可复用面向对象软件的基础 设计模式(Design pattern)是一套被反复使用.多数人知晓的.经过分类编目的.代码设计经验的总结.使用设计模式是为了 ...
- c语言练习 二维数组 年平均降水量 月平均降水量
#define YEARS 5#define MONTHES 12 int main(void) { const float rain[YEARS][MONTHES] = { {4.3,4.3,4.3 ...
- 自定义JTabbedPane 标签风格
自定义JTabbedPane 标签风格 2012年03月09日 20:33:12 阅读数:2435 demo 下载地址:http://download.csdn.net/detail/jinannan ...
- IDEA小技巧:添加代码快捷方式
非常怀恋eclipse的的代码快捷方式tryc,今天给IDEA也添加了一个
- hdu 4278 Faulty Odometer(进制转换)
十进制转八进制的变形: #include<stdio.h> int main() { int n; while(scanf("%d",&n)!=EOF& ...
- C#中的异步编程Async 和 Await
谈到C#中的异步编程,离不开Async和Await关键字 谈到异步编程,首先我们就要明白到底什么是异步编程. 平时我们的编程一般都是同步编程,所谓同步编程的意思,和我们平时说的同时做几件事情完全不同. ...
- C/C++预处理指令常见的预处理指令
C/C++预处理指令常见的预处理指令如下: #空指令,无任何效果 #include包含一个源代码文件 #define定义宏 #undef取消已定义的宏 #if如果给定条件为真,则编译下面代码 #ifd ...
- 【Azure Active Directory】单一登录 (SAML 协议)
Azure Active Directory 支持 SAML 2.0 Web 浏览器单一登录 (SSO) 配置文件. 若要请求 Azure Active Directory 对用户进行身份验证时,云服 ...
- UVa 10245 The Closest Pair Problem (分治)
题意:给定 n 个点,求最近两个点的距离. 析:直接求肯定要超时的,利用分治法,先把点分成两大类,答案要么在左边,要么在右边,要么一个点在左边一个点在右边,然后在左边或右边的好求,那么对于一个在左边一 ...