权限题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的更多相关文章

  1. bzoj 3118: Orz the MST(单纯形)

    题目链接:http://www.lydsy.com:808/JudgeOnline/problem.php?id=3118 题意:给出一个图以及图中指定的n-1条边组成的生成树.每条边权值加1或者减去 ...

  2. BZOJ 2654 & 玄学二分+MST

    题意: 给一张图,边带权且带颜色黑白,求出一棵至少包含k条白边的MST SOL: 正常人都想优先加黑边或者是白边,我也是这么想的...你看先用白边搞一棵k条边的MST...然后维护比较黑边跟白边像堆一 ...

  3. BZOJ3118 : Orz the MST

    对于树边显然只需要减少权值,对于非树边显然只需要增加权值 设i不为树边,j为树边 X[i]:i增加量 X[j]:j减少量 C[i]:修改1单位i的代价 对于每条非树边i(u,v),在树上u到v路径上的 ...

  4. BZOJ 2654: tree( 二分 + MST )

    我们给白色的边增加权值 , 则选到的白色边就会变多 , 因此可以二分一下. 不过这道题有点小坑... ------------------------------------------------- ...

  5. BZOJ3118 Orz the MST 【单纯形 + 生成树】

    题目链接 BZOJ3118 题解 少有的单纯形好题啊 我们先抽离出生成树 生成树中的边只可能减,其它边只可能加 对于不在生成树的边,其权值一定要比生成树中其端点之间的路径上所有的边都大 然后就是一个最 ...

  6. bzoj3118: Orz the MST(线性规划+单纯形法)

    传送门 不难发现,对于每一条树边肯定要减小它的权值,对于每一条非树边要增加它的权值 对于每一条非树边\(j\),他肯定与某些树边构成了一个环,那么它的边权必须大于等于这个环上的所有边 设其中一条边为\ ...

  7. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

  8. BZOJ4652 [Noi2016]循环之美 【数论 + 莫比乌斯反演 + 杜教筛】

    题目链接 BZOJ 题解 orz 此题太优美了 我们令\(\frac{x}{y}\)为最简分数,则\(x \perp y\)即,\(gcd(x,y) = 1\) 先不管\(k\)进制,我们知道\(10 ...

  9. luogu4169 [Violet]天使玩偶/SJY摆棋子 / bzoj2648 SJY摆棋子 k-d tree

    k-d tree + 重构的思想,就能卡过luogu和bzoj啦orz #include <algorithm> #include <iostream> #include &l ...

随机推荐

  1. 深入理解Vuex 模块化(module)

    todo https://www.jb51.net/article/124618.htm

  2. LeetCode124----二叉树中最大路径和

    给定一个非空二叉树,返回其最大路径和. 本题中,路径被定义为一条从树中任意节点出发,达到任意节点的序列.该路径至少包含一个节点,且不需要经过根节点. 示例 1: 输入: [1,2,3] 1 / \ 2 ...

  3. LeetCode 18. 四数之和(4Sum)

    题目描述 给定一个包含 n 个整数的数组 nums 和一个目标值 target,判断 nums 中是否存在四个元素 a,b,c 和 d ,使得 a + b + c + d 的值与 target 相等? ...

  4. LeetCode 19. 删除链表的倒数第N个节点(Remove Nth Node From End Of List)

    题目描述 给定一个链表,删除链表的倒数第 n 个节点,并且返回链表的头结点. 示例: 给定一个链表: 1->2->3->4->5, 和 n = 2. 当删除了倒数第二个节点后, ...

  5. SpringBoot启动加载yml配置文件出现编码格式错误

    Caused by: org.yaml.snakeyaml.error.YAMLException: java.nio.charset.MalformedInputException: Input l ...

  6. Hidden的应用

    在写jsp中如果一个 请求的参数(例如:paramTypeCode)不能在另一个请求中使用,我们为了能让他在请求中使用可以利用隐藏域来表示,下面介绍他的用法: 1    <input type= ...

  7. 【9】letter-spacing / box-shadow

    1.letter-spacing :增加或减少字符间的空白(字符间距),如:h1 {letter-spacing:2px} 2.box-shadow : box-shadow: 10px 10px 5 ...

  8. bulk_create(lst) 批量创建数据

    # 批量创建数据 # Create your views here. from django.db import models from django.shortcuts import HttpRes ...

  9. Jmeter(七)参数化

    初识Jmeter的时候, 除了感觉安装和配置都很轻量以外, 还有一个最大的感触就是, 翻译真硬啊, 真的够够的! 和他磨合了挺长一段时间之后, 终于开悟了, 这些硬硬的翻译, 其实还是基本靠谱的, 看 ...

  10. 分期花呗 账户交易通知:尾号6932客户,您的申请已通过,账户余额38139元,无手续费,点t.cn/Aijsx9vq取款,回T退订。

    10692285499 分期花呗 账户变动通知:尾号6932客户,您的申请已通过,账户余额5000元,请及时点击t.cn/AiOMsNAm取款,回T退订. 106935276259002分期花呗 账户 ...