4568: [Scoi2016]幸运数字

Time Limit: 60 Sec  Memory Limit: 256 MB
Submit: 848  Solved: 336
[Submit][Status][Discuss]

Description

A 国共有 n 座城市，这些城市由 n-1 条道路相连，使得任意两座城市可以互达，且路径唯一。每座城市都有一个
幸运数字，以纪念碑的形式矗立在这座城市的正中心，作为城市的象征。一些旅行者希望游览 A 国。旅行者计划
乘飞机降落在 x 号城市，沿着 x 号城市到 y 号城市之间那条唯一的路径游览，最终从 y 城市起飞离开 A 国。
在经过每一座城市时，游览者就会有机会与这座城市的幸运数字拍照，从而将这份幸运保存到自己身上。然而，幸
运是不能简单叠加的，这一点游览者也十分清楚。他们迷信着幸运数字是以异或的方式保留在自己身上的。例如，
游览者拍了 3 张照片，幸运值分别是 5，7，11，那么最终保留在自己身上的幸运值就是 9（5 xor 7 xor 11）。
有些聪明的游览者发现，只要选择性地进行拍照，便能获得更大的幸运值。例如在上述三个幸运值中，只选择 5 
和 11 ，可以保留的幸运值为 14 。现在，一些游览者找到了聪明的你，希望你帮他们计算出在他们的行程安排中
可以保留的最大幸运值是多少。
 
lyd说这道题维护倍增数组后暴力合并就可以。
于是。。。。
 
我不是有意要卡评测啊啊啊啊

然后发现果然合并有更NB的方法。

不是高斯消元(原谅我zz)，而是应该一个一个插。

终于19sA掉了。

52s

 #include<iostream>

 #include<cstdio>

 #include<algorithm>

 #include<cstring>

 #define ll long long

 #define N 20005

 using namespace std;

 int read()

 {

     int p=;char c=getchar();

     while(c<''||c>'')c=getchar();

     while(c>=''&&c<='')p=p*+c-'',c=getchar();

     return p;

 }

 int n,m;

 int head[N],ver[N*],nxt[N*],tot;

 ll mi[N];

 void add(int a,int b)

 {

     tot++;nxt[tot]=head[a];head[a]=tot;ver[tot]=b;return ;

 }

 ll ji[N][][];

 int fa[N][],dep[N];

 ll a[N];

 void dfs(int x,int f)

 {

     for(int i=head[x];i;i=nxt[i])

     {

         if(ver[i]==f)continue;

         fa[ver[i]][]=x;

         dep[ver[i]]=dep[x]+;

         ji[ver[i]][][]=;

         ji[ver[i]][][]=a[ver[i]];

         dfs(ver[i],x);

     }

     return ;

 }

 ll tmp[];

 void guess()

 {

     int k=;int now=tmp[];

     for(int i=;i>=;i--)

     {

         int p=;ll o=mi[i];

         for(int j=k;j<=now;j++)

         {

             if((tmp[j]&o))

             {

                 p=j;break;

             }

         }

         if(p)

         {

             swap(tmp[p],tmp[k]);

             for(int j=;j<=now;j++)

             {

                 if(j!=k&&(tmp[j]&o))

                 {

                     tmp[j]^=tmp[k];

                 }

             }k++;

         }

     }

     tmp[]=k-;

 }

 void yu()

 {

     for(int i=;i<=;i++)

     {

         for(int j=;j<=n;j++)

         {

             tmp[]=ji[j][i-][]+ji[fa[j][i-]][i-][];

             for(int k=;k<=ji[j][i-][];k++)tmp[k]=ji[j][i-][k];

             for(int k=;k<=ji[fa[j][i-]][i-][];k++)tmp[k+ji[j][i-][]]=ji[fa[j][i-]][i-][k];

             guess();

             for(int k=;k<=tmp[];k++)ji[j][i][k]=tmp[k];

             fa[j][i]=fa[fa[j][i-]][i-];

         }

     }return ;

 }

 void solve(int x,int y)

 {

     if(dep[x]<dep[y])swap(x,y);

     for(int i=;i>=;i--)

     {

         if(dep[fa[x][i]]>=dep[y])

         {

             for(int j=;j<=ji[x][i][];j++)tmp[]++,tmp[tmp[]]=ji[x][i][j];

             guess();

             x=fa[x][i];

         }

     }

     if(x!=y)

     {

         for(int i=;i>=;i--)

         {

             if(fa[x][i]!=fa[y][i])

             {

                 for(int k=;k<=ji[x][i][];k++)tmp[]++,tmp[tmp[]]=ji[x][i][k];

                 guess();

                 for(int k=;k<=ji[y][i][];k++)tmp[]++,tmp[tmp[]]=ji[y][i][k];

                 guess();

                 x=fa[x][i];y=fa[y][i];

             }

         }

         tmp[]++;tmp[tmp[]]=a[fa[x][]];

         tmp[]++;tmp[tmp[]]=a[x];

         tmp[]++;tmp[tmp[]]=a[y];

     }

     else

     {

         tmp[]++,tmp[tmp[]]=a[x];

     }

     guess();

     ll now=;

     for(int i=tmp[];i>=;i--)

     {

         now^=tmp[i];

     }

     printf("%lld\n",now);

     return ;

 }

 int main()

 {

     int q;

     mi[]=;

     for(int i=;i<=;i++)mi[i]=mi[i-]*;

     scanf("%d%d",&n,&q);

     for(int i=;i<=n;i++)

     {

         scanf("%lld",&a[i]);

     }

     int t1,t2,t3;

     for(int i=;i<n;i++)

     {

         t1=read();t2=read();

         add(t1,t2);add(t2,t1);

     }

     dep[]=;

     ji[][][]=;ji[][][]=a[];

     dfs(,-);

     yu();

     for(int i=;i<=q;i++)

     {

         t1=read();t2=read();

         memset(tmp,,sizeof(tmp));

         solve(t1,t2);

     }

     return ;

 }

