题目:http://www.lydsy.com/JudgeOnline/problem.php?id=2001

cdq分治+重建图。

可以保留当前一定会被选的非修改边然后把点缩起来。这样的话每次点数至多只有r-l+1个,边数只有2*(r-l+1)-1个(生成树+修改边)。

时间复杂度O(nlog^2(n))

#include<cstring>
#include<iostream>
#include<cstdio>
#include<algorithm>
#define rep(i,l,r) for (int i=l;i<=r;i++)
#define down(i,l,r) for (int i=l;i>=r;i--)
#define clr(x,y) memset(x,y,sizeof(x))
#define maxn 50050
#define ll long long
#define inf 1000000000
using namespace std;
struct data{int x,y;ll z;int k,id;
}a[][maxn];
struct node{int k;ll w;int id;
}b[maxn];
int fa[maxn],pos[maxn],p[maxn],n,m,Q;
ll dis[maxn];
int read(){
int x=,f=; char ch=getchar();
while (!isdigit(ch)) {if (ch=='-') f=-; ch=getchar();}
while (isdigit(ch)) {x=x*+ch-''; ch=getchar();}
return x*f;
}
int find(int x){
if (fa[x]==x) return x;
return fa[x]=find(fa[x]);
}
bool cmp(data a,data b){
return a.z<b.z;
}
void cdq(int dep,int l,int r,int n,int m,ll ans){
rep(i,,m){
data &c=a[dep][i],&c2=a[dep-][i];
a[dep][i]=(data){c2.x,c2.y,dis[c2.id],,c2.id};
pos[c.id]=i;
}
//当l=r,计算答案,直接赋值当前修改边。
if (l==r){
data &c=a[dep][pos[b[l].k]]; c.z=dis[c.id]=b[l].w;
rep(i,,n) fa[i]=i;
sort(a[dep]+,a[dep]++m,cmp);
rep(i,,m){
int x=find(a[dep][i].x),y=find(a[dep][i].y);
if (x!=y) fa[x]=y,ans+=a[dep][i].z;
}
printf("%lld\n",ans);
return;
}
//删边,如果k=0,那么这是一条无用的边。
rep(i,l,r) a[dep][pos[b[i].k]].k=;
rep(i,,n) fa[i]=i;
sort(a[dep]+,a[dep]++m,cmp);
rep(i,,m) if (a[dep][i].k==){
int x=find(a[dep][i].x),y=find(a[dep][i].y);
if (x!=y) fa[x]=y,a[dep][i].k=;
}
//删点,找出必连边(k=3)
rep(i,,m) if (a[dep][i].k==) a[dep][i].z=-inf;
rep(i,,n) fa[i]=i;
sort(a[dep]+,a[dep]++m,cmp);
rep(i,,m) if (a[dep][i].k!=){
int x=find(a[dep][i].x),y=find(a[dep][i].y);
if (x!=y) {
fa[x]=y;
if (a[dep][i].k!=) a[dep][i].k=;
}
}
//重建,把必连边(k=3)取为答案,将两端点缩为一个点。
rep(i,,n) fa[i]=i;
rep(i,,m) if (a[dep][i].k==){
data &c=a[dep][i];
int x=find(c.x),y=find(c.y);
if (x!=y) fa[x]=y, ans+=c.z;
}
int nn=,nm=;
rep(i,,n) p[i]=;
rep(i,,n) if (!p[find(i)]) p[find(i)]=++nn;
rep(i,,m) if (a[dep][i].k!=){
int x=find(a[dep][i].x),y=find(a[dep][i].y);
data &c=a[dep][i];
if (x!=y) a[dep][++nm]=(data){p[x],p[y],c.z,c.k,c.id};
}
int mid=(l+r)/;
cdq(dep+,l,mid,nn,nm,ans);
cdq(dep+,mid+,r,nn,nm,ans);
}
int main(){
n=read(); m=read(); Q=read();
int x,y; ll z;
rep(i,,m){
x=read(); y=read(); z=read();
a[][i]=(data){x,y,z,,i};
}
rep(i,,Q){
x=read(); y=read();
b[i]=(node){x,y,i};
}
rep(i,,m) dis[a[][i].id]=a[][i].z;
cdq(,,Q,n,m,);
return ;
}

