【ZJOI2011】【BZOJ2229】最小割
来源:互联网 发布:java网上商城视频教程 编辑:程序博客网 时间:2024/04/29 13:57
Description
小白在图论课上学到了一个新的概念——最小割,下课后小白在笔记本上写下了如下这段话: “对于一个图,某个对图中结点的划分将图中所有结点分成两个部分,如果结点s,t不在同一个部分中,则称这个划分是关于s,t的割。 对于带权图来说,将所有顶点处在不同部分的边的权值相加所得到的值定义为这个割的容量,而s,t的最小割指的是在关于s,t的割中容量最小的割” 现给定一张无向图,小白有若干个形如“图中有多少对点它们的最小割的容量不超过x呢”的疑问,小蓝虽然很想回答这些问题,但小蓝最近忙着挖木块,于是作为仍然是小蓝的好友,你又有任务了。
Input
输入文件第一行有且只有一个正整数T,表示测试数据的组数。 对于每组测试数据, 第一行包含两个整数n,m,表示图的点数和边数。 下面m行,每行3个正整数u,v,c(1<=u,v<=n,0<=c<=106),表示有一条权为c的无向边(u,v) 接下来一行,包含一个整数q,表示询问的个数 下面q行,每行一个整数x,其含义同题目描述。
Output
对于每组测试数据,输出应包括q行,第i行表示第i个问题的答案。对于点对(p,q)和(q,p),只统计一次(见样例)。
两组测试数据之间用空行隔开。
Sample Input
1 5 0 1 0
Sample Output
10
【数据范围】
对于100%的数据 T<=10,n<=150,m<=3000,q<=30,x在32位有符号整数类型范围内。
图中两个点之间可能有多条边
HINT
Source
Day1
学了一下最小割树,做了个裸题.
#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<algorithm>#include<queue>#define MAXN 200#define MAXINT 1ll<<62#define GET (ch>='0'&&ch<='9')#define LL long longusing namespace std;int Tst,n,m,S,T,Q,top;int vis[MAXN];int num[MAXN],tmp[MAXN];LL mncut[MAXN][MAXN];int dis[MAXN],cnt[MAXN];LL C[MAXN][MAXN];struct edge { int st,to;LL c; edge *next,*rev; }e[MAXN*MAXN],*prev[MAXN];inline void insert(int u,int v,LL c) { e[++top].to=v;e[top].st=u;e[top].next=prev[u];prev[u]=&e[top];e[top].c=c; }inline void add(int u,int v,LL c) { insert(u,v,c);insert(v,u,c);prev[u]->rev=prev[v];prev[v]->rev=prev[u]; }LL ISAP(){ edge *E[MAXN],*rep[MAXN];int now=S;LL ret=0; memset(cnt,0,sizeof(cnt));memset(dis,0,sizeof(dis)); for (int i=1;i<=n;i++) E[i]=prev[i];cnt[0]=n; while (dis[S]<=n) { edge *i;bool t=0; for (i=E[now];i;i=i->next) if (i->c>0&&dis[i->to]+1==dis[now]) { t=1;E[now]=i;break; } if (t) { rep[now=i->to]=i; if (now==T) { LL minn=MAXINT; for (int i=T;i!=S;i=rep[i]->st) minn=min(minn,rep[i]->c); for (int i=T;i!=S;i=rep[i]->st) rep[i]->c-=minn,rep[i]->rev->c+=minn; ret+=minn;now=S; } } else { if (!(--cnt[dis[now]])) break; int mind=n+1;E[now]=prev[now]; for (edge *i=prev[now];i;i=i->next) if (i->c>0) mind=min(mind,dis[i->to]); dis[now]=mind+1;++cnt[dis[now]]; if (now!=S) now=rep[now]->st; } } return ret;}void dfs(int x){ vis[x]=1; for (edge *i=prev[x];i;i=i->next) if (i->c>0&&!vis[i->to]) dfs(i->to);}void solve(int l,int r){ if (l==r) return; for (int i=1;i<=top;i+=2) e[i].c=e[i].rev->c=((e[i].c+e[i].rev->c)>>1); S=num[l];T=num[r];LL t=ISAP(); memset(vis,0,sizeof(vis));dfs(S); for (int i=1;i<=n;i++) if (vis[i]) for (int j=1;j<=n;j++) if (!vis[j]) mncut[i][j]=min(mncut[i][j],t),mncut[j][i]=min(mncut[j][i],t); int tp=0,len,L=l-1; for (int i=l;i<=r;i++) if (vis[num[i]]) tmp[++tp]=num[i];len=tp; for (int i=l;i<=r;i++) if (!vis[num[i]]) tmp[++tp]=num[i]; for (int i=1;i<=tp;i++) num[++L]=tmp[i]; solve(l,l+len-1);solve(l+len,r);}inline void in(int &x){ char ch=getchar();x=0; while (!GET) ch=getchar(); while (GET) x=x*10+ch-'0',ch=getchar();}int main(){ for (in(Tst);Tst;--Tst) { in(n);in(m);top=0;int u,v,w,k; memset(mncut,0x3f,sizeof(mncut));memset(prev,0,sizeof(prev));memset(C,0,sizeof(C)); for (int i=1;i<=n;i++) num[i]=i; for (int i=1;i<=m;i++) in(u),in(v),in(w),C[u][v]+=w,C[v][u]+=w; for (int i=1;i<=n;i++) for (int j=i+1;j<=n;j++) if (C[i][j]) add(i,j,C[i][j]); solve(1,n); for (in(Q);Q;Q--) { in(k);int ans=0; for (int i=1;i<=n;i++) for (int j=i+1;j<=n;j++) ans+=(mncut[i][j]<=k); printf("%d\n",ans); } puts(""); }}
1 0
- 【ZJOI2011】【BZOJ2229】最小割
- bzoj2229: [Zjoi2011]最小割
- Bzoj2229:[Zjoi2011]最小割:分治最小割
- BZOJ2229: [Zjoi2011]最小割 分治最小割
- 【BZOJ2229】[Zjoi2011]最小割【Gomory-Hu树】
- BZOJ2229—— [Zjoi2011]最小割
- 【BZOJ2229】【ZJOI2011】最小割 {没有错,这道题的算法跟题帽是一样的!!!}
- 【bzoj2229】【ZJOI2229】【最小割】【最小割】
- 2229: [Zjoi2011]最小割
- P3329 [ZJOI2011]最小割
- 2229: [Zjoi2011]最小割 最小割
- ZJOI2011最小割 最小割树
- bzoj 2229: [Zjoi2011]最小割
- BZOJ 2229:[Zjoi2011]最小割
- BZOJ 2229: [Zjoi2011]最小割
- BZOJ 2229 [Zjoi2011]最小割
- BZOJ 2229 [Zjoi2011] 最小割
- 【BZOJ2229】最小割【BZOJ4519】不同的最小割
- Android之XML序列化和解析
- 排序算法--冒泡排序(一)
- oracle 使用sql获取数据库表,表的字段
- Smile - Statistical Machine Intelligence and Learning Engine
- ubuntu14.04系统加固(4)
- 【ZJOI2011】【BZOJ2229】最小割
- LayoutInflater.inflate的用法总结
- [4] 数组和指针的前世今生 - 指针篇
- 华为机试题ex002
- java之重载和重写
- 批量漏洞检测框架
- Java中的Map汇总——思维导图
- Latex 算法Algorithm
- Spark运行任务