【bzoj5452】[Hnoi2016]大数（莫队）

　　题目传送门：https://www.lydsy.com/JudgeOnline/problem.php?id=4542

　　首先若p=2,5则这题就是道傻逼题，前缀和搞一下没了。如果p为其他质数，那么可以这么处理：

　　我们先预处理出数组num[i]表示原串第i~n位表示的数模p的余数，那么第l~r位表示的数模p的余数为(num[l]-num[r+1])/10^(n-r)，因为10^(n-r)与p互质，所以若num[l]=num[r+1]，则第l~r位表示的数是p的倍数。于是莫队一下就好了。

　　代码：

#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<string>
#include<iostream>
#include<algorithm>
#include<queue>
#include<vector>
#include<map>
#define ll long long
#define ull unsigned long long
#define max(a,b) (a>b?a:b)
#define min(a,b) (a<b?a:b)
#define lowbit(x) (x& -x)
#define mod 1000000007
#define inf 0x3f3f3f3f
#define eps 1e-18
#define maxn 2000010
inline ll read(){ll tmp=; char c=getchar(),f=; for(;c<''||''<c;c=getchar())if(c=='-')f=-; for(;''<=c&&c<='';c=getchar())tmp=(tmp<<)+(tmp<<)+c-''; return tmp*f;}
inline ll power(ll a,ll b){ll ans=; for(;b;b>>=){if(b&)ans=ans*a%mod; a=a*a%mod;} return ans;}
inline ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
inline void swap(int &a,int &b){int tmp=a; a=b; b=tmp;}
using namespace std;
struct data{
    int l,r,id;
}a[];
struct data2{
    ll val;
    int id;
}x[];
char s[];
ll rk[],cnt[],ans[],sum1[],sum2[];
int n,m,size;
ll p,tot;
bool cmp(data a,data b){return a.l/size!=b.l/size?a.l/size<b.l/size:a.r<b.r;}
bool cmp2(data2 a,data2 b){return a.val<b.val;}
int main()
{
    p=read();
    scanf("%s",s); n=strlen(s); size=sqrt(n);
    m=read();
    if(p==||p==){
        sum1[]=sum2[]=;
        for(int i=;i<=n;i++){
            sum1[i]=sum1[i-]; sum2[i]=sum2[i-];
            if((s[i-]-'')%p==)++sum1[i],sum2[i]+=i;
        }
        for(int i=;i<=m;i++){
            int l=read(),r=read();
            printf("%lld\n",sum2[r]-sum2[l-]-(sum1[r]-sum1[l-])*(l-));
        }
        fclose(stdin); fclose(stdout);
        return ;
    }
    for(int i=;i<=m;i++)
        a[i].l=read()-,a[i].r=read(),a[i].id=i;
    sort(a+,a+m+,cmp);
    ll tmp=; x[n].val=; x[n].id=n;
    for(int i=n-;i>=;i--,tmp=tmp*%p)x[i].val=(x[i+].val+(s[i]-'')*tmp)%p,x[i].id=i;
    sort(x,x+n+,cmp2);
    rk[x[].id]=;
    for(int i=;i<=n;i++)
        if(x[i].val==x[i-].val)rk[x[i].id]=rk[x[i-].id];
        else rk[x[i].id]=i;
    tot=;
    for(int i=a[].l;i<=a[].r;i++)
        tot+=cnt[rk[i]]++;
    ans[a[].id]=tot;
    for(int i=;i<=m;i++){
        if(a[i-].l<a[i].l){
            for(int j=a[i-].l;j<a[i].l;j++)
                tot-=--cnt[rk[j]];
        }
        else{
            for(int j=a[i].l;j<a[i-].l;j++)
                tot+=cnt[rk[j]]++;
        }
        if(a[i-].r<a[i].r){
            for(int j=a[i-].r+;j<=a[i].r;j++)
                tot+=cnt[rk[j]]++;
        }
        else{
            for(int j=a[i].r+;j<=a[i-].r;j++)
                tot-=--cnt[rk[j]];
        }
        ans[a[i].id]=tot;
    }
    for(int i=;i<=m;i++)
        printf("%lld\n",ans[i]);
}

bzoj4542