19s

 #include<iostream>

 #include<cstdio>

 #include<algorithm>

 #include<cstring>

 #define ll long long

 #define N 20005

 using namespace std;

 int read()

 {

     int p=;char c=getchar();

     while(c<''||c>'')c=getchar();

     while(c>=''&&c<='')p=p*+c-'',c=getchar();

     return p;

 }

 int n,m;

 int head[N],ver[N*],nxt[N*],tot;

 ll mi[N];

 void add(int a,int b)

 {

     tot++;nxt[tot]=head[a];head[a]=tot;ver[tot]=b;return ;

 }

 ll ji[N][][];

 int fa[N][],dep[N];

 ll a[N];

 void dfs(int x,int f)

 {

     for(int i=head[x];i;i=nxt[i])

     {

         if(ver[i]==f)continue;

         fa[ver[i]][]=x;

         dep[ver[i]]=dep[x]+;

         for(int j=;j>=;j--)

         {

             if(a[ver[i]]&mi[j])

             {

                 ji[ver[i]][][j]=a[ver[i]];

                 break;

             }

         }

         dfs(ver[i],x);

     }

     return ;

 }

 ll tmp[];

 void insert(ll x)

 {

     for(int i=;i>=;i--)

     {

         if(x&mi[i])

         {

             if(!tmp[i])

             {

                 tmp[i]=x;

                 break;

             }

             else x^=tmp[i];

         }

     }return ;

 }

 void yu()

 {

     for(int i=;i<=;i++)

     {

         for(int j=;j<=n;j++)

         {

             for(int k=;k<=;k++)tmp[k]=ji[j][i-][k];

             for(int k=;k<=;k++)

             {

                 if(ji[fa[j][i-]][i-][k])insert(ji[fa[j][i-]][i-][k]);

             }

             for(int k=;k<=;k++)ji[j][i][k]=tmp[k];

             fa[j][i]=fa[fa[j][i-]][i-];

         }

     }return ;

 }

 void solve(int x,int y)

 {

     if(dep[x]<dep[y])swap(x,y);

     for(int i=;i>=;i--)

     {

         if(dep[fa[x][i]]>=dep[y])

         {

             for(int j=;j<=;j++)if(ji[x][i][j])insert(ji[x][i][j]);

             x=fa[x][i];

         }

     }

     if(x!=y)

     {

         for(int i=;i>=;i--)

         {

             if(fa[x][i]!=fa[y][i])

             {

                 for(int k=;k<=;k++)if(ji[x][i][k])insert(ji[x][i][k]);

                 for(int k=;k<=;k++)if(ji[y][i][k])insert(ji[y][i][k]);

                 x=fa[x][i];y=fa[y][i];

             }

         }

         insert(a[fa[x][]]);insert(a[x]);insert(a[y]);

     }

     else

     {

         insert(a[x]);

     }

     ll now=;

     for(int i=;i>=;i--)

     {

         now=max(now,now^(tmp[i]));

     }

     printf("%lld\n",now);

     return ;

 }

 int main()

 {

     int q;

     mi[]=;

     for(int i=;i<=;i++)mi[i]=mi[i-]*;

     scanf("%d%d",&n,&q);

     for(int i=;i<=n;i++)

     {

         scanf("%lld",&a[i]);

     }

     int t1,t2,t3;

     for(int i=;i<n;i++)

     {

         t1=read();t2=read();

         add(t1,t2);add(t2,t1);

     }

     dep[]=;

     for(int i=;i>=;i--)

     {

         if(a[]&mi[i])

         {

             ji[][][i]=a[];

             break;

         }

     }

     dfs(,-);

     yu();

     for(int i=;i<=q;i++)

     {

         t1=read();t2=read();

         memset(tmp,,sizeof(tmp));

         solve(t1,t2);

     }

     return ;

 }

