这题一场模拟赛我们出了弱化版(n<=1e6),抄题面给的程序能拿到71分的好成绩
其实后面的29分是加了几个1e9的数据卡人
这糟老头子真是坏得很
正解我们机房看了三天

在这里感谢这篇题解的作者,代码解释得很清晰~

经过打表观察,可以发现:当\(1\le x \le k\)时

如果 \(k\) 为奇数,\(x*2^{y}\) \(mod\) \(k\)的值成环

如果 \(k\) 为偶数,质因数分解\(x\),如果所含因子\(2\)的次数大于\(k\)所含因子\(2\)的次数,那\(x*2^{y}\) \(mod\) \(k\)的值也成环

然而考试时永远发现不了

于是可以对每条链开一棵动态开点线段树,再开一棵线段树处理不在单点上的点,一起用来模拟题目所给的树状数组(\(N\le10^{9}\),直接开树状树组会炸)。

于是可以这么做修改操作:

对于不在链上的点,最多只会跳\(log_{2}^{k}*log_{k}^{n}=log_{2}^{n}\)次就会跳到\(n\)之外或跳到链上,暴力更新就完事了。

对于在链上的点,我们可以找到这条链的链头(找链头的方法详见代码),求出这个点到链头的距离\(dis\),区间更新\(dis\)和链尾之间的点。

查询也差不多的。

更详细的解释见代码:

#include <bits/stdc++.h>
#define N 200010
using namespace std;

int k;
inline void rd(int &x){
    int y=0,f=1;char c=getchar();
    while(c<'0' || c>'9' ){if(c=='-')f=-1;c=getchar();}
    while(c>='0' && c<='9')y=y*10+(c^48),c=getchar();
    x=y*f;
}

int num[N],jw[N],go[N][20],to[N][20],pr[N],a[N],fr,dis,u;
bool vis[N];

inline int lowbit(int x){
    while(x%k==0) x/=k;
    return x%k;
}
inline int lowbitv(int x){
    register int lg=0;
    while (x%k==0) x/=k,lg++;
    x%=k;
    while (lg) x*=k,lg--;
    return x;
}
void init()//预处理
{
    //jw[i]:i所在链上的进位次数
    register int cnt,i,j,w,l;
    for(i=1;i<=k;++i) for(j=i;j%2==0;j/=2,++pr[i]);
    for(i=1;i<k;++i){
        cnt=w=0;
        if(vis[i] || pr[i]<pr[k]) continue;
        int x=i*2%k;a[++cnt]=i,vis[i]=1;
        while(x!=i) a[++cnt]=x,vis[x]=1,x=x*2%k;
        for(j=1;j<=cnt;++j) w+=(a[j]*2>=k);
        for(j=1;j<=cnt;++j){
            jw[a[j]]=w,num[a[j]]=cnt;
            (j==1)?x=cnt:x=j-1;
            to[a[j]][0]=a[x];
            go[a[j]][0]=(a[x]*2>=k);
        }
        //to[i][j]:i向前跳j步跳到的点
        //go[i][j]:i向前跳j步的进位次数
        for(j=1;j<=18;++j)
            for(l=1;l<=cnt;++l){
                to[a[l]][j]=to[to[a[l]][j-1] ][j-1];
                go[a[l]][j]=go[a[l]][j-1]+go[to[a[l]][j-1] ][j-1];
            }
    }
}
void GetF(int x)//求链头
{
    int lg=0;
    while(x%k==0) x/=k,lg++;
    fr=lowbit(x),dis=1;int u=(x-fr)/k,i;//u:x的进位次数
    //dis=跳的环数*环的长度+剩余长度(没跳完整个环)
    dis+=u/jw[fr]*num[fr];
    u%=jw[fr];
    for(i=18;i>=0;--i){
        if(go[fr][i]<=u){
            u-=go[fr][i];
            dis+=(1<<i);
            fr=to[fr][i];
        }
    }
    while(lg) fr*=k,lg--;//由lowbit(x)的链头还原回去x的链头
}

