[HNOI 2016]最小公倍数
Description
给定一张 \(N\) 个顶点 \(M\) 条边的无向图(顶点编号为 \(1,2,\cdots,n\) ),每条边上带有权值。所有权值都可以分解成 \(2^a\times 3^b\) 的形式。 \(q\) 个询问,每次询问给定四个参数 \(u,v,a,b\) ,请你求出是否存在一条顶点 \(u\) 到 \(v\) 之间的路径,使得路径依次经过的边上的权值的最小公倍数为 \(2^a\times 3^b\) 。
\(1\leq n,q\leq 50000,1\leq m\leq 100000\)
Solution
容易发现题中要求的就是找一条可以是非简单路径的路径连接 \(u,v\) 并且边权 \(a'\) 最大值为 \(a\) ,边权 \(b'\) 最大值为 \(b\) 。
考虑暴力就是对于每组询问,将边的 \(a\) 值小于等于询问的 \(a\) 值且 \(b\) 值小于等于询问的 \(b\) 值的边加入图中。判断 \(u,v\) 是否联通,且联通块内 \(a_{max},b_{max}\) 是否与询问的值相等。
显然这样的复杂度是假的。优化暴力,考虑分块。
边按 \(a\) 排序,然后分块。对于每一块 \(i\) ,处理 \(a'\) 在这一块中的询问。这时候之前块的 \(a<a'\) 这一个关系一定满足,按 \(b\) 排序后也满足 \(b<b'\) 了。
块 \(i\) 中还有一些满足的,最多 \(\sqrt m\) 。
暴力加入然后撤销就可以了。不路径压缩的并查集,按秩合并的话复杂度 \(log_2 n\) 。
总复杂度 \(O(m\sqrt m log_2 n)\) 。
Code
//It is made by Awson on 2018.3.3
#include <bits/stdc++.h>
#define LL long long
#define dob complex<double>
#define Abs(a) ((a) < 0 ? (-(a)) : (a))
#define Max(a, b) ((a) > (b) ? (a) : (b))
#define Min(a, b) ((a) < (b) ? (a) : (b))
#define Swap(a, b) ((a) ^= (b), (b) ^= (a), (a) ^= (b))
#define writeln(x) (write(x), putchar('\n'))
#define lowbit(x) ((x)&(-(x)))
using namespace std;
const int N = 50000, M = 100000;
void read(int &x) {
char ch; bool flag = 0;
for (ch = getchar(); !isdigit(ch) && ((flag |= (ch == '-')) || 1); ch = getchar());
for (x = 0; isdigit(ch); x = (x<<1)+(x<<3)+ch-48, ch = getchar());
x *= 1-2*flag;
}
void print(int x) {if (x > 9) print(x/10); putchar(x%10+48); }
void write(int x) {if (x < 0) putchar('-'); print(Abs(x)); }
int n, m, q, lim, ans[N+5];
struct tt {int u, v, a, b, id; }edge[M+5], query[N+5], tmp[N+5], pre[N+5];
bool comp1(const tt &a, const tt &b) {return a.a == b.a ? a.b < b.b : a.a < b.a; }
bool comp2(const tt &a, const tt &b) {return a.b == b.b ? a.a < b.a : a.b < b.b; }
struct Set {
int fa[N+5], ma[N+5], mb[N+5], rank[N+5], pos;
void clear() {for (int i = 1; i <= n; i++) fa[i] = rank[i] = 0, ma[i] = mb[i] = -1; }
int find(int u) {return fa[u] ? find(fa[u]) : u; }
void merge(int u, int v, int a, int b) {
u = find(u), v = find(v); if (rank[u] > rank[v]) Swap(u, v);
ma[v] = max(ma[v], max(ma[u], a)), mb[v] = max(mb[v], max(mb[u], b));
if (u != v) fa[u] = v, rank[v] += (rank[v] == rank[u]);
}
void unin(int u, int v, int a, int b) {
u = find(u), v = find(v); if (rank[u] > rank[v]) Swap(u, v); ++pos;
pre[pos].u = u, pre[pos].v = v, pre[pos].a = ma[v], pre[pos].b = mb[v], pre[pos].id = rank[v];
ma[v] = max(ma[v], max(ma[u], a)), mb[v] = max(mb[v], max(mb[u], b));
if (u != v) fa[u] = v, rank[v] += (rank[v] == rank[u]);
}
void undo() {
while (pos) {
int u = pre[pos].u, v = pre[pos].v, a = pre[pos].a, b = pre[pos].b, rk = pre[pos].id;
if (u != v) fa[u] = 0, rank[v] = rk;
ma[v] = a, mb[v] = b; --pos;
}
}
}T;
void work() {
read(n), read(m); lim = sqrt(m);
for (int i = 1; i <= m; i++) read(edge[i].u), read(edge[i].v), read(edge[i].a), read(edge[i].b);
read(q);
for (int i = 1; i <= q; i++) read(query[i].u), read(query[i].v), read(query[i].a), read(query[i].b), query[i].id = i;
sort(edge+1, edge+1+m, comp1); sort(query+1, query+1+q, comp2);
for (int i = 1; i <= m; i = i+lim) {
int cnt = 0, k = 1, ed = Min(i+lim-1, m);
for (int j = 1; j <= q; j++) if (query[j].a >= edge[i].a && (i+lim > m || query[j].a < edge[i+lim].a)) tmp[++cnt] = query[j];
sort(edge+1, edge+i, comp2); T.clear();
for (int j = 1; j <= cnt; j++) {
while (k < i && edge[k].b <= tmp[j].b) T.merge(edge[k].u, edge[k].v, edge[k].a, edge[k].b), ++k;
for (int q = i; q <= ed; q++) if (edge[q].a <= tmp[j].a && edge[q].b <= tmp[j].b) T.unin(edge[q].u, edge[q].v, edge[q].a, edge[q].b);
int u = T.find(tmp[j].u), v = T.find(tmp[j].v);
if (u == v && T.ma[u] == tmp[j].a && T.mb[u] == tmp[j].b) ans[tmp[j].id] = 1;
T.undo();
}
}
for (int i = 1; i <= q; i++) puts(ans[i] ? "Yes" : "No");
}
int main() {
work(); return 0;
}
[HNOI 2016]最小公倍数的更多相关文章
- [BZOJ 4537][Hnoi 2016]最小公倍数
传送门 并查集+分块 看到题目可以想到暴力做法, 对于每个询问, 将所有a和b小于等于询问值的的边加入图中(用并查集), 如果询问的u和v在一个联通块中, 且该联通块的maxa和maxb均等与询问的a ...
- [HNOI 2016]树
Description 题库链接 给你一棵 \(N\) 个节点根节点为 \(1\) 的有根树,结点的编号为 \(1\sim N\) :我们称这颗树为模板树.需要通过这棵模板树来构建一颗大树.构建过程如 ...
- 【BZOJ 4539】【HNOI 2016】树
http://www.lydsy.com/JudgeOnline/problem.php?id=4539 今天测试唯一会做的一道题. 按题目要求,如果暴力的把模板树往大树上仍,最后得到的大树是$O(n ...
- hnoi 2016 省选总结
感觉省选好难的说...反正我数据结构太垃圾正解想到了也打不出来打一打暴力就滚粗了! DAY1 0+20+30 DAY2 60-20+0+60 最后170-20分,暴力分还是没有拿全! 然而这次是给了5 ...
- HNOI 2016 省队集训日记
第一天 DeepDarkFantasy 从东京出发,不久便到一处驿站,写道:日暮里. ——鲁迅<藤野先生> 定义一个置换的平方为对1~n的序列做两次该置换得到的序列.已知一个置换的平方, ...
- 数据结构(树链剖分,堆):HNOI 2016 network
2215. [HNOI2016]网络 ★★★☆ 输入文件:network_tenderRun.in 输出文件:network_tenderRun.out 简单对比时间限制:2 s 内存 ...
- [HNOI 2016]序列
Description 题库链接 给你一个长度为 \(n\) 的序列 \(A\) ,给出 \(q\) 组询问.每次询问 \([l,r]\) ,求该区间内所有的子序列中最小值的和. \(1\leq n, ...
- [HNOI 2016]大数
Description 题库链接 给你一个长度为 \(n\) ,可含前导零的大数,以及一个质数 \(p\) . \(m\) 次询问,每次询问你一个大数的子区间 \([l,r]\) ,求出子区间中有多少 ...
- [HNOI 2016]网络
Description 一个简单的网络系统可以被描述成一棵无根树.每个节点为一个服务器.连接服务器与服务器的数据线则看做 一条树边.两个服务器进行数据的交互时,数据会经过连接这两个服务器的路径上的所有 ...
随机推荐
- Java基础学习笔记十三 常用API之正则表达式、Date、DateFormat、Calendar
正则表达式 正则表达式(英语:Regular Expression,在代码中常简写为regex).正则表达式是一个字符串,使用单个字符串来描述.用来定义匹配规则,匹配一系列符合某个句法规则的字符串.在 ...
- 云计算--网络原理与应用--20171123--网络地址转换NAT
NAT的概述 NAT的配置 实验 一. NAT的概述 NAT(Network address translation,网络地址转换)通过将内部网络的的私有地址翻译成全球唯一的共有网络IP地址,是内部网 ...
- img之间的间隙问题
前言:关于基线(base line),中线(middle line),行高(line height)的了解还是比较浅的,所以引用前辈的成果,稍带解释下 1)行高:两行文字之间"基线" ...
- 位图(BitMap)
位图 引入:2.给40亿个不重复的无符号整数,没排过序.给一个无符号整数,如何快速判断一个数是否在这40亿个数中. [腾讯] 思路:如果内存够的话,40亿个整型使用位图存储需要500M左右的空间. 分 ...
- continue和break的特殊用法。
break在程序中一般来说的作用就是跳出当前循环,然后再据需执行循环外的语句.continue也是对当前循环来说直接进入到下一次循环.其实我们在程序中有时候循环体嵌套太多,进行到某一步是希望直接bre ...
- css变化代码2
<!DOCTYPE html><html> <head> <meta charset="utf-8" /> ...
- Ubuntu安装使用latex
TeX Live is a TeX distribution to get up and running with the TeX document production system. To ins ...
- Linux的rsync 配置,用于服务器之间远程传大量的数据
[教程主题]:rsync [课程录制]: 创E [主要内容] [1] rsync介绍 Rsync(Remote Synchronize) 是一个远程资料同步工具,可通过LAN/WAN快速同步多台主机, ...
- MongoDb进阶实践之五 MongoDB修改命令详述
一.引言 上一篇文章我们已经详细介绍了MongoDB数据库的有关查询的内容,但是这只是所有查询命令的冰山一角.所有查询命令都写完也没有必要,我只是写了一些常用的命令,对MongoDB的 ...
- Python内置函数(13)——bytearray
英文文档: class bytearray([source[, encoding[, errors]]]) Return a new array of bytes. The bytearray cla ...