链接

loj

一道阔爱的分块

题意

边权是二元组(A, B),每次询问u, v, a, b,求u到v是否存在一条简单路径,使得各边权上\(A_{max} = a, B_{max} = b\)

分析

对于这种有两种限制的题目

一般的套路就是条件按照第一种权值为关键字排序,询问按照第二种关键字排序

然后给条件分块,然后对于一个块只把第一关键字符合条件的询问放进去

在把当前块前面的整块里的点按照第二关键字排序

这样当前块前面的点都是符合当前询问点对于第一关建字条件的

而且第二关键字都是单调的,所以扫一下

然后对于每个询问,暴力处理一下当前块的贡献

[参考BeNoble_] (https://blog.csdn.net/benoble_/article/details/79777757)

对于这道题 如果暴力怎么做?

对于询问u, v, a, b

把所有满足A<= a, B <= b的边加进来

因为只要最大值,所以可以维护一个带权并查集(find的时候不更新father哦)

(小声:其实这个带权并查集就像一个树一样

然后查询一下是否连通,连通的话所在并查集最大权满不满足条件(即Amax == a && Bmax == b)

所以说,分块的本质都是暴力

然后复杂度就在这个加边上了

就像最前面说的那样加 加完把不整块的删了 okk

框架就像这样

询问按b排序 边按a排序

for(对于每一个块){
收集a大小在该块范围内的询问 按b排前面整块的点(这样后面就单调了 初始化并查集 for(对于每一个询问){
加前面整块里 b满足条件的边(a必然满足条件) 加不整块的a,b都满足条件的边 判定是否联通且满足条件 还原不整块的边
}
}

最后附上代码

#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <vector>
#include <set>
using namespace std;
const int N = 1e5 + 5; int n, m, qsize;
struct E{
int u, v, a, b, id;
}e[N], q[N];
struct Opt{
int u, v, a, b, fa, size;
}opt[N];
bool ans[N];
int blsize;
int b[N], bcnt, tim;
int fa[N], mxa[N], mxb[N], size[N]; bool rulea(E x, E y){return x.a < y.a;}
bool ruleb(E x, E y){return x.b < y.b;} int find(int x){
return x == fa[x] ? x : find(fa[x]);
} inline void merge(E x, bool type){
x.u = find(x.u); x.v = find(x.v);
if(size[x.u] > size[x.v]) swap(x.u, x.v);
//printf("do %d %d %d\n", x.u, x.v, type);
if(type) {
++tim;
opt[tim].u = x.u, opt[tim].v = x.v; opt[tim].fa = fa[x.u];
opt[tim].a = mxa[x.v], opt[tim].b = mxb[x.v];
opt[tim].size = size[x.v];
} if(find(x.u) == find(x.v)){
mxa[x.v] = max(mxa[x.v], x.a);
mxb[x.v] = max(mxb[x.v], x.b);
}
else {// u -> v
fa[x.u] = x.v;
size[x.v] += size[x.u];
mxa[x.v] = max(mxa[x.v], x.a);
mxb[x.v] = max(mxb[x.v], x.b);
mxa[x.v] = max(mxa[x.v], mxa[x.u]);
mxb[x.v] = max(mxb[x.v], mxb[x.u]);
}
} inline void undo(){
while(tim){
// printf("undo %d %d\n", opt[tim].u, opt[tim].v);
fa[opt[tim].u] = opt[tim].fa;
mxa[opt[tim].v] = opt[tim].a;
mxb[opt[tim].v] = opt[tim].b;
size[opt[tim].v] = opt[tim].size;
--tim;
}
} int main(){
scanf("%d%d", &n, &m);
blsize = sqrt(15 * m);//
for(int i = 1; i <= m; ++i){
scanf("%d%d%d%d", &e[i].u, &e[i].v, &e[i].a, &e[i].b);
}
sort(e + 1, e + m + 1, rulea);
scanf("%d", &qsize);
for(int i = 1; i <= qsize; ++i){
scanf("%d%d%d%d", &q[i].u, &q[i].v, &q[i].a, &q[i].b);
q[i].id = i;
}
sort(q + 1, q + qsize + 1, ruleb); // for(int i = 1; i <= m; ++i)
// printf("%d %d %d %d\n", e[i].u, e[i].v, e[i].a, e[i].b); for(int i = 1, lim; i <= m; i += blsize){
bcnt = 0;
lim = min(m, i + blsize - 1);
for(int j = 1; j <= qsize; ++j)
if(e[i].a <= q[j].a
&& (i + blsize > m || e[i + blsize].a > q[j].a))
b[++bcnt] = j;
sort(e + 1, e + i, ruleb);//这里排的是前面整块的点!
for(int j = 1; j <= n; ++j){
fa[j] = j, size[j] = 1, mxa[j] = mxb[j] = -1;
}
for(int j = 1, top = 1; j <= bcnt; ++j){
while(top < i && q[b[j]].b >= e[top].b){
//printf("mer %d 0\n");
merge(e[top], 0);
++top;
}
for(int k = i; k <= lim; ++k){
if(q[b[j]].a >= e[k].a && q[b[j]].b >= e[k].b){
merge(e[k], 1);
}
} int x = find(q[b[j]].u), y = find(q[b[j]].v);
//printf("%d %d %d %d %d %d %d\n", q[b[j]].u, q[b[j]].v, q[b[j]].id, mxa[x], mxb[x], x, y);
ans[q[b[j]].id] = ((x == y) && (mxa[x] == q[b[j]].a) && (mxb[x] == q[b[j]].b)); undo();
}
}
for(int i = 1; i <= qsize; ++i)
if(ans[i]) printf("Yes\n");
else printf("No\n");
return 0;
}

「HNOI2016」最小公倍数的更多相关文章

  1. 「HNOI2016」最小公倍数 解题报告

    「HNOI2016」最小公倍数 考虑暴力,对每个询问,处理出\(\le a,\le b\)的与询问点在一起的联通块,然后判断是否是一个联通块,且联通块\(a,b\)最大值是否满足要求. 然后很显然需要 ...

  2. loj2048 「HNOI2016」最小公倍数

    link 题意: 给定一张$N$个顶点$M$条边的无向图(顶点编号为$1,2,...,n$),每条边上带有权值.所有权值都可以分解成$2^a \cdot 3^b$的形式. 现在有$q$个询问,每次询问 ...

  3. LOJ #2048. 「HNOI2016」最小公倍数

    题意 有 \(n\) 个点,\(m\) 条边,每条边连接 \(u \Leftrightarrow v\) 且权值为 \((a, b)\) . 共有 \(q\) 次询问,每次询问给出 \(u, v, q ...

  4. 「HNOI2016」数据结构大毒瘤

    真是 \(6\) 道数据结构毒瘤... 开始口胡各种做法... 「HNOI2016」网络 整体二分+树状数组. 开始想了一个大常数 \(O(n\log^2 n)\) 做法,然后就被卡掉了... 发现直 ...

  5. 「HNOI2016」树 解题报告

    「HNOI2016」树 事毒瘤题... 我一开始以为每次把大树的子树再接给大树,然后死活不知道咋做,心想怕不是个神仙题哦 然后看题解后才发现是把模板树的子树给大树,虽然思维上难度没啥了,但是还是很难写 ...

  6. 「HNOI2016」序列 解题报告

    「HNOI2016」序列 有一些高妙的做法,懒得看 考虑莫队,考虑莫队咋移动区间 然后你在区间内部找一个最小值的位置,假设现在从右边加 最小值左边区间显然可以\(O(1)\),最小值右边的区间是断掉的 ...

  7. 「HNOI2016」网络 解题报告

    「HNOI2016」网络 我有一个绝妙的可持久化树套树思路,可惜的是,它的空间是\(n\log^2 n\)的... 注意到对一个询问,我们可以二分答案 然后统计经过这个点大于当前答案的路径条数,如果这 ...

  8. loj #2051. 「HNOI2016」序列

    #2051. 「HNOI2016」序列 题目描述 给定长度为 n nn 的序列:a1,a2,⋯,an a_1, a_2, \cdots , a_na​1​​,a​2​​,⋯,a​n​​,记为 a[1: ...

  9. 【LOJ】#2052. 「HNOI2016」矿区

    题解 之前尝试HNOI2016的时候弃坑的一道,然后给补回来 (为啥我一些计算几何就写得好长,不过我写啥都长orz) 我们尝试给这个平面图分域,好把这个平面图转成对偶图 怎么分呢,我今天也是第一次会 ...

随机推荐

  1. 1171: lfx捧杯稳啦!

    escription Lfx在复习离散的时候突然想到了一个算法题,毕竟是lfx, 算法题如下: 他想知道这样的问题,先定义1~n中即是3的倍数,又是11的倍数的那些数的和sum, 他想知道sum有多少 ...

  2. p211有界自共轭算子T是实数集合的子集

    对条件 取非   是 ∉谱集合的实数 才对  现在是  入  属于正则点集 他  然后  又说T 的谱是实数 这不矛盾吗 这里根据   必要性    推出 蓝色和红色矛盾    矛盾就是 这是谱点 然 ...

  3. latex中插入eps文件

    \documentclass{article} \usepackage{graphicx}\usepackage{epstopdf} \begin{document}\begin{figure}  \ ...

  4. MySQL之优化

    当 MySQL 单表记录数过大时,增删改查性能都会急剧下降,本文会提供一些优化参考,大家可以参考以下步骤来优化. 一. 单表优化 除非单表数据未来会一直不断上涨,否则不要一开始就考虑拆分,拆分会带来逻 ...

  5. PHP的优化建议(仅借鉴)

    转载: https://www.awaimai.com/1050.html 1 字符串 1.1 少用正则表达式 能用PHP内部字符串操作函数的情况下,尽量用他们,不要用正则表达式, 因为其效率高于正则 ...

  6. React-Native之截图组件view-shot的介绍与使用

    React-Native之截图组件view-shot的介绍与使用 一,需求分析 1,需要将分享页生成图片,并分享到微信好友与朋友圈. 二,react-native-view-shot介绍 1,可以截取 ...

  7. web项目中的监听器,过滤器以及自定义servlet的执行顺序

    可以看到web容器一启动就会实例化监听器的contextInitialized(ServletContextEvent event)方法,然后是过滤器的init()方法,最后在用户访问web应用的 时 ...

  8. Node & CLI

    Node & CLI cli 生成文件的原理是什么 https://nodejs.org/api/cli.html http://nodejs.cn/api/cli.html CLI & ...

  9. Scrapy网络爬虫框架的开发使用

    1.安装 2.使用scrapy startproject  project_name 命令创建scrapy项目 如图: 3.根据提示使用scrapy genspider spider_name dom ...

  10. 解析xml文件 selectSingleNode取不到节点

    今天在做批量生成XML的时候,碰到一个情况 解析xml文件 selectSingleNode一直返回NULL. XML的格式开头有一句这个<CE401Message xmlns="ht ...