bzoj 4568: [Scoi2016]幸运数字的更多相关文章

  1. BZOJ 4568: [Scoi2016]幸运数字 [线性基 倍增]

    4568: [Scoi2016]幸运数字 题意:一颗带点权的树,求树上两点间异或值最大子集的异或值 显然要用线性基 可以用倍增的思想,维护每个点向上\(2^j\)个祖先这些点的线性基,求lca的时候合 ...

  2. BZOJ 4568 [Scoi2016]幸运数字（树链剖分 + 异或线性基）

    题目链接  BZOJ 4568 考虑树链剖分+线段树维护每一段区域的异或线性基 对于每个询问,求出该点集的异或线性基.然后求一下这个线性基里面能异或出的最大值即可. #include <bits ...

  3. bzoj 4568: [Scoi2016]幸运数字【树链剖分+线段树+线性基】

    一眼做法,好处是好想好写坏处是常数大,容易被卡(bzoj loj 洛谷开O2 能AC,不开有90分-- 大概就是树剖之后维护线段树,在线段树的每个节点上上维护一个线性基,暴力\( 60^2 \)的合并 ...

  4. BZOJ 4568 [Scoi2016]幸运数字 ——线性基 倍增

    [题目分析] 考虑异或的最大值,维护线性基就可以了. 但是有多次的询问,树剖或者倍增都可以. 想了想树剖动辄数百行的代码. 算了,我还是写倍增吧. 注:被位运算和大于号的优先级坑了一次,QaQ [代码 ...

  5. BZOJ 4568: [Scoi2016]幸运数字(倍增+线性基)

    传送门 解题思路 异或最大值肯定线性基了,树上两点那么就倍增搞一搞,就维护每个点到各级祖先的线性基,时间复杂度\(O(nlog^3n)\),并不知道咋过去的. 代码 #include<iostr ...

  6. 【BZOJ 4568】 4568: [Scoi2016]幸运数字 （线性基+树链剖分+线段树）

    4568: [Scoi2016]幸运数字 Description A 国共有 n 座城市,这些城市由 n-1 条道路相连,使得任意两座城市可以互达,且路径唯一.每座城市都有一个 幸运数字,以纪念碑的形 ...

  7. BZOJ.4516.[SCOI2016]幸运数字(线性基 点分治)

    题目链接 线性基可以\(O(log^2)\)暴力合并.又是树上路径问题,考虑点分治. 对于每个点i求解 LCA(u,v)==i 时的询问(u,v),只需求出这个点到其它点的线性基后,暴力合并. LCA ...

  8. [BZOJ4568][Scoi2016]幸运数字 倍增+线性基

    4568: [Scoi2016]幸运数字 Time Limit: 60 Sec  Memory Limit: 256 MBSubmit: 1791  Solved: 685[Submit][Statu ...

  9. [BZOJ4568][SCOI2016]幸运数字(倍增LCA,点分治+线性基)

    4568: [Scoi2016]幸运数字 Time Limit: 60 Sec  Memory Limit: 256 MBSubmit: 2131  Solved: 865[Submit][Statu ...

随机推荐

  1. 【视频编解码·学习笔记】13. 提取PPS信息程序

    PPS结构解析 与之前解析SPS方式类似 一.定义PPS类: 在3.NAL Unit目录下,新建PicParamSet.cpp和PicParamSet.h,在这两个文件中写入类的定义和函数实现. 类定 ...

  2. MariaDB远程连接问题

    MariaDB在设置完通过Navicat Premium远程连接账号验证通过,但是无法正常使用工具的功能,只能使用sql语句查询,但是通过控制台命令功能正常. 经过修改账号权限,添加新用户等功能都无法 ...

  3. dirname命令详解

    基础命令学习目录首页 原文链接:https://blog.csdn.net/xiaofei125145/article/details/50620281 示例一 来自手册页的例子 $ dirname ...

  4. [BUAA OO]第三次博客作业

    OO第三次博客作业 1. 规格化设计的发展 我认为,规格化设计主要源自于软件设计的两次危机.第一次是由于大量存在的goto语句,让当时被广泛应用的面向过程式的编程语言臃肿不堪,在逻辑性上与工程规模上鱼 ...

  5. CS小分队第二阶段冲刺站立会议（6月4日）

    昨日成果:昨天一直在对主界面进行修改,遇到问题没有进展 遇到的问题:我代码写的不够缜密,各按钮信息添加的删除的时候总是有重名或者覆盖现象,需要有一次大的检查 今日计划:冲刺已经结束,项目的难度超过了预 ...

  6. iOS学习资源搜集

    swift 2.0 新的开始 iOS7初学者入门 斯坦福大学公开课:iOS 8开发 苹果官方开发 中文 iOS/Mac 开发博客列表    git

  7. 冲刺One之站立会议8 /2015-5-21

    今天我们把聊天界面做了优化和改进,主要实现了聊天的功能.显示了正在进行通信的成员列表,和当前状态,是否连通和正常通信,大体完成了预期的目标. 燃尽图8

  8. 初学Cocos2dx

    初学cocos2dx Cocos2dx 中的主要概念包括:应用.导演.场景.层.精灵.动画.动作. Cocos2dx里面的主要类 1.CCObject Object Object Object 是co ...

  9. 45度炸队Alpha冲刺博客集

    博客集如下: Alpha冲刺Day1:第一天冲刺记录 Alpha冲刺Day2:第二天冲刺记录 Alpha冲刺Day3:第三天冲刺记录 Alpha冲刺Day4:第四天冲刺记录 Alpha冲刺Day5:第 ...

  10. 团队作业7——第二次项目冲刺（Beta版本12.04——12.07）

    1.当天站立式会议照片 本次会议在5号公寓3楼召开,本次会议内容:①:熟悉每个人想做的模块.②:根据项目要求还没做的完成. 2.每个人的工作 经过会议讨论后确定了每个人的分工 组员 任务 陈福鹏 实现 ...