P5072 [Ynoi2015] 盼君勿忘

于 2024-04-11 22:10:47 发布
阅读量236 收藏 4
点赞数 4
分类专栏： ynoi_ds题 文章标签： 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Lanthamum/article/details/137653464
版权
文章讨论了解决一个关于数对贡献的问题，利用双向链表存储出现次数，unordered_set存储特定次数的数，采用光速幂优化查询时间复杂度，实现区间内特定数出现次数的高效计算。
摘要由CSDN通过智能技术生成

如果有来世希望还能与你相遇

分析一下询问

考虑计算每个数对答案所作的贡献

假设存在一个x，它在[l,r]区间出现了k次，我们想计算它的贡献

正向考虑似乎很困难

反向考虑一下，我们知道子字串的数量，没有x的子字串的数量

因此贡献应该是(2^{r-l+1}-2^{r-l+1-k})*x

因此我们只需要维护每个数的贡献即可->普通莫队

出现的次数可以用c[N]计数，出现该次数的数可以用sum[N]记录，有几种k可以用双向链表记录或者set,（选择unordered_set常数比较小）

感觉完事了，pow(2,i)快速幂不就完了吗

分析一下

k最多有sqrt(n)个 如果用快速幂预处理的话时间复杂度是O(sqrt(n)*log(sn(sqrt(2e5))))

sn是等差求和 k=1,k=2,k=3...理论上来说能过？

但是这是ynoi的题

考虑用光速幂 预处理O(sqrt(n)),查询O(1)

a^{b}=a^{\left \lfloor b/sqrt(n) \right\rfloor*sqrt(n)}*a^{bmodsqrt(n)}

因此我们只需要处理a^{sqrt(n)*i}以及a^{i}即可

inline void init(int mod){//处理sqrt(n)
    t1[0]=1;
    for(int i=1;i<=mi;i++){
        t1[i]=t1[i-1]*2%mod;
    }
    t2[0]=1;
    for(int i=1;i<=mi;i++){
        t2[i]=t2[i-1]*t1[mi]%mod;
    }
}
inline ll cqmi(int x,int mod){//O(1)询问
    return t1[x%mi]*t2[x/mi]%mod;
}

完整代码 unordered_set