BZOJ2001: [Hnoi2010]City 城市建设的更多相关文章

  1. BZOJ2001 [Hnoi2010]City 城市建设 CDQ分治

    2001: [Hnoi2010]City 城市建设 Time Limit: 20 Sec  Memory Limit: 162 MB Description PS国是一个拥有诸多城市的大国,国王Lou ...

  2. BZOJ2001 [Hnoi2010]City 城市建设 【CDQ分治 + kruskal】

    题目链接 BZOJ2001 题解 CDQ分治神题... 难想难写.. 比较朴素的思想是对于每个询问都求一遍\(BST\),这样做显然会爆 考虑一下时间都浪费在了什么地方 我们每次求\(BST\)实际上 ...

  3. 【BZOJ2001】 [Hnoi2010]City 城市建设

    BZOJ2001 [Hnoi2010]City 城市建设 Solution 我们考虑一下这个东西怎么求解? 思考无果...... 咦? 好像可以离线cdq,每一次判断一下如果这条边如果不选就直接删除, ...

  4. BZOJ 2001: [Hnoi2010]City 城市建设

    2001: [Hnoi2010]City 城市建设 Time Limit: 20 Sec  Memory Limit: 162 MBSubmit: 1132  Solved: 555[Submit][ ...

  5. 【bzoj2001】 Hnoi2010—City 城市建设

    http://www.lydsy.com/JudgeOnline/problem.php?id=2001 (题目链接) 题意 给出一张无向图,$m$组操作,每次修改一条边的权值,对于每次操作输出修改之 ...

  6. 2001: [Hnoi2010]City 城市建设 - BZOJ

    DescriptionPS国是一个拥有诸多城市的大国,国王Louis为城市的交通建设可谓绞尽脑汁.Louis可以在某些城市之间修建道路,在不同的城市之间修建道路需要不同的花费.Louis希望建造最少的 ...

  7. 【刷题】BZOJ 2001 [Hnoi2010]City 城市建设

    Description PS国是一个拥有诸多城市的大国,国王Louis为城市的交通建设可谓绞尽脑汁.Louis可以在某些城市之间修建道路,在不同的城市之间修建道路需要不同的花费.Louis希望建造最少 ...

  8. [HNOI2010]CITY 城市建设

    问题: 给一张图,支持边长度修改,求MST 题解: 自己想就想不到了.. 考虑cdq分治 1.首先求出一定有用的边 对于未处理的边,全部设为-INF,求一次MST,出现在MST上的边一定最终出现在后面 ...

  9. 【HNOI2010】城市建设(对时间分治 & Kruskal)

    Description \(n\) 个点 \(m\) 条边的带边权无向图.\(q\) 次操作,每次修改一条边的权值. 求每次修改后的最小生成树的边权和. Hint \(1\le n\le 2\time ...

随机推荐

  1. iOS APP内购

    看到网上文章一大把,看了这个觉得挺不错的,谢谢 iOS大全 公众平台; 原文:http://mp.weixin.qq.com/s?__biz=MzAxMzE2Mjc2Ng==&mid=2652 ...

  2. 从头开始搭建一个Spring boot+ActiveMQ高可用分布式环境

    *:first-child { margin-top: 0 !important; } body > *:last-child { margin-bottom: 0 !important; } ...

  3. intelij Android 搭建 java 项目

    1.打开ide创建页面,点击创建项目 2.选择Android,发现下面有一段红字,解释很清楚,由于第一使用为设置Android的开发环境sdk,去设置一下就好,点击取消回到上一界面 3 配置 4.当然 ...

  4. cp 命令详解

    作用:  cp 指令用于复制文件或目录,如同时指定两个以上的文件或目录,且最后的目的地是一个已经存在的目录, 则它会把前面指定的所有文件或目录复制到此目录下, 若同时指定多个文件或目录, 而最后的目的 ...

  5. 深入理解用户权限rwx

    其实在UNIX的实现中,文件权限用12个二进制位表示,如果该位置上的值是1,表示有相应的权限,如果是0则没有相应权限第11位为SUID位,第10位为SGID位,第9位为sticky位,第8-0位对应于 ...

  6. 手撕vue-cli配置文件——config篇

    最近一直在研究webpack,突然想看看vue-cli中的webpack是如何配置,查阅了很多相关的文章,所以也想出几篇关于vue-cli配置的东西.正所谓"工欲善其事必先利其器" ...

  7. Git知识总览(一) 从 git clone 和 git status 谈起

    本篇博客是整理git相关知识的第一篇,因为之前一直是用SourceTree对Git的命令行操作用的不是特别熟,于是乎过了一遍ProGit(链接:https://git-scm.com/book/zh/ ...

  8. Oracle12c_安装1——准备工作

    1.建议用户和组 su root #切换到root groupadd oinstall #创建用户组oinstall groupadd dba #创建用户组dba useradd -g oinstal ...

  9. tomcat中使用mysql连接池的配置

    1.下载相应的jar包,添加到工程中 需要下载的包主要有commons-pool2-2.2 commons-dbcp2-2.0.1-src commons-dbcp2-2.0.1  commons-c ...

  10. redis实现分布式可重入锁

    利用redis可以实现分布式锁,demo如下: /** * 保存每个线程独有的token */ private static ThreadLocal<String> tokenMap = ...