bzoj 3118: Orz the MST(单纯形)
题目链接:http://www.lydsy.com:808/JudgeOnline/problem.php?id=3118
题意:给出一个图以及图中指定的n-1条边组成的生成树。每条边权值加1或者减去1都有相应的代价。求一个最小代价使得给出的边是最小生成树。
思路:对于每条非树边,必与某些树边形成环。设该非树边的权值为w2,某树边的权值为 w1。最后非树边增加x2,树边减少x1,那么w1-x1<=w2+x2。这样我们可以得到一些式子。代价也知道,这样就转化成线性规划问题。题目求的是最小值,我们可以将目标方程的系数取反求最大值。
单纯形的步骤:
(1)求出一个初始解;
(2)迭代。
(1)这个题的系数矩阵A是全么模:1、元素都是0,-1,1;2、任意子方阵的行列式为0,-1,1。
(2)据说A是全么模时解是整数解,因此此题可直接用单纯形。
const int COL=1005; const int ROW=30005; int n,m,B[ROW],N[COL]; double A[ROW][COL],b[ROW],c[COL],v; double ans[COL]; int sgn(double x) { if(x>1e-8) return 1; if(x<-1e-8) return -1; return 0; } //B中第l个替换N中第e个 void pivot(int l,int e) { int i,j; double temp=A[l][e]; b[l]/=temp; A[l][e]=1/temp; for(i=1;i<=n;i++) if(i!=e) A[l][i]/=temp; for(i=1;i<=m;i++) if(i!=l) { b[i]-=A[i][e]*b[l]; for(j=1;j<=n;j++) if(j!=e) A[i][j]-=A[i][e]*A[l][j]; A[i][e]=-A[i][e]/temp; } v+=b[l]*c[e]; for(i=1;i<=n;i++) if(i!=e) c[i]-=c[e]*A[l][i]; c[e]*=-A[l][e]; swap(B[l],N[e]); } void simplex() { int i,j,k,x; int l,s; double temp,temp1,temp2,temp3; while(1) { temp2=-dinf; s=-1; for(i=1;i<=n;i++) if(sgn(c[i])>0) { temp=dinf; for(k=1;k<=m;k++) if(sgn(A[k][i])>0) { temp3=b[k]/A[k][i]; if(temp3<temp) temp=temp3,x=k; } if(temp2<temp*c[i]) { s=i,l=x,temp2=temp*c[i]; } } if(s==-1) break; pivot(l,s); } for(i=1;i<=n;i++) { for(j=1;j<=m;j++) if(B[j]==i) break; if(j<=m) ans[i]=b[j]; else ans[i]=0; } } void print() { int i,j; printf("v: %.3lf\n",v); printf("c:\n"); for(i=1;i<=n;i++) printf("%.3lf ",c[i]); puts(""); printf("A:\n"); for(i=1;i<=m;i++) { for(j=1;j<=n;j++) printf("%.3lf ",A[i][j]); puts(""); } printf("b:\n"); for(i=1;i<=m;i++) printf("%.3lf ",b[i]); puts(""); printf("B:\n"); for(i=1;i<=m;i++) printf("%d ",B[i]); puts(""); printf("N:\n"); for(i=1;i<=n;i++) printf("%d ",N[i]); puts(""); } int init() { int i,j; int k=1; for(i=1;i<=m;i++) if(b[i]<b[k]) k=i; if(sgn(b[k])>=0) { for(i=1;i<=n;i++) N[i]=i; for(i=1;i<=m;i++) B[i]=n+i; v=0; simplex(); return 1; } static double tmpC[COL]; for(i=1;i<=n;i++) tmpC[i]=c[i]; tmpC[n+1]=0; n++; for(i=1;i<=m;i++) A[i][n]=-1; for(i=1;i<=n;i++) N[i]=i; for(i=1;i<=m;i++) B[i]=n+i; v=0; for(i=1;i<=n;i++) { if(i<n) c[i]=0; else c[i]=-1; } pivot(k,n); simplex(); if(sgn(ans[n])!=0) return 0; static int belongB[COL]; clr(belongB,0); for(i=1;i<=m;i++) { if(B[i]>n) continue; belongB[B[i]]=i; } map<int,int> mp; for(i=1;i<=n;i++) mp[N[i]]=i; clr(c,0); v=0; for(i=1;i<=n;i++) { if(!belongB[i]) { c[mp[i]]+=tmpC[i]; } else { v+=tmpC[i]*b[belongB[i]]; int j; for(j=1;j<=n;j++) { c[j]+=tmpC[i]*(-A[belongB[i]][j]); } } } c[mp[n]]=0; for(i=1;i<=m;i++) A[i][mp[n]]=0; simplex(); n--; return 1; } struct node { int u,v,id,w,next; }; node edges[ROW]; int e; int head[COL]; void add(int u,int v,int w,int id) { e++; edges[e].u=u; edges[e].v=v; edges[e].w=w; edges[e].id=id; edges[e].next=head[u]; head[u]=e; } int h[COL],up[COL],down[COL]; int eNum; int inq[COL],KK; int pre[COL]; void build(int s,int t,int p) { int i; for(i=pre[t];i!=-1;) { int w1=edges[i].w; int w2=edges[p*2].w; m++; A[m][edges[i].id]=-1; A[m][p]=-1; b[m]=w2-w1; int u=edges[i].u; i=pre[u]; } } void bfs(int s,int t,int p) { queue<int> Q; Q.push(s); KK++; inq[s]=KK; pre[s]=-1; while(!Q.empty()) { int u=Q.front(); Q.pop(); int i; for(i=head[u];i!=-1;i=edges[i].next) { int v=edges[i].v; int id=edges[i].id; if(!h[id]||KK==inq[v]) continue; pre[v]=i; if(v==t) { build(s,t,p); return; } Q.push(v); inq[v]=KK; } } } int main() { n=myInt(); eNum=myInt(); clr(head,-1); int i; for(i=1;i<=eNum;i++) { int u,v,w; scanf("%d%d%d%d%d%d",&u,&v,&w,&h[i],&up[i],&down[i]); add(u,v,w,i); add(v,u,w,i); } for(i=1;i<=eNum;i++) if(!h[i]) { int t=i*2; bfs(edges[t].u,edges[t].v,i); } for(i=1;i<=eNum;i++) c[i]=h[i]?-down[i]:-up[i]; n=eNum; if(!init()) { puts("no solution"); } double res=-v; if(res<1e-10) res=0; printf("%.0lf\n",-v); }
bzoj 3118: Orz the MST(单纯形)的更多相关文章
- BZOJ 3118 Orz the MST
权限题qwq 如果我们要使得某棵生成树为最小生成树,那么上面的边都不能被替代,具体的,对于一个非树边,它的权值要\(\ge\)它两端点在树上的路径上的所以边的权值,所以对于每个非树边就可以对一些树边列 ...
- BZOJ3118 Orz the MST 【单纯形 + 生成树】
题目链接 BZOJ3118 题解 少有的单纯形好题啊 我们先抽离出生成树 生成树中的边只可能减,其它边只可能加 对于不在生成树的边,其权值一定要比生成树中其端点之间的路径上所有的边都大 然后就是一个最 ...
- BZOJ 2654 & 玄学二分+MST
题意: 给一张图,边带权且带颜色黑白,求出一棵至少包含k条白边的MST SOL: 正常人都想优先加黑边或者是白边,我也是这么想的...你看先用白边搞一棵k条边的MST...然后维护比较黑边跟白边像堆一 ...
- BZOJ 3112 Zjoi2013 防守战线 单纯形
题目大意: 单纯形*2.. . #include <cmath> #include <cstdio> #include <cstring> #include < ...
- BZOJ 3265 志愿者招募加强版(单纯形)
3265: 志愿者招募加强版 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 848 Solved: 436[Submit][Status][Disc ...
- 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 )
我们给白色的边增加权值 , 则选到的白色边就会变多 , 因此可以二分一下. 不过这道题有点小坑... ------------------------------------------------- ...
- bzoj [Noi2008] 1061 志愿者招募 单纯形
[Noi2008]志愿者招募 Time Limit: 20 Sec Memory Limit: 162 MBSubmit: 5437 Solved: 3267[Submit][Status][Di ...
- bzoj3118: Orz the MST(线性规划+单纯形法)
传送门 不难发现,对于每一条树边肯定要减小它的权值,对于每一条非树边要增加它的权值 对于每一条非树边\(j\),他肯定与某些树边构成了一个环,那么它的边权必须大于等于这个环上的所有边 设其中一条边为\ ...
随机推荐
- [OrangePi] Building the system
You can try to build Debian/Ubuntu for OrangePI yourself. Clone my GitHub repository. You will need ...
- wc移植sae笔记
1.wc移植到sae---上传图片 ①先看profile.ptl.html中的ajax代码修改functions.js中G_BASE_URL的值.在这里我先写死成->'http://2.idan ...
- Android压力测试快速入门教程(图解)——Monkey工具
文章目录: 一.Monkey简介 二.Monkey的基本用法 三.Monkey测试示例图解 四.Monkey命令参数介绍 五.Monkey log分析 一.Monkey简介 Monkey:Androi ...
- nginx负载均衡的实现
将一台nginx主机当作前端负载均衡服务器,后面通过交换机链接多台web服务器,提供html和php的web界面服务.通过配置前端负载均衡服务器,可以实现将html界面和php界面的分开访问,即htm ...
- asp显示记录条数
Sql = "select * from xin126 where ID=" & id Rs.Open Sql,Conn,1,1 %> 共有<strong st ...
- 记录一下:chrome上,把网页保存为文件的插件
插件地址: https://chrome.google.com/webstore/detail/full-page-screen-capture/fdpohaocaechififmbbbbbknoal ...
- C#:关联程序和文件
一.关联代码 /// <summary> /// 关联程序和类型 /// </summary> private void RegFileExt() { try { string ...
- ACM题目————A simple problem
Description Zty很痴迷数学问题..一天,yifenfei出了个数学题想难倒他,让他回答1 / n.但Zty却回答不了^_^. 请大家编程帮助他. Input 第一行整数T,表示测试组 ...
- linux 时钟源初步分析linux kernel 时钟框架详细介绍
初步概念: 看datasheet的关于时钟与定时器的部分, FCLK供给cpu, HCLK供给AHB总线设备(存储器控制器,中断控制器.LCD控制器.DMA.USB主机控制器等), PCLK供给APB ...
- HDU 自动刷题机 Auto AC (轻轻松松进入HDU首页)
前言: 在写这篇文章之前,首先感谢给我思路以及帮助过我的学长们 以下4篇博客都是学长原创,其中有很多有用的,值得学习的东西,希望能够帮到大家! 1.手把手教你用C++ 写ACM自动刷题神器(冲入HDU ...