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\)联通. ...
随机推荐
- 利用yield 实现Xrange功能
def xrange(n): start = 0 while True: if start>n: return yield start start+=1 obj = xrange(5) n1 = ...
- JAVA和Android的回调机制
本文出自xiaanming的博客(http://blog.csdn.net/xiaanming/article/details/17483273),请尊重他人的辛勤劳动成果,谢谢 以 前不理解什么叫回 ...
- ThreadLocal源码分析:(二)get()方法
在ThreadLocal的get(),set()的时候都会清除线程ThreadLocalMap里所有key为null的value. 而ThreadLocal的remove()方法会先将Entry中对k ...
- 老帖收藏,留供参考:SpringMvc2.5+Mybatis3.2.7
一.项目背景 SpringMvc+Mybatis 数据库连接池是阿里巴巴的druid.日志框架式logback 二.配置文件 1.SpringMvc-servlet.xml <?xml vers ...
- Linux实战案例(7)安装jdk
一.文件准备 1.1 文件名称 jdk-8u121-linux-x64.tar.gz 1.2 下载地址 http://www.oracle.com/technetwork/java/javase/do ...
- HashMap就是这么简单【源码剖析】
前言 声明,本文用得是jdk1.8 前面已经讲了Collection的总览和剖析List集合以及散列表.Map集合.红黑树的基础了: Collection总览 List集合就这么简单[源码剖析] Ma ...
- python开发:python字符串操作方法
name = "my \tname is {name} and i am {year} old" capitalize:第一个单词的首字母大写的方法 print(name.capi ...
- Web 项目报错No suitable driver found for jdbc:mysql://localhost:3306/book 的一个解决办法
确认jar包加入到了build path中,然后注意版本是否与数据库相配,还要留意将jar包放入WEB-INF下的lib文件夹中
- hive:默认允许动态分区个数为100,超出抛出异常:
在创建好一个分区表后,执行动态分区插入数据,抛出了错误: Caused by: org.apache.hadoop.hive.ql.metadata.HiveFatalException: [Erro ...
- 1.0 添加WEB API项目并按注释生成文档(多项目结构)
1.新建ASP.NET 项目,模板选择如图 2.选择Web API,并选择不进行身份验证方式 成功后我们看到这个结果. 至于其它三种身份验证方式,不太适合我的使用.而且这种方式也可以在代码里去实现身份 ...