#include<iostream>
#include<algorithm>
#include<cmath>
#include<unordered_set>
#include<vector>
#define INF (1ll<<60)
using namespace std;
typedef long long ll;
namespace Lan {
    inline string sread() {
        string s=" ";char e=getchar();
        while(e==' '||e=='\n')e=getchar();
        while(e!=' '&&e!='\n')s+=e,e=getchar();
        return s;
    }
    inline void swrite(string s){
        for(char e:s)putchar(e);
        printf("\n");
    }
    inline ll read() {
        ll x=0,y=1;char c=getchar();
        while(!isdigit(c)){if(c=='-')y=-1;c=getchar();}
        while(isdigit(c)){x=(x<<3)+(x<<1)+(c^48);c=getchar();}
        return x*=y;
    }
    inline void write(ll x) {
        if(x<0){x=-x,putchar('-');}ll sta[35],top=0;
        do sta[top++]=x%10,x/=10;while(x);
        while(top)putchar(sta[--top]+'0');
    }
}using namespace Lan;
const int N=1e5+9;
const int B=1e3+9;
int a[N],c[N],sum[N];//数，数量，出现k次的和
int L[B],R[B],pos[N];
ll t1[N],t2[N];
ll ans[N];
int cur,mi;
unordered_set<int> st;
struct Q{
    int l,r,p,id;
    friend bool operator < (const Q &a,const Q &b){
        return pos[a.l]^pos[b.l]?pos[a.l]<pos[b.l]:pos[a.l]&1?a.r<b.r:a.r>b.r;
    }
}que[N];
inline void init(int mod){//处理sqrt(n)
    t1[0]=1;
    for(int i=1;i<=mi;i++){
        t1[i]=t1[i-1]*2%mod;
    }
    t2[0]=1;
    for(int i=1;i<=mi;i++){
        t2[i]=t2[i-1]*t1[mi]%mod;
    }
}
inline ll cqmi(int x,int mod){//O(1)询问
    return t1[x%mi]*t2[x/mi]%mod;
}
inline ll query(int l,int r,int p){
    ll res=0;
    for(auto &i : st){
        res+=sum[i]*(cqmi(r-l+1,p)-cqmi(r-l+1-i,p)+p)%p;
    }
    res%=p;
    return res;
}
inline void add(int pos){
    sum[c[a[pos]]]-=a[pos];
    if(!sum[c[a[pos]]]){
        st.erase(c[a[pos]]);
    }
    c[a[pos]]++;
    if(!sum[c[a[pos]]]){
        st.insert(c[a[pos]]);
    }
    sum[c[a[pos]]]+=a[pos];
}
inline void del(int pos){
    sum[c[a[pos]]]-=a[pos];
    if(!sum[c[a[pos]]]){
        st.erase(c[a[pos]]);
    }
    c[a[pos]]--;
    if(!sum[c[a[pos]]]){
        st.insert(c[a[pos]]);
    }
    sum[c[a[pos]]]+=a[pos]; 
}
int main(){
    ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);
    int n,m;
    n=read();
    m=read();
    for(int i=1;i<=n;i++){
        a[i]=read();
    }
    for(int i=1;i<=m;i++){
        que[i].l=read(),que[i].r=read(),que[i].p=read();
        que[i].id=i;
    }
    mi=sqrt(n)+1;
    int blo=n/sqrt(m);
    int t=ceil(1.0*n/blo);
    for(int i=1;i<=t;i++){
        L[i]=(i-1)*blo+1;
        R[i]=i*blo;
    }
    R[t]=n;
    for(int i=1;i<=t;i++){
        for(int j=L[i];j<=R[i];j++){
            pos[j]=i;
        }
    }
    sort(que+1,que+1+m);
    int l=1,r=0;
    for(int i=1;i<=m;i++){
        while(que[i].l<l){
            add(--l);
        }
        while(que[i].l>l){
            del(l++);
        }
        while(que[i].r>r){
            add(++r);
        }
        while(que[i].r<r){
            del(r--);
        }
        init(que[i].p);
        ans[que[i].id]=query(que[i].l,que[i].r,que[i].p);
    }
    for(int i=1;i<=m;i++){
        cout<<ans[i]<<'\n';
    }
    return 0;
}

双向链表

