bzoj 4568 [SCOI 2016] 幸运数字

最新推荐文章于 2018-04-22 16:59:26 发布
最新推荐文章于 2018-04-22 16:59:26 发布
阅读量91 收藏
点赞数
原文链接：http://www.cnblogs.com/acha/p/7145254.html
版权

题目大意

给定一棵\(n\)个点的树，每个点有权值
\(q\)次询问树上路径中
每个点权值可选可不选的最大异或和
\(n\le 2*10^4,q\le 2*10^5,val[i]\le 2^{60}\)

分析

线性基
但数据范围不太对啊woc
两线性基合并\(O(60^2)\)
如果 树剖+线性基 则\(O(q\log^2n ~60^2)\)
如果 树上倍增+线性基 则是一个常数较大的\(O(q\log n~60^2)\)
不是很妙啊

做法1

注意到，这不是常规的树上询问，因为是异或(线性基)
线性基同一个值重复插入时没有问题的
即树上某一段算重是可以的
那么我们可以这样算
1086046-20170710110919728-505044825.jpg
其中橙色部分时一个长度\(2^k\)的段
这样就只用统计\(4\)次了
复杂度\(O(4*60^2~q)\)

做法2

点分治
求出重心到每个点的线性基\(O(60n)\)
然后解决询问
如果询问与该重心有关，单次\(O(60^2)\)
否则将该询问传到儿子
点分每层扫过的询问数\(O(q)\)的
总复杂度
\(O(60n\log n~+~q \log n~+~60^2 q)\)

solution 1

#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cctype>
#include <cmath>
#include <algorithm>
using namespace std;
const int M=2e4+7;
const int B=60;
typedef long long LL;
 
inline int rd(){
    int x=0;bool f=1;char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c=='-') f=0;
    for(;isdigit(c);c=getchar()) x=x*10+c-48;
    return f?x:-x;
}
 
inline LL lrd(){
    LL x=0;bool f=1;char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c=='-') f=0;
    for(;isdigit(c);c=getchar()) x=x*10+c-48;
    return f?x:-x;
}
 
struct vec{
    int g[M],te;
    struct edge{
        int y,nxt;
        edge(int _y=0,int _nxt=0){
            y=_y,nxt=_nxt;
        }
    }e[M<<1];
    vec(){memset(g,0,sizeof g);te=0;}
    inline void push(int x,int y){e[++te]=edge(y,g[x]);g[x]=te;}
    inline void push2(int x,int y){push(x,y);push(y,x);}
    inline int& operator () (int x){return g[x];}
    inline edge& operator [] (int x){return e[x];}
}e;
 
struct Base{
    LL a[B+3];
    Base(){memset(a,0,sizeof a);}
    void ins(LL x){
        for(int i=B;i>=0;i--) if(x>>i&1){
            if(a[i]) x^=a[i];
            else {a[i]=x;return;}
        }
    }
 
    LL getmx(){
        LL res=0; int i;
        for(i=B;i>=0;i--) if((res^a[i])>res) res^=a[i];
        return res;
    }
 
    friend Base merge(const Base &x,const Base &y){
        Base res; int i;
        for(i=B;i>=0;i--) if(x.a[i]) res.ins(x.a[i]);
        for(i=B;i>=0;i--) if(y.a[i]) res.ins(y.a[i]);
        return res;
    }
 
    friend Base merge(const Base &x,const LL &y){
        Base res=x;
        res.ins(y);
        return res;
    }
};
 
int n,m,D;
LL val[M];
int dep[M];
int pre[M][15];
Base f[M][15];
 
void dfs(int x){
    int p,y;
    for(p=e(x);p;p=e[p].nxt)
        if((y=e[p].y)!=pre[x][0]){
            dep[y]=dep[x]+1;
            pre[y][0]=x;
            f[y][0].ins(val[y]);
            dfs(y);
        }
}
 
void init(){
    D=(int)log2(n);
    int i,j;
    for(j=1;j<=D;j++)
        for(i=1;i<=n;i++){
            pre[i][j]=pre[pre[i][j-1]][j-1];
            f[i][j]=merge(f[i][j-1],f[pre[i][j-1]][j-1]);
        }
}
 
