[BZOJ4542] [JZYZOJ2014][Hnoi2016] 大数(莫队+离散化)

正经题解在最下面

http://blog.csdn.net/qq_32739495/article/details/51286548

写的时候看了大神的题解[就是上面那个网址],看到下面这段话

观察题目，发现一串数s(l~r)整除p满足s(l~n-1)%p==s(r+1~n-1)%p 
但p值为2或5不满足这个性质需要特判（不过数据中好像没有，于是笔者没写，有兴趣的可以自己去写写。。。。。。）

然后问题转化为求一段区间中有几对相等的f值。

看到这里,我感觉豁然开朗,完全忽视了离散化的要求,我以为把余数值存起来扫一遍就行了离散个p啊..

写着写着完全参透这道题之后发现离散化的是余数啊,你不离散化怎么存数量啊,不存某个余数数量硬扫肯定超时啊....

然后我暴力硬扫果然[dian]超时了.........

 然后老老实实写离散化..........

 最重要的:2和5要特判

 

更重要的:离散化的时候要注意判定0的情况...即等于0时特判,不等0时离散化的赋值不应该从0开始,看加注释的那一段即可,不然会像我一样不停错两个点.....

我的程序200+大牛程序100-行..被吊着打.......我觉得我写的还挺清晰的...虽然完全不简洁

#include<iostream>
#include<cstring>
#include<cmath>
#include<cstdio>
#include<algorithm>
using namespace std;
int p,m,s=0,sz;
int a[100010]={};
int bel[100010]={};
int ans[100010]={};
long long mo[100010]={};
long long b[100010]={};
long long re[10][100010]={};
int vis[100010]={};
int tot[10]={};
struct nod{
    int x,y;
    int id;
}e[100010];
void readin(){
    char c=getchar();
    while(c<'0'||c>'9'){
        c=getchar();
    }
    while(c>='0'&&c<='9'){
        a[++s]=(int)(c-'0');
        c=getchar();
    }
}
bool mmp(nod aa,nod bb){
    if(bel[aa.x]==bel[bb.x]){
        if(aa.y==bb.y){
            return aa.x<bb.x;
        }
        return aa.y<bb.y;
    }
    return bel[aa.x]<bel[bb.x];
}
void work(){
    int l=1,r=0;
    int an=0;
    for(int i=1;i<=m;i++){
        while(l>e[i].x){
            l--;
            an+=vis[mo[l]];
            if(mo[r+1]==mo[l]){
                an+=1;
            }
            vis[mo[l]]++;
        }
        while(r<e[i].y){
            r++;
            vis[mo[r]]++;
            an+=vis[mo[r+1]];
        }
        while(l<e[i].x){
            vis[mo[l]]--;
            an-=vis[mo[l]];
            if(mo[r+1]==mo[l]){
                an-=1;
            }
            l++;
        }
        while(r>e[i].y){
            an-=vis[mo[r+1]];
            vis[mo[r]]--;
            r--;
        }
        ans[e[i].id]=an;
    }
    for(int i=1;i<=m;i++){
        printf("%d\n",ans[i]);
    }
}
void work5(){
    int l=1,r=0;
    int an=0;
    for(int i=1;i<=m;i++){
        while(l>e[i].x){
            l--;
            vis[a[l]%p]++;
            an+=vis[0];
        }
        while(r<e[i].y){
            r++;
            vis[a[r]%p]++;
            if(a[r]%p==0){
                an+=r-l+1;
            }
        }
        while(l<e[i].x){
            an-=vis[0];
            vis[a[l%p]]--;
            l++;
        }
        while(r>e[i].y){
            if(a[r]%p==0){
                an-=r-l+1;
            }
            vis[a[r]%p]--;
            r--;
        }
        ans[e[i].id]=an;
    }
    for(int i=1;i<=m;i++){
        printf("%d\n",ans[i]);
    }
}
int main(){
    //freopen("wtf.in","r",stdin);
    scanf("%d",&p);
    readin();
    sz=(int)sqrt((double)s);
    scanf("%d",&m);
    for(int i=1;i<=m;i++){
        scanf("%d%d",&e[i].x,&e[i].y);
        if(e[i].x>=s){
            e[i].x=s;
        }
        if(e[i].y>s){
            e[i].y=s;
        }
        e[i].id=i;
    }
    for(int i=1;i<=s;i++){
        bel[i]=(i-1)/sz+1;
    }
    sort(e+1,e+1+m,mmp);
    if(p==2||p==5){
        work5();
        return 0;
    }
    for(int i=0;i<=9;i++){
        re[i][1]=i%p;
        tot[i]=1;
    }
    for(int i=s,w=1;i>=1;i--){
        int x=a[i];
        while(tot[x]<w){
            tot[x]++;
            re[x][tot[x]]=re[x][tot[x]-1]*10%p;
        }
        mo[i]=(re[x][w]+mo[i+1])%p;
        b[i]=mo[i];
        w++;
    }
    sort(b+1,b+1+s);
    int size=unique(b+1,b+1+s)-b-1;
    for(int i=1;i<=s;i++){//离散化部分....注意一定要特判..
        if(mo[i]==0){
            mo[i]==0;
        }
        else{
            mo[i]=lower_bound(b+1,b+1+s,mo[i])-b;//这里-1且输入字符串中没有0时,1会被离散化为0
        }
    }
    work();
    return 0;
}