#include<iostream>
#include<algorithm>
#include<cmath>
#include<vector>
#define INF (1ll<<60)
using namespace std;
typedef long long ll;
namespace Lan {
    inline string sread() {
        string s=" ";char e=getchar();
        while(e==' '||e=='\n')e=getchar();
        while(e!=' '&&e!='\n')s+=e,e=getchar();
        return s;
    }
    inline void swrite(string s){
        for(char e:s)putchar(e);
        printf("\n");
    }
    inline ll read() {
        ll x=0,y=1;char c=getchar();
        while(!isdigit(c)){if(c=='-')y=-1;c=getchar();}
        while(isdigit(c)){x=(x<<3)+(x<<1)+(c^48);c=getchar();}
        return x*=y;
    }
    inline void write(ll x) {
        if(x<0){x=-x,putchar('-');}ll sta[35],top=0;
        do sta[top++]=x%10,x/=10;while(x);
        while(top)putchar(sta[--top]+'0');
    }
}using namespace Lan;
const int N=1e6+9;
const int B=1e3+9;
int a[N],c[N],sum[N];//数，数量，出现k次的和
int L[B],R[B],pos[N];
int pre[N],nxt[N];//双向链表
ll t1[N],t2[N];
ll ans[N];
int t,cur,mi;
struct Q{
    int l,r,p,id;
    friend bool operator < (const Q &a,const Q &b){
        return pos[a.l]^pos[b.l]?pos[a.l]<pos[b.l]:pos[a.l]&1?a.r<b.r:a.r>b.r;
    }
}que[N];
inline void init(int mod){//处理sqrt(n)
    t1[0]=1;
    for(int i=1;i<=mi;i++){
        t1[i]=t1[i-1]*2%mod;
    }
    t2[0]=1;
    for(int i=1;i<=mi;i++){
    	t2[i]=t2[i-1]*t1[mi]%mod;
    }
}
inline ll cqmi(int x,int mod){//O(1)询问
    return (t1[x%mi]*t2[x/mi])%mod;
}
inline void insert(int k){
    nxt[cur]=k;//cur的后面是k
    pre[k]=cur;//k的前面是cur
    cur=k;//现在最后一个是k
}
inline void erase(int k){
    if(cur==k){//在最后一个直接删
        nxt[pre[k]]=0;
        cur=pre[k];	
    }else{
        nxt[pre[k]]=nxt[k];//中间断开
        pre[nxt[k]]=pre[k];
    }
    pre[k]=0;
    nxt[k]=0;
}
inline ll query(int l,int r,int p){
    ll res=0;
    for(int i=nxt[0];i;i=nxt[i]){
        res+=sum[i]*(cqmi(r-l+1,p)-cqmi(r-l+1-i,p)+p)%p;
    }
    res%=p;
    return res;
}
inline void add(int pos){
    sum[c[a[pos]]]-=a[pos];
    if(!sum[c[a[pos]]]){
        erase(c[a[pos]]);
    }
    c[a[pos]]++;
    if(!sum[c[a[pos]]]){
        insert(c[a[pos]]);
    }
    sum[c[a[pos]]]+=a[pos];
}
inline void del(int pos){
    sum[c[a[pos]]]-=a[pos];
    if(!sum[c[a[pos]]]){
        erase(c[a[pos]]);
    }
    c[a[pos]]--;
    if(!sum[c[a[pos]]]){
        insert(c[a[pos]]);
    }
    sum[c[a[pos]]]+=a[pos]; 
}
int main(){
    ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);
    int n,m;
    n=read();
    m=read();
    for(int i=1;i<=n;i++){
        a[i]=read();
    }
    for(int i=1;i<=m;i++){
        que[i].l=read(),que[i].r=read(),que[i].p=read();
        que[i].id=i;
    }
    mi=sqrt(n)+1;
    int blo=n/sqrt(m);
    t=ceil(1.0*n/blo);
    for(int i=1;i<=t;i++){
        L[i]=(i-1)*blo+1;
        R[i]=i*blo;
    }
    R[t]=n;
    for(int i=1;i<=t;i++){
        for(int j=L[i];j<=R[i];j++){
            pos[j]=i;
        }
    }
    sort(que+1,que+1+m);
    int l=1,r=0;
    for(int i=1;i<=m;i++){
        while(que[i].l<l){
            add(--l);
        }
        while(que[i].l>l){
            del(l++);
        }
        while(que[i].r>r){
            add(++r);
        }
        while(que[i].r<r){
            del(r--);
        }
        init(que[i].p);
        ans[que[i].id]=query(que[i].l,que[i].r,que[i].p);
    }
    for(int i=1;i<=m;i++){
        cout<<ans[i]<<'\n';
    }
    return 0;
}

学了一下双向链表,unordered_set, 光速幂感觉比做思维题有收获

我是数据结构魔怔人

