博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
BZOJ 2561 最小生成树
阅读量:4970 次
发布时间:2019-06-12

本文共 2159 字,大约阅读时间需要 7 分钟。

第一眼瞎那啥贪心,然后觉得不太对劲,就滚去看题解,发现是网络流OTZ

模拟Kruskal的过程发现,若<u,v>要在最小生成树中出现,权值则小于<u,v>的边不能让u,v联通,转换成最小割模型,最大生成树同理。

跑两遍最大流。

注意边要建双向的啊,被这点坑死了。然后板子不要瞎那啥乱打。

顺便发现果然ISAP跑得很快。

//Twenty#include
#include
#include
#include
#include
#include
#include
#include
#define INF 0x3f3f3f3ftypedef long long LL;using namespace std;const int maxn=1e6+299;int ecnt=1,ans,n,m,u,v,w,fir[maxn],pre[maxn];struct Edge { int u,v,w,no,vv; Edge(){} Edge(int u,int v,int w):u(u),v(v),w(w){}}ee[maxn];struct edge { int from,to,flow,cap,nxt; edge(){} edge(int from,int to,int cap,int flow,int nxt):from(from),to(to),cap(cap),flow(flow),nxt(nxt){}}e[maxn];int d[maxn],c[maxn],cur[maxn];void add(int from,int to,int cap) { e[++ecnt]=edge(from,to,cap,0,fir[from]); fir[from]=ecnt; e[++ecnt]=edge(to,from,0,0,fir[to]); fir[to]=ecnt;}queue
que;void bfs(int s,int t) { for(int i=1;i<=n;i++) d[i]=n; d[t]=0; que.push(t); while(!que.empty()) { int x=que.front(); que.pop(); for(int i=fir[x];i;i=e[i].nxt) { edge tp=e[i]; if(d[e[i].to]==n&&e[i].flow==e[i].cap) { int to=e[i].to; d[to]=d[x]+1; que.push(to); } } }}int cal(int s,int t) { int fl=1e9; for(int i=t;i!=s;i=e[pre[i]].from) fl=min(fl,e[pre[i]].cap-e[pre[i]].flow); for(int i=t;i!=s;i=e[pre[i]].from) { e[pre[i]].flow+=fl; e[pre[i]^1].flow-=fl; } return fl;} int maxflow(int s,int t) { int res=0; bfs(s,t); for(int i=1;i<=n;i++) cur[i]=fir[i],c[d[i]]++; for(int x=s;d[x]
e[i].flow) { pre[x=e[i].to]=i; ok=1; break; } if(!ok) { cur[x]=fir[x]; int M=n; for(int i=fir[x];i;i=e[i].nxt) { if(e[i].cap>e[i].flow&&(M>d[e[i].to]+1)) M=d[e[i].to]+1; } if(!(--c[d[x]])) break; c[d[x]=M]++; if(x!=s) x=e[pre[x]].from; } } return res;}int main(){ //freopen(".in","r",stdin); //freopen(".out","w",stdout); scanf("%d%d",&n,&m); for(int i=1;i<=m;i++) { scanf("%d%d%d",&u,&v,&w); ee[i]=Edge(u,v,w); } scanf("%d%d%d",&u,&v,&w); for(int i=1;i<=m;i++) { if(ee[i].w
w) { add(ee[i].u,ee[i].v,1); add(ee[i].v,ee[i].u,1); } } ans+=maxflow(u,v); printf("%d\n",ans); return 0;}
View Code

 

转载于:https://www.cnblogs.com/Achenchen/p/7604812.html

你可能感兴趣的文章
Learning Cocos2d-x for XNA(2)——深入剖析Hello World
查看>>
软件建模——第9章 毕业论文管理系统—面向对象方法
查看>>
Http协议
查看>>
手机端web开发必备代码
查看>>
[SDOI2008]洞穴勘测
查看>>
NOI2014 购票
查看>>
Difference between Linearizability and Serializability
查看>>
电影《绿皮书》
查看>>
IDEA使用操作文档
查看>>
如何对网课、游戏直播等进行录屏
查看>>
UIView
查看>>
有关去掉谷歌及火狐浏览器文本框 数字类型 上下箭头的方法
查看>>
MySQL数据迁移到SQL Server
查看>>
复杂链表的复制(python)
查看>>
添加日期选择控件
查看>>
jquery.cookie.js操作cookie
查看>>
javascript遍历数组
查看>>
bzoj4765: 普通计算姬 (分块 && BIT)
查看>>
thinkphp5-----模板中函数的使用
查看>>
POJ-3211 Washing Clothes[01背包问题]
查看>>