CF 704 D. Captain America
CF 704 D. Captain America
题目大意:给出\(n\)个点的坐标,你要将每个点染成红色或者蓝色。染一个红色要付出\(r\)的代价,染一个蓝色要付出\(b\)的代价。有\(m\)个限制,每个限制为“第\(i\)行(或者第\(i\)列)两种颜色的球数量差不能超过\(d\)”。
对于第\(i\)行,有\(x\)个点,数量差不能超过\(d\),假设最终红球的数量为\(R\),则:
(x-R)-R\leq d\\
\]
得到:
\]
这就是个上下界网络流。将每一行,每一列分别建成一个点。源点向行点连边,列点向汇点连边。对于给定的交叉点上的行列之间连边,容量为\(1\)。先跑有源汇可行流,如果\(r<b\),就跑最大可行流,否则跑最小可行流。
记得当前弧优化。
代码:
#include<bits/stdc++.h>
#define ll long long
#define N 200005
using namespace std;
inline int Get() {int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9') {if(ch=='-') f=-1;ch=getchar();}while('0'<=ch&&ch<='9') {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}return x*f;}
int n,m;
ll r,b;
int x[N],y[N];
int t[N],l[N],d[N];
struct road {
int to,nxt;
int f;
}s[N<<4];
int h[N<<1],cnt=1;
void add(int i,int j,int f) {
if(!f) return ;
s[++cnt]=(road) {j,h[i],f};h[i]=cnt;
s[++cnt]=(road) {i,h[j],0};h[j]=cnt;
}
vector<int>px,py;
map<int,int>mpx,mpy;
int S,T;
int SS,TT;
ll ans;
int mnx[N],mny[N];
int totx[N],toty[N];
int nx,ny;
int ux[N],uy[N];
int dx[N],dy[N];
int dis[N];
queue<int>q;
bool dinic_bfs(int S,int T) {
while(!q.empty()) q.pop();
memset(dis,0x3f,sizeof(dis));
q.push(S);
dis[S]=0;
while(!q.empty()) {
int v=q.front();q.pop();
for(int i=h[v];i;i=s[i].nxt) {
int to=s[i].to;
if(s[i].f&&dis[to]>dis[v]+1) {
dis[to]=dis[v]+1;
if(to==T) return 1;
q.push(to);
}
}
}
return 0;
}
int dep;
bool ins[N];
int cur[N<<1];
int dfs(int v,int T,int maxf) {
if(v==T) return maxf;
ins[v]=1;
int ret=0;
for(int &i=h[v];i;i=s[i].nxt) {
int to=s[i].to;
if(s[i].f&&dis[to]==dis[v]+1&&!ins[to]) {
int dlt=dfs(to,T,min(maxf,s[i].f));
if(!dlt) dis[to]=0;
s[i].f-=dlt;
s[i^1].f+=dlt;
ret+=dlt;
maxf-=dlt;
if(!maxf) break;
}
}
ins[v]=0;
return ret;
}
int dinic(int S,int T) {
int ans=0;
int dep=0;
while(dinic_bfs(S,T)) {
memcpy(cur,h,sizeof(h));
ans+=dfs(S,T,1e9);
memcpy(h,cur,sizeof(h));
}
return ans;
}
int e[N];
bool tag[N<<1];
void dfs(int v) {
tag[v]=1;
for(int i=h[v];i;i=s[i].nxt) {
int to=s[i].to;
if(s[i].f&&!tag[to]) dfs(to);
}
}
int main() {
n=Get(),m=Get();
r=Get(),b=Get();
for(int i=1;i<=n;i++) {
x[i]=Get(),y[i]=Get();
px.push_back(x[i]);
py.push_back(y[i]);
mpx[x[i]]=1;
mpy[y[i]]=1;
}
sort(px.begin(),px.end());
sort(py.begin(),py.end());
px.resize(unique(px.begin(),px.end())-px.begin());
py.resize(unique(py.begin(),py.end())-py.begin());
for(int i=1;i<=n;i++) {
x[i]=lower_bound(px.begin(),px.end(),x[i])-px.begin();
y[i]=lower_bound(py.begin(),py.end(),y[i])-py.begin();
totx[x[i]]++;
toty[y[i]]++;
}
nx=px.size(),ny=py.size();
for(int i=0;i<nx;i++) mnx[i]=totx[i];
for(int i=0;i<ny;i++) mny[i]=toty[i];
for(int i=1;i<=m;i++) {
t[i]=Get(),l[i]=Get(),d[i]=Get();
if(t[i]==1) {
if(mpx.find(l[i])==mpx.end()) continue ;
l[i]=lower_bound(px.begin(),px.end(),l[i])-px.begin();
mnx[l[i]]=min(mnx[l[i]],d[i]);
} else {
if(mpy.find(l[i])==mpy.end()) continue ;
l[i]=lower_bound(py.begin(),py.end(),l[i])-py.begin();
mny[l[i]]=min(mny[l[i]],d[i]);
}
}
S=nx+ny,T=nx+ny+1;
SS=nx+ny+2,TT=nx+ny+3;
int E=cnt+1;
add(T,S,1e9);
int tot1=0,tot2=0;
for(int i=0;i<nx;i++) {
dx[i]=totx[i]-mnx[i]+1>>1;
ux[i]=totx[i]+mnx[i]>>1;
tot1+=dx[i];
add(S,i,ux[i]-dx[i]);
if(dx[i]) add(SS,i,dx[i]);
if(dx[i]>ux[i]) {
cout<<-1;
return 0;
}
}
for(int i=0;i<ny;i++) {
dy[i]=toty[i]-mny[i]+1>>1;
uy[i]=toty[i]+mny[i]>>1;
tot2+=dy[i];
add(i+nx,T,uy[i]-dy[i]);
if(dy[i]) add(i+nx,TT,dy[i]);
if(dy[i]>uy[i]) {
cout<<-1;
return 0;
}
}
for(int i=1;i<=n;i++) {
e[i]=cnt+1;
add(x[i],y[i]+nx,1);
}
add(S,TT,tot1),add(SS,T,tot2);
int x=dinic(SS,TT);
if(x!=tot1+tot2) {
cout<<-1;
return 0;
}
s[E].f=s[E^1].f=0;
int x2;
if(r>b) {
x2=dinic(S,T);
} else {
dinic(T,S);
}
for(int i=1;i<=n;i++) {
if(s[e[i]].f==0) {
ans+=b;
} else {
ans+=r;
}
}
cout<<ans<<"\n";
for(int i=1;i<=n;i++) {
if(s[e[i]].f==0) {
cout<<"b";
} else {
cout<<"r";
}
}
return 0;
}
CF 704 D. Captain America的更多相关文章
- CF#366 704D Captain America 上下界网络流
CF上的题,就不放链接了,打开太慢,直接上题面吧: 平面上有n个点, 第 i 个点的坐标为 ($X_i ,Y_i$), 你需要把每个点染成红色或者蓝色, 染成红色的花费为 r , 染成蓝色的花费为 b ...
- 【CF704D】Captain America(上下界网络流)
[CF704D]Captain America(上下界网络流) 题面 CF 洛谷 题解 如果没有限制,似乎就不用做了...因为我们只需要贪心的选择代价较小的颜色就行了. 那么我们不妨假设染红色的代价较 ...
- CF704D Captain America
http://codeforces.com/problemset/problem/704/D 题解 对于两种颜色的染色,我们可以把它看做选择问题. 比如说红色的代价小,所以我们尽可能多的染红色. 然后 ...
- Codeforces 704D Captain America
题意:平面上有n个点,每个点必须涂成红色和蓝色中的一种,花费各为r和b(对所有的点花费都一样).m条限制,每条限制形如"y=b这条直线上两种颜色的点的数目之差的绝对值不能超过c"或 ...
- codeforces704D Captain America【上下界最大流】
分别给行和列hash建两排点,对(x,y)坐标连x行y列的点 设红色价格低,那么就要尽量多选红色 设一个点出度为s,要求最小的最大差值为d,又,假设有流量表示选红没流量表示选蓝,那么要求就变成了这个点 ...
- CF704D Captain America(上下界网络流)
传送门 题意: 二维平面给出\(n\)个点,现在可以给每个点进行染色,染红色的代价为\(r\),染蓝色的代价为\(b\). 之后会有\(m\)个限制,形式如:\(t_i\ l_i\ d_i\),当\( ...
- CF704D Captain America 上下界网络流
传送门 现在相当于说每一个条件都有一个染成红色的盾牌的数量限制\([l,r]\),需要满足所有限制且染成红色的盾牌数量最小/最大. 注意到一个盾牌染成红色对于一行和一列都会产生影响.如果选中一个物品对 ...
- Codeforces VP/补题小记 (持续填坑)
Codeforces VP/补题小记 1149 C. Tree Generator 给你一棵树的括号序列,每次交换两个括号,维护每次交换之后的直径. 考虑括号序列维护树的路径信息和,是将左括号看做 ...
- 五、Pandas玩转数据
Series的简单运算 import numpy as np import pandas as pd s1=pd.Series([1,2,3],index=['A','B','C']) print(s ...
随机推荐
- VS2017创建的单元测试不支持顺序测试
问题:使用IDE创建的单元测试项目,标准引用是,导致不能添加顺序测试,复制其它项目的顺序测试文件进行编辑时,也会提示基于MSTest V2的测试不能用于顺序测试 解决办法: 移除自带的NU ...
- [ERR] Node 172.16.6.154:7002 is not empty. Either the node already knows other nodes (check with CLUSTER NODES) or contains some key in database 0.
关于启动redis集群时: [ERR] Node 172.168.63.202:7001 is not empty. Either the nodealready knows other nodes ...
- pytorch 中的Variable一般常用的使用方法
Variable一般的初始化方法,默认是不求梯度的 import torch from torch.autograd import Variable x_tensor = torch.randn(2, ...
- 【转】开发一个这样的 APP 要多长时间?
作者:蒋国刚 www.cnblogs.com/guogangj/p/4676836.html 呵呵. 这是一个“如有雷同,纯属巧合”的故事,外加一些废话,大家请勿对号入座.开始了…… 我有些尴尬地拿着 ...
- vue 中的小知识点
1)使用is解决小bug <!DOCTYPE html><html lang="en"> <head> <meta charset=&qu ...
- 剑指Offer-38.平衡二叉树(C++/Java)
题目: 输入一棵二叉树,判断该二叉树是否是平衡二叉树. 分析: 可以从根节点开始遍历每一个节点,求得节点左右子树的最大高度,判断是不是平衡二叉树.这样做的问题在于会重复遍历节点,造成不必要的浪费. 所 ...
- 安装swoole出现make报错的原因与解决方法
安装swoole报错,错误信息如下: 报错原因 报错很明显,找不到 openssl/ssl.h ,首先要确认电脑上有没有安装 openssl Mac安装openssl 既然它找不到,那么就需要手动的指 ...
- leetcode-字符串篇
Implement strStr() /** * Implement strStr(). * * Return the index of the first occurrence of needle ...
- Python 周刊第 418 期
新闻 PyCon US 2020 开始接受财务赞助! https://pycon.blogspot.com/2019/10/financial-aid-launches-for-pycon-us-20 ...
- Glide缓存流程
本文首发于 vivo互联网技术 微信公众号 链接:https://mp.weixin.qq.com/s/cPLkefpEb3w12-uoiqzTig作者:连凌能 Android上图片加载的解决方案有多 ...