Lanthanmum
关注
专栏目录
[Ynoi2018]未来日记
StaroForgin的博客
08-15 260
所以说卡半天块长不如加一个没啥用的特判。再update上一个「望月悲叹的最初分块」。
P7126 [Ynoi2008] rdCcot
Z_Y_S_的博客
03-03 375
点分治套路
[Ynoi2015]盼君勿忘
12-11 310
题目大意： 给定一个序列，每次查询一个区间\([l,r]\)中所有子序列分别去重后的和\(\bmod p\)（每次询问模数不同）。 解题思路： 在太阳西斜的这个世界里，置身天上之森。等这场战争结束之后，不归之人与望眼欲穿的众人， 人人本着正义之名，长存不灭的过去、逐渐消逝的未来。我回来了，纵使日薄西山，即便看不到未来，此时此刻的光辉，盼君勿忘。————世界上最幸福的女孩 珂朵莉，...
【LG-P5072 [Ynoi2015]】盼君勿忘
531 の 流 水 线
03-17 3194
做的第一道莫队题 QwQ。 P5072 [Ynoi2015] 盼君勿忘 给定一个序列，每次查询一个区间 [l,r][l,r][l,r] 中所有子序列分别去重后的和  mod  p\bmod\ pmod p。 思路 莫队 + 手写双向链表 1 题意有点绕，但注意是“分别去重”。 那么假设一个数 xxx，它在长为 nnn 的序列中出现了 mmm 次，即一共有 2n2^n2n 个子序列，不含 xxx 的子序列有 2n−m2^{n-m}2n−m 个。 则 xxx 在所有去重后的子序列中出现了 (
Ynoi(5/35)
weixin_30664539的博客
06-15 1156
luoguP5309 [Ynoi2012]D1T1 luoguP5310 [Ynoi2012]D1T2 luoguP5311 [Ynoi2012]D1T3 luoguP5312 [Ynoi2012]D2T1 luoguP5313 [Ynoi2012]D2T2 luoguP5314 [Ynoi2012]D2T3 luoguP5062 [Ynoi2014]在太阳西斜的这个世界里 luog...
题解 P4117 【[Ynoi2018]五彩斑斓的世界】
09-12 455
题目链接 我觉得AVX2指令集不够爽,于是写了AVX512.到官网查了一天手册，直接拿下最优解 Solution [Ynoi2018]五彩斑斓的世界 题目大意:给定一个长度为\(n\)的序列.每次将\([l,r]\)内大于\(x\)的数减\(x\).询问\([l,r]\)内\(x\)出现了多少次 分析:还能有啥分析,直接暴力,指令集优化一下就可以了 所有函数都可以在Intel手册查到 首...
P4690 [Ynoi2016] 镜中的昆虫
06-01
TokenType.SYMBOL_LEFT_PAREN) { match(TokenType.SYMBOL_LEFT_PAREN); expression(); match(TokenType.SYMBOL_RIGHT_PAREN); } else { throw new Exception("Syntax error"); ... private void match(TokenType ...
Ynoi 做题记录
最新发布
FuDuoDuo_qwq的博客
10-25 117
为了方便处理，我们先将相邻两个数对应区间的答案作为线段树叶子结点对应的区间。然后在线段树上维护一个类似 dp 的数组。，表示区间最左端和最右端的数 没有被覆盖 / 被覆盖 / 是极大值。所以只需要关注操作中心的集合就可以了。然后我们考虑用线段树维护。然后类似 dp 地在线段树上合并。减一的操作一直进行到。，不会改变最终答案。
[Ynoi2018]天降之物
StaroForgin的博客
08-19 198
「 弑尽破净的第四分块」，看来出题人越来越中二了。
【[Ynoi2018]五彩斑斓的世界】
dichanpan2699的博客
02-10 373
分块毒瘤题。如果真的想好好练习思维以及代码的话就请不要使用\(\text{WC}\)讲过的黑科技指令集。 正文部分： 由乃题必定分块。 我们可以将同一个块内值相同的用并查集维护起来，因为如果是一个块内相同的值打标机怎么改都一样。 接着谈修改，我们设块内最大值为\(\text{mx}\),修改的值为\(\text{v}\) 如果\(v*2>mx\)的话就改大于\(\text{v}\)...
bzoj5145 [Ynoi2018]未来日记 (多校第4场1013 Yuno and Claris)
yfzcsc的博客
01-23 2429
2018目标：让nzhtl1477nzhtl1477打通罚抄！这题其实不难写，就是我写完后MLEMLE了。。。对权值和序列都分块，假设nn和值域同阶记b[i][j]b[i][j]表示前ii块，权值在第jj块的有几个记c[i][j]c[i][j]表示前ii块，权值为jj的有几个预处理这个是O(nn√)O(n\sqrt{n})对于查询：直接暴力每块，O(n√)O(\sqrt{n})对于修改：两边的块暴力
[Ynoi2014]不归之人与望眼欲穿的人们
02-17 423
题目大意： 给定一个序列，每次单点修改一个数，或给定$x$，询问最短的or起来大于等于$x$的区间的长度（不存在输出-1）。 解题思路： 在太阳西斜的这个世界里，置身天上之森。等这场战争结束之后，不归之人与望眼欲穿的众人， 人人本着正义之名，长存不灭的过去、逐渐消逝的未来。我回来了，纵使日薄西山，即便看不到未来，此时此刻的光辉，盼君勿忘。————世界上最幸福的女孩 盼你归来，珂朵...
[Ynoi2012]D1T3
weixin_34246551的博客
05-28 272
https://www.luogu.org/problemnew/show/P5311 题解 先把点分树建出来。 对于吗，每一个询问\((l,r,x)\),我们对于x要找到它在点分树上最靠上的父亲节点使得两点之间的点在\(l \sim r\)中。 然后问题就变成了从一个根出发，可以经过\(l \sim r\)的点，能访问的颜色个数。 那么因为有了点分树，所以我们对于每个点\(dfs\)它的子树的复...
[Ynoi2014]人人本着正义之名（题解）
蒟蒻QMQ
10-23 718
我们不打代码，我们只是思想的搬运工。 题意 [Ynoi2014]人人本着正义之名 背景 略 题目描述 你需要帮珂朵莉维护一个长为n的01序列a，有m个操作： 1 l r : 把区间[l,r]的数变成0 2 l r : 把区间[l,r]的数变成1 3 l r : [l,r-1]内所有数a[i],变为a[i]与a[i+1]按位或的值,这些数同时进行这个操作 4 l r : [l+1,r]内所有数a[i...
【Ynoi2011】D1T3（点分树）（树状数组）
zxyoi_dreamer的博客（不定期诈尸）
10-01 401
传送门 题解： 考虑保留[l,r][l,r][l,r]的点的时候xxx所在的联通块怎么表示。 考虑点分治的时候第一次把分治中心放到这个联通块里面的那次分治。 取出这个分治中心，则统计颜色就变成了在它子树内部，只通过编号在[l,r][l,r][l,r]内的点能够到达多少种不同颜色。 离线，二维数点，树状数组维护一下就没了。 代码： #include<bits/stdc++.h> #d...
[YNOI2017]bzoj 4811 由乃的OJ - 树剖 - 线段树
Mys_C_K的博客
08-03 271
树剖，线段树维护每一位经过这个区间从左到右和从右到左0和1分别变成了啥，三个log稳T。发现其实只需要维护00000…000和1111…111会变成啥即可，这样就两个log了。 #include&amp;lt;iostream&amp;gt; #include&amp;lt;cstring&amp;gt; #include&amp;lt;cstdio&amp;gt; #include&amp;lt;algorithm&amp;gt; #inc
Ynoi2019模拟赛划水记
weixin_30896825的博客
03-19 162
几天前被各种多项式搞自闭了，就去\(Ynoi2019模拟赛\)划了一波水，简单记一下（不过代码是不会有的，因为我一个题都没A而且看了讲评完全没有改题的欲望） T1 正解码量巨大而我只写了暴力，但是居然过了\(50pts\)，为什么呢……回头审视代码发现原来是我习惯性卡常而且自己的代码框架里有这么一句： #pragma GCC optimize(3, "inline", "Ofast") 于是多水了...
[Ynoi2016]炸脖龙I/Nephren Ruq Insania
NephrenRuqInsania
05-13 323
Ynoi竟然出了个奈芙莲的题，少见啊 Ynoi竟然有个不卡常的题，更少见了 这道题其实思路非常好想，本质上就是一个暴力，但是要用到拓展欧拉定理，是什么呢？ ac≡{ac,c<ϕ(p)ac mod ϕ(p)+ϕ(p),c≥ϕ(p) mod pa^c\equiv\begin{cases}a^c,c<\phi(p) \\ a^{c\bmod \phi(p)+\phi(p)},c\geq\phi(p)\end{cases} \bmod pac≡{ac,c<ϕ(p)acmodϕ(p)+ϕ(p),c
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值