int LCA(int x,int y){
    int i;
    if(dep[x]<dep[y]) swap(x,y);
    for(i=D;i>=0;i--)
        if(dep[pre[x][i]]>=dep[y]) x=pre[x][i];
    if(x==y) return x;
    for(i=D;i>=0;i--)
        if(pre[x][i]!=pre[y][i]) x=pre[x][i],y=pre[y][i];
    return pre[x][0];
}
 
int jp(int x,int kth){
    for(int i=D;i>=0;i--) if(kth>>i&1) x=pre[x][i];
    return x;
}
 
Base calc(int x,int to){
    for(int i=D;i>=0;i--) if(dep[pre[x][i]]>=dep[to]){
        int y=jp(x,dep[pre[x][i]]-dep[to]);
        return merge(f[x][i],f[y][i]);
    }
    return Base();
}
 
LL get(int x,int y){
    Base res;
    int lca=LCA(x,y),i;
    res=merge(merge(calc(x,lca),calc(y,lca)),val[lca]);
    return res.getmx();
}
 
int main(){
 
    int i,x,y;
 
    n=rd(),m=rd();
    for(i=1;i<=n;i++) val[i]=lrd();
 
    for(i=1;i<n;i++) e.push2(rd(),rd());
 
    pre[1][0]=0; dep[1]=1;
    dfs(1);
    init();
 
    for(i=1;i<=m;i++){
        x=rd(),y=rd();
        printf("%lld\n",get(x,y));
    }
 
    return 0;
}

solution 2

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <cmath>
#include <algorithm>
using namespace std;
const int N=2e5+7;
const int M=2e4+7;
const int B=60;
typedef long long LL;
 
inline int ri(){
    int x=0;bool f=1;char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c=='-') f=0;
    for(;isdigit(c);c=getchar()) x=x*10+c-48;
    return f?x:-x;
}
 
inline LL rl(){
    LL x=0;bool f=1;char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c=='-') f=0;
    for(;isdigit(c);c=getchar()) x=x*10+c-48;
    return f?x:-x;
}
 
struct vec{
    int g[M],te;
    struct edge{
        int y,nxt;
        edge(int _y=0,int _nxt=0){y=_y,nxt=_nxt;}
    }e[M<<1];
    vec(){memset(g,0,sizeof g);te=0;}
    inline void push(int x,int y){e[++te]=edge(y,g[x]);g[x]=te;}
    inline void push2(int x,int y){push(x,y);push(y,x);}
    inline int& operator () (int x){return g[x];}
    inline edge& operator [] (int x){return e[x];}
}e;
 
struct vec2{
    int g[M],te;
    struct ques{
        int x,y,id,nxt;
        ques(int _x=0,int _y=0,int _id=0,int _nxt=0){x=_x,y=_y,id=_id,nxt=_nxt;}
    }e[N*17];
    vec2(){memset(g,0,sizeof g);te=0;}
    inline void push(int u,int x,int y,int id){e[++te]=ques(x,y,id,g[u]);g[u]=te;}
    inline int& operator () (int x){return g[x];}
    inline ques& operator [] (int x){return e[x];}
}ask;
 
struct Base{
    LL a[B+3];
 
    Base(){memset(a,0,sizeof a);}
    void clear(){memset(a,0,sizeof a);}
 
    void ins(LL x){
        for(int i=B;i>=0;i--) if(x>>i&1){
            if(a[i]) x^=a[i];
            else {a[i]=x;break;}
        }
    }
 
    LL getmx(){
        LL res=0;
        for(int i=B;i>=0;i--) if((res^a[i])>res) res^=a[i];
        return res;
    }
 
    friend Base merge(const Base &x,const Base &y){
        Base res;
        for(int i=B;i>=0;i--) if(x.a[i]) res.ins(x.a[i]);
        for(int i=B;i>=0;i--) if(y.a[i]) res.ins(y.a[i]);
        return res;
    }
 
    friend Base merge(const Base &x,const LL &y){
        Base res=x;
        res.ins(y);
        return res;
    }
};
 
