BZOJ 3118 Orz the MST
权限题qwq
如果我们要使得某棵生成树为最小生成树,那么上面的边都不能被替代,具体的,对于一个非树边,它的权值要\(\ge\)它两端点在树上的路径上的所以边的权值,所以对于每个非树边就可以对一些树边列出不等关系,并且显然要把树边减小边权,非树边增加边权才会更优,所以记树边减少量为\(dt_a\),非树边增加量为\(dt_b\),就可以改写成不等式\(dt_a+dt_b\ge \max(w_a-w_b,0)\).那么问题变成若干变量,每个变量每增加\(1\)要付出一定代价,求最小代价满足所有不等式
这个不等式的形式有点不好下手,考虑这个问题的对偶问题,这等价于对于每个不等式\(i\)确定一个非负权值,每有一个权值获得\(max(w_a-w_b,0)\)的收益,同时满足权值\(\le\)对应树边和非树边修改一次的代价,求最大收益.那么可以网络流解决,\(S\)向所有树边连容量/费用为\((c_a,0)\)的边,非树边向\(T\)连\((c_b,0)\)的边,有不等式关系的一组边连\((+\infty,max(w_a-w_b,0))\)的边,然后最大费用流即可
#include<bits/stdc++.h>
#define LL long long
#define uLL unsigned long long
#define db double
using namespace std;
const int N=2000+10,M=N*100,inf=1<<30;
int rd()
{
int x=0,w=1;char ch=0;
while(ch<'0'||ch>'9'){if(ch=='-') w=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
return x*w;
}
int to[M],nt[M],c[M],w[M],hd[N],tot=1;
void add(int x,int y,int z,int zz)
{
++tot,to[tot]=y,nt[tot]=hd[x],c[tot]=z,w[tot]=zz,hd[x]=tot;
++tot,to[tot]=x,nt[tot]=hd[y],c[tot]=0,w[tot]=-zz,hd[y]=tot;
}
int ps,pt,di[N],fw[N],pr[N],cst;
bool v[N];
queue<int> q;
bool csfl()
{
for(int i=0;i<=pt;++i) di[i]=inf;
di[ps]=0,fw[ps]=inf,fw[pt]=0,v[ps]=1,q.push(ps);
while(!q.empty())
{
int x=q.front();
q.pop();
for(int i=hd[x];i;i=nt[i])
{
int y=to[i];
if(c[i]>0&&di[y]>di[x]+w[i])
{
di[y]=di[x]+w[i];
fw[y]=min(fw[x],c[i]),pr[y]=i;
if(!v[y]) v[y]=1,q.push(y);
}
}
v[x]=0;
}
if(!fw[pt]||di[pt]>=0) return 0; //边权取反做最大费用流,增广的权值不优就停止
cst+=di[pt]*fw[pt];
int x=pt;
while(x!=ps)
{
int i=pr[x];
c[i]-=fw[pt],c[i^1]+=fw[pt];
x=to[i^1];
}
return 1;
}
int n,m,ee[N][2],ew[N],cs[N],fa[N],pre[N],de[N];
bool fg[N];
void dd(int x)
{
for(int i=hd[x];i;i=nt[i])
{
int y=to[i];
if(y==fa[x]) continue;
if(fg[i>>1]) fa[y]=x,pre[y]=i>>1,de[y]=de[x]+1,dd(y);
}
}
void link(int i,int j)
{
int ss=max(0,ew[i]-ew[j]);
add(i,j,inf,-ss);
}
int main()
{
freopen("3118.in","r",stdin);
freopen("3118.out","w",stdout);
n=rd(),m=rd();
for(int i=1;i<=m;++i)
{
int x=rd(),y=rd();
add(ee[i][0]=x,ee[i][1]=y,0,0);
ew[i]=rd(),fg[i]=rd();
if(fg[i]) rd(),cs[i]=rd();
else cs[i]=rd(),rd();
}
de[1]=1,dd(1);
memset(hd,0,sizeof(int)*(n+3)),tot=1;
ps=0,pt=m+2;
for(int i=1;i<=m;++i)
{
if(fg[i]) add(ps,i,cs[i],0);
else
{
add(i,pt,cs[i],0);
int x=ee[i][0],y=ee[i][1];
if(de[x]<de[y]) swap(x,y);
while(de[x]>de[y]) link(pre[x],i),x=fa[x];
while(x!=y)
{
link(pre[x],i),x=fa[x];
link(pre[y],i),y=fa[y];
}
}
}
while(csfl());
printf("%d\n",-cst);
return 0;
}
BZOJ 3118 Orz the MST的更多相关文章
- bzoj 3118: Orz the MST(单纯形)
题目链接:http://www.lydsy.com:808/JudgeOnline/problem.php?id=3118 题意:给出一个图以及图中指定的n-1条边组成的生成树.每条边权值加1或者减去 ...
- BZOJ 2654 & 玄学二分+MST
题意: 给一张图,边带权且带颜色黑白,求出一棵至少包含k条白边的MST SOL: 正常人都想优先加黑边或者是白边,我也是这么想的...你看先用白边搞一棵k条边的MST...然后维护比较黑边跟白边像堆一 ...
- BZOJ3118 : Orz the MST
对于树边显然只需要减少权值,对于非树边显然只需要增加权值 设i不为树边,j为树边 X[i]:i增加量 X[j]:j减少量 C[i]:修改1单位i的代价 对于每条非树边i(u,v),在树上u到v路径上的 ...
- BZOJ 2654: tree( 二分 + MST )
我们给白色的边增加权值 , 则选到的白色边就会变多 , 因此可以二分一下. 不过这道题有点小坑... ------------------------------------------------- ...
- BZOJ3118 Orz the MST 【单纯形 + 生成树】
题目链接 BZOJ3118 题解 少有的单纯形好题啊 我们先抽离出生成树 生成树中的边只可能减,其它边只可能加 对于不在生成树的边,其权值一定要比生成树中其端点之间的路径上所有的边都大 然后就是一个最 ...
- bzoj3118: Orz the MST(线性规划+单纯形法)
传送门 不难发现,对于每一条树边肯定要减小它的权值,对于每一条非树边要增加它的权值 对于每一条非树边\(j\),他肯定与某些树边构成了一个环,那么它的边权必须大于等于这个环上的所有边 设其中一条边为\ ...
- bzoj AC倒序
Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...
- BZOJ4652 [Noi2016]循环之美 【数论 + 莫比乌斯反演 + 杜教筛】
题目链接 BZOJ 题解 orz 此题太优美了 我们令\(\frac{x}{y}\)为最简分数,则\(x \perp y\)即,\(gcd(x,y) = 1\) 先不管\(k\)进制,我们知道\(10 ...
- luogu4169 [Violet]天使玩偶/SJY摆棋子 / bzoj2648 SJY摆棋子 k-d tree
k-d tree + 重构的思想,就能卡过luogu和bzoj啦orz #include <algorithm> #include <iostream> #include &l ...
随机推荐
- R_Studio(癌症)数据连续属性离散化处理
对“癌症.csv”中的肾细胞癌组织内微血管数进行连续属性的等宽离散化处理(分为3类),并用宽值找替原来的值 癌症.csv setwd('D:\\data') list.files() dat=read ...
- Unity3D_(游戏)跳一跳超简单制作过程
跳一跳 工程文件界面 游戏界面 脚本 using DG.Tweening; using System.Collections; using System.Collections.Generic; us ...
- JS框架_(JQuery.js)网页文字评论弹幕
百度云盘 传送门 密码:3azl jQuery网页右下角文字评论弹幕效果 <!DOCTYPE html> <html> <head> <title>jQ ...
- windows 安装使用 Memcached
Windows无官方版本:下载地址http://static.runoob.com/download/memcached-win64-1.4.4-14.zip 安装: 1.解压下载的压缩包2.命令行模 ...
- 套接字之 getsockname && getpeername
getsockname-获取本地地址:比如,在绑定的时候设置端口号为0由系统自动选择端口绑定,或者使用了INADDR_ANY通配所有地址的情况下,后面需要用到具体的地址和端口,就可以用getsockn ...
- JAVA-ThreadPoolExecutor 线程池
一.创建线程池 /** * @param corePoolSize 核心线程池大小 * 当提交一个任务到线程池时,如果当前 poolSize < corePoolSize 时,线程池会创建一个线 ...
- [HTML辅助方法-Html.Raw()的简单应用]
Html.Raw(); 当我们使用 文本编辑器,存入到数据库中的数据会带 html 标签,如果我们需要在前台显示存入时的相同样式,不输出为带有html标签的字符串 ,不通过富文本显示的话,可以通过ht ...
- xcode dyld: Library not loaded: @rpath/libswiftCore.dylib问题解决
app安装好了之后就报这个错误,这个时候可以将xcode工程clear一下,删除已经安装好的app,再重新安装即可
- WPF VLC客户端和SDK的简单应用
VLC_SDK编程指南 VLC 是一款自由.开源的跨平台多媒体播放器及框架,可播放大多数多媒体文件,以及 DVD.音频 CD.VCD 及各类流媒体协议.它可以支持目前市面上大多数的视频解码,除了Rea ...
- STM32 PWM注意事项
频率: f = 时钟周期 / Prescaler / Period: 占空比 = ( Pulse / Period )* 100 %: 需要注意的是,Pulse 不能大于 Period