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. 半年收入超2亿RMB 独立游戏开发者的艰苦创业路

    一款叫做<监狱建筑师>的模拟经营游戏,目前在Steam平台获得了3000万美元(近2亿元)以上的收入.这款游戏由英国独立工作室Introversion Software发布,而团队最困难的 ...

  2. 【机器学习】无监督学习Autoencoder和VAE

    众所周知,机器学习的训练数据之所以非常昂贵,是因为需要大量人工标注数据. autoencoder可以输入数据和输出数据维度相同,这样测试数据匹配时和训练数据的输出端直接匹配,从而实现无监督训练的效果. ...

  3. PASSWORD MySQL 5.6.21-1ubuntu14.04_amd64

    /***************************************************************************** The main idea is that ...

  4. 1.Python3.6环境部署

    标题:Python3.6环境部署文档 作者:刘耀 内容 Linux部署Python3.6环境 Mac部署Python3.6环境 Window10部署Python3.6环境 Pycharm安装 1. L ...

  5. Node.js开发入门—套接字(socket)编程

    Node.js的net模块提供了socket编程接口,方便我们利用较为底层的套接字接口来实现应用协议.这次我们看一个简单的回显服务器示例,包括服务端和客户端的代码. 代码 分服务器和客户端两部分来说吧 ...

  6. 关于React面试题汇总

    1.redux中间件 中间件提供第三方插件的模式,自定义拦截 action -> reducer 的过程.变为 action -> middlewares -> reducer .这 ...

  7. java的第一个实验

    实验一 Java开发环境的熟悉 北京电子科技学院(BESTI) 实     验    报     告 课程:Java程序设计 班级:1352 姓名:林涵锦 学号:20135213 成绩:        ...

  8. Sprint8

    进展:添加事件主界面实现之后,实现事件添加部分代码的编写,进行设置事件提醒,选择时间.

  9. Arduino Leonardo读取DHT22温湿度传感器

    首先在该地址下载库:https://codeload.github.com/nethoncho/Arduino-DHT22/zip/master 使用以下代码测试: /**************** ...

  10. Java通用oracle和mysql数据库连接

    Java中oracle数据库连接写一个通用类UBUtil(){} import java.io.InputStream; import java.sql.*; import java.util.Pro ...