int n,m;
LL val[M];
int sz[M];
int mi,size,rt;
bool vis[M];
int bl[M];
LL ans[N];
Base f[M];
 
LL calc(int x,int y){
    Base res=merge(f[x],f[y]);
    return res.getmx();
}
 
void getsz(int x,int fa){
    sz[x]=1;
    int p,y;
    for(p=e(x);p;p=e[p].nxt)
    if((y=e[p].y)!=fa&&!vis[y]){
        getsz(y,x);
        sz[x]+=sz[y];
    }
}
 
void getrt(int x,int fa){
    int f=size-sz[x],p,y;
    for(p=e(x);p;p=e[p].nxt)
    if((y=e[p].y)!=fa&&!vis[y]){
        getrt(y,x);
        f=max(f,sz[y]);
    }
    if(f<mi) mi=f,rt=x;
}
 
void dfs(int x,int fa){
    f[x]=merge(f[fa],val[x]);
    int p,y;
    for(p=e(x);p;p=e[p].nxt)
    if((y=e[p].y)!=fa&&!vis[y]) bl[y]=bl[x],dfs(y,x);
}
 
void work(int fr){
    getsz(fr,0);
    mi=size=sz[fr];
    getrt(fr,0);
    int x=rt,i,p,y;
    vis[x]=1;
 
    bl[x]=x;
    f[x].clear(); f[x].ins(val[x]);
    if(x!=fr) {ask(x)=ask(fr); ask(fr)=0;}//
 
    for(p=e(x);p;p=e[p].nxt)
    if(!vis[y=e[p].y]){
        bl[y]=y;
        dfs(y,x);
    }
 
    for(p=ask(x);p;p=ask[p].nxt){
        if(bl[ask[p].x]==bl[ask[p].y]) ask.push(bl[ask[p].x],ask[p].x,ask[p].y,ask[p].id);
        else ans[ask[p].id]=calc(ask[p].x,ask[p].y);
    }
     
    for(p=e(x);p;p=e[p].nxt)
    if(!vis[y=e[p].y]) work(y);
}
 
int main(){
 
    int i,x,y;
    n=ri(),m=ri();
    for(i=1;i<=n;i++) val[i]=rl();
 
    for(i=1;i<n;i++) e.push2(ri(),ri());
 
    for(i=1;i<=m;i++){
        x=ri(), y=ri();
        if(x==y) ans[i]=val[x];
        else ask.push(1,x,y,i);
    }
 
    work(1);
 
    for(int i=1;i<=m;i++) printf("%lld\n",ans[i]);
 
    return 0;
}

转载于:https://www.cnblogs.com/acha/p/7145254.html