#define mid ((l+r)>>1)
int root1,root2,root3;
struct tr{ int ls,rs,x; };
struct qwq{
    tr t[N*60];
    int sz=0;
    void update(int &o,int l,int r,int L,int R,int v){
        if(!o) o=++sz;
        if(L<=l && r<=R){
            t[o].x^=v;
            return;
        }
        if(L<=mid) update(t[o].ls,l,mid,L,R,v);
        if(R>mid) update(t[o].rs,mid+1,r,L,R,v);
    }
    int query(int o,int l,int r,int x){
        if(l==r || !o) return t[o].x;
        if(x<=mid) return t[o].x^query(t[o].ls,l,mid,x);//不用pushup的原因(异或只用异或一遍)
        else return t[o].x^query(t[o].rs,mid+1,r,x);
    }
} Point,Chain,Front;
//Point:不在链上的点
//Chain:链上的点的更新(以链头的点值为根)
//Front:链头在Chain上的位置
int main()
{
    int n,q,op,x,v,ans;
    rd(n),rd(q),rd(k);
    init();
    while(q--){
        rd(op),rd(x);
        if(op==1){
            rd(v);
            while(x<=n && pr[lowbit(x)]<pr[k]){
                Point.update(root2,1,n,x,x,v);
                x+=lowbitv(x);
            } //暴力跳
            if(x>n) continue;
            GetF(x);
            root1=Front.query(root3,1,n,fr);
            int tmp=root1;Chain.update(root1,1,n,dis,n,v);
            //Chain相当于重构了编号,比如说链头为x,则线段树里编号为2的元素实际上是原树状数组编号为x+lowbit(x)的元素
            if(!tmp) Front.update(root3,1,n,fr,fr,root1);
        } else{
            ans=0;
            while(x){
                if(pr[lowbit(x)]<pr[k]) ans^=Point.query(root2,1,n,x);
                else{
                    GetF(x);
                    root1=Front.query(root3,1,n,fr);
                    if(root1) ans^=Chain.query(root1,1,n,dis);
                }
                x-=lowbitv(x);
            }
            printf("%d\n",ans);
        }
    }
}

LOJ#510 北校门外的回忆(找性质+倍增+线段树)的更多相关文章

  1. LOJ 510: 「LibreOJ NOI Round #1」北校门外的回忆

    题目传送门:LOJ #510. 题意简述: 给出一个在 \(K\) 进制下的树状数组,但是它的实现有问题. 形式化地说,令 \(\mathrm{lowbit}(x)\) 为在 \(K\) 进制下的 \ ...

  2. LOJ #2359. 「NOIP2016」天天爱跑步(倍增+线段树合并)

    题意 LOJ #2359. 「NOIP2016」天天爱跑步 题解 考虑把一个玩家的路径 \((x, y)\) 拆成两条,一条是 \(x\) 到 \(lca\) ( \(x, y\) 最近公共祖先) 的 ...

  3. LOJ#510. 「LibreOJ NOI Round #1」北校门外的回忆(线段树)

    题面 传送门 题解 感谢\(@M\_sea\)的代码我总算看懂题解了-- 这个操作的本质就是每次把\(x\)的\(k\)进制最低位乘\(2\)并进位,根据基本同余芝士如果\(k\)是奇数那么最低位永远 ...

  4. LOJ 北校门外的回忆 倍增+线段树

    正解:倍增+线段树 解题报告: 传送门! $umm$这题有个对正解毫无启发的部分分还有个正解,都挺神仙的所以我都写了趴$QAQ$ 先说部分分 可以考虑把$x$向$x+lowbit(x)$连边,然后当$ ...

  5. 51nod 1463 找朋友 (扫描线+线段树)

    1463 找朋友  基准时间限制:1.5 秒 空间限制:262144 KB 分值: 80 难度:5级算法题  收藏  关注 给定: 两个长度为n的数列A .B 一个有m个元素的集合K 询问Q次 每次询 ...

  6. 【NOI P模拟赛】校门外歪脖树上的鸽子(树链剖分)

    题面 2 ≤ n ≤ 2 × 1 0 5 , 1 ≤ m ≤ 2 × 1 0 5 , 1 ≤ l ≤ r ≤ n , 1 ≤ d ≤ 1 0 8 2 ≤ n ≤ 2 × 10^5,1 ≤ m ≤ 2 ...

  7. loj#6029. 「雅礼集训 2017 Day1」市场(线段树)

    题意 链接 Sol 势能分析. 除法是不能打标记的,所以只能暴力递归.这里我们加一个剪枝:如果区间内最大最小值的改变量都相同的话,就变成区间减. 这样复杂度是\((n + mlogn) logV\)的 ...

  8. LOJ.6062.[2017山东一轮集训]Pair(Hall定理 线段树)

    题目链接 首先Bi之间的大小关系没用,先对它排序,假设从小到大排 那么每个Ai所能匹配的Bi就是一个B[]的后缀 把一个B[]后缀的匹配看做一条边的覆盖,设Xi为Bi被覆盖的次数 容易想到 对于每个i ...

  9. 51nod 1463 找朋友(线段树+离线处理)

    http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1463 题意: 思路: 好题! 先对所有查询进行离线处理,按照右区间排序, ...

