BZOJ3218: a + b Problem
题解:
先做60分。。。
考虑最小割,连边容量为需要付出的代价。不妨设在s割为黑色,t割为白色。
(s,i,b[i])(i,t,w[i])
关于奇怪,因为不是按份数来的。所以我们这样建图:
(i,i+n,p[i])(i+n,j,inf) l[i]<=a[j]<=r[i]
代表只要有一个j属于t割,那么i+n就会属于t割,而如果i属于s割,就会付出p[i]的代价。
注意:(x,y,inf)表示x在s割,那么y一定在s割 或者说 y在t割,那么x一定在t割。而有可能出现x在t割,而y在s割的情况。这证实了上面算法的正确性。
代码:
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<string>
#define inf 1000000000
#define maxn 1000000+5
#define maxm 1000000+5
#define eps 1e-10
#define ll long long
#define pa pair<int,int>
#define for0(i,n) for(int i=0;i<=(n);i++)
#define for1(i,n) for(int i=1;i<=(n);i++)
#define for2(i,x,y) for(int i=(x);i<=(y);i++)
#define for3(i,x,y) for(int i=(x);i>=(y);i--)
#define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go)
#define mod 1000000007
using namespace std;
inline int read()
{
int x=,f=;char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}
return x*f;
}
int a[maxn],b[maxn],w[maxn],l[maxn],r[maxn],p[maxn];
int n,m,s,t,sum,maxflow,tot=,head[maxn],cur[maxn],h[maxn];
queue<int>q;
struct edge{int go,next,v;}e[maxm];
void add(int x,int y,int v)
{
e[++tot]=(edge){y,head[x],v};head[x]=tot;
e[++tot]=(edge){x,head[y],};head[y]=tot;
}
bool bfs()
{
for(int i=s;i<=t;i++)h[i]=-;
q.push(s);h[s]=;
while(!q.empty())
{
int x=q.front();q.pop();
for(int i=head[x];i;i=e[i].next)
if(e[i].v&&h[e[i].go]==-)
{
h[e[i].go]=h[x]+;q.push(e[i].go);
}
}
return h[t]!=-;
}
int dfs(int x,int f)
{
if(x==t) return f;
int tmp,used=;
for(int i=cur[x];i;i=e[i].next)
if(e[i].v&&h[e[i].go]==h[x]+)
{
tmp=dfs(e[i].go,min(e[i].v,f-used));
e[i].v-=tmp;if(e[i].v)cur[x]=i;
e[i^].v+=tmp;used+=tmp;
if(used==f)return f;
}
if(!used) h[x]=-;
return used;
}
void dinic()
{
maxflow=;
while(bfs())
{
for (int i=s;i<=t;i++)cur[i]=head[i];maxflow+=dfs(s,inf);
}
}
int main()
{
freopen("input.txt","r",stdin);
freopen("output.txt","w",stdout);
n=read();s=;t=*n+;
for1(i,n)a[i]=read(),b[i]=read(),w[i]=read(),l[i]=read(),r[i]=read(),p[i]=read(),sum+=b[i]+w[i];
for1(i,n)
{
add(s,i,b[i]);add(i,t,w[i]);add(i,i+n,p[i]);
for1(j,i-)if(a[j]>=l[i]&&a[j]<=r[i])add(i+n,j,inf);
}
dinic();
cout<<sum-maxflow<<endl;
return ;
}
现在考虑满分做法:
奇怪的格子的约束条件是:存在j<i,且l[i]<=a[j]<=r[i]。
首先我们先忽略j<i。
那我们可以先把所有的a[i]插入一棵线段树中,由每个i+n像[l[i],r[i]]所包含的区间连边inf。
然后每个包含i的区间向 i 连边。
这样就实现了 s-> i -> i+n -> 线段树的节点(表示区间) -> j -> t
只不过在原来的算法上多转了几个点。
然后考虑限制:j<i
我们想到对每个i,建一棵1-i的线段树,然后执行上面的算法。
我们想到了可持久化线段树。
这样就可以解决了。
需要注意的细节:
1.离散化
2.
如果出现a[i]相同的情况,我们需不需要向每个i连边?
不需要,如果在前一个版本a[i]出现过,记为last.对于当前版本出现的a[i],记为now
(now,last,INF)
这样子,我们可以实现i’通过线段树对所有满足0<j<i,l <= a[j] <= r的点的控制。--谢图图
代码:真是道好题+神题!!!
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<string>
#define inf 2000000000
#define maxn 500000+5
#define maxm 500000+5
#define eps 1e-10
#define ll long long
#define pa pair<int,int>
#define for0(i,n) for(int i=0;i<=(n);i++)
#define for1(i,n) for(int i=1;i<=(n);i++)
#define for2(i,x,y) for(int i=(x);i<=(y);i++)
#define for3(i,x,y) for(int i=(x);i>=(y);i--)
#define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go)
#define mod 1000000007
using namespace std;
inline int read()
{
int x=,f=;char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}
return x*f;
}
int a[maxn],b[maxn],w[maxn],l[maxn],r[maxn],p[maxn],rt[maxn],ss[maxm],ls[maxm],rs[maxm];
int n,m,s,t,cnt,maxflow,sum,tot=,head[maxn],cur[maxn],h[maxn];
queue<int>q;
struct edge{int go,next,v;}e[maxm];
void add(int x,int y,int v)
{
e[++tot]=(edge){y,head[x],v};head[x]=tot;
e[++tot]=(edge){x,head[y],};head[y]=tot;
}
bool bfs()
{
for(int i=;i<=cnt;i++)h[i]=-;
q.push(s);h[s]=;
while(!q.empty())
{
int x=q.front();q.pop();
for(int i=head[x];i;i=e[i].next)
if(e[i].v&&h[e[i].go]==-)
{
h[e[i].go]=h[x]+;q.push(e[i].go);
}
}
return h[t]!=-;
}
int dfs(int x,int f)
{
if(x==t) return f;
int tmp,used=;
for(int i=cur[x];i;i=e[i].next)
if(e[i].v&&h[e[i].go]==h[x]+)
{
tmp=dfs(e[i].go,min(e[i].v,f-used));
e[i].v-=tmp;if(e[i].v)cur[x]=i;
e[i^].v+=tmp;used+=tmp;
if(used==f)return f;
}
if(!used) h[x]=-;
return used;
}
void dinic()
{
maxflow=;
while(bfs())
{
for (int i=;i<=cnt;i++)cur[i]=head[i];maxflow+=dfs(s,inf);
}
}
inline void update(int l,int r,int x,int &y,int z,int p)
{
y=++cnt;
ss[y]=ss[x]+;
add(y,p,inf);
if(x)add(y,x,inf);
if(l==r)return;
ls[y]=ls[x];rs[y]=rs[x];
int mid=(l+r)>>;
if(z<=mid)update(l,mid,ls[x],ls[y],z,p);
else update(mid+,r,rs[x],rs[y],z,p);
}
inline void query(int k,int l,int r,int x,int y,int z)
{
if(!k)return;
if(l==x&&r==y){add(z,k,inf);return;}
int mid=(l+r)>>;
if(y<=mid)query(ls[k],l,mid,x,y,z);
else if(x>mid)query(rs[k],mid+,r,x,y,z);
else query(ls[k],l,mid,x,mid,z),query(rs[k],mid+,r,mid+,y,z);
}
int main()
{
freopen("input.txt","r",stdin);
freopen("output.txt","w",stdout);
n=read();s=;t=*n+;cnt=t;
for1(i,n)
{
a[i]=b[i]=read();
int x=read(),y=read();
l[i]=read();r[i]=read();
int z=read();
sum+=x+y;
add(s,i,x);add(i,t,y);
add(i,i+n,z);
}
sort(b+,b+n+);
for1(i,n)
{
int x=lower_bound(b+,b+n+,l[i])-b;
int y=upper_bound(b+,b+n+,r[i])-b-;
int z=lower_bound(b+,b+n+,a[i])-b;
if(x<=y)query(rt[i-],,n,x,y,i+n);
update(,n,rt[i-],rt[i],z,i);
}
dinic();
cout<<sum-maxflow<<endl;
return ;
}
3218: a + b Problem
Time Limit: 20 Sec Memory Limit: 40 MB
Submit: 426 Solved: 168
[Submit][Status]
Description
BZOJ3218: a + b Problem的更多相关文章
- [bzoj3218]a + b Problem 网络流+主席树优化建图
3218: a + b Problem Time Limit: 20 Sec Memory Limit: 40 MBSubmit: 2229 Solved: 836[Submit][Status] ...
- [bzoj3218] a+b problem [最小割+数据结构优化建图]
题面 传送门 思路 最小割 我们首先忽略掉那个奇♂怪的限制,就有一个比较显然的最小割模型: 建立源点$S$和汇点$T$ 对于每个元素$i$建立一个点$i$,连边$<S,i,w[i]>$和$ ...
- bzoj3218 a + b Problem(网络流+主席树)
$ans=\sum_{color_i=black}\ b_i+\sum_{color_i=white}\ w_i-\sum_{i=abnormal}\ p_i$ 把它转化一下 $ans=\sum_{i ...
- bzoj3218 a+b Problem(最小割+主席树优化建边)
由于6.22博主要学测,大半时间学文化课,近期刷题量&写题解的数量会急剧下降. 这题出得挺经典的,首先一眼最小割,考虑朴素的做法:与S联通表示白色,与T联通表示黑色,S向i连流量为w[i]的边 ...
- OI刷题录——hahalidaxin
16-3-25 —— bzoj 2049 [Sdoi2008]Cave 洞穴勘测:LCT入门 bzoj 2002 [Hnoi2010]Bounce 弹飞绵羊:LCT Tsinsen A1303. t ...
- 【bzoj3218】 a + b Problem
http://www.lydsy.com/JudgeOnline/problem.php?id=3218 (题目链接) 题意 给${n}$个格子涂白或黑色,白则${w_i}$,黑则${b_i}$的好看 ...
- 【BZOJ3218】a + b Problem 可持久化线段树优化建图
[BZOJ3218]a + b Problem 题解:思路很简单,直接最小割.S->i,容量为Bi:i->T,容量为Wi:所有符合条件的j->new,容量inf:new->i, ...
- 【BZOJ-3218】a+b Problem 最小割 + 可持久化线段树
3218: a + b Problem Time Limit: 20 Sec Memory Limit: 40 MBSubmit: 1320 Solved: 498[Submit][Status] ...
- BZOJ3218 UOJ#77 A+B Problem(最小割+主席树)
竟然在BZOJ上拿了Rank1太给力啦. p.s.:汗,一发这个就被一堆人在2月27号强势打脸-- 传送门(BZOJ) 传送门(UOJ) 说说这道题目吧: 首先是说说这个构图吧.因为有选择关系,我们很 ...
随机推荐
- [转载]非常完善的Log4net详细说明
前言 此篇文章是我见过写得最好的一片关于Log4Net的文章,内容由简入难,而且面面俱到,堪称入门和精通的佳作,特从懒惰的肥兔的转载过来. 1.概述 log4net是.Net下一个非常优秀的开源日志记 ...
- 公众号开发学习Day01
登录https://mp.weixin.qq.com/后点击注册进入注册微信公众号界面,使用邮箱进行注册 注意点1:一个邮箱只能创建一个公众号,并且公众号的三种类型(订阅号,服务号,企业号)只能选一种 ...
- JS读取UserAgent信息并做判断
JS读取UserAgent信息并做判断 userAgent信息可以由navigator.userAgent拿到.例子: <script type="text/javascript&qu ...
- JavaScript高级---桥模式设计
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01//EN" "http://www.w3.org/TR/html4/stri ...
- [C/CPP系列知识] 那些程序C语言可以编译通过但C++无法编译成功 Write a C program that won’t compile in C++
http://www.geeksforgeeks.org/write-c-program-wont-compiler-c/ 1) C++中在函数声明之前调用一个函数会引发错误,但是在C中有可能可以. ...
- C++代码反汇编后的堆栈寄存器EBP和ESP
最近在分析一个进程崩溃的严重问题,其中有些过程分析需要对ebp, esp 有清晰的理解,对于ebp 和esp 相信大家都很熟悉了,但是为了使本文自成体系,我还是解释一下. ebp--栈底指针 esp- ...
- 原 Linux搭建SVN 服务器
原 Linux搭建SVN 服务器 发表于1年前(2014-08-05 17:55) 阅读(12257) | 评论(3) 31人收藏此文章, 我要收藏 赞3 摘要 Linux搭建SVN 服务器 目录 ...
- selenium2 安装、简单使用及浏览器启动问题解决汇总
一.selenium2 安装 1.安装jdk并配置环境变量 jdk需要1.6版本及以上的,这个从网上下载就可以,环境变量的配置前边的随笔整理过就不多说了. 2.安装Firefox,Selenium I ...
- eclipse怎么集成配置JDK和Tomcat
eclipse怎么集成配置JDK和Tomcat_电脑软件_百度经验http://jingyan.baidu.com/album/db55b609ab6a1f4ba30a2f2c.html?picind ...
- GuessNum
import java.util.Scanner;import java.util.Random;/***猜数字,使用随机一个0-100的数字,然后用户猜,猜中了就提示猜中了,*否则提示猜大了还是猜小 ...