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]最小割的更多相关文章

  1. bzoj 2229 [Zjoi2011]最小割(分治+最小割)

    [题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=2229 [题意] 回答若干个关于割不超过x的点对数目的询问. [思路] [最小割最多有n ...

  2. BZOJ.2229.[ZJOI2011]最小割(最小割树)

    题目链接 题意:给定一张无向图,求任意两点之间的最小割. 在所有点中任选两个点作为源点\(S\).汇点\(T\),求它们之间的最小割\(ans\),并把原图分成两个点集\(S',T'\),用\(ans ...

  3. bzoj 2229: [Zjoi2011]最小割【Gomory–Hu tree最小割树】

    这个算法详见http://www.cnblogs.com/lokiii/p/8191573.html 求出两两之间最小割之后暴力统计即可 #include<iostream> #inclu ...

  4. ●BOZJ 2229 [Zjoi2011]最小割

    题链: http://www.lydsy.com/JudgeOnline/problem.php?id=2229 题解: 首先先去看看这个博客:http://blog.csdn.net/jyxjyx2 ...

  5. 2229: [Zjoi2011]最小割(最小割树)

    Description 小白在图论课上学到了一个新的概念——最小割,下课后小白在笔记本上写下了如下这段话: “对于一个图,某个对图中结点的划分将图中所有结点分成两个部分,如果结点s,t不在同一个部分中 ...

  6. bzoj千题计划139:bzoj2229: [Zjoi2011]最小割

    http://www.lydsy.com/JudgeOnline/problem.php?id=2229 最小割树介绍:http://blog.csdn.net/jyxjyx27/article/de ...

  7. BZOJ2229: [Zjoi2011]最小割

    题解: 真是一道神题!!! 大家还是围观JZP的题解吧(网址找不到了...) 代码: #include<cstdio> #include<cstdlib> #include&l ...

  8. bzoj2229: [Zjoi2011]最小割(分治最小割+最小割树思想)

    2229: [Zjoi2011]最小割 题目:传送门 题解: 一道非常好的题目啊!!! 蒟蒻的想法:暴力枚举点对跑最小割记录...绝对爆炸啊.... 开始怀疑是不是题目骗人...难道根本不用网络流?? ...

  9. 【BZOJ2229】[ZJOI2011]最小割(网络流,最小割树)

    [BZOJ2229][ZJOI2011]最小割(网络流,最小割树) 题面 BZOJ 洛谷 题解 戳这里 那么实现过程就是任选两点跑最小割更新答案,然后把点集划分为和\(S\)联通以及与\(T\)联通. ...

随机推荐

  1. c# windows service 实现监控其他程序是否被关闭,关闭则报警

    namespace MonitorService { public partial class MonitorSv : ServiceBase { string AppName = "&qu ...

  2. python 关键字的操作

    声明:本文章默认使用的是python 3.6.1 1.要想当个牛逼的程序员,就要精通各种hello world的写法,当然,我不牛逼,只能用python去写^..^! print("Hell ...

  3. thinkphp中的常见静态常亮

    thinkphp __PUBLIC__的定义 __ROOT__等常量的定义 1 2 3 4 5 6 7 8 9 '__TMPL__'      =>  APP_TMPL_PATH,  // 项目 ...

  4. Java.nio-随机读写汉字

    笔者最近在用多线程来计算中文文本的标点符号数目,遇到了以下问题: 在Windows下,文本中汉字通常采用Unicode编码,这就导致需要随机(RandomAccessFile)读取文本时,产生乱码现象 ...

  5. J2ee入门:servlet-mapping的映射配置

    <servlet-mapping>元素在Servlet和URL样式之间定义一个映射.它包含了两个子元素<servlet- name>和<url-pattern> & ...

  6. 新概念英语(1-25)Mrs. Smith's Kitchen

    What colour is the electric cooker? Mrs. Smith's Kitchen is small. There is a refrigerator in the ki ...

  7. Spring Security 入门(1-3-3)Spring Security - logout 退出登录

    要实现退出登录的功能我们需要在 http 元素下定义 logout 元素,这样 Spring Security 将自动为我们添加用于处理退出登录的过滤器 LogoutFilter 到 FilterCh ...

  8. Spring Security 入门(1-6-2)Spring Security - 内置的filter顺序、自定义filter、http元素和对应的filterChain

    Spring Security 的底层是通过一系列的 Filter 来管理的,每个 Filter 都有其自身的功能,而且各个 Filter 在功能上还有关联关系,所以它们的顺序也是非常重要的. 1.S ...

  9. Python之线程

    操作系统线程理论 线程概念的引入背景 进程 之前我们已经了解了操作系统中进程的概念,程序并不能单独运行,只有将程序装载到内存中,系统为它分配资源才能运行,而这种执行的程序就称之为进程.程序和进程的区别 ...

  10. 表单提交中的input、button、submit的区别

    1.input[type=submit] 我们直接来看例子: 代码如下: <form> <input name="name"> <input type ...