随机推荐

  1. FAIL - Deploy Upload Failed, Exception: [org.apache.tomcat.util.http.fileupload.FileUploadBase$SizeLimitExceededException: the request was rejected because its size (112503036) exceeds the configured

    Message:  FAIL - Deploy Upload Failed, Exception: [org.apache.tomcat.util.http.fileupload.FileUpload ...

  2. mysql之整型数据int

    mysql数据库设计,其中,对于数据性能优化,字段类型考虑很重要,mysql整型bigint.int.mediumint.smallint 和 tinyint的语法介绍,如下:1.bigint 从 - ...

  3. Javascript与C#对变量的处理方式

      先来看一下Javascript的情况(下面所说的基本类型和简单类型是一个意思): Javascript中变量会存在两种情况,一种是基本类型的,一共有五种,有null.Bollean.undefin ...

  4. 转 JQuery:常用方法一览

    出处 :http://www.cnblogs.com/Fooo/archive/2010/02/01/1661157.html 代码 Attribute:$(”p”).addClass(css中定义的 ...

  5. ASP.NET Core Building chat room using WebSocket

    Creating “Login form” We use here simple form where user can insert his or her preferred nick name f ...

  6. Modification Default Identity Table

    Step 1. Open IdentityModel.cs to add following code. protected override void OnModelCreating(DbModel ...

  7. HTML5 ----- deviceorientation API

    当我们把设备举到面前,就是坐公交车时玩手机的状态,手机的宽度(水平方向)就是X轴,从左向右依次变大:手机的高度(垂直方向)就是y轴,从下到上依次变大: 眼光盯着屏幕,我们的视线就是Z轴,离我们的眼睛越 ...

  8. 训练赛-Building Numbers

    题意:首先告诉你,一个数字从1开始有两种变换方式:1.当前数字的值加1 2.当前的数字值乘2: 思路:首先把数组里的数字需要的变换次数算出来,然后用前缀和解决: 代码: #include<ios ...

  9. JarvisOJ Basic veryeasyRSA

    已知RSA公钥生成参数: p = 3487583947589437589237958723892346254777 q = 8767867843568934765983476584376578389 ...

  10. Django中Form组件的使用

    Form介绍 HTML页面中利用form表单向后端提交数据时,都会写一些获取用户输入的标签并且用form标签把它们包起来. 与此同时我们在好多场景下都需要对用户的输入做校验,比如校验用户是否输入,输入 ...