weixin_30810239
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
BZOJ全代码
02-05
BZOJ平台全部代码，解压到一个文件夹在打开使用。BZOJ平台全部代码，解压到一个文件夹在打开使用。
BZOJ第一部分
11-03
BZOJ原题-BZOJP1000-P2000的题目，下载后可以离线做题。
bzoj4568: [Scoi2016]幸运数字
miaom的博客
05-24 1158
一句话题意：树上两点间路径子集最大异或和。 看了一下线性基发现好强，但我也不太会，先贴代码吧。。 #include #include #include #define N 20005 #define M 200005 #define Mlg 6000005 #define ll long long using namespace std; bool ok[N<<1]; int n,m,x,y,
BZOJ 4568 [Scoi2016]幸运数字 【倍增线性基
Flaze_的博客
07-16 492
嗯今年的省选题……我才不是从省选颓到现在才把它写了呢哼唧 嗯……其实这题真的就是标准的模板套模板啊┑(°Д °)┍ 倍增模板里套个线性基……都不难写【但是一般都用链剖写lca的我毫不犹豫写跪了倍增lca【助けて下さいQuQ 嗯……脑子是个好东西……感觉写了这道题智力--一定是我太辣鸡【望天 嗯对……调了好几天的原因啊 EXCITING -LCA写挂了智力-=INF -在线性基的
BZOJ4568 [Scoi2016]幸运数字
dgoh41514的博客
03-25 140
本文版权归ljh2000和博客园共有，欢迎转载，但须保留此声明，并给出原文链接，谢谢合作。 本文作者：ljh2000 作者博客：http://www.cnblogs.com/ljh2000-jump/转载请注明出处，侵权必究，保留最终解释权！ 题目链接：BZOJ4568 正解：线性基+链剖+线段树 解题报告： 　　考虑直接用线段树维护线性基，构出来之后...
【BZOJ4568】幸运数字（SCOI2016）-树上倍增+异或线性基合并
Maxwei_wzj的OI世界
03-22 335
测试地址：幸运数字 做法：本题需要用到树上倍增+异或线性基合并。 我们发现，每次询问都是求从一条路径上选一个子集，使得这个子集的异或和最大。如果只有一个询问，那敢开这题的大家肯定都会O(60n)O(60n)O(60n)的线性基+贪心做法，问题是有qqq个询问，暴力做显然会炸，所以需要找到优化的方法。 我们发现只要求出了对应向量空间的线性基，那么求最大异或和就是O(60)O(60)O(60)的...
BZOJ 4568 [Scoi2016]幸运数字
Poi的博客
05-20 731
题目大意：给你一颗树，多个询问，问你树上任意两点的路径上选任意几个点使得异或和最大。我是参考的claris大神的代码（%%%），点分治，对于询问在两个子树间或者有一个在重心上的进行回答，否则把问题用链表接到询问点所在的子树上。 具体方法可以选中重心都对每个子树染色，染为这个子树的根节点。在子树处理问题之前一定要记住先把询问拖出来，然后把子树的问题清零。因为你选的是这个子树的重心开始分治，而不是根，
bzoj 4568 [Scoi2016]幸运数字
Elijahqi
04-22 182
http://www.elijahqi.win/archives/3161 Description A 国共有 n 座城市，这些城市由 n-1 条道路相连，使得任意两座城市可以互达，且路径唯一。每座城市都有一个 幸运数字，以纪念碑的形式矗立在这座城市的正中心，作为城市的象征。一些旅行者希望游览 A 国。旅行者计划 乘飞机降落在 x 号城市，沿着 x 号城市到 y 号城市之间那条唯一的路径...
BZOJ4568 SCOI2016 幸运数字 倍增的思想维护线性基（线性基详解）
YxuanwKeith的博客
07-15 2075
题目大意给你一颗NN个节点的树，每个节点都有一个权值AiA_i，现在有MM组询问，每组询问有3个数u,v,u,v,，要求你输出在树上节点uu到节点vv的路径上，每个节点的权值可以选或不选，求选出的点的权值的异或最大值。N≤2∗104N \leq 2 * 10^4 M≤2∗105M \leq 2 * 10^5 Ai≤260A_i \leq 2^{60}解题思路看到这种区间异或最大值的题，那么我们就
【bzoj4568 scoi2016】幸运数字
weixin_33834075的博客
03-26 67
题目描述 A 国共有 n 座城市，这些城市由 n-1 条道路相连，使得任意两座城市可以互达，且路径唯一。每座城市都有一个幸运数字，以纪念碑的形式矗立在这座城市的正中心，作为城市的象征。 一些旅行者希望游览 A 国。旅行者计划乘飞机降落在 x 号城市，沿着 x 号城市到 y 号城市之间那条唯一的路径游览，最终从 y 城市起飞离开 A 国。在经过每一座城市时，游览者就会有机会与这座城市的幸运数字拍...
BZOJ 4568: [Scoi2016]幸运数字
limboman的博客
12-08 309
暴力合并线性基可以过？
BZOJ3230相似子串测试数据
08-19
BZOJ3230相似子串的测试数据，希望能够帮到大家。
bzoj problem
08-21
八中OJ所有题目
BZOJ上面的题目以及数据
10-13
BZOJ题目以及数据
resnet模型-基于图像分类算法对汉字写的是否工整识别-不含数据集图片-含逐行注释和说明文档.zip
05-23
resnet模型_基于图像分类算法对汉字写的是否工整识别-不含数据集图片-含逐行注释和说明文档 本代码是基于python pytorch环境安装的。 下载本代码后，有个环境安装的requirement.txt文本 如果有环境安装不会的，可自行网上搜索如何安装python和pytorch,这些环境安装都是有很多教程的，简单的 环境需要自行安装,推荐安装anaconda然后再里面推荐安装python3.7或3.8的版本，pytorch推荐安装1.7.1或1.8.1版本 首先是代码的整体介绍 总共是3个py文件，十分的简便 且代码里面的每一行都是含有中文注释的,小白也能看懂代码 然后是关于数据集的介绍。 本代码是不含数据集图片的，下载本代码后需要自行搜集图片放到对应的文件夹下即可 在数据集文件夹下是我们的各个类别,这个类别不是固定的，可自行创建文件夹增加分类数据集 需要我们往每个文件夹下搜集来图片放到对应文件夹下，每个对应的文件夹里面也有一张提示图，提示图片放的位置 然后我们需要将搜集来的图片，直接放到对应的文件夹下，就可以对代码进行训练了。 运行01生成txt.py,
计算机毕业设计选题精品毕设分享+源码+论文+PPT+asp.net0班级网站的设计与实现演示录像.rar
最新发布
05-23
博主给大家详细整理了计算机毕业设计最新项目,对项目有任何疑问(部署跟文档),都可以问博主哦~ 一、JavaWeb管理系统毕设项目【计算机毕设选题】计算机毕业设计选题，500个热门选题推荐,更多作品展示 计算机毕业设计|PHP毕业设计|JSP毕业程序设计|Android毕业设计|Python设计论文|微信小程序设计 更多作品展示 +微亻言 biyesheji02
Reborn 是使用 Go 开发的，基于 Redis 存储的配置库，简单配置，易于使用
05-23
Go语言（也称为Golang）是由Google开发的一种静态强类型、编译型的编程语言。它旨在成为一门简单、高效、安全和并发的编程语言，特别适用于构建高性能的服务器和分布式系统。以下是Go语言的一些主要特点和优势： 简洁性：Go语言的语法简单直观，易于学习和使用。它避免了复杂的语法特性，如继承、重载等，转而采用组合和接口来实现代码的复用和扩展。 高性能：Go语言具有出色的性能，可以媲美C和C++。它使用静态类型系统和编译型语言的优势，能够生成高效的机器码。 并发性：Go语言内置了对并发的支持，通过轻量级的goroutine和channel机制，可以轻松实现并发编程。这使得Go语言在构建高性能的服务器和分布式系统时具有天然的优势。 安全性：Go语言具有强大的类型系统和内存管理机制，能够减少运行时错误和内存泄漏等问题。它还支持编译时检查，可以在编译阶段就发现潜在的问题。 标准库：Go语言的标准库非常丰富，包含了大量的实用功能和工具，如网络编程、文件操作、加密解密等。这使得开发者可以更加专注于业务逻辑的实现，而无需花费太多时间在底层功能的实现上。 跨平台：Go语言支持多种操作系统和平台，包括Windows、Linux、macOS等。它使用统一的构建系统（如Go Modules），可以轻松地跨平台编译和运行代码。 开源和社区支持：Go语言是开源的，具有庞大的社区支持和丰富的资源。开发者可以通过社区获取帮助、分享经验和学习资料。 总之，Go语言是一种简单、高效、安全、并发的编程语言，特别适用于构建高性能的服务器和分布式系统。如果你正在寻找一种易于学习和使用的编程语言，并且需要处理大量的并发请求和数据，那么Go语言可能是一个不错的选择。
其他类别Jsp考试系统-jspks.7z
05-23
[其他类别]Jsp考试系统_jspks.7z
洛谷 P1640 BZOJ 1854 [SCOI2010]连续攻击游戏
05-14
这道题目需要进行博弈论的思考。我们可以使用动态规划来解决这道题目。 首先，我们可以定义 $f_{i,j}$ 表示当前剩下 $i$ 到 $j$ 这些数时，当前玩家最多能获得的分数。我们可以发现，在这个状态中，当前玩家只能从 ...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值