bzoj 2229: [Zjoi2011]最小割
Description
小白在图论课上学到了一个新的概念——最小割,下课后小白在笔记本上写下了如下这段话: “对于一个图,某个对图中结点的划分将图中所有结点分成两个部分,如果结点s,t不在同一个部分中,则称这个划分是关于s,t的割。 对于带权图来说,将所有顶点处在不同部分的边的权值相加所得到的值定义为这个割的容量,而s,t的最小割指的是在关于s,t的割中容量最小的割” 现给定一张无向图,小白有若干个形如“图中有多少对点它们的最小割的容量不超过x呢”的疑问,小蓝虽然很想回答这些问题,但小蓝最近忙着挖木块,于是作为仍然是小蓝的好友,你又有任务了。
解题报告
对于最小割,其实本质不同的只有n-1个,考虑怎么找出这n-1个
可以简单证明:最小割不会相交,可以简单的用反证法证明
考虑一个四格矩形1,2,3,4四个区域(象限),分割线分别为A,B,C,D。
假设1,4的最小割为(A,C) 2,3的最小割为(B,D),那么可以推出 \(A<B\),\(C<D\),所以2,3的最小割不为 (A,C) 产生矛盾,固不相交.
既然不相交,考虑每次找出本质不同的,因为不相交,所以可以考虑分治.
找出当前层的最小割,化为S,T两个集合,然后再选择S,T集合分别做同样的操作,因为最小割不相交,所以找出的一定不同,对于一个点对,在每一个最小割里取Min即可
#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#define RG register
#define il inline
#define iter iterator
#define Max(a,b) ((a)>(b)?(a):(b))
#define Min(a,b) ((a)<(b)?(a):(b))
using namespace std;
const int N=155,inf=2e9,M=3005;
int id[N],n,m,b[N],head[N],nxt[M<<1],to[M<<1],dis[M<<1],num=1,a[N][N],T,S;
void link(int x,int y,int z){nxt[++num]=head[x];to[num]=y;dis[num]=z;head[x]=num;}
void addedge(int x,int y,int z){link(x,y,z);link(y,x,z);}
struct edge{int x,y,dis;}e[M];
int q[N],dep[N];
il bool bfs(){
int x,u,t=0,sum=1;
memset(dep,0,sizeof(dep));
dep[S]=1;q[1]=S;
while(t!=sum){
x=q[++t];
for(int i=head[x];i;i=nxt[i]){
u=to[i];if(dep[u] || dis[i]<=0)continue;
dep[u]=dep[x]+1;q[++sum]=u;
}
}
return dep[T];
}
il int dfs(int x,int flow){
if(x==T || !flow)return flow;
RG int u,tot=0,tmp;
for(int i=head[x];i;i=nxt[i]){
u=to[i];
if(dis[i]<=0 || dep[u]!=dep[x]+1)continue;
tmp=dfs(u,Min(flow,dis[i]));
dis[i]-=tmp;dis[i^1]+=tmp;
tot+=tmp;flow-=tmp;
if(!flow)break;
}
if(!tot)dep[x]=-1;
return tot;
}
il int maxflow(int ss,int tt){
S=ss;T=tt;int tot=0,tmp;
while(bfs()){
tmp=dfs(S,inf);
while(tmp)tot+=tmp,tmp=dfs(S,inf);
}
return tot;
}
void Clear(){memset(head,0,sizeof(head));num=1;}
il void solve(int l,int r){
if(!(l^r))return ;
int L=l,R=r;
Clear();
for(RG int i=1;i<=m;i++)
addedge(e[i].x,e[i].y,e[i].dis);
int tmp=maxflow(id[l],id[r]);
for(int i=1;i<=n;i++)
if(dep[i])for(RG int j=1;j<=n;j++)if(!dep[j])a[i][j]=a[j][i]=Min(a[i][j],tmp);
for(int i=l;i<=r;i++){
if(dep[id[i]])b[L++]=id[i];
else b[R--]=id[i];
}
for(int i=l;i<=r;i++)id[i]=b[i];
solve(l,R);solve(L,r);
}
void work()
{
Clear();
int x,y,z;
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++){
scanf("%d%d%d",&x,&y,&z);
e[i].x=x;e[i].y=y;e[i].dis=z;
}
for(int i=1;i<=n;i++)
for(RG int j=i+1;j<=n;j++)
a[i][j]=a[j][i]=inf;
for(int i=1;i<=n;i++)id[i]=i;
solve(1,n);
int Q;cin>>Q;
while(Q--){
scanf("%d",&x);
int ans=0;
for(int i=1;i<=n;i++)
for(RG int j=i+1;j<=n;j++)
if(a[i][j]<=x)ans++;
printf("%d\n",ans);
}
}
int main()
{
int T;cin>>T;
while(T--){
work();
if(T)puts("");
}
return 0;
}
bzoj 2229: [Zjoi2011]最小割的更多相关文章
- bzoj 2229 [Zjoi2011]最小割(分治+最小割)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=2229 [题意] 回答若干个关于割不超过x的点对数目的询问. [思路] [最小割最多有n ...
- BZOJ.2229.[ZJOI2011]最小割(最小割树)
题目链接 题意:给定一张无向图,求任意两点之间的最小割. 在所有点中任选两个点作为源点\(S\).汇点\(T\),求它们之间的最小割\(ans\),并把原图分成两个点集\(S',T'\),用\(ans ...
- bzoj 2229: [Zjoi2011]最小割【Gomory–Hu tree最小割树】
这个算法详见http://www.cnblogs.com/lokiii/p/8191573.html 求出两两之间最小割之后暴力统计即可 #include<iostream> #inclu ...
- ●BOZJ 2229 [Zjoi2011]最小割
题链: http://www.lydsy.com/JudgeOnline/problem.php?id=2229 题解: 首先先去看看这个博客:http://blog.csdn.net/jyxjyx2 ...
- 2229: [Zjoi2011]最小割(最小割树)
Description 小白在图论课上学到了一个新的概念——最小割,下课后小白在笔记本上写下了如下这段话: “对于一个图,某个对图中结点的划分将图中所有结点分成两个部分,如果结点s,t不在同一个部分中 ...
- bzoj千题计划139:bzoj2229: [Zjoi2011]最小割
http://www.lydsy.com/JudgeOnline/problem.php?id=2229 最小割树介绍:http://blog.csdn.net/jyxjyx27/article/de ...
- BZOJ2229: [Zjoi2011]最小割
题解: 真是一道神题!!! 大家还是围观JZP的题解吧(网址找不到了...) 代码: #include<cstdio> #include<cstdlib> #include&l ...
- bzoj2229: [Zjoi2011]最小割(分治最小割+最小割树思想)
2229: [Zjoi2011]最小割 题目:传送门 题解: 一道非常好的题目啊!!! 蒟蒻的想法:暴力枚举点对跑最小割记录...绝对爆炸啊.... 开始怀疑是不是题目骗人...难道根本不用网络流?? ...
- 【BZOJ2229】[ZJOI2011]最小割(网络流,最小割树)
[BZOJ2229][ZJOI2011]最小割(网络流,最小割树) 题面 BZOJ 洛谷 题解 戳这里 那么实现过程就是任选两点跑最小割更新答案,然后把点集划分为和\(S\)联通以及与\(T\)联通. ...
随机推荐
- djangoueditor 集成xadmin
1.安装Python3兼容版本 https://github.com/twz915/DjangoUeditor3/ 2.model加入字段 from DjangoUeditor.models impo ...
- highcharts 具体参数详解
<script type="text/javascript" src="js/jquery.min.js"></script> < ...
- Winserver+Apache+django部署
废话不多说,干活直接上. winserver2012 + django2.0.1 + apache 部署过程 python ==> 3.4 64位 https://www.python.org/ ...
- 关于移动web教程免费发布
各位老铁大家好,最近经历了太多太多,精力一直不能集中做自己愿意做的事情. 移动Web课程一开始设置收费10块,其实本意是让大家感觉有支出,就会相对珍惜好好学习,但是发现收费把大部分人挡在门外,现在恢复 ...
- JAVAEE——BOS物流项目09:业务受理需求分析、创建表、实现自动分单、数据表格编辑功能使用方法和工作单快速录入
1 学习计划 1.业务受理需求分析 n 业务通知单 n 工单 n 工作单 2.创建业务受理环节的数据表 n 业务通知单 n 工单 n 工作单 3.实现业务受理自动分单 n 在CRM服务端扩展方法根据手 ...
- JAVA_SE基础——60.初识Object
java是面向对象的语言,核心思想:找适合 的对象做适合 的事情:方式一:自定义类,然后通过自定义的类创建对象.方式二:sun提供了很多的类给我使用,我们只需要认识这些类,我们就可以通过这些类创建对象 ...
- GIT入门笔记(10)- 多种撤销修改场景和对策
场景1:当你改乱了工作区某个文件的内容,想直接丢弃工作区的修改时,用命令git checkout -- file. 场景2:当你不但改乱了工作区某个文件的内容,还添加到了暂存区时,想丢弃修改,分两步, ...
- yum 安装Apache
1.查看是否安装Apache,命令: rpm -qa httpd 2.yum install httpd ,yum安装Apache 3.chkconfig httpd on s ...
- Python/MySQL(四、MySQL数据库操作)
Python/MySQL(四.MySQL数据库操作) 一.数据库条件语句: case when id>9 then ture else false 二.三元运算: if(isnull(xx)0, ...
- Hive&SqlServerql:inner join on条件中如果两边都是空值的情况下,关联结果中会把数据给过滤掉。
今天遇到的一个大坑,话不多少,看sql和下边的查询结果: --问题:恰好把buildingid is null的记录给过滤掉 ),buildingid ),)); ); ); ); ); ); ); ...