【BZOJ3218】a+b problem (最小割 + 主席树)

传送门
在这里插入图片描述
继续优化:把a[ ]离散化

#include<bits/stdc++.h>
using namespace std;
const int inf=1000000007;
const int N=200010;
const int M=1000010;struct Edge{int u,v,f,next;
}edge[M];
int head[N],cnt;
int s,t,flow,level[N];struct Node{int x,id;
}e[5010];
int L[5010],R[5010],num,key[5010],seq[5010];int tot,ls[N],rs[N],rt[N],Cnt[N];int n,ans;void add(int u,int v,int f){edge[cnt].u=u;edge[cnt].v=v;edge[cnt].f=f;edge[cnt].next=head[u];head[u]=cnt++;edge[cnt].u=v;edge[cnt].v=u;edge[cnt].f=0;edge[cnt].next=head[v];head[v]=cnt++;
} 
bool bfs(){queue<int> q;memset(level,-1,sizeof(level));level[s]=0;q.push(s);while(!q.empty()){int u=q.front();q.pop();for(int i=head[u];i!=-1;i=edge[i].next){int v=edge[i].v;int f=edge[i].f;if(f>0&&level[v]<0){level[v]=level[u]+1;q.push(v);}}}return level[t]>0;
}
int dfs(int u,int mx){if(u==t||!mx) return mx;int tmp,res=0;for(int i=head[u];i!=-1&&mx;i=edge[i].next){int v=edge[i].v;int f=edge[i].f;if(level[v]==level[u]+1&&(tmp=dfs(v,min(mx,f)))){edge[i].f-=tmp;edge[i^1].f+=tmp;res+=tmp;mx-=tmp;}}if(!res)level[u]=-1;return res;
}
int dinic(){int maxflow=0;while(bfs()){maxflow+=dfs(s,inf);}return maxflow;
}bool cmp(Node x,Node y){return x.x<y.x;
}
int findl(int x){int l=1,r=num,ans=num;while(l<=r){int mid=(l+r)/2;if(seq[mid]>=x) ans=mid,r=mid-1; else l=mid+1;}return ans;
}
int findr(int x){int l=1,r=num,ans=1;while(l<=r){int mid=(l+r)/2;if(seq[mid]<=x) ans=mid,l=mid+1; else r=mid-1;}return ans;
}
int modify(int pre,int l,int r,int i){int now=++tot;if(l==r){ls[now]=rs[now]=0;Cnt[now]=Cnt[pre]+1;}else{int mid=(l+r)/2;if(key[i]<=mid){ls[now]=modify(ls[pre],l,mid,i);rs[now]=rs[pre];}else{rs[now]=modify(rs[pre],mid+1,r,i);ls[now]=ls[pre];}Cnt[now]=Cnt[ls[now]]+Cnt[rs[now]];}add(i,2*n+1+now,inf);add(2*n+1+pre,2*n+1+now,inf);return now;
}
void query(int o,int l,int r,int i){if (!Cnt[o]) return;if (L[i]<=l&&r<=R[i]){add(2*n+1+o,n+i,inf);return;}int mid=(l+r)/2;if(L[i]<=mid) query(ls[o],l,mid,i);if(mid<R[i]) query(rs[o],mid+1,r,i);
}
int main(){cnt=0;memset(head,-1,sizeof(head));scanf("%d",&n);s=0;t=100000;for(int i=1;i<=n;i++){int b,w,p;scanf("%d%d%d%d%d%d",&e[i].x,&b,&w,&L[i],&R[i],&p);e[i].id=i;ans+=b+w;add(s,i,w);add(i,t,b);add(n+i,i,p);}sort(e+1,e+n+1,cmp);num=1;key[e[1].id]=1;seq[1]=e[1].x;for(int i=2;i<=n;i++)if(e[i].x==e[i-1].x) key[e[i].id]=num;else key[e[i].id]=++num,seq[num]=e[i].x;tot=0;rt[0]=ls[0]=rs[0]=0;for(int i=1;i<=n;i++){L[i]=findl(L[i]);R[i]=findr(R[i]);rt[i]=modify(rt[i-1],1,num,i);query(rt[i-1],1,num,i);     }ans-=dinic();printf("%d\n",ans); return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/320325.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

ST表讲解

ST表主要用于解决RMQ问题&#xff08;区间最值问题&#xff09; 当然你可以用线段树等&#xff0c;但今天用一种ST表&#xff08;倍增算法&#xff09; ST表是倍增算法的一个典型应用 暴力做RMQ问题&#xff0c;往往会超时&#xff0c;ST表利用对其进行优化 给定一段序列A&am…

P3346-[ZJOI2015]诸神眷顾的幻想乡【广义SAM】

正题 题目链接:https://www.luogu.com.cn/problem/P3346 题目大意 一棵树&#xff0c;求树上所有路径构成的字符串有多少种。&#xff08;叶子不超过303030个&#xff09; 解题思路 如果是根节点到一些节点的路径的话很好做&#xff0c;直接建广义SAMSAMSAM即可&#xff0c;…

纪中B组模拟赛总结(2020.2.1)

成绩 rankrankranknamenamenamescorescorescoreT1T1T1T2T2T2T3T3T3T4T4T4414141lyflyflyf909090000909090000000 总结 T1T1T1直接爆炸&#xff0c;没想到是数学问题 T2T2T2想了出来&#xff0c;但忘记开mapmapmap了 T3T3T3直接懵逼 T4T4T4没有太了解欧拉函数导致翻车 题解 T1…

微软技术直通车(第三期)现场实录

微软技术直通车&#xff08;第三期&#xff09;于本月18日14时&#xff0c;在中国微软总部圆满举行。本次活动得到了微软最佳有价值专家&#xff08;Microsoft MVP&#xff09;中国区项目组的鼎力支持。同时&#xff0c;Microsoft MVP中国区项目组负责人也亲临现场致辞&#xf…

[NOI2019] 序列(模拟费用流)

原先自己想的建图&#xff1a; 正确建图&#xff1a; 但是 n 太大了&#xff0c;所以考虑模拟费用流&#xff1a; 注意&#xff1a; 在1中&#xff0c; 若选的两个位置相同&#xff0c;则为情况2&#xff0c;不用减 f&#xff1b; 若选的位置在另一序列中已被选&#xff0c;…

【结论】游戏(jzoj 5536)

游戏 jzoj 5536 题目大意&#xff1a; 给一个范围和一个矩形&#xff0c;让你在这个范围内放若干个这样的矩形&#xff08;不能改变方向&#xff0c;不能重叠&#xff09;&#xff0c;让你求出最少放多个矩形可以使范围内无法再放矩形 输入样例#1 11 4 3 2输入样例#2 10 …

[2020.11.26NOIP模拟赛]询问【字符串hash】

正题 题目链接:https://www.luogu.com.cn/problem/U142342?contestId37784 题目大意 一个字符串&#xff0c;定义两个字符串相似为用一些字母代替相同的字母后可以相同。 如urbbrgroorurbbrgroorurbbrgroor&#xff0c;apple≠abcdeapple\neq abcdeapple​abcde 要求支持…

LCA总结

文章目录LCA介绍解决方法概括&#xff1a;倍增法&#xff1a;Tarjan欧拉序树剖解法&#xff1a;看了很多关于LCA的文章&#xff0c;这里是一个总结我们通过这个题洛谷P3379 【模板】最近公共祖先来讲LCA LCA介绍 lca是啥&#xff1f;最近公共祖先 就是&#xff1a;两个点在这…

ASP.NET Core Web API + Ng6 实战视频 Day 2

第一天课程&#xff1a;ASP.NET Core Web API Identity Server 4 Angular 6 实战小项目视频Day 2 第一部分:Day 2 第二部分:视频专辑持续更新中....地址请点击原文链接.原文地址: http://v.qq.com/vplus/4cfb00af75c16eb8d198c58fb86eb4dc/foldervideos/8hk0029019k2fft.NET…

浮沉子

浮沉子制作及其原理 浮沉子的下沉过程&#xff1a;

纪中B组模拟赛总结(2020.2.3)

成绩 rankrankranknamenamenamescorescorescoreT1T1T1T2T2T2T3T3T3333lyflyflyf252.6252.6252.610010010010010010052.652.652.6 总结 T1以前做过&#xff0c;十分嗨皮,就直接切掉了 T2发现是一道傻逼题,不多说 T3比赛时理解题意错了&#xff0c;以为一个字符串属于一个组别就…

[2020.11.26NOIP模拟赛]勇者的后缀【SA,RMQ,主席树,二分】

正题 题目链接:https://www.luogu.com.cn/problem/U142356?contestId37784 题目大意 一个字符串&#xff0c;询问给出(x,l,r)(x,l,r)(x,l,r)表示询问在[l,r][l,r][l,r]中作为起点找一个后缀它与xxx作为起点的后缀的LCPLCPLCP最长&#xff0c;且满足最长的情况下字典序最小。…

.NET Core 项目指定SDK版本

一. 版本里的坑自从 .NET Core 2.1.0版本发布以后&#xff0c;近几个月微软又进行了几次小版本的发布&#xff0c;可见 .NET Core 是一门生命力非常活跃的技术。经过一段时间的实践&#xff0c;目前做 ASP.NET Core 开发时&#xff0c;使用的 Nuget 包&#xff0c;比如 Microso…

牛客网【每日一题】4月17日题目精讲 华华给月月准备礼物

文章目录题目描述题解&#xff1a;代码&#xff1a;推荐例题&#xff1a;试题链接时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 32768K&#xff0c;其他语言65536K 64bit IO Format: %lld 题目描述 二月中旬虐狗节前夕&#xff0c;华华决定给…

[BZOJ 3173] [TJOI 2013] 最长上升子序列(splay)

根据题意很容易得到最原始的做法&#xff1a; #include<iostream> #include<cstdio> #include<cstring> #include<vector> using namespace std; const int N100005; int n,m,f[N]; int main(){scanf("%d",&n);int x;for(int i1;i<n…

P4300-[AHOI2006]上学路线【网络流,最短路】

正题 题目链接:https://www.luogu.com.cn/problem/P4300 题目大意 nnn个点mmm条边的无向图。求1∼n1\sim n1∼n的最短路和删除cic_ici​和最小的边使得最短路变长。 解题思路 显然我们需要跑一次最短路。 之后考虑如何求第二问&#xff0c;我们发现我们要割掉最短路上的边&a…

Mathematics(jzoj 1747)

Mathematics jzoj 1747 题目大意 有n堆石子&#xff0c;总和为2k2^k2k&#xff0c;现在对于两堆石子&#xff0c;你可以从a中取b的分量到b&#xff08;a要大于b&#xff09;&#xff0c;问合成一堆大小为2k2^k2k的石子要怎么做&#xff08;输出a,b&#xff09; 输入样例 2…

Bing.com在.NET Core 2.1上运行

相关知识请参考.netCore开发团队博客 https://blogs.msdn.microsoft.com/dotnet/ Bing.com是一种云服务&#xff0c;运行在遍布全球许多数据中心的数千台服务器上。Bing服务器每秒处理来自全球消费者的数千个用户查询&#xff0c;通过他们的浏览器&#xff0c;使用Microsoft认…

Jelly

链接&#xff1a; 时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 131072K&#xff0c;其他语言262144K 64bit IO Format:%lld 题目描述 Nancy喜欢吃果冻&#xff01; Nancy钻进了一个n \times n \times nnnn的果冻里&#xff0c;她想从(1,1,1)一…

[bzoj 4811] 由乃的OJ(贪心 + 树链剖分)

前置技能&#xff1a;[Noi2014]起床困难综合症。 不难看出&#xff0c;这道题其实就是上一道题的加强版 在上一道题中&#xff0c;因为位运算时位与位之间互不干扰 所以从高位到低位枚举初始值二进制上的每一位为0和为1时&#xff0c;经过n次计算后这一位的结果&#xff0c;贪…