$分块$

$一个很有趣的技巧$

$在树上选sqrt(n)个关键点,每两个关键点之间的距离<=sqrt(n),每个关键点属于一条链$

$预处理出每两个关键点的bitset$

$每次询问就暴力向上爬,合并bitset$

$由于要查询,所以要手写bitset$

#include<bits/stdc++.h>
using namespace std;
typedef unsigned long long ull;
const int bl = , N = 1e5 + ;
ull Num[];
vector<int> G[N];
int n, m, f, bbl, last;
ull Max;
int fa[N][], dep[N], mx[N], w[N], mark[N], top[N], mir[N];
int read() {
int x = , f = ;
char c = getchar();
while(!isdigit(c)) {
if(c == '-') {
f = -;
}
c = getchar();
}
while(isdigit(c)) {
x = x * + c - '';
c = getchar();
}
return x * f;
}
struct Bitset {
int len;
ull s[bl];
void clear() {
for(int i = ; i <= len; ++i) {
s[i] = ;
}
len = ;
}
void operator |= (const Bitset &a) {
len = max(len, a.len);
for(int i = ; i <= len; ++i) {
s[i] |= a.s[i];
}
}
void operator |= (int b) {
len = max(len, b >> );
s[b >> ] |= 1LL << (b & );
}
int num() {
int ret = ;
for(int i = ; i <= len; ++i) {
ret += Num[s[i] >> ] + Num[(s[i] >> ) & ] + Num[(s[i] >> ) & ] + Num[s[i] & ];
}
return ret;
}
int mex() {
for(int i = ; i <= len; ++i) {
if(s[i] != Max) {
for(int j = ; j < ; ++j) {
if(!(s[i] & (1LL << j))) {
return i * + j;
}
}
}
}
}
} b[][], ans;
void dfs(int u) {
mx[u] = dep[u];
for(int i = ; i < G[u].size(); ++i) {
int v = G[u][i];
if(v == fa[u][]) {
continue;
}
fa[v][] = u;
dep[v] = dep[u] + ;
dfs(v);
mx[u] = max(mx[u], mx[v]);
}
if(mx[u] - dep[u] >= bbl || u == ) {
mark[u] = ++mark[];
mir[mark[u]] = u;
mx[u] = ;
}
}
int lca(int u, int v) {
if(dep[u] < dep[v]) {
swap(u, v);
}
int d = dep[u] - dep[v];
for(int i = ; i < ; ++i) {
if(d & ( << i)) {
u = fa[u][i];
}
}
if(u == v) {
return u;
}
for(int i = ; i >= ; --i) {
if(fa[u][i] != fa[v][i]) {
u = fa[u][i];
v = fa[v][i];
}
}
return fa[u][];
}
int main() {
n = read();
m = read();
f = read();
bbl = sqrt(n);
for(int i = ; i < ; ++i) {
Max |= 1LL << i;
}
for(int i = ; i < ; ++i) {
for(int j = ; j < ; ++j) {
if(i & ( << j)) {
++Num[i];
}
}
}
for(int i = ; i <= n; ++i) {
w[i] = read();
}
for(int i = ; i < n; ++i) {
int u = read(), v = read();
G[u].push_back(v);
G[v].push_back(u);
}
dep[] = ;
dfs();
for(int j = ; j <= ; ++j) {
for(int i = ; i <= n; ++i) {
fa[i][j] = fa[fa[i][j - ]][j - ];
}
}
for(int i = ; i <= mark[]; ++i) {
ans.clear();
ans |= w[mir[i]];
int x;
b[i][i] |= w[mir[i]];
for(x = fa[mir[i]][]; x; x = fa[x][]) {
ans |= w[x];
if(mark[x]) {
b[i][mark[x]] = ans;
if(!top[mir[i]]) {
top[mir[i]] = x;
}
}
}
}
while(m--) {
int a = read(), x, y, z;
ans.clear();
while(a--) {
x = read();
y = read();
if(f) {
x = x ^ last;
y = y ^ last;
}
ans |= w[x];
ans |= w[y];
z = lca(x, y);
for(; x && !mark[x] && dep[x] > z; x = fa[x][]) {
ans |= w[x];
}
for(; y && !mark[y] && dep[y] > z; y = fa[y][]) {
ans |= w[y];
}
int xx, yy;
for(xx = x; dep[top[xx]] >= dep[z]; xx = top[xx]);
for(yy = y; dep[top[yy]] >= dep[z]; yy = top[yy]);
ans |= b[mark[x]][mark[xx]];
ans |= b[mark[y]][mark[yy]];
for(x = xx; dep[x] >= dep[z]; x = fa[x][]) {
ans |= w[x];
}
for(y = yy; dep[y] >= dep[z]; y = fa[y][]) {
ans |= w[y];
}
}
int ta = ans.num(), tb = ans.mex();
last = ta + tb;
printf("%d %d\n", ta, tb);
}
return ;
}

最新文章

  1. JS中的常量(基本数据类型)和内置对象
  2. 洛谷P3379 【模板】最近公共祖先(LCA)
  3. 怎样用命令行编译C#程序
  4. Hibernate包及相关工具包下载地址
  5. shell常识总结
  6. Android 纯代码加入点击效果
  7. 13.如何生成订单号,用uuid
  8. (转发)RequestDispatcher的include()方法和forward()方法的区别
  9. 浏览器加载 CommonJS 模块的原理与实现
  10. localhost和127.0.0.1的区别
  11. 數據監控与診斷--環形緩沖區(RING BUFFER)
  12. 使用korofileheader插件vs code添加文件头注释和函数注释
  13. 移动端reset
  14. Stanford机器学习---第十一讲.异常检测
  15. Xshell中vim退出内容仍停留在屏幕的问题
  16. [React] Validate React Forms with Formik and Yup
  17. Goroutine(协程)为何能处理大并发?
  18. 使用 MSBuild 响应文件 (rsp) 来指定 dotnet build 命令行编译时的大量参数
  19. application/x-www-form-urlencoded接口响应报文中文乱码
  20. Python 函数 min()

热门文章

  1. 获取网站的BaseURL
  2. [python学习] 简单爬取图片站点图库中图片
  3. 关于Gradle和Gradle插件的问题
  4. Elasticsearch集群UNASSIGNED
  5. async &amp; await (转载)
  6. 九度OJ 1037:Powerful Calculator(强大的计算器) (大整数运算)
  7. Brotli
  8. cocos2d-js添加360广告联盟插屏(通过jsb反射机制)
  9. ES5中的类与继承
  10. 后